麻烦计算方差

Question

我正在尝试手动计算R中矩阵的方差，但是遇到一些问题，因为我得到的维数不正确。

我有x，一个12x250矩阵，Temp.means，一个12x1矩阵，然后将其馈送到内部带有for外观的函数中。

我基本上希望x的方差，所以我正在使用apply函数尝试将我的方差函数应用于我的代码...但是出现了问题。

这是我的代码：

Variance.fun <- function(x,Temp.means) {
for (i in 1:250) {
     vari <- (matrix(x[i,]-Temp.means)^2)/250
 }
 }

Answer 1

如果通过“方差”表示与“标准行”（如for循环所建议）的行差异，则一种获取方法是使用apply

res <- apply( x, 1, "-", Temp.means)

那可以做成一个函数：

 row.variation.from.std <- function (mat, std){
                            apply( mat, 1, "-", std) }

我怀疑您对“方差”一词的使用会误导一些潜在的响应者，他们认为您的意思是该概念在统计实践中已定义。