R索引数字（0）而不是0

Question

这就是全部，但它仍然不起作用。我猜问题是第二个 Amat 条件和 if else 结构：

Dmat <- matrix(0,nrow= J, ncol=J)
diag(Dmat) <- 1
dvec <- rep(0,J)
dvec <- -hsmooth
Aeq <- matrix(nrow= J, ncol=0)
beq <- rep(0)
Amat <- matrix(0,J,2*J-3)
bvec <- vector(mode="numeric", length=2*J-3)

for(j in 1:J)
{
Amat[j-1,j-1] = -1
Amat[j,j-1]   = 1
bvec[j-1]     = Delta1
}

for(j in 2:J){
Amat[j,J-1+j-2]   = -1/ (Q[j] - Q[j-1])  
if (j==2)
{Amat[j-1,J-1+j-2] = 1/(Q[j] - Q[j-1]) + 1/(Q[j-1] - 0)
}
else
{Amat[j-1,J-1+j-2] = 1/(Q[j] - Q[j-1]) + 1/(Q[j-1] - Q[j-2])
}
Amat[j-2,J-1+j-2] = -1/ (Q[j-1] - Q[j-2])
bvec[J-1+j-2]     = Delta2
}

solution1 <- solve.QP(Dmat, dvec, Amat, bvec, meq=0, factorized=FALSE)

数学形式如下：

min sum(mj- m^)，其中 j 从 1 到 J。

圣。

mj-1

1/(Qj-1 -Qj-2)(mj-2 -mj-1)

如果能在这里得到帮助，我将不胜感激。

Answer 1

R 不是C。 R 中的索引以 1 开头，而不是从零开始。所以Q 中的第一个元素是Q[1]。这与许多其他 (C-inspired) 语言不同，但更像是通常的数学符号。

所以您可能想用Q[j-1] 替换Q[j-2]。

编辑：

自从我提出答案以来，问题已经发生了很大变化。问题依然存在：

for(j in 1:J)
{
    Amat[j-1,j-1] = -1
    

当索引应该从 1 开始时，Wenn j 是 1，然后 j-1 为零，Amat[0,0] 不是一个有用的表达式。