如何迭代熊猫数据框的列以运行回归

Posted

技术标签:

【中文标题】如何迭代熊猫数据框的列以运行回归【英文标题】:How to iterate over columns of pandas dataframe to run regression 【发布时间】:2015-03-28 21:56:08 【问题描述】:

我确信这很简单,但作为一个完整的 python 新手,我无法弄清楚如何迭代 pandas 数据帧中的变量并对每个变量运行回归。

这就是我正在做的事情:

all_data = 
for ticker in ['FIUIX', 'FSAIX', 'FSAVX', 'FSTMX']:
    all_data[ticker] = web.get_data_yahoo(ticker, '1/1/2010', '1/1/2015')

prices = DataFrame(tic: data['Adj Close'] for tic, data in all_data.iteritems())  
returns = prices.pct_change()

我知道我可以像这样运行回归:

regs = sm.OLS(returns.FIUIX,returns.FSTMX).fit()

但假设我想对数据框中的每一列执行此操作。特别是,我想在 FSTMX 上回归 FIUIX,然后在 FSTMX 上回归 FSAIX,然后在 FSTMX 上回归 FSAVX。每次回归后,我都想存储残差。

我尝试了以下各种版本,但我一定是语法错误:

resids = 
for k in returns.keys():
    reg = sm.OLS(returns[k],returns.FSTMX).fit()
    resids[k] = reg.resid

我认为问题是我不知道如何按键引用returns列,所以returns[k]可能是错误的。

任何关于最佳方式的指导将不胜感激。也许我缺少一种常见的 pandas 方法。

【问题讨论】:

你可以像这样给列下标:for i in len(df): if i + 1 != len(df): # sm.OLS(returns[returns.coloumns[i]], returns[returns.columns[ i+1]]), fit()os similar 【参考方案1】:

假设 X 因子,y 标签(多列):

columns = [c for c in _df.columns if c in ['col1', 'col2','col3']]  #or '..c not in..'
_df.set_index(columns, inplace=True)
print( _df.index)

X, y =  _df.iloc[:,:4].values, _df.index.values

【讨论】:

【参考方案2】:

我只是在寻找一个干净的列迭代器(Series,没有名字)时遇到了这个问题。

除非我弄错了,否则没有这样的事情,如果是真的,那就有点烦人了。特别是,有时希望将几个单独的列(系列)分配给变量,例如:

x, y = df[['x', 'y']]  # does not work

有一个接近的df.items(),但它给出了一个元组迭代器(column_name, column_series)。有趣的是,有一个对应的df.keys() 返回df.columns,即列名称为Index,因此a, b = df[['x', 'y']].keys() 正确分配a='x'b='y'。但是没有对应的df.values(),并且有充分的理由,因为df.values 是一个属性并返回底层numpy 数组。

一种(不优雅的)方法是:

x, y = (v for _, v in df[['x', 'y']].items())

但它不像我想要的那样pythonic。

【讨论】:

嘿@Pierre D 我遇到了你的答案并且正在寻找类似的东西。我不知道this link 是否有帮助,但它可能值得一看。【参考方案3】:

你可以使用iteritems():

for name, values in df.iteritems():
    print('name: value'.format(name=name, value=values[0]))

【讨论】:

很好的答案。顺便说一句,df.iteritems() 也可以写成df.items() 给出相同的结果。【参考方案4】:

这个答案是迭代 选定的列 以及 DF 中的所有列。

df.columns 给出一个包含 DF 中所有列名的列表。现在,如果您想遍历所有列,这不是很有帮助。但是当您只想迭代您选择的列时,它会派上用场。

我们可以很容易地使用 Python 的列表切片来根据我们的需要对 df.columns 进行切片。例如,要遍历除第一列之外的所有列,我们可以这样做:

for column in df.columns[1:]:
    print(df[column])

类似地以相反的顺序遍历所有列,我们可以这样做:

for column in df.columns[::-1]:
    print(df[column])

