Displaying 2 results from an estimated 2 matches for "arlv2".
Did you mean:
arlv
2013 Feb 18
2
error: Error in if (is.na(f0$objective)) { : argument is of length zero
...*delta-(1-x[2])*z_i)/x[2]))-s
down<-(((-x[1]+(j-1)*delta-(1-x[2])*z_i)/x[2]))-s
H[i,j]<-pnorm(up,mean=0,sd=1)-pnorm(down,mean=0,sd=1)}
}
K<-solve(I-H)
one <- matrix(1,nrow=t,ncol=1)
u <- matrix(0,nrow=m,ncol=1)
y<-c(u, 1, u)
z<-t(y)
ARLV<-K%*%one
ARLV2<-t(ARLV)
return((ARLV2%*%y)/3)
}
# constraint function
eval_g0 <- function(x) {m1<-100
t1<-2*m1+1
H1<-matrix(data=NA,nrow=t1,ncol=t1)
I1<-diag(t1)
delta1<-2*x[1]/t1
z1<--x[1]+.5*delta1
range11<-t1
for (i1 in 1:range11){
z_i1=-x[1]+(i1-.5)*delta1
for (j1 in 1:range11){...
2013 Feb 27
0
A program running for a too long time
...delta-(1-x[2])*z_i)/x[2]))-y
down<-(((-x[1]+(j-1)*delta-(1-x[2])*z_i)/x[2]))-y
H[i,j]<-pnorm(up,mean=0,sd=1)-pnorm(down,mean=0,sd=1)}
}
K<-solve(I-H)
one <- matrix(1,nrow=t,ncol=1)
u <- matrix(0,nrow=m,ncol=1)
y<-c(u, 1, u)
z<-t(y)
ARLV<-K%*%one
ARLV2<-t(ARLV)
return(ARLV2%*%y)
}
# constraint function
eval_g0 <- function(x) abs(f(x,0)-500)-1
#objective function
eval_f0 <- function(x) romberg(function(y) ((1+y*y)*f(x,y)/6),0,6, tol =
.Machine$double.eps^(2/3))$value
obj <- nloptr( x0=c(0.653,0.09),
eval_f=eval_f...