Question

我正在努力解决一些非常简单的问题。我试图找到矩阵子集的均值。例如，我想找到矩阵中每行的最后196个元素的平均值。

我的代码如下：

for (i in 1:length(ratios))  
{  
    mean(ratios[length(ratios[,1])-196:length(ratios[,1]),i])  
}

任何帮助将不胜感激！

Answer 1

这是一种更直接的方法

>rowMeans(tail(ratios, 196))

tail(ratios, 196)返回ratios每行的最后196个元素，rowMeans计算每行的平均值。

在你的例子中，你正在采取col的意思，行测量，所以尝试：

>colMeans(tail(ratios,10))

这是一个例子

> ratios <- iris[, sapply(iris, is.numeric)]
> Means <- numeric(ncol(ratios))

> for (i in 1:length(ratios))
  {
    Means[i] <- (mean(ratios[(length(ratios[,1])-10):length(ratios[,1]),i]))
 }

> COLMeans <- colMeans(tail(ratios,10))

> Means  ## Using `for` loop
[1] 6.490909 3.036364 5.336364 2.163636

> COLMeans # using `tail` and `colMeans`
Sepal.Length  Sepal.Width Petal.Length  Petal.Width 
        6.45         3.03         5.33         2.17

Answer 2

我实际上已经解决了。出于某种原因，您需要在比率矩阵的行部分周围插入括号，如下所示：

for (i in 1:length(ratios)){  
  mean(ratios[(length(ratios[,1])-196):length(ratios[,1]),i])  
}

希望这可以帮助任何有类似问题的人

找到矩阵子集的均值

2 个答案: