我有一个向量'gprop $ par',所以他的长度是1072。 我正在尝试计算'gprop $ par'的元素总和,从0到每次都会有所不同的值。 每次求和的元素将有所不同,例如第一次,求和将基于“ tt”向量中的数字,从0:0开始,然后是0:1,依此类推。
这是我要做的代码:
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])
}
但是它也不起作用,因为我没有合适的尺寸。
任何人都有解决此错误的想法吗?
谢谢。