[R] extracting significance test for individual lm() parameters after using by

Michael Dewey info at aghmed.fsnet.co.uk
Sun Jan 25 13:35:52 CET 2015


Dear Jerad

I may have completely misunderstood your question but you do know that 
you can write your own function and use it in sapply where you have 
summary? You could incorporate calls to summary or to coef or somet 
other extractor or you could use the $ tool.

On 25/01/2015 02:01, Moxley, Jerad wrote:
> I’m trying to test what growth functions best fit individual subjects. I’m wanting compare linear, quadratic, cubic etc.  Here is the example from the cubic curve.
>
> b3a<-by(c,id,function(x) lm(w~agec+ageq+agecub,data=x))
>
> I can get quiet a bit of information out of sapply(b3a,summary) but it reports something like this for each person
>                37
> call          Expression
> terms         Expression
> residuals     Numeric,62
> coefficients  Numeric,16
> aliased           Logical,4  4
> sigma         67.05895
> df            Integer,3
> r.squared     0.9822921
> adj.r.squared 0.9813762
> fstatistic    Numeric,3
> cov.unscaled  Numeric,16
>
> I could obviously compute by hand the r squared change and then compute a p value based for what the partial r is for the variable but I’d like to have a simpler solution to it. Can I get the information or is there an package that I should download that will do what I’m trying to do much simpler?
>
> 	[[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> 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.
>
> -----
> No virus found in this message.
> Checked by AVG - www.avg.com
> Version: 2015.0.5645 / Virus Database: 4273/8997 - Release Date: 01/25/15
>

-- 
Michael
http://www.dewey.myzen.co.uk



More information about the R-help mailing list