Я хочу получить различные значения прибл. N в зависимости от значения n.grid. Код следующий:
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)])
Я хочу 4 * ч / 3 для четных значений в n.value и 2 * ч / 3 для нечетных значений, скорректированных для каждого значения n
Извините за усложнение. спасибо