我们可以使用这种技术以很酷的方式遍历所有列。另请记住,您可以使用以下方法轻松获取所有列的索引:

for ind, column in enumerate(df.columns):
    print(ind, column)

【讨论】:

【参考方案5】:

基于the accepted answer,如果每列对应的index也是desired

for i, column in enumerate(df):
    print i, df[column]

上面的df[column]类型是Series,可以简单转换成numpyndarrays:

for i, column in enumerate(df):
    print i, np.asarray(df[column])

【讨论】:

【参考方案6】:

我有点晚了,但我是这样做的。步骤:

    创建所有列的列表 使用 itertools 取 x 个组合 将每个结果 R 平方值与排除的列列表一起附加到结果数据框 按照 R 平方的降序对结果 DF 进行排序,看看哪个最合适。

这是我在 DataFrame 上使用的代码aft_tmt。随意推断您的用例..

import pandas as pd
# setting options to print without truncating output
pd.set_option('display.max_columns', None)
pd.set_option('display.max_colwidth', None)

import statsmodels.formula.api as smf
import itertools

# This section gets the column names of the DF and removes some columns which I don't want to use as predictors.
itercols = aft_tmt.columns.tolist()
itercols.remove("sc97")
itercols.remove("sc")
itercols.remove("grc")
itercols.remove("grc97")
print itercols
len(itercols)

# results DF
regression_res = pd.DataFrame(columns = ["Rsq", "predictors", "excluded"])

# excluded cols
exc = []

# change 9 to the number of columns you want to combine from N columns.
#Possibly run an outer loop from 0 to N/2?
for x in itertools.combinations(itercols, 9):
    lmstr = "+".join(x)
    m = smf.ols(formula = "sc ~ " + lmstr, data = aft_tmt)
    f = m.fit()
    exc = [item for item in x if item not in itercols]
    regression_res = regression_res.append(pd.DataFrame([[f.rsquared, lmstr, "+".join([y for y in itercols if y not in list(x)])]], columns = ["Rsq", "predictors", "excluded"]))

regression_res.sort_values(by="Rsq", ascending = False)

【讨论】:

【参考方案7】:

使用列表推导,您可以获得所有列名(标题):

[column for column in df]

【讨论】:

短版:list(df.columns)[c for c in df]【参考方案8】:
for column in df:
    print(df[column])

【讨论】:

我似乎只有在使用这种方法时才能取回列标题。例如: print(df) 向我显示数据框列中的数据,但对于 df 中的 c: print(c) 仅打印标题而不是数据。 好吧忽略我——我在做 print(column) 而不是 print (df[column]) 注意同名列! 简洁明了。不过,我希望for x in df 能够遍历行。 :-/ for idx, row in df.iterrows() 遍历行。由于基于 col 的操作是矢量化的,因此主迭代自然是在列上 :)【参考方案9】:

一种解决方法是转置DataFrame 并遍历行。

for column_name, column in df.transpose().iterrows():
    print column_name

【讨论】:

换位相当昂贵:) 可能很昂贵,但对于相对较小的数据帧来说,这是一个很好的解决方案。谢谢 kdauria!【参考方案10】:

您可以使用 ix 按位置索引数据框列。

df1.ix[:,1]

例如,这将返回第一列。 (0 是索引)

df1.ix[0,]

这将返回第一行。

df1.ix[:,1]

这将是第 0 行和第 1 列交叉处的值:

df1.ix[0,1]

等等。所以你可以enumerate()returns.keys(): 并使用数字来索引数据框。

【讨论】:

ix 已弃用,请使用 iloc

以上是关于如何迭代熊猫数据框的列以运行回归的主要内容,如果未能解决你的问题,请参考以下文章

规范化熊猫数据框的列

如何在以下数据框的列之间进行迭代? [复制]

如何组合熊猫数据透视表中的列?

如何删除熊猫数据框的最后一列数据

如何加快熊猫数据框迭代

如何在python中获取熊猫数据框的行列表? [复制]