现有一m×n矩阵，如何对每三行三列组成的矩阵求取均值，并将该均值赋给NA

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

现有一m×n矩阵，如何对每三行三列组成的矩阵求取均值，并将该均值赋给NA

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

不清楚有没有现成的函数
如果是自己写的话，能想到的是你可以用Corbi包的submatrix，或者自己写一个简单函数提取子矩阵。然后求平均值就好。至于你说的赋值给NA，你意思是本身就知道卷积后的矩阵维度了吗？如果是的话，那你用循环每一步把计算的结果赋值就好了。如果不知道维度的话，你通过你的m,n和给定的矩阵可以计算出来。

楼主的想法，类似卷积神经网络中的部分方法，草草编个函数，请楼主和owenqi老师验证一下，我得给领导们做菜去了。

1. subsetMatMean <- function (M,#原始矩阵
   
2.                           m = 3,n = 3,#子矩阵维度
   
3.                           #jump为子矩阵移动间隔
   
4.                           jump = 1) {
   
5.   if (jump > min(m, n))
   
6.     return("jump值须小于子集维度")
   
7.   nr <- ceiling((dim(M)[1] - m) / jump + 1)
   
8.   nc <- ceiling((dim(M)[2] - n) / jump + 1)
   
9.   new <- matrix(data = NA, nr, nc)
   
10.   for (i in 1:nr) {
    
11.     for (j in 1:nc) {
    
12.       new[i, j] <-
    
13.         mean(M[((i - 1) * jump + 1):min(((i - 1) * jump + m), dim(M)[1]),
    
14.                ((j - 1) * jump + 1):min(((j - 1) * jump + n), dim(M)[2])])
    
15.     }
    
16.   }
    
17.   return(new)
    
18. }
    
19. 
    
20. 
    
21. X<-matrix(1:120,12,10)#测试矩阵
    
22. 
    
23. subsetMatMean(X, m=3, n=3, jump=3)
    
24. subsetMatMean(X, 3, 3, jump=1)
    
25. subsetMatMean(X, 5, 5, 2)
    
26. subsetMatMean(X, 5, 3, 3)

复制代码

owenqi 发表于 2022-4-8 23:53
不清楚有没有现成的函数
如果是自己写的话，能想到的是你可以用Corbi包的submatrix，或者自己写一个简单函 ...

谢谢老师，是已经知道矩阵维度的

llb_321 发表于 2022-4-9 17:18
楼主的想法，类似卷积神经网络中的部分方法，草草编个函数，请楼主和owenqi老师验证一下，我得给领导们做菜 ...

谢谢老师，我运行试试看！谢谢您的帮助！