矩阵运算求助

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

adjacency = matrix(1,nrow = 3600,ncol = 3600)
Sdatexpr = as.numeric(Sdatexpr)
sc = as.matrix(as.numeric(Sdatexpr))
for(i in 1:3600){
  for(j in 1:3600){
    sd = as.matrix(as.numeric(round(sc[i,j],2)))
    adjacency[i,j] = 1/(1+ml_fun(0.7,1,-0.3*abs(sd)^0.7,3,0.1))
  }
}

Sdatexpr是一个3600*3600的实对称矩阵。我想把Sdatexpr中m行n列处的元素通过运算之后放到新矩阵(adjacency矩阵)中的m行n列处。ml_fun函数如下：

library(expm)
ml_fun <- function(a,b,x,n,eps0){
  if (nargs() < 5) eps0 <- .Machine$double.eps
  if (nargs() < 4) n <- 0
  am = matrix(1,nrow = 1,ncol = 1)
  f <- 0 ; fa <- am ; j <-0
  while(norm(fa,"1")>=eps0){
    # fa <- gamma(j+n+1)%/%gamma(j+1)%/%gamma(a%*%j+a%*%n+b)%*%x%^%j
    fa <- gamma(j+n+1)/gamma(j+1)/gamma(a*j+a*n+b)*x%^%j
    f <- f+fa ; j<-j+1
  }
}
运行之后提示：
Error in adjacency[i, j] <- 1/(1 + ml_fun(0.7, 1, -0.3 * abs(sd)^0.7,  :
  replacement has length zero
求同志们帮助，是什么原因导致的这个错误呢。

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

顶顶，求助。。

顶顶~~

针对矩阵元素的运算，大多数不需要循环吧，直接对矩阵用做运算函数，得到对应的值矩阵

你得把计算公式和问题贴出来，R使用的是向量运算，感觉代码不需要这么复杂。

zhumengjin 发表于 2017-12-5 00:02
针对矩阵元素的运算，大多数不需要循环吧，直接对矩阵用做运算函数，得到对应的值矩阵

感谢！

zhumengjin 发表于 2017-12-5 00:02
针对矩阵元素的运算，大多数不需要循环吧，直接对矩阵用做运算函数，得到对应的值矩阵

感谢！

zhumengjin 发表于 2017-12-5 00:02
针对矩阵元素的运算，大多数不需要循环吧，直接对矩阵用做运算函数，得到对应的值矩阵

感谢！

zhumengjin 发表于 2017-12-5 00:02
针对矩阵元素的运算，大多数不需要循环吧，直接对矩阵用做运算函数，得到对应的值矩阵

感谢！