从rpy2中的R lme模型中提取系数答案

【问题标题】：extract coefficients from R lme model in rpy2从rpy2中的R lme模型中提取系数
【发布时间】：2016-02-14 18:29:23
【问题描述】：

我需要使用 rpy2 提取在 python 中调用的 R lme 模型的结果。用 rx2 提取系数如下：

model = nlme.lme(fixed=fixed, data=dfr, random=random, method="REML")
print model.rx2("coefficients")

产生所有系数如下：

$fixed
(Intercept)   log.var1 
  12.571692   -2.929928 

$random
$random$item
                (Intercept)    log.var1
12545646546 -5.189606e-16 8.276929e-16

$random$category
                 (Intercept)  log.var1
0001544848484/DLMX  -3.1909917  2.3938670

我想在模型的固定部分提取log.var1 的系数。我尝试了以下，但我得到了 NULL

print model.rx2("coefficients$fixed[2]")
#gives NULL

如何获得 log.var1 的系数？

【问题讨论】：

标签： python r rpy2

【解决方案1】：

方法rx2对应R的[[，据我了解和$是一样的。

考虑到这一点，您可以像访问“系数”一样访问“系数”中称为“固定”的项目：

# Python sequences are zero-based
model.rx2("coefficients").rx2("fixed")[1]

或

# R vectors are one-based (so 1+1=2)
model.rx2("coefficients").rx2("fixed").rx(2)

【讨论】：