কিউকিউ লাইনের জন্য আত্মবিশ্বাস ব্যান্ড


14

এই প্রশ্নটি বিশেষভাবে সম্পর্কিত নয় 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)) গণনা করার জন্য ব্যবহৃত সূত্রটির ন্যায়সঙ্গততা কী ?

এফডাব্লুআইডাব্লিউআইডাব্লু লিনিয়ার রিগ্রেশন-এ ব্যবহৃত সাধারণ আত্মবিশ্বাস ব্যান্ডের সূত্র থেকে এই সূত্রটি খুব আলাদা


2
আমি আশা করি এটি আদেশের পরিসংখ্যান f x ( k ) ( x ) = n এর বিতরণের সাথে করবে !এবং আরও বিশেষতঅ্যাসিম্পোটিক ফলাফল:এক্স(nপি)এন(এফ-1(পি)
এক্স()(এক্স)=এন!(-1)!(এন-)![এফএক্স(এক্স)]-1[1-এফএক্স(এক্স)]এন-এক্স(এক্স)
এক্স(এনপি)~একজনএন(এফ-1(পি),পি(1-পি)এন[(এফ-1(পি))]2)
গ্লেন_বি -রিনস্টেট মনিকা

4
@ Glen_b ঠিক আছে। জন ফক্স পৃষ্ঠাগুলি 35-36 লিখেছেনঃ "অর্ডার পরিসংখ্যাত আদর্শ ত্রুটি হয় এস ( এক্স ( আমি ) ) = σএক্স(আমি)
এস(এক্স(আমি))=σ^পি(z- রআমি)পিআমি(1-পিআমি)এন
পি(z- র)পি(z- র)এক্স^(আমি)=μ^+ +σ^z- রআমিএক্স^(আমি)±2×এস(এক্স(আমি))

2
(এফ-1(পি))(পি(z- রআমি)/σ^)

উত্তর:


6

এক্স()(এক্স)=এন!(-1)!(এন-)![এফএক্স(এক্স)]-1[1-এফএক্স(এক্স)]এন-এক্স(এক্স)
এক্স(এনপি)~একজনএন(এফ-1(পি),পি(1-পি)এন[(এফ-1(পি))]2)

COOLSerdash মন্তব্যে যেমন উল্লেখ করেছে, জন ফক্স [১] ৩৫-৩6 পৃষ্ঠায় লিখেছেন:

এক্স(আমি)

এস(এক্স(আমি))=σ^পি(z- রআমি)পিআমি(1-পিআমি)এন
পি(z- র)পি(z- র)এক্স^(আমি)=μ^+ +σ^z- রআমিএক্স^(আমি)±2×এস(এক্স(আমি))

(এফ-1(পি)) দ্বারা অনুমান করা হয় (পি(z- রআমি)/σ^)

[1] ফক্স, জে। (২০০৮),
প্রয়োগিত রিগ্রেশন অ্যানালাইসিস এবং জেনারালাইজড লিনিয়ার মডেলস, ২ য় এড। ,
সেজ পাবলিকেশনস, ইনক

আমাদের সাইট ব্যবহার করে, আপনি স্বীকার করেছেন যে আপনি আমাদের কুকি নীতি এবং গোপনীয়তা নীতিটি পড়েছেন এবং বুঝতে পেরেছেন ।
Licensed under cc by-sa 3.0 with attribution required.