आर: आत्मविश्वास दृश्य पैकेज ग्रहण

नमस्ते

आर-हब से अंतिम पोस्ट में "एक विमान पर दो आयामी गॉसियन का दृश्य" , एक सहसंयोजक मैट्रिक्स से एक आत्मविश्वास दीर्घवृत्त के निर्माण के लिए एक एल्गोरिथ्म का वर्णन किया गया था। एल्गोरिथ्म एक उदाहरण और एक आर-स्क्रिप्ट के साथ था।

शायद Gaussian mephistopheies के विज़ुअलाइज़ेशन पर पोस्ट के लेखक और आर-हब के पाठकों को निम्नलिखित जानकारी उपयोगी मिलेगी। आर रिपॉजिटरी में एक दीर्घवृत्त पैकेज है। इस पैकेज में विश्वास क्षेत्रों के दीर्घवृत्तों के निर्माण की विभिन्न प्रक्रियाएँ हैं।

एक उदाहरण पर विचार करें।

एक दो आयामी सामान्य वितरण से नमूना

दो-आयामी सामान्य वितरण से एक नमूना उत्पन्न करने के लिए, हम mvtnorm पैकेज का उपयोग करते हैं।

हम म्यू और एक सहसंयोजक मैट्रिक्स सिग्मा के औसत मूल्यों के वेक्टर के साथ 1000 तत्वों के डेटा नमूने का निर्माण करते हैं:

require(mvtnorm) mu <- c(1,2) sigma <- matrix(c(4,2,2,3), ncol=2) set.seed(100) data <- rmvnorm(n=1000, mean=mu, sigma=sigma) #      

एलीपस 95% कॉन्फिडेंस एरिया

सिग्मा सहसंयोजक मैट्रिक्स के लिए, हमें 95% विश्वास क्षेत्र के दीर्घवृत्त के 100 बिंदुओं के निर्देशांक मिलते हैं:

 require(ellipse) confidence.ellipse <- ellipse(sigma,centre=mu,level=0.95,npoints=100) 

दृश्य

विमान पर डेटा बिंदुओं और आत्मविश्वास को प्रदर्शित करें

 plot(data,pch=19,col=rgb(0, 0.5, 1, 0.2), xlab="x", ylab="y", xlim=range(data[,1]), ylim=range(data[,2])) par(new=TRUE) plot(confidence.ellipse,type="l", xlab="",ylab="", xlim=range(data[,1]),ylim=range(data[,2])) 

परिणाम: