Re: [R] Computational Probability

2008-12-26 Thread Rory Winston
Fantastic - thanks for all the helpful replies. Best Rory On Fri, Dec 26, 2008 at 6:51 PM, Matthias Kohl wrote: > indeed one could use package distr ... > > library(distr) > X <- Unif(Min = 0, Max = 1) > Y <- convpow(X, 10) > p(Y)(6) - p(Y)(4) > > Best > Matthias > > Prof Brian Ripley wrote: > >

Re: [R] Computational Probability

2008-12-26 Thread Matthias Kohl
indeed one could use package distr ... library(distr) X <- Unif(Min = 0, Max = 1) Y <- convpow(X, 10) p(Y)(6) - p(Y)(4) Best Matthias Prof Brian Ripley wrote: Look at packages distr* : they can do your example and might do what your real applications. On Fri, 26 Dec 2008, Rory Winston wrote:

Re: [R] Computational Probability

2008-12-26 Thread Prof Brian Ripley
Look at packages distr* : they can do your example and might do what your real applications. On Fri, 26 Dec 2008, Rory Winston wrote: Hi Firstly , happy Christmas to R-Help! Secondly, I wonder if anyone can help me with the following query: I am trying to reproduce some explicit probability c

Re: [R] Computational Probability

2008-12-26 Thread Gabor Grothendieck
Actually the last line could be simplified to just: > s > 4 & s < 6 mean sd sims [1] 0.72 0.45 2500 On Fri, Dec 26, 2008 at 8:33 AM, Gabor Grothendieck wrote: > Try a simulation approach. vignette("rv") for more info. > >> set.seed(1) >> library(rv) >> x <- rvunif(10) >> s <- simapply(x,

Re: [R] Computational Probability

2008-12-26 Thread Gabor Grothendieck
Try a simulation approach. vignette("rv") for more info. > set.seed(1) > library(rv) > x <- rvunif(10) > s <- simapply(x, sum) > mean(s > 4 & s < 6) mean sd 1% 2.5% 25% 50% 75% 97.5% 99% sims [1] 0.72 0.45 00 0 1 1 1 1 2500 On Fri, Dec 26, 2008 at 7:42 AM, Rory Winston

[R] Computational Probability

2008-12-26 Thread Rory Winston
Hi Firstly , happy Christmas to R-Help! Secondly, I wonder if anyone can help me with the following query: I am trying to reproduce some explicit probability calculations performed in APPL (a Maple extension for computational probability). For instance, in APPL, to compute the probability that the