ポワソン回帰で求めた傾きと切片のSE(標準誤差)の求め方が分かりません。Rのglm(poisson, log)を使用しました。パラメータのSEは、共分散行列から求められることは分かりましたが、どうしてそのようになるのか、また、共分散行列からではなく、手計算で求めるにはどのような式になるのか分かりません。どなたかアドバイスをお願いいたします。
x<-c(2,2,3,3,3,4,4) y<-c(15,9,13,10,7,11,5) model<-glm(y~x,family=poisson(link=log)) (result<-summary(model)) # Estimate Std. Error z value Pr(>|z|) #(Intercept) 2.8940 0.4749 6.094 1.1e-09 *** #x -0.2010 0.1593 -1.262 0.207 (SEvcov<-diag(sqrt(vcov(model))))#parameter SE #(Intercept) x # 0.4748514 0.1592545
あなたの回答
tips
プレビュー