Hi Greg
On 21/06/13 08:02, Greg Dropkin wrote:
hi
I'm trying to understand (a little) the code behind summary.gam, and have
the Biometrika article referred to in the Help. But am stuck early on.
In the code, starting at line 167:
if (est.disp)
rdf - residual.df
else rdf - -1
res - testStat(p, Xt, V, df[i], type = p.type,
res.df = rdf)
1) shouldn't there be some brackets in the if else statement?
- I don't think so. rdf is either the residual degrees of freedom used
to compute the scale parameter estimate or a signal of fixed scale
parameter, to be passed to testStat.
2) what is testStat?
mgcv:::testStat
- probably worth looking at the original source code which includes some
comments. It does the actual computation.
best,
Simon
thanks
Greg
--
Simon Wood, Mathematical Science, University of Bath BA2 7AY UK
+44 (0)1225 386603 http://people.bath.ac.uk/sw283 On
--
Simon Wood, Mathematical Science, University of Bath BA2 7AY UK
+44 (0)1225 386603 http://people.bath.ac.uk/sw283
__
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.
