Just when I think I'm starting to learn .... Statement z1 works, statement z doesn't. Why doesn't z work and what do I do to fix it ? Clearly the problem is with the first NA, but I would think it's handled through the loop vectorization.
y1 <- rnorm(20, 0, .013) y1 [1] -0.0068630836 -0.0101106230 -0.0169663344 -0.0066314769 0.0075063818 [6] -0.0033548024 0.0015647863 0.0119815982 -0.0021430336 0.0044617167 [11] 0.0053447708 -0.0005590323 0.0063195781 0.0073059640 -0.0181872678 [16] -0.0098094568 0.0013679040 -0.0028490887 -0.0131129191 0.0126610358 z1 <- ifelse(is.na(y1), 10000, 10000*cumprod(1+y1)) z1 [1] 9931.369 9830.957 9664.162 9600.074 9672.136 9639.688 9654.772 9770.451 [9] 9749.513 9793.012 9845.354 9839.850 9902.034 9974.378 9792.971 9696.907 [17] 9710.172 9682.506 9555.541 9676.524 y <-c(NA, rnorm(19,0, .013)) y [1] NA 0.0056258152 -0.0117690116 0.0163961630 0.0007818773 [6] 0.0007761957 0.0139769376 0.0041086982 -0.0049545337 0.0059587216 [11] -0.0079022056 0.0083076357 -0.0075823658 0.0173806814 -0.0034915869 [16] -0.0045480358 0.0168642491 0.0038681635 -0.0123010077 0.0087494624 z <-ifelse(is.na(y), 10000, 10000*cumprod(1+y)) z [1] 10000 NA NA NA NA NA NA NA NA NA NA NA [13] NA NA NA NA NA NA NA NA -- View this message in context: http://r.789695.n4.nabble.com/How-do-I-fix-this-tp3239239p3239239.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ 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.
