आप क्या कर सकते हैं यह इस तरह मानक आर कार्यों का उपयोग कर:
c <- 1
d <- 2
a <- -2
b <- 3.5
ll <- pnorm(a, c, d)
ul <- pnorm(b, c, d)
x <- qnorm(runif(3000, ll, ul), c, d)
hist(x)
range(x)
mean(x)
sd(x)
plot(x, type='l')
वर्दी विचलन के लिए उपयोग करने के लिए सीमाओं को खोजने के लिए निमंत्रण समारोह का उपयोग किया जाता है, तब डेटा को वर्दी से उत्पन्न किया जाता है और फिर सामान्य में बदल जाता है।
यह distr पैकेज का उपयोग भी सरल है:
library(distr)
N <- Norm(c,d)
N2 <- Truncate(N, lower=a, upper=b)
plot(N2)
x <- r(N2)(3000)
hist(x)
range(x)
mean(x)
sd(x)
plot(x, type='l')
ध्यान दें कि दोनों ही मामलों में मतलब ग नहीं है और एसडी घ नहीं है। यदि आप परिणामी छंटनी वाले डेटा का मतलब और एसडी चाहते हैं तो सी और डी हो, तो आपको अलग-अलग मूल्यों (उच्च एसडी, मतलब छंटनी मूल्यों पर निर्भर करता है) के लिए अभिभावक वितरण (छंटनी से पहले) की आवश्यकता होती है, यह मान एक मूल्य होगा गणित/स्टेट सिद्धांत पाठ्यक्रम के लिए अच्छी होमवर्क समस्या। यदि आपको वही ज़रूरत है तो एक टिप्पणी जोड़ें या विशेष रूप से कहने के लिए प्रश्न संपादित करें।
आप untruncated सामान्य से डेटा उत्पन्न करने के लिए है, लेकिन केवल रेंज [क, ख] तो बस ylim तर्क प्लॉट करने के लिए उपयोग करने के भीतर डेटा प्लॉट करना चाहते हैं:
plot(rnorm(3000, c, d), ylim=c(a,b))
स्रोत
2010-08-26 14:46:22
यह वही नहीं करता है ओपी के लिए पूछा। वह सामान्य वितरण के बाद 3000 अंक चाहता है, इसलिए 'रेनॉर्म' – nico
जाने का तरीका है यह एक खराब गठित प्रश्न है ... मैंने अनुमान लगाया। यदि यह वास्तव में वाई-अक्ष पर यादृच्छिक मानों की साजिश कर रहा था, तो इसका एक हिस्सा यह बताता है कि एक्स-अक्ष सीमा अनुरोध क्या है? इसका कोई मतलब नही बनता। मेरा अनुमान है कि यादृच्छिक मानों का अनुरोध घनत्व फ़ंक्शन उत्पन्न करने के लिए किया जाता है जैसे @joFrhld और @Joris Mays। घनत्व समारोह उत्पन्न करने का यह सही तरीका है। – John