【发布时间】:2021-01-08 05:52:35
【问题描述】:
我想用 t_ij = a_ij / sum a_ij for i 从 1 到 n 制作一个矩阵 T。我知道矩阵 A,它是一个 100x100 的矩阵。
这是我迄今为止尝试过的:
n2<-10000
mynumbers<-c(1,0)
myprobs<-c(0.2,0.8)
trials<-sample(mynumbers, n2, replace=TRUE, prob=myprobs)
matrixA<-matrix(trials, 100, 100)
matrixA
n=100
matrixT=matrix(0,100,100)
for(i in 1:n){
for (j in 1:n){
matrixT[i,j]=matrixA[i,j]/sum(matrixA[1:n,j])
}
}
matrixT
它没有给出任何错误,但我怀疑它是否正确,尤其是因为我的输出有很多零。但这可能是合乎逻辑的,我不确定。我特别为 i 从 1 到 n 的 a_ij 求和部分而苦苦挣扎。提前谢谢!
【问题讨论】: