矩阵避免循环

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

问题是这样的：
BL<-diag(w/abs(lm1))
其中w,lm1均为向量。比如说六维
这样的话，做出BL为6*6矩阵。
但是实际数据中，lm1中的某些值有可能是0，当存在0的时候，希望BL[i,i]<-0
如果不用循环的话可以做吗？

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

究竟是循环还是训话？

> lm1 <- c(0, 2, 0, 5, 2, 1)
> w <- 1:6
> BL<-diag(w/abs(lm1))
> BL
     [,1] [,2] [,3] [,4] [,5] [,6]
[1,]  Inf    0    0  0.0  0.0    0
[2,]    0    1    0  0.0  0.0    0
[3,]    0    0  Inf  0.0  0.0    0
[4,]    0    0    0  0.8  0.0    0
[5,]    0    0    0  0.0  2.5    0
[6,]    0    0    0  0.0  0.0    6
> BL[BL == Inf] <- 0
> BL
     [,1] [,2] [,3] [,4] [,5] [,6]
[1,]    0    0    0  0.0  0.0    0
[2,]    0    1    0  0.0  0.0    0
[3,]    0    0    0  0.0  0.0    0
[4,]    0    0    0  0.8  0.0    0
[5,]    0    0    0  0.0  2.5    0
[6,]    0    0    0  0.0  0.0    6

qoiqpwqr 发表于 2012-12-25 13:09
> lm1  w  BL BL
     [,1] [,2] [,3] [,4] [,5] [,6]
[1,]  Inf    0    0  0.0  0.0    0

感谢版主的热心解答，可以了。

lm1 <- c(0,2,0,5,2,1)
w <- 1:6
bl <- diag(ifelse(lm1!=0,w/abs(lm1),0))

也可以这样试试

qoiqpwqr 发表于 2012-12-25 13:09
> lm1  w  BL BL
     [,1] [,2] [,3] [,4] [,5] [,6]
[1,]  Inf    0    0  0.0  0.0    0

学习了