我正在处理矩阵中的计算，但遇到困难，非常感谢您的时间和帮助。 我正在尝试使用850行300列的数据化合物的不同阈值来估计频率。我要估计每一行的频率，因为每一行都有300个城市的洪水值。
library(matrixStats)
library(data.table)
df <- fread("res5.out",sep="\t",head=T)
#Example res5.out it has 850 rows and 300 columns, here an extract of what I am trying:
H_x1y2 H_x2y2 H_x3y3 ...
0.5 1.2 3.0
0.2 2.3 2.1
0.6 0.9 0.2
.. .. ..
num <- seq(0, 15, by = 0.1)
dfff <- t(df)
probs <- c(0.05,0.1,0.16,0.25,0.5,0.75,0.84,0.90,0.95,0.99)
quantCalc <- rowQuantiles(dfff, probs=probs)
timer=10000
## frequency
H0.3Z <- (rowSums(dfff >= 0.3)/timer)
我可以估计一个值的频率，但是如何使用“ num”序列来估计每个值0.1、0.2、0.3 ... 15以上的频率。我是否应该使用sapply或data.table来估计频率，如H0.3Z所示，以每一列到每个num序列的频率避免典型循环？
谢谢您的帮助。