[R] Why there is no p-value from likelihood ratio test using anova in GAM model fitting?

Dieter Menne dieter.menne at menne-biomed.de
Tue Apr 28 15:08:08 CEST 2009


willow1980 <jianghua.liu <at> shef.ac.uk> writes:

> However, when I run LRT to compare
> them, the test did not return F value and p-value for me. What's the reason?
> 
> ####################################################
> Analysis of Deviance Table
> 
> Model 1: sum_surv15 ~ s(FLBS) + s(byear) + s(FLBS, byear)
> Model 2: sum_surv15 ~ s(FLBS) + SES + s(byear) + s(FLBS, byear)
>    Resid. Df Resid. Dev         Df Deviance F Pr(>F)
> 1 1202.21094     601.27                             
> 2 1201.43848     601.29    0.77246    -0.02
> ####################################################

-- Provide at least the command line you used
-- Wild guess: Try to add test="F" (which is indicated in the docs)

Dieter




More information about the R-help mailing list