所以我有两个四方矩阵,每个4x4,我试图创建一个矩阵C,其元素为(a11,b11),(a12,b12),(a13,b13),...,(a44,b44) ),共16对。我试图在R中编码。我有我的初始矩阵a_ij和b_ij,我想要矩阵C。有人可以帮我这个吗?
这是我的矩阵:
mu_ijA <- (matrix(c(seq(4, 16, by=4), seq(10, 22, by=4), seq(16, 28, by=4), seq(22, 34, by=4)), nrow= 4)/100)
a_ijA <- (4* mu_ijA)
b_ijA <- (4* (1- mu_ijA))
我希望C =((a11,b11),(a12,b12),...,(a44,b44))
感谢您的帮助!
答案 0 :(得分:1)
矩阵只能保存原子对象,因此无法保持坐标对。我相信你想要一个数组:
C <- array(c(a_ijA,b_ijA),dim=c(dim(a_ijA),2))
C[1,1,] # notice that comma at the end
[1] 0.16 3.84