Re: [Rd] Residual DF of NLS models (PR#14194)

From: Prof Brian Ripley <ripley_at_stats.ox.ac.uk>
Date: Tue, 26 Jan 2010 12:04:23 +0000 (GMT)

Can we please have a reproducible example (as we did ask in the FAQ, the posting guide ...).

On Tue, 26 Jan 2010, han_at_enfor.dk wrote:

> Full_Name: Henrik Aalborg Nielsen
> Version: 2.10
> OS: Linux (SLES 10 / openSUSE 11.1)
> Submission from: (NULL) (77.66.63.89)
>
>
> There seems to be a bug in df.residual.nls which is triggered when nls is called
> with argument na.action = na.exclude; in that case 'resid(object)' will contain
> NA-values which should be disregarded when counting the number of residuals:
>
> df.residual.nls <- function(object, ...) {
> w <- object$weights
> n <- if(!is.null(w)) sum(w != 0) else length(resid(object))
> n - length(coef(object))
> }
>
> The bug cause the F-test of anova.nls to be wrong.
>
> Replace 'length(resid(object))' with 'sum(!is.na(resid(object)))' ?
>
> ... and thank you for producing this fantastic software!
>
> BR
> Henrik
>
> ______________________________________________
> R-devel_at_r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel
>

-- 
Brian D. Ripley,                  ripley_at_stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595

______________________________________________
R-devel_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
Received on Tue 26 Jan 2010 - 12:06:07 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Tue 26 Jan 2010 - 22:30:19 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-devel. Please read the posting guide before posting to the list.

list of date sections of archive