Displaying 1 result from an estimated 1 matches for "fit_beta".
Did you mean:
fi_beta
2005 Jul 27
1
Problem specifying "function" for "mle" operation
...re 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)
CASE 2
The error message is as follows => Error in validObject(.Object) : invalid class "mle" object: invalid object for slot "fullcoef" in class "mle": got class "list", should be or extend class "numeric"
x,s,n, and alpha ar...