2017-05-16 8 views
2

. 나는 임의 효과없이 LM를 사용하는 경우그람의 줄거리는? 나는 패키지 NLME에서 LME 기능을 사용하고</p> <pre><code>m <- lme(y ~ Time, random = ~1|Subject) </code></pre> <p>를 작성하고 다음</p> <pre><code>Variogram(m, form = ~Time|Subject) </code></pre> <p>는 베리오 그램 아무 문제를 생산을 작성하는 경우

그러나

m <- lm(y ~ Time) 

과는 무슨 일이야

Variogram(m, form = ~Time) 

Error in Variogram.default(m, form = ~Time) : 
    argument "distance" is missing, with no default 

을 생산 쓰기? 나는 lm을 사용하기 전에 그것을 필요로하지 않을 때, 내가 영화에 적합 할 때 왜 거리가 필요한가?

그런데 어떻게 "거리"를 지정하지 않고도 Variogram을 그립니다. 나는 다른 모델링 방법을 사용하여 동일한 문제가 : GLM, 할머니, GAMM 등을

편집 :

당신은 예를 들어, 사용하여 자신이 모든 것을 확인할 수 있습니다 nlme의 BodyWeight 데이터 nlme에서

> m <- lm(weight ~ Time, data = BodyWeight) 
> Variogram(m, form =~Time) 
Error in Variogram.default(m, form = ~Time) : 
    argument "distance" is missing, with no default 
+0

예, 그게 내가 무슨 짓을했는지. 나는 실수로 제목을 썼다. – Maron

답변

1

lme 맞춤하는 방법 Variogram.lme 함수가 있지만 lm 모델 동등한 방법은 없다. 다음과 같이

당신은 Variogram.default을 사용할 수 있습니다 : 우리는 또한 베리오 그램을 플롯 할 수 있습니다

library(nlme) 
mod1 <- lm(weight ~ Time, data = BodyWeight) 
n <- nrow(BodyWeight) 
variog <- Variogram(resid(mod1), distance=dist(1:n)) 
head(variog) 

############ 
     variog dist 
1 17.4062805 1 
2 23.1229516 2 
3 29.6500135 3 
4 15.6848617 4 
5 3.1222878 5 
6 0.9818238 6 

:

plot(variog) 

enter image description here