4

np パッケージ (npreg) を使用してノンパラメトリック回帰を実行し、他のすべての変数を平均値/モードで保持している対象変数 x1 の結果をプロットしようとしました。

library("np")

y  <- rnorm(100)
x1 <- rnorm(100,10,30)
x2 <- rbinom(100,1,0.5)
x3 <- rbinom(100,1,0.5)

model.np <- npreg(y ~ x1 + x2 + x3)
plot(model.np)

プロットはまさに私が望むものですが、それらを「手動で」個別に生成する方法がわかりません。特に、(3 つのうちの) 最初の出力プロットのみが必要です。

4

1 に答える 1

2

明らかに、詳細な回答は npplot-routine のヘルプ ファイルにあり、重要な引数は plot.behavior です。

私の例では、x1-graph のみをプロットするには、次の方法で実行できます。

nlmodel.plot <- plot(model.np, plot.behavior="data")

y.eval <- fitted(nlmodel.plot$r1) #fit partial regression model for r1=airnoise
y.se <- se(nlmodel.plot$r1) #grab SE from botstrap 
y.lower.ci <- y.eval + logp.se[,1] #lower CI
y.upper.ci <- y.eval + logp.se[,2] #upper CI

x1.eval <- nlmodel.plot$r1$eval[,1]  #grab x1 values saved in plot$r1

plot(x1,y)
 lines(x1.eval,y.eval)
 lines(x1.eval,y.lower.ci,lty=3)
 lines(x1.eval,y.upper.ci,lty=3)
于 2012-11-15T16:22:21.273 に答える