我有一个物种矩阵的网站.尺寸为375 x 360.每个值代表该场地样本中物种的频率.
我试图将这个矩阵从频率转换为每个站点的相对丰度.
我已经尝试了几种方法来实现这一点,唯一有效的方法是使用for循环.然而,这需要相当长的时间或完全没有完成.
是否有功能或矢量化方法来实现这一目标?我把我的for循环作为我想要做的一个例子.
relative_abundance <- matrix(0, nrow= nrow(data_wide),
ncol=ncol(data), dimnames = dimnames(data))
i=0
j=0
for(i in 1:nrow(relative_abundance)){
for(j in 1:ncol(relative_abundance)){
species_freq <- data[i,j]
row_sum <- sum(data[i,])
relative_abundance[i,j] <- species_freq/row_sum
}
}
Run Code Online (Sandbox Code Playgroud)