Message-ID: <49999922.7000901@pburns.seanet.com>
Date: 2009-02-16T16:49:38Z
From: Patrick Burns
Subject: Alternate to for-loop
In-Reply-To: <49997751.4000902@statistik.tu-dortmund.de>
If the goal is to "look" professional, then
'replicate' probably suits. If the goal is to
compute as fast as possible, then that isn't
the case because 'replicate' is really a 'for'
loop in disguise and there are other ways.
Here's one other way:
function (size, replicates, distfun, ...)
{
colMeans(array(distfun(size * replicates, ...), c(size,
replicates)))
}
Patrick Burns
patrick at burns-stat.com
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of "The R Inferno" and "A Guide for the Unwilling S User")
Uwe Ligges wrote:
>
>
> megh wrote:
>> No, it is not homework. I obviously
>
>
> For some value of "obvious" as you has not given a single line of code
> as the posting guide suggests.
>
> You probably want:
>
> replicate(10, mean(rnorm(100)))
>
> Uwe Ligges
>
>
>> could do that using a for-loop, and that
>> I already did. However I thought whether there could be a better
>> approach as
>> it was looking very messy and unprofessional.
>>
>>
>>
>> Uwe Ligges-3 wrote:
>>>
>>>
>>> megh wrote:
>>>> Hi, I am trying to create a vector of length 10 (say), wherein each
>>>> element
>>>> will be average of random sample of size 100, from a distribution, say
>>>> Normal. Can anyone please tell me without creating a "for" loop, how I
>>>> can
>>>> do that?
>>>
>>> Homework? Then please ask you course material or teacher.
>>>
>>> PLEASE do read the posting guide
>>> http://www.R-project.org/posting-guide.html
>>> and provide commented, minimal, self-contained, reproducible code.
>>>
>>> Uwe Ligges
>>>
>>>> Regards,
>>>>
>>>>
>>> ______________________________________________
>>> R-help at 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.
>>>
>>>
>>
>
> ______________________________________________
> R-help at 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.
>
>