Displaying 1 result from an estimated 1 matches for "alfz2".
Did you mean:
alf2
2009 Mar 23
0
Problems with adapt
...----------------------------------------
DATA<-read.table("Data.csv",header=TRUE,sep=",")
library(adapt)
mnl.lik<-function(theta,y){
th1<-theta[1]
th2<-theta[2]
tha<-theta[3]
thb<-theta[4]
thc<-theta[5]
thp<-theta[6]
thmu<-theta[7]
alfz1<-theta[8]
alfz2<-theta[9]
alfc<-theta[10]
alf<-theta[11]
r<- 1
s<- 1
n<-2000
lik<-numeric(n)
int<-numeric(n)
v<- numeric(2)
for (i in 1:n) { #Beggin Loop
lstarpre<- function(v){ #This is the fuction to be integrated
e1<-y$p1[i]-alfz1*y$z_a1[i] - alfz2*y$z_b1[i] -a...