0

私はそれらのデータを使用しています:

'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 ...

「NDVI」を連続説明変数として glm を実行してから、データの散布図を使用してモデルをプロットしています (R の本、Crawley、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);

私の問題は、グラフに線が表示されないことです...

何が問題なのですか？

4

1 に答える 1

1

Gavinの洞察に続いて、ここに提案があります：

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

1L係数は、デフォルトでラベルのアルファ順の割り当てで始まる整数ベクトルです。したがって、「no」<「yes」でOKになり、Noが1になり、Yesが2になり、正しいスケール[0,1]にシフトダウンします。あなたも見る必要があるかもしれませんstr(yv)

于 2012-07-05T15:33:47.310 に答える