Displaying 1 result from an estimated 1 matches for "arlv21".
Did you mean:
arlv2
2013 Feb 18
2
error: Error in if (is.na(f0$objective)) { : argument is of length zero
...i1)/x[2])
down1<-((-x[1]+(j1-1)*delta1-(1-x[2])*z_i1)/x[2])
H1[i1,j1]<-pnorm(up1,mean=5,sd=1)-pnorm(down1,mean=5,sd=1)}
}
K1<-solve(I1-H1)
one1 <- matrix(1,nrow=t1,ncol=1)
u1 <- matrix(0,nrow=m1,ncol=1)
y1<-c(u1, 1, u1)
z1<-t(y1)
ARLV1<-K1%*%one1
ARLV21<-t(ARLV1)
return( abs((ARLV21%*%y1)-500)-1)
}
eval_f0 <- function(x) romberg(function(s) f(x, s), 0, 6)
of <- nloptr( x0=c(0.653,0.09),
eval_f=eval_f0,
lb = c(0,0),
ub = c(6,1),
eval_g_ineq = eval_g0,
opts...