Displaying 1 result from an estimated 1 matches for "z_b3".
Did you mean:
z8b3
2009 Mar 23
0
Problems with adapt
...umeric(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] -alfc*y$c1[i] -alf
e2<-y$p2[i]-alfz1*y$z_a2[i] - alfz2*y$z_b2[i] -alfc*y$c2[i] -alf
e3<-y$p3[i]-alfz1*y$z_a3[i] - alfz2*y$z_b3[i] -alfc*y$c3[i] -alf
U1<- (th1 +tha*y$a1[i] +thb*y$b1[i] +thc*y$c1[i] +thp*y$p1[i]
+thmu*(e1[i]+v[1]))
U2<- (th2 +tha*y$a2[i] +thb*y$b2[i] +thc*y$c2[i] +thp*y$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...