for循环计算R中每个列的返回

我将价格保存在7行（7个月）和8列（8个城市）的矩阵中。我想使用for循环来计算每列的几何平均回报率

     [,1]  [,2]  [,3]  [,4]  [,5] [,6] [,7] [,8]
[1,] 22940 25206 35206 52104 63716 5992 7228 7005
[2,] 22019 25271 35160 52549 67951 5953 7172 6869
[3,] 21743 25730 35138 53087 66165 6061 7245 6846
[4,] 20941 25549 35291 55779 66428 6319 7315 6953
[5,] 20786 25500 36221 58128 66250 6562 7196 6990
[6,] 21177 25812 36735 60738 63204 6783 7155 6968
[7,] 20684 25911 37354 63716 62389 6942 7194 6923

我的代码如下：

for (i in 1 : 8) {
policyPricesReturns[i] <- 1 + diff(log(policyPrices[,i]))[-1]
meanPolicyPricesReturns[i] <- geoMean(policyPricesReturns[i])
print(meanPolicyPricesReturns[i])
}

for（1:8中的i）{
policyPricesReturns[i]
geometric.mean（）
from
psych
在这里很有用。请尝试：

library(psych)
geometric.mean(mat)
# output
       V1        V2        V3        V4        V5        V6        V7        V8 
21456.963 25567.233 35862.503 56437.768 65129.888  6362.601  7214.829  6936.061

# data
mat <- structure(c(22940L, 22019L, 21743L, 20941L, 20786L, 21177L, 20684L, 
25206L, 25271L, 25730L, 25549L, 25500L, 25812L, 25911L, 35206L, 
35160L, 35138L, 35291L, 36221L, 36735L, 37354L, 52104L, 52549L, 
53087L, 55779L, 58128L, 60738L, 63716L, 63716L, 67951L, 66165L, 
66428L, 66250L, 63204L, 62389L, 5992L, 5953L, 6061L, 6319L, 6562L, 
6783L, 6942L, 7228L, 7172L, 7245L, 7315L, 7196L, 7155L, 7194L, 
7005L, 6869L, 6846L, 6953L, 6990L, 6968L, 6923L), .Dim = 7:8, .Dimnames = list(
    NULL, c("V1", "V2", "V3", "V4", "V5", "V6", "V7", "V8")))

图书馆（心理学）
几何平均值（mat）
#输出
V1 V2 V3 V4 V5 V6 V7 V8
21456.963 25567.233 35862.503 56437.768 65129.888  6362.601  7214.829  6936.061
#资料
mat你可能想使用@ANG的建议（如果你的矩阵只包含正数），但是如果你在寻找一个基R解，你可以直接用

(-1) ^ colSums(policyPricesReturns < 0) * exp(colMeans(log(abs(policyPricesReturns))))

这里不需要循环，因为这些函数在整个矩阵上运行。
您可以使用
colMeans
直接获得每列的平均值。非常感谢您的回答！您知道“diff”和“log”是否也是逐列计算的吗？在使用之前，我还需要先计算每列的返回值“geometric.mean”。谢谢！！谢谢你，ANG，geometric.mean有效！我已经查找到“diff”也可以应用于矩阵。不客气。你收到错误消息“object'policyPricesReturns”找不到，因为你没有在循环之前定义
policyPricesReturns
。请尝试设置
policyPricesReturns
exp(colMeans(log(policyPricesReturns)))