Displaying 1 result from an estimated 1 matches for "prelik".
Did you mean:
prelink
2009 Mar 23
0
Problems with adapt
...p2[i]
+thmu*(e2[i]+v[2]))
U3<- ( +tha*y$a3[i] +thb*y$b3[i] +thc*y$c3[i] +thp*y$p3[i] +thmu*(e3[i]))
Pni<- (y$Ch1[i]*exp(U1) + y$Ch2[i]*exp(U2)+
y$Ch3[i]*exp(U3))/(exp(U1)+exp(U2)+exp(U3))
fe<- exp((-e1[i]*e1[i]-e2[i]*e2[i]-e3[i]*e3[i])/(2*r))
fv<- exp((-v[1]*v[1]-v[2]*v[2])/(2*s))
prelik<- Pni*fe*fv
return(prelik)
}
int[i]<- adapt(2,lo = c(-4,-4), up = c(4,4), functn= lstarpre, eps=0.01)$value
print(int[i])
lik[i]<- log(max(1E-20,int))
print(lik[i])
} #End Loop
logl<-sum(lik)
return(-logl)
}
p<-optim(c(0,0,0,0,0,0,0,0,0,0,0),mnl.lik,y=DATA,method="BFGS"...