R语言修改矩阵元素

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

有个4*52的零矩阵A，有个矩阵b=【3,1；23,50；42,2；10,11】，维度4*2. 现想把A矩阵中第一行的第3个和第1个元素修改为1，第二行的第23和第50个元素修改为1，依次类推。。。。。也就是矩阵A中的每行，以矩阵b中每行元素为对应位置，将其修改为1. 如果用循环，能够写出来，但是程序太大，想避免循环，请教大神如何不用循环，一步实现上述过程？在此先谢谢了。

扫码加我 拉你入群

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

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

分享0 收藏3 回帖

关键词：矩阵元素 R语言 程序 如何 元素

ding yi xia tiezi

1. A=matrix(0,nrow=4,ncol=52)
   
2. b=matrix(c(3,1,23,50,42,2,10,11),nrow=4,byrow=TRUE)
   
3. ##可以实现,temp是中间结果，不用关注。
   
4. temp<-sapply(1:nrow(A),function(i) A[i,b[i,]]<<-1)
   
5. ##这种方法也可以
   
6. b=4*(b-1)+1:nrow(b)
   
7. A[t(b)]<-1

复制代码

十里春风 发表于 2015-8-8 17:42

第二种方法降一维index很巧妙，赞；第一种temp后面是不是还有其他code？只是中间结果吧

rosenbloog 发表于 2015-8-8 23:34
第二种方法降一维index很巧妙，赞；第一种temp后面是不是还有其他code？只是中间结果吧

用sapply来对A元素进行修改，修改只是中间过程，temp是sapply输出的结果，不过这个结果并不是我们需要的。

十里春风 发表于 2015-8-8 17:42

厉害，方法非常巧妙。非常感谢。不过，两个程序都没有看懂，请再细细指教一下。第一个程序，function（i) A[i,b[i,]]<<-1, 特别是这个<<-，期初我以为你打错了，应该是<-,我用了<-运行完之后居然没有任何变化，又用了<<-，果然实现了，请问这两个赋值号之间的区别是什么?如何使用？
第二个程序，这里A[t(b)]<-1为什么要给b转置，我尝试了不转置b，显示错误，下表出界。A[]这个括号内，如果是个矩阵，它是怎么引用的？为什么不转置会出界？
谢谢了，我是菜鸟一个，希望指点一下。

salesone 发表于 2015-8-10 14:06
厉害，方法非常巧妙。非常感谢。不过，两个程序都没有看懂，请再细细指教一下。第一个程序，function（i) ...

sapply计算中间结果是局部变量，<-是局部赋值，<<-可以对全局变量进行修改。
A是按矩阵index取值，b是4x2，每一行是一个位置。A[t(b)]是把A当做向量index取值。也可以用另一种，更直观，A[c(b)]，把b转换为向量，然后赋值。

十里春风 发表于 2015-8-10 16:53
sapply计算中间结果是局部变量，

非常感谢，学习了。。。

其实我在好奇sapply函数和循环会差多少

谢谢分享