我想根据n.grid的值获得不同的大约n值 代码如下：
ll <- 0
ul <- 0.5
mean <- 0
sdsq <- 1
n.value <- 100
n.value <- c(1:n.value)
xi <- lapply(n.value, function(i) seq(ll,ul,length.out = i+1))
xi <- xi[-1]
h <- (ul-ll)/n.value
fx <- sapply(xi, dnorm, mean= mean,sd=sqrt(sdsq))
indv.exact <- pnorm(ul, mean= mean,sd=sqrt(sdsq)) - pnorm(ll, mean= mean,sd=sqrt(sdsq))
approx.n <- ( h/3 ) * ( fx [1] + fx[n.value] ) + ( 4*h/3 ) * sum(fx[seq(2,n.value-2,2)]) + ( 2*h/3 ) * sum(fx[seq(3,n.value -1,2)])
搜索文件会出现错误
Error in fx[1] + fx[n.grid] : non-numeric argument to binary operator
此后也要注意以下内容
( 4*h/3 ) * sum(fx[seq(2,n.value-2,2)]) + ( 2*h/3 ) * sum(fx[seq(3,n.value -1,2)])
我希望n.value中的偶数为4 * h / 3，为n的每个值调整后的奇数为2 * h / 3
很抱歉给您带来的麻烦。谢谢