两个问题：
1.你不展示数据。
1.不显示输出
我想在控制台对话框的这个小片段中可以看到可能发生的情况：

> 1+numeric()
 numeric(0)

所以结果实际上不是零，而是一个空的（长度= 0）数值向量，对空加一些东西是可以论证的未定义的，所以人们可能期望NA作为结果，但在这个例子中肯定不是1。
您可以尝试将mse设置为0，然后重新运行代码。

我不回答你的问题，因为我想上面有人做了。
但是考虑到你的代码，这里有一个建议。你应该在循环外创建向量，矩阵等，然后在循环内填充它们。这是为了避免在启动一些很长的循环时计算时间太长。
请看下面我的意思作为一个简单的例子。

#your method
    x <- 0
    system.time(
    for(i in 2:10000)
    {
    x <- c(x,x[i-1]+i)
    }
    )
utilisateur     système      écoulé 
       0.10        0.07        0.16 

    #better 
    z <- rep(NA,10000)
    z[1] <- 0
    system.time(
    for(i in 2:10000)
    {
    z[i] <- z[i-1]+i
    }
    )
utilisateur     système      écoulé 
          0           0           0 
    all(x == z)
    TRUE

您的方法需要一点时间，而它不应该。