आर

2015-09-11 18 views
5

में एक डबल योग लिखना मैं आर के लिए नया हूं, और दो इनपुट पी और क्यू के साथ एक डबल योग फ़ंक्शन लिखने का प्रयास कर रहा हूं। उदाहरण के लिए, वॉलफ्रेम अल्फा पर मेरे राशि की गणना करने के, मैंआर

Sum[p^25 * (1-q)^j * Sum[(25 choose r) * ((j-1) choose (r-1)) * ((q * (1-p))/(p * (1 - q)))^r,{r,1,25}],{j,0,23}] 

p = .6 के लिए Here's an example और q = .5

यह आर में लिखने के लिए आसान होना चाहिए का उपयोग करें, लेकिन मैं ऐसा करने के लिए कैसे को समझ नहीं सकता। (मुझे एहसास है कि चुनिंदा फ़ंक्शन choose(n,k) बन जाता है, लेकिन रकम अचूक साबित हो रही हैं।)

किसी भी मदद के लिए धन्यवाद।

+1

SO में आपका स्वागत है। WolframAlpha पर उदाहरण प्रदान करने के लिए +1। – jlhoward

उत्तर

3

यहाँ एक तरीका है:

p <- 0.6 
q <- 0.5 
qp <- q*(1-p)/(p*(1-q)) 
f <- function(r,j) choose(25,r)*choose(j-1,r-1)*qp^r 
g <- function(j) p^25*(1-q)^j*sum(sapply(1:25,f,j)) 
sum(sapply(1:23,g)) 
# [1] 0.721113 

आप f और g के बजाय गुमनाम कार्यों का उपयोग और एक लाइन पर यह सब डाल सकता है, लेकिन IMO यह थोड़ा कम भ्रामक है।