A coding question involving variable assignments
Ted and List, What I need is I need to know what max of rnorm(rpois(1,10)) is before R does sum(), replicate(10, ...) and replicate(5, ...). The fact that you have set.seed(99) twice, does that mean, say, entry [1,1] 0.4896243 in 'mx' is one of the z number of values generated by rnorm(rpois(1,10)) that add up to [1,1] -2.0071674 in 'x'? Another way to ask the question, I guess, is, by doing set.seed(99) twice are the values generated by rnorm(rpois(1,10)) for 'x' same as those for 'mx'?
On 8/11/05, ecatchpole <e.catchpole at adfa.edu.au> wrote:
Sorry, I don't follow. What's wrong with this? Ted.
> set.seed(99) > x <- replicate(5,replicate(10,sum(rnorm(rpois(1,10))))) > set.seed(99) > mx <- replicate(5,replicate(10,max(rnorm(rpois(1,10))))) > x
[,1] [,2] [,3] [,4] [,5] [1,] -2.0071674 -7.0883335 -0.03766649 2.6859415 -5.4685172 [2,] -4.7211799 2.4177121 -0.92575948 -7.2201952 -2.5969177 [3,] 0.2325584 -0.1790635 -3.17988580 -2.0249829 1.6994276 [4,] -1.7738725 1.5836438 5.06193854 -5.4798269 -2.4363479 [5,] -3.0394562 2.8859440 2.67993750 5.0534413 0.6560877 [6,] 6.2436591 -4.0226431 1.97545757 -1.5641548 4.0443831 [7,] 0.4641453 -10.4417831 1.08048629 2.4675178 -5.5114109 [8,] 1.1570728 -3.6081361 1.37858782 0.3534015 1.8282236 [9,] 4.3988625 3.0692562 -0.69898483 -1.6882952 -1.1548913 [10,] -3.7288105 -1.4455309 5.80146323 -6.1962790 -1.3698381
> mx
[,1] [,2] [,3] [,4] [,5] [1,] 0.4896243 1.4110132 0.7329387 3.2700493 1.861891 [2,] 1.0989215 1.6215407 1.5980779 2.2921963 1.417614 [3,] 1.4000518 1.4867612 1.0130372 0.5686142 1.442630 [4,] 0.5981696 2.0916016 2.0894395 0.9760749 1.419316 [5,] 1.0066032 1.8084703 1.6556502 2.1431458 2.393037 [6,] 2.7329641 1.8689793 1.1494738 1.2899945 1.702919 [7,] 0.5851713 0.6224785 2.4466643 1.1955567 0.951106 [8,] 1.3850466 0.9305735 1.4003689 1.5209779 1.864211 [9,] 1.8645760 1.1958389 0.9270208 0.4971312 1.576020 [10,] 1.5889265 1.0253490 4.6865908 0.9852816 1.032410 On 12/08/05 12:17, xpRt.wannabe wrote,:
Ted and List, In your code that produced 'mx', you dropped sum() from my original code though. As a result, the 10 x 5 max's are of the same value. Unfortunately, that's not what I need. On 8/11/05, ecatchpole <e.catchpole at adfa.edu.au> wrote:
On 12/08/05 05:55, xpRt.wannabe wrote,:
Dear List, I have the following code that does what I want: x <- replicate(5,replicate(10,sum(rnorm(rpois(1,10))))) How might one change it such that the maximum value generated by rnorm(rpois(1,10)) can be retrieved for later use?
set.seed(99) x <- replicate(5,replicate(10,sum(rnorm(rpois(1,10))))) set.seed(99) mx <- replicate(5,replicate(10,max(rnorm(rpois(1,10))))) should work? Ted. -- Dr E.A. Catchpole Visiting Fellow Univ of New South Wales at ADFA, Canberra, Australia and University of Kent, Canterbury, England - www.ma.adfa.edu.au/~eac - fax: +61 2 6268 8786 - ph: +61 2 6268 8895
-- Dr E.A. Catchpole Visiting Fellow Univ of New South Wales at ADFA, Canberra, Australia and University of Kent, Canterbury, England - www.ma.adfa.edu.au/~eac - fax: +61 2 6268 8786 - ph: +61 2 6268 8895