У меня есть вектор 'gprop $ par' такой, что его длина равна 1072. Я пытаюсь вычислить сумму элементов 'gprop $ par' от 0 до значения, которое будет отличаться в каждый раз.Каждый раз элементы для суммы будут отличаться, например, в первый раз сумма будет от 0: 0, затем 0: 1 и так далее, в зависимости от числа в векторе 'tt'.
Вот код, который я сделал, чтобы иметь это:
tt = seq(1,1073)
fctG <- function(t,gamma,l0,beta){
log(gamma)-beta*log(l0+t^2)
}
regG <- function(x, data=lg) {
sqrt(10000*mean( log(sum(gprop$par[0:tt-1]))-fctG(tt,x[1],x[2],x[3]) ) ^ 2)
}
propG <- nlminb(start=rep(0.1,3), objective=regG, lower =rep(0,3), control=list(iter.max=10000))
Я получил эти сообщения об ошибках:
Warning messages:
1: In 0:tt : numerical expression has 1073 elements: only the first used
2: In 0:tt : numerical expression has 1073 elements: only the first used
3: In 0:tt : numerical expression has 1073 elements: only the first used
4: In 0:tt : numerical expression has 1073 elements: only the first used
.
.
.
50: In 0:tt : numerical expression has 1073 elements: only the first used
Я также пытался использовать код, которыйЯ нашел вопрос, который совпадает с моим в Stackoverflow:
for (rowIdx in 1:nrow(gprop$par)){
startCol <- tt[rowIdx, "b"]
gprop$par[rowIdx, "sum"] <- sum(gprop$par[rowIdx, 0:startCol])
}
Но он тоже не сработал, потому что у меня не было хороших измерений.
У кого-нибудь есть идеярешить эту ошибку?
Спасибо.