[R] confint() error

Prof Brian D Ripley ripley at stats.ox.ac.uk
Fri Jul 21 15:21:34 CEST 2000


On Fri, 21 Jul 2000, Tord Snall wrote:

> Dear all,
> I have run the confint() function according to below and I get the
> following error message:
> 
> > confint(stepAIC.glm.spe.var.konn2.abund, level=0.95)
> Waiting for profiling to be done...
> Error: missing value where logical needed
> In addition: Warning message: 
> NaNs produced in: sqrt((fm$deviance - OriginalDeviance)/DispersionParameter)
> 
> or
> 
> > confint(stepAIC.glm.spe.var.konn2.abund, parm="DIA", level=0.95)
> Waiting for profiling to be done...
> Error in object[[pm]] : subscript out of bounds
> 
> if only one covariate is analysed.
> 
> Profiling seem to be the problem. I would too very much like to have know
> what the warning and error messages mean. 

(Try reading the code which contains the messages.)

> All clues are welcome.

First clue: confint is in package MASS (not MASS) in the VR bundle, and you
will get better service if you give some credit.  If you do not have the
latest version, you will need to upgrade.  This was fixed some while ago,
then broken again, and AFAIK the second point at least is fixed in version
6.1-9.  The first point may well be a profiling failure due to an
inadequately fitting model.  You can u	se traceback() and debugger() to
find out what is wrong.

In any case, I can't check as you have not provided a reproducible example
so please read the section on bug reporting in the FAQ.

-- 
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 272860 (secr)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list