[R] invalid function value in 'nlm' optimizer

Hi, R users, I am using nlm function to get the MLE of parameter alpha and lambda from a parametric survival model (Weibull distribution). However, this message always came out: ' invalid function value in 'nlm' optimizer'. Could anyone help me? Code is project<-read.table(file="C://data.txt",

Re: [R] invalid function value in 'nlm' optimizer

function(theta,x1,x2 ) { -sum(delta)*log(theta[1]*theta[2])-sum(delta)*(theta[1]-1)*log(x[delta==1])+theta[2]*sum(x^theta[1]) } #MLE nlm(ln,theta<-c(1,1),x1=x, x2=delta, hessian=TRUE) On Tue, Apr 28, 2015 at 6:40 AM, Duncan Murdoch wrote: > On 28/04/2015 2:43 AM, Hanze Zhang wrote: &g

Re: [R] invalid function value in 'nlm' optimizer

t; > > Bill Dunlap > TIBCO Software > wdunlap tibco.com > > On Tue, Apr 28, 2015 at 6:40 PM, Hanze Zhang > wrote: > >> I still cannot solve the problem: 'invalid function value in 'nlm' >> optimizer' >> >> I want to get the MLE for t

[R] R2WINBUGS Error message

Hi, guys, I am a new user for package R2winbugs. When I run the code a=bugs(...), an error message always comes out, see below: Error in file(con, "wb") : cannot open the connection In addition: Warning messages: 1: In file.create(to[okay]) : cannot create file 'c:/Program Files/WinBUGS14//Syst

[R] how to do poisson regression for a complex sample design data by svyglm

Hi, I am in trouble on doing poisson regression for a complex sample design data by svyglm (survey package). I have a dataset contained these variables: caseumber, y(which is count), x1, x2, stratum, weight. It is a complex sample design with equal probability without replacement. How can I fit a

[R] problem on package "survey" , function svyglm,

Hi, I want to do logistic regression based on a complex sample design. I used package survey, but when I ran svyglm, an error message came out: Error in onestrat(x[index, , drop = FALSE], clusters[index], nPSU[index][1], : Stratum (16) has only one PSU at stage 1 My code is below: a.design<-

Re: [R] problem on package "survey" , function svyglm,

Thank you, guys. I merged strata and fitted my model. On Wed, Apr 16, 2014 at 5:33 PM, Thomas Lumley wrote: > > On Tue, Apr 15, 2014 at 7:07 AM, Milan Bouchet-Valat wrote: > >> Le lundi 14 avril 2014 à 13:59 -0400, Hanze Zhang a écrit : >> > Hi, >> > >