R में द्विपद बंटन कैसे आलेखित करें
आर में द्विपद वितरण के लिए संभाव्यता द्रव्यमान फ़ंक्शन को प्लॉट करने के लिए, हम निम्नलिखित फ़ंक्शन का उपयोग कर सकते हैं:
- संभाव्यता द्रव्यमान फ़ंक्शन बनाने के लिए dbinom(x, size, prob)।
- प्लॉट (x, y, टाइप = ‘h’) प्रायिकता द्रव्यमान फ़ंक्शन को प्लॉट करने के लिए, यह निर्दिष्ट करते हुए कि प्लॉट एक हिस्टोग्राम है (टाइप = ‘h’)
संभाव्यता द्रव्यमान फ़ंक्शन को प्लॉट करने के लिए, बस आकार निर्दिष्ट करें (उदाहरण के लिए परीक्षणों की संख्या) और dbinom() फ़ंक्शन में संभावना (उदाहरण के लिए किसी दिए गए परीक्षण में सफलता की संभावना)।
उदाहरण के लिए, निम्नलिखित कोड दिखाता है कि आकार = 20 और संभावना = 0.3 के साथ द्विपद वितरण के लिए संभाव्यता द्रव्यमान फ़ंक्शन को कैसे प्लॉट किया जाए:
success <- 0:20 plot(success, dbinom(success, size=20, prob=.3), type='h')
x-अक्ष सफलताओं की संख्या दर्शाता है और y-अक्ष 20 परीक्षणों में सफलताओं की उतनी संख्या प्राप्त करने की संभावना दर्शाता है।
कथानक को सौंदर्य की दृष्टि से अधिक आकर्षक बनाने के लिए हम एक शीर्षक जोड़ सकते हैं, अक्ष लेबल बदल सकते हैं और पंक्तियों की चौड़ाई बढ़ा सकते हैं:
success <- 0:20 plot(success,dbinom(success,size=20,prob=.3), type='h', main='Binomial Distribution (n=20, p=0.3)', ylab='Probability', xlab ='# Successes', lwd=3)
आप चार्ट में प्रदर्शित प्रत्येक सफलता की संख्या की वास्तविक संभावनाएँ प्राप्त करने के लिए निम्नलिखित कोड का उपयोग कर सकते हैं:
#prevent R from displaying numbers in scientific notation options(scipen=999) #define range of successes success <- 0:20 #display probability of success for each number of trials dbinom(success, size=20, prob=.3) [1] 0.00079792266297612 0.00683933711122388 0.02784587252426865 [4] 0.07160367220526231 0.13042097437387065 0.17886305056987975 [7] 0.19163898275344257 0.16426198521723651 0.11439673970486122 [10] 0.06536956554563482 0.03081708090008504 0.01200665489613703 [13] 0.00385928193090119 0.00101783259716075 0.00021810698510587 [16] 0.00003738976887529 0.00000500755833151 0.00000050496386536 [19] 0.00000003606884753 0.00000000162716605 0.00000000003486784