Re: [Rd] Windows crash in confint() with nls fit (PR#8428)

From: <murdoch_at_stats.uwo.ca>
Date: Thu 22 Dec 2005 - 19:19:27 GMT


On 12/22/2005 1:23 PM, bolker@zoo.ufl.edu wrote:
> Full_Name: Ben Bolker
> Version: 2.2.1
> OS: Windows XP and 2000
> Submission from: (NULL) (128.227.60.124)
>
>
> The following code, using confint() to try
> to get confidence intervals on an nls object
> that has been fitted with algorithm="port"
> reliably crashes R 2.2.0 and 2.2.1 with the
> latest version of MASS on a Windows 2000 and
> a Windows XP machine here. I *think* earlier
> versions of MASS may not have crashed, but
> I have now updated the versions on most of
> my machines and don't have a test version handy.
>
> On XP, I get an error-report-do-you-want-
> to-send-this-to-Microsoft dialog box. On 2000
> I get a "Rgui has generated errors ..." dialog box.
>
> I can get this far:
>
> debug: profiledModel <- .Call("nls_iter", fittedModel, ctrl, trace,
> PACKAGE = "stats")
>
> Browse[1]> where
> where 1: prof$getProfile()
> where 2: profile.nls(object, which = parm, alphamax = (1 - level)/4)
> where 3: profile(object, which = parm, alphamax = (1 - level)/4)
> where 4: confint.nls(n1)
> where 5: confint(n1)
>
> I'm not set up to debug compiled code on Windows, and I haven't
> been able to reproduce the problem on Linux.

I see it in R-devel too. It's in a strcmp; I'll see if I can track down the cause.

Duncan Murdoch

>
> set.seed(1001)
> x = runif(200)
> a =1
> b = 1
> c = -0.1
> y = a+b*x+c*x^2+rnorm(200,sd=0.05)
> plot(x,y)
> curve(a+b*x+c*x^2,add=TRUE)
> n1 = nls(y~a+b*x+c*I(x^2),start=c(a=1,b=1,c=0.1),algorithm="port")
> confint(n1) ## boom
>
> N.B.: It only crashes if algorithm="port" is specified.
>
> ## [1] "Windows 2000 Professional (build 2195) Service Pack 4.0"
> ##
> ## > R.version
> ## _
> ## platform i386-pc-mingw32
> ## arch i386
> ## os mingw32
> ## system i386, mingw32
> ## status
> ## major 2
> ## minor 2.0
> ## year 2005
> ## month 10
> ## day 06
> ## svn rev 35749
> ## language R
> ##
> ## MASS
> ##
> ## Version: 7.2-23
> ## Date: 2005-12-08
> ##
> ##
>
> ## "Windows XP Professional (build 2600) Service Pack 2.0"
> ## _
> ##platform i386-pc-mingw32
> ##arch i386
> ##os mingw32
> ##system i386, mingw32
> ##status
> ##major 2
> ##minor 2.0
> ##year 2005
> ##month 10
> ##day 06
> ##svn rev 35749
> ##language R
>
> ##
> ##Version: 7.2-23
> ##Date: 2005-12-08 ## CRASH
>
> ##> R.version
> ## _
> ##platform i386-pc-mingw32
> ##arch i386
> ##os mingw32
> ##system i386, mingw32
> ##status
> ##major 2
> ##minor 2.1
> ##year 2005
> ##month 12
> ##day 20
> ##svn rev 36812
> ##language R
> ## CRASH
>
> ______________________________________________
> R-devel@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel



R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel Received on Fri Dec 23 08:10:05 2005

This archive was generated by hypermail 2.1.8 : Fri 23 Dec 2005 - 02:22:20 GMT