I am trying to compare two different GAM fits.
I have something like
Course.bam20 <-bam(zGrade ~ Rep + ISE + White + Female + Years + AP_TOTAL
+ MATH + HSGPA+ EOP + factor(P7APrior, ordered = FALSE)+s(Yfrm7A,k=20),
data= Course, na.action = na.exclude,samfrac =0.1)
Course.bam4 <-bam(zGrade ~ Rep + ISE + White + Female + Years + AP_TOTAL
+ MATH + HSGPA+ EOP + factor(P7APrior, ordered = FALSE)+s(Yfrm7A,k=4),
data= Course, na.action = na.exclude,samfrac =0.1)
anova(Course.bam20, Course.bam4)
Model 1: zGrade ~ Rep + ISE + White + Female + Years + AP_TOTAL + MATH +
HSGPA + EOP + factor(P7APrior, ordered = FALSE) + s(Yfrm7A,
k = 20)
Model 2: zGrade ~ Rep + ISE + White + Female + Years + AP_TOTAL + MATH +
HSGPA + EOP + factor(P7APrior, ordered = FALSE) + s(Yfrm7A,
k = 4)
Resid. Df Resid. Dev Df Deviance
1 4721.7 1907.0
2 4724.5 1913.5 -2.7919 -6.4986
How can I get a p-value out of the anova?
[[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.
