2012-07-05 2 views
0

선을 예측 플롯 없습니다제가 그 데이터를 사용하고 R

'data.frame': 1584 obs. of 3 variables: 
$ Individual: Factor w/ 3 levels "AG201","AG202",..: 1 1 1 1 1 1 1 1 1 1 ... 
$ Used  : Factor w/ 2 levels "no","yes": 2 2 2 2 2 2 2 2 2 2 ... 
$ NDVI  : int 4724 4576 4894 4297 4670 4932 4346 3810 3481 4058 ... 

I 연속 설명 변수로서 "NDVI"와 GLM하고있어 한 다음 I 통해 모델을 플로팅있어 데이터의 산점도 (필자는 R 책과 같은 스크립트를 재현하고있어, 크롤리, p.596)

model<-glm(Used~NDVI,binomial); 
xv<-seq(0,10000,0.2); 
yv<-predict(model,list(NDVI=xv),type="response"); 
plot(NDVI,Used); 
lines(xv,yv); 

내 문제 ... 더 라인 내 그래프에 표시 없다는 것입니다

모든 무슨 일이야? 개빈의 통찰력에 따라

+0

yv를 보셨습니까? 축 범위 안에 범위 (xv)와 범위 (yv)가 있습니까? – Roland

+1

짐작할 수 있듯이 y 축에 대한 요소가 있습니다 ('str()'의 출력에서'$ Used' 참조). 이것은 당신이 생각하는 줄거리를 그려지지 않았습니다. 플롯을 그리는 사용자가'par ("usr")'**로 플롯의 좌표를 보았는지 확인한 다음'range (yv) '를 보면 그 점이 잘 알려줄 것입니다. 범위 0, ..., 1을 포함하는 연속적인 그림이 필요하며 Y 축에 계수를 그릴 때 그 값을 얻지 못할 수도 있습니다. –

+0

당신은 분명히'attach (dfrm)'도 실행했는데, 일반적으로 권하지는 않습니다 만,'attach'를 사용하려면 게시 된 코드와 함께 포함시켜야합니다. –

답변

1

, 여기 제안입니다 :

plot(NDVI, as.numeric(Used)-1) 
lines(xv,yv) 

요인 라벨의 알파 위해 기본적으로 할당과 1L에서 시작 정수 벡터입니다. 따라서 "아니오"< "예"로 설정하면 No가 1로되고 hte Yes가 2가되고 올바른 배율 [0,1]로 이동합니다. 봐야 할 수도 있습니다 str(yv)

관련 문제