Я занимаюсь вычислениями в матрице, но застрял, и я очень ценю ваше время и помощь.Я пытаюсь оценить частоты с использованием различных порогов для соединения данных по 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, чтобы избежать типичного зацикливания?
Спасибо за помощь.