GAM p-values

[kate.bassil@utoronto.ca]

I am using S-plus 7.0 to run GAM models in my analysis.
An example of the model is:
model<-gam(Def1~T.Max+day+lo(RH)+lo(time)+offset(log(All.Calls)),family=poisson(link=log),na.action=na.omit).
To get the coefficients I enter:
coefficients(summary.glm(model,c=F,dispersion=0).

However, this only produces the value, std error, and t-value. I would  
also like p-values for the GAM models. Any suggestions on how to do  
this?
Thanks,
Kate Bassil