You are using a random function to generate this values(runif).
So, in each replication, the value is not equal.
On Wed, Mar 3, 2010 at 4:50 PM, Val <valkr...@gmail.com> wrote:
> Thank you Henrique,
>
> However,
>
> F1<-sum(!findInterval(colMeans(replicate(100, z1(100, 4))), 0.2:0.3))
> F2<-sum(findInterval(colMeans(replicate(100, z1(100, 4))), 0.2:0.3))
>
> the sum of the two (F1+F2) = number of replicates (in this case 100).
> Sometimes I do not get that sum. Do you know why?
>
> Val
>
>
>
>
>
>
>
>
> On Wed, Mar 3, 2010 at 2:33 PM, Henrique Dallazuanna <www...@gmail.com> wrote:
>> Try this:
>>
>> sum(!findInterval(colMeans(replicate(1000, z1(100, 4))), 0.2:0.3))
>>
>> On Wed, Mar 3, 2010 at 4:15 PM, Val <valkr...@gmail.com> wrote:
>>> Hi all,
>>>
>>> Assume the following function that generate a random number.
>>>
>>> z1<-function (n, eta)
>>> {
>>> wv <- runif(n)
>>> wz <- (-1/eta) * log(wv)
>>> wz
>>> }
>>> y <- z1(100,4)
>>> mean(y)
>>>
>>>
>>> I want to run this function say 1000 times and I want to count if
>>> the mean(y) outside the following range 0.20 to 0.30.
>>>
>>> How do I do it in R?
>>>
>>> Thanks in advance
>>>
>>> ______________________________________________
>>> R-help@r-project.org mailing list
>>> https://stat.ethz.ch/mailman/listinfo/r-help
>>> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
>>> and provide commented, minimal, self-contained, reproducible code.
>>>
>>
>>
>>
>> --
>> Henrique Dallazuanna
>> Curitiba-Paraná-Brasil
>> 25° 25' 40" S 49° 16' 22" O
>>
>
--
Henrique Dallazuanna
Curitiba-Paraná-Brasil
25° 25' 40" S 49° 16' 22" O
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
