[R] extracting significance test for individual lm() parameters after using by
info at aghmed.fsnet.co.uk
Sun Jan 25 13:35:52 CET 2015
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
> 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
> 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
More information about the R-help