Prof Brian Ripley <[EMAIL PROTECTED]> writes: > The short answer is no, as there is no way to recover the fitted values > and residuals so you can't get a proper fit object of class "lm" (and > hence get `summaries and all'). > > Your pseudo-data method needs to fix the u_i to be mean zero, variance > one in the sample.
That's what he did, no? A slightly more readable version is m + s*scale(rnorm(n)) > That is probably the quickest method. The elegant one > is to create a new class "groupedlm" and write a constructor etc for it > .... However, it begs the question whether it wouldn't have been better to design the RSS into the lm class rather than computing it from residuals in summary.lm and anova.lm and predict.lm and... -- O__ ---- Peter Dalgaard Blegdamsvej 3 c/ /'_ --- Dept. of Biostatistics 2200 Cph. N (*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918 ~~~~~~~~~~ - ([EMAIL PROTECTED]) FAX: (+45) 35327907 ______________________________________________ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
