Re: [R] for loop help

From: jim holtman <jholtman_at_gmail.com>
Date: Tue, 29 Jan 2008 17:54:07 -0500

Is this what you want?

> numSim <- 15
>
> genx<-double(numSim)
> N <- rep(20, numSim)
> A <- F <- M <- numeric(numSim)
>
> result <- lapply(1:5, function(.x){

+     for (i in 1:numSim) {
+         PN<-(runif(N[i], 0, 1))
+         A[i]<-sum(ifelse(PN>0.2, 1, 0))
+         PF<- runif((A[i]*0.5), 0, 1)
+         F[i]<-sum(ifelse(PF>0.2, 1, 0))
+         PM<- runif((A[i]*0.5), 0, 1)
+         M[i]<-sum(ifelse(PM>0.2, 1, 0))
+         genx[i]<- (min(F[i], M[i])*12)
+     }
+
+     data <- data.frame(N, A, F, M, genx)
+     N <<- genx   # assign at the top level
+     data   # return the data frame
+ })

> result

[[1]]

    N A F M genx
1 20 16 6 7 72
2 20 12 6 4 48
3 20 16 5 8 60
4 20 18 7 6 72
5 20 14 6 7 72
6 20 18 8 8 96
7 20 15 3 7 36
8 20 13 5 6 60
9 20 15 7 5 60
10 20 15 6 7 72
11 20 16 6 5 60
12 20 13 5 6 60
13 20 15 6 6 72
14 20 15 6 7 72
15 20 20 8 8 96

[[2]]

    N A F M genx
1 72 60 21 25 252
2 48 39 17 15 180
3 60 50 18 22 216
4 72 55 20 18 216
5 72 57 22 25 264
6 96 72 29 25 300
7 36 28 12 8 96
8 60 51 18 22 216
9 60 49 19 20 228
10 72 56 24 21 252
11 60 47 19 17 204
12 60 44 18 20 216
13 72 59 21 26 252
14 72 56 24 21 252
15 96 75 33 32 384

[[3]]

     N A F M genx
1 252 209 81 84 972
2 180 137 54 59 648
3 216 172 66 69 792
4 216 166 59 65 708
5 264 226 92 93 1104
6 300 246 106 93 1116
7 96 74 29 28 336
8 216 171 68 70 816
9 228 183 74 75 888
10 252 207 83 77 924
11 204 166 65 63 756
12 216 169 69 67 804
13 252 195 74 79 888
14 252 200 88 89 1056
15 384 297 118 119 1416

[[4]]

      N A F M genx
1 972 799 313 316 3756
2 648 523 213 217 2556
3 792 646 259 256 3072
4 708 575 242 231 2772
5 1104 888 372 351 4212
6 1116 913 365 353 4236
7 336 272 114 100 1200
8 816 650 250 262 3000
9 888 718 285 296 3420
10 924 741 302 291 3492
11 756 603 260 245 2940
12 804 637 249 245 2940
13 888 720 288 281 3372
14 1056 847 343 353 4116
15 1416 1111 438 446 5256

[[5]]

      N A F M genx
1 3756 3021 1229 1232 14748
2 2556 2057 830 793 9516
3 3072 2463 1009 982 11784
4 2772 2214 884 884 10608
5 4212 3385 1355 1361 16260
6 4236 3380 1329 1368 15948
7 1200 986 376 400 4512
8 3000 2407 966 960 11520
9 3420 2702 1091 1090 13080
10 3492 2778 1130 1133 13560
11 2940 2361 932 952 11184
12 2940 2289 932 903 10836
13 3372 2680 1045 1054 12540
14 4116 3279 1305 1304 15648
15 5256 4233 1713 1715 20556

>

On Jan 29, 2008 2:23 PM, Adriana Bejarano <bejaranoac_at_gmail.com> wrote:
> Hi,
>
>
>
> I have written the following code which works fine
>
>
>
> step<-5
>
> numSim<-15
>
>
>
> N<-double(numSim)
>
> A<-double(numSim)
>
> F<-double(numSim)
>
> M<-double(numSim)
>
> genx<-double(numSim)
>
>
>
> for (i in 1:numSim) {
>
>
>
> N[i]<-20
>
> PN<-(runif(N[i], 0, 1))
>
> A[i]<-sum(ifelse(PN>0.2, 1, 0))
>
>
>
> PF<- runif((A[i]*0.5), 0, 1)
>
> F[i]<-sum(ifelse(PF>0.2, 1, 0))
>
>
>
> PM<- runif((A[i]*0.5), 0, 1)
>
> M[i]<-sum(ifelse(PM>0.2, 1, 0))
>
>
>
> genx[i]<- (min(F[i], M[i])*12)
>
>
>
> }
>
>
>
> data<-data.frame(N, A, F, M, genx)
>
>
>
> What I would like to do next is to make 'genx' from step '1' the seed for N
> in step '2'- run the above code 15 times (numSim), make 'genx' from step '2'
> the seed for N in step '3'- run the above code 15 times (numSim) etc…
>
>
>
> How do I do that?. Any advice would be greatly appreciated
>
>
>
> Adriana (R-beginner)
>
> [[alternative HTML version deleted]]
>
>
> ______________________________________________
> 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.
>
>

-- 
Jim Holtman
Cincinnati, OH
+1 513 646 9390

What is the problem you are trying to solve?

______________________________________________
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.
Received on Tue 29 Jan 2008 - 22:57:23 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Tue 29 Jan 2008 - 23:30:09 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive