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 <tianwen...@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@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. 538...@gmail.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.