rpy2中rpy2模型系数的提取

Question

我需要使用rpy2提取python中调用的rpy2模型的结果。用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的系数？

Answer 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)