【发布时间】:2020-03-09 08:29:26
【问题描述】:
我正在尝试隔离数据集的哪些特征(即 Pandas DataFrame 的列)用于线性回归,并且我想选择那些不强相关的特征(假设自变量需要彼此不相关,因此我们要删除任何看似强相关的)。
我已经隔离了与我的目标变量相关的初始特征列表,如下所示:
# get the absolute correlations for the target variable
correlations_target = abs(df.corr()[target_variable_name])
# filter out those that are below our threshold
correlated_features = correlations_target[correlations_target >= correlation_threshold]
# drop the target variable's column
correlated_features.drop(target_variable_name, inplace=True)
# get the column names for later use
correlated_feature_variable_names = correlated_features.index
我现在想检查这些相关的特征变量中的每一个,并确保它们中没有一个具有很强的相关性,如果确实存在,则将与目标变量相关性最弱的那个丢弃。这是我为此准备的:
# the collection of feature variable names we'll drop due to their being correlated to other features
correlated_feature_variable_names_to_drop = []
# loop over the feature combinations
for name_1 in correlated_feature_variable_names:
for name_2 in correlated_feature_variable_names:
# only look at correlations between separate feature variables
if name_1 != name_2:
# drop one of the feature variables if there's a strong correlation
if abs(df[[name_1, name_2]].corr()[name_1][name_2]) > 0.6:
# only worry about it if neither of the variables have been added to the drop list
if (name_1 not in correlated_feature_variable_names_to_drop) and \
(name_2 not in correlated_feature_variable_names_to_drop):
# drop the one which has the least correlation to the target variable
if correlated_features[name_1] >= correlated_features[name_2]:
correlated_feature_variable_names_to_drop.append(name_2)
else:
correlated_feature_variable_names_to_drop.append(name_1)
# drop the variables we've found that qualify
correlated_features.drop(correlated_feature_variable_names_to_drop, inplace=True)
# get the remaining variables' column names for later use
filtered_feature_variable_names = correlated_features.index
过滤后的特征集将用作简单回归模型的输入。例如:
# fit a simple ordinary least squares model to the features
X = df[filtered_feature_variable_names]
y = df[target_variable_name]
estimate = sm.OLS(y, np.asarray(X)).fit()
# display the regression results
estimate.summary()
因为这是我第一次尝试这个,所以我不确定这是否是正确的方法,如果是,那么可能有一种更聪明或更有效(或“Pythonic”)的方法比我上面使用的循环方法执行过滤。
【问题讨论】:
标签: python pandas linear-regression