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 <[EMAIL PROTECTED]> 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 <[EMAIL PROTECTED]> 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 > ______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
