不同行的频率估算

时间:2018-10-28 19:56:02

标签: r dataframe sapply frequency-analysis

我正在处理矩阵中的计算,但遇到困难,非常感谢您的时间和帮助。 我正在尝试使用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序列的频率避免典型循环?

谢谢您的帮助。

0 个答案:

没有答案