QQラインの信頼帯
この質問は特にに関係するものRではありませんがR、説明のために使用することにしました。 (通常の)qqラインの周りに信頼帯を生成するコードを考えます。 library(car) library(MASS) b0<-lm(deaths~.,data=road) qqPlot(b0$resid,pch=16,line="robust") 私はこれらの信頼帯がどのように構築されるかの説明(または代替の紙/オンライン文書へのリンク)を探しています(RのヘルプファイルでFox 2002への参照を見ましたが、悲しいことに私はこれを持っていません便利な本)。 私の質問は例を使用してより正確になります。Rこれらの特定のCIの計算方法は次のとおりです(で使用するコードを短縮/簡略化しましたcar::qqPlot)。 x<-b0$resid good<-!is.na(x) ord<-order(x[good]) ord.x<-x[good][ord] n<-length(ord.x) P<-ppoints(n) z<-qnorm(P) plot(z,ord.x,type="n") coef<-coef(rlm(ord.x~z)) a<-coef[1] b<-coef[2] abline(a,b,col="red",lwd=2) conf<-0.95 zz<-qnorm(1-(1-conf)/2) SE<-(b/dnorm(z))*sqrt(P*(1-P)/n) #[WHY?] fit.value<-a+b*z upper<-fit.value+zz*SE lower<-fit.value-zz*SE lines(z,upper,lty=2,lwd=2,col="red") lines(z,lower,lty=2,lwd=2,col="red") 問題は、これらのSEを計算するために使用される式の正当化とは何ですか(例:line SE<-(b/dnorm(z))*sqrt(P*(1-P)/n))。 FWIWこの式は、線形回帰で使用される通常の信頼帯の式とは大きく異なります