Python中使用Scipy的方差分析，无需手动输入每个组的nam

Question

我有一个像这样的测向仪

           Ens_prot_ID      Ens_gene_ID Sample     TPM      ppm   ppm/TPM
0      ENSP00000416240  ENSG00000109072  liver  2540.4    0.003  0.000001
21597  ENSP00000226218  ENSG00000109072  liver  2540.4  110.000  0.043300
...

下面的代码不工作并给出“KeyError:0”

^{pr2}$

我跑起来没问题

from scipy import stats
proteins=df['Ens_prot_ID'].unique()
for prot in proteins:
    df[df['Ens_prot_ID'] == prot]['ppm/TPM']

所以问题似乎是f帴oneway（）不喜欢我的输入形式。有没有办法让函数接受不是手工输入的组名？在

Answer 1

你可以这样做：

stats.f_oneway(*(df[df['Ens_prot_ID'] == prot]['ppm/TPM'] for prot in proteins))

*允许您将多个参数传递给函数：https://pythontips.com/2013/08/04/args-and-kwargs-in-python-explained/

我还将列表理解改为generator comprehension，以防你的DataFrame很大。在