<Bill.Venables <at> csiro.au> writes: > > Finally, I'm a bit puzzled why you use glm() when the simpler lm() would > have done the job. You are fitting a linear model and do not need the > extra paraphernaila that generalized linear models require. > > Bill Venables. >
Perhaps the original poster is confused about the difference between general (a la PROC GLM) and generalized (glm) linear models? The code is also a little puzzling because the same tests seem to be run whether p>0.05 or not. Perhaps the code will eventually be written to log-transform the data if it fails the normality test? [ hint: ?boxcox in the MASS package might be a better way to go ] Ben Bolker ______________________________________________ [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 and provide commented, minimal, self-contained, reproducible code.
