Displaying 1 result from an estimated 1 matches for "fit_alpha".
Did you mean:
dst_alpha
2005 Jul 27
1
Problem specifying "function" for "mle" operation
...mle": got class "list", should be or extend class "numeric"
x,s,n, and alpha are vectors of the same length (5 dimensions in this specific case).
ll_alpha<-function(alpha=c(0.5,0.5,0.5,0.5,0.5)) -sum(s*b0+s*b1*x+s*alpha-n*log(1+exp(b0+b1*x+alpha))-0.5*(alpha/sigma)^2)
fit_alpha<-mle(ll_alpha)
MY THEORY:
I think it has something to do with the fact that the function argument is a vector in CASE 2. Hope this doesn't mean I have to re-write the function in a way that it doesn't require vector inputs.....
Any suggestions on how to make CASE 2 work would be appre...