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.