From: Narcyz Ghinea <Narcyz.Ghinea_at_swsahs.nsw.gov.au>

Date: Wed 27 Jul 2005 - 13:01:26 EST

R-help@stat.math.ethz.ch mailing list

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Wed Jul 27 13:06:13 2005

Date: Wed 27 Jul 2005 - 13:01:26 EST

Hello fellow R users,

Below are two cases using the "mle" operation from the stats4 package. In CASE 1 the code runs fine, in CASE 2 errors occur:

**CASE 1
**
x, alpha_current, s, and n are vectors of the same length.

ll_beta<-function(b0=0,b1=0) -sum(s*b0+s*b1*x+s*alpha_current-n*log(1+exp(b0+b1*x+alpha_current)))

fit_beta<-mle(ll_beta)

x,s,n, and alpha are vectors of the same length (5 dimensions in this specific case).

fit_alpha<-mle(ll_alpha)

Any suggestions on how to make CASE 2 work would be appreciated.

peace

Narcyz.

This message is intended for the addressee named and may con...{{dropped}}

R-help@stat.math.ethz.ch mailing list

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Wed Jul 27 13:06:13 2005

*
This archive was generated by hypermail 2.1.8
: Fri 03 Mar 2006 - 03:34:02 EST
*