3

多変量 gls モデルを実行しています:

m <- gls(y ~ x + factor1 + factor2, cor = corPagel(1,phylogeny), weight= ~1/log(n))

結果をプロットしたいのですが、次のような予測値を取得できます。

newdata <- data.frame(expand.grid( x = mean(x), factor1= unique(factor1), factor2 = unique(factor2)))

predvals <- predict(m,newdata)

プロットしたい 2 つの変数は因子ですが、これらの予測値の信頼区間または標準誤差を取得したいと考えています。

誰もそれを行う方法を知っていますか?

よろしくお願いします。乾杯、ジャスミン

4

0 に答える 0