এই প্রশ্নটি বিশেষভাবে সম্পর্কিত নয় R
, তবে আমি R
এটি চিত্রিত করার জন্য ব্যবহার করতে পছন্দ করেছি।
একটি (সাধারণ) কিউকি-লাইনের চারপাশে আত্মবিশ্বাস ব্যান্ড উত্পাদন করার কোডটি বিবেচনা করুন:
library(car)
library(MASS)
b0<-lm(deaths~.,data=road)
qqPlot(b0$resid,pch=16,line="robust")
আমি কীভাবে এই আত্মবিশ্বাস ব্যান্ডগুলি তৈরি করা হয় (বা আরক্সের ফাইলগুলিতে ফক্স 2002 এর একটি রেফারেন্স দেখেছি তবে দুঃখের সাথে আমার কাছে এটি নেই) এর ব্যাখ্যা (বা একটি কাগজ / অনলাইন ডকুমেন্টের সাথে একটি লিঙ্কের বিকল্পের ব্যাখ্যা করে) খুঁজছি বই সহজ)।
আমার প্রশ্ন একটি উদাহরণ সঙ্গে আরও সুনির্দিষ্ট করা হবে। R
এই নির্দিষ্ট সিআইয়ের গণনাগুলি কীভাবে করা হয়েছে (এতে ব্যবহৃত কোডটি আমি সংক্ষিপ্ত / সংক্ষিপ্ত করে রেখেছি 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<-(b/dnorm(z))*sqrt(P*(1-P)/n)
) গণনা করার জন্য ব্যবহৃত সূত্রটির ন্যায়সঙ্গততা কী ?
এফডাব্লুআইডাব্লিউআইডাব্লু লিনিয়ার রিগ্রেশন-এ ব্যবহৃত সাধারণ আত্মবিশ্বাস ব্যান্ডের সূত্র থেকে এই সূত্রটি খুব আলাদা