有些复杂的抽样+循环，求问

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

x <- rnorm(40000)
W <- matrix(x,200,200)
P <- sum(W*W)
S <- apply(W,2,sum)
h<-S/P
H=W[,rank(h) > length(h) - 10]#抽出来10列
Q=svd(H)
U=Q$u#提取左奇异向量
T=W-U %*% t(U)%*%W
N = ncol(W)
j = 1:N
s=ncol(H)
i=1:s
1.对矩阵W更新被选择的概率:
如果j∈[1:N]，重复从Q到T步骤，如果j包含在s中（第j列已被抽出去），令被选择的概率P_j=0,否则，P_j等于T_j的谱范数（P_j=sum(E*E)）。
2.按照更新后的被选择概率再从W中抽取10列，并到H中，整个过程重复5次，共抽取50列。
对于这两步，抽样+循环，有没有大佬可以帮忙指导一下代码呀？

扫码加我 拉你入群

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

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

分享0 收藏0 回帖