[R] save results in a loop

[Greg Snow]

Generally if you want to save the results of a loop then it is time to
learn to use the lapply and sapply functions instead.

Try something like:

Tukey <- lapply( 6:n, function(i) Tukey1 = HSD.test(lm(sdata_mg[,i] ~
sdata_mg$Medium+sdata_mg$color+sdata_mg$type+sdata_mg$Micro),
'sdata_mg$Micro') )

though that is still ugly with all those subscripting `$`s.  A little
nicer might be:

tmpfun <- function(vname) {
  tmp.f <- as.formula( paste( vname, '~ Medium + color + type + Micro') )
  HSD.test( lm(tmp.f, data=sdata_mg), 'Micro' )
}

Tukey <- lapply( names(sdata_mg)[6:n], tmpfun )

or

tmpfun <- function(vname) {
  vname <- as.name(vname)
  HSD.test(
    eval(substitute(
       lm(vname ~ Medium + color + type + Micro, data=sdata_mg)
    )), 'Micro')
}

Tukey <- lapply(names(sdata_mg)[6:n], tmpfun)



On Wed, Aug 13, 2014 at 1:43 PM, Wenlan Tian <tianwenlan at gmail.com> wrote:
> Hi, i'm new to R. I have a question about how to save results in a loop to
> a file. Here is an example:
>
> for (i in 6:n){
> Tukey1 = HSD.test(lm(sdata_mg[,n] ~
> sdata_mg$Medium+sdata_mg$color+sdata_mg$type+sdata_mg$Micro),
> 'sdata_mg$Micro')
> }
>
>
> I don't know how to do it with the loop for all data, so i just tried one
> of them,
> Tukey1 = HSD.test(lm(sdata_mg[,1] ~
> sdata_mg$Medium+sdata_mg$color+sdata_mg$type+sdata_mg$Micro),
> 'sdata_mg$Micro')
> Tukey1
>
> The results look like this:
>
> $statistics
>       Mean       CV   MSerror       HSD r.harmonic
>   11.87421 3.102479 0.1357148 0.5288771   7.384615
>
> $parameters
>   Df ntr StudentizedRange
>   24   4         3.901262
>
> $means
>        sdata_mg[, 6]       std  r      Min      Max
> 111d                  11.86369 0.5317421  6 11.08623 12.45651
> 125d                  11.74433 0.1663130  6 11.53504 12.02412
> 14d                   11.54073 0.3877921  8 10.80300 11.96797
> Ground                12.16673 0.3391952 12 11.56278 12.86199
>
> $comparison
> NULL
>
> $groups
>      trt    means  M
> 1 Ground 12.16673  a
> 2 111d   11.86369 ab
> 3 125d   11.74433 ab
> 4 14d    11.54073  b
>
>
> How could i get all the results in the loop?
>
>         [[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.



-- 
Gregory (Greg) L. Snow Ph.D.
538280 at gmail.com