rpy2의 상관 구조 corAR1() "정의되지 않음"일반화 최소 자승 모델 gls

Question

파이썬에서 rpy2 모듈을 사용하여 두 timeseries (ts1과 ts2) 간의 회귀를 계산합니다. 잔차는 자기 상관 (autocorrelated)되어 있으므로, 작품보다는 gl 모델을 사용해야합니다. 나는 아래의 코드에서와 같이 corAR1()으로 내 상관 구조를 설정할 수 있어야한다.

import rpy2.robjects as robjects 
from rpy2.robjects import FloatVector 
from rpy2.robjects.packages import importr 

nlme = importr('nlme') 

y = FloatVector(ts1) 
x = FloatVector(ts2) 
fmla = robjects.Formula('y ~ x - 1') 
env = fmla.environment 
env['x'] = x 
env['y'] = y 
fit = nlme.gls(fmla, cor=corAR1(value=c(0.5))) 

는 그러나, 나는 오류를 얻을

내가 R 사용자와 RPY2 아주 새로운 아니에요으로

*** NameError: name 'corAR1' is not defined                          

, 내가 여기에 무슨 일이 일어나고 있는지 확실하지 않다! 어떤 아이디어라도 높이 평가할 수 있습니다.

건배, 경사는

Answer 1

importr('nlme')를 호출, 반환 된 객체는 파이썬 패키지/네임 스페이스 같다. corAR1()가 nlme 패키지에 정의되어 경우, 당신은이 있음을 파이썬을 말해야한다 :

fit = nlme.gls(fmla, cor=nlme.corAR1(value=c(0.5)))