Я хотел бы построить взвешенный CDF, используя ggplot. Некоторые старые обсуждения, не относящиеся к SO (например, это от 2012) предполагают, что это невозможно, но решил сделать ре-рейз.
Например, рассмотрим эти данные:
df <- data.frame(x=sort(runif(100)), w=1:100)
Я могу показать невзвешенный CDF с
ggplot(df, aes(x)) + stat_ecdf()
Как мне взвесить это на w? Для этого примера я бы ожидал, что функция будет выглядеть как x^2, так как большие числа имеют больший вес.
В вашем ответе ошибка.
Это правильный код для вычисления взвешенного ECDF:
df <- df[order(df$x), ] # Won't change anything since it was created sorted
df$cum.pct <- with(df, cumsum(w) / sum(w))
ggplot(df, aes(x, cum.pct)) + geom_line()
ECDF представляет собой функцию F(a), равную сумме весов (вероятностей) наблюдений, где x<a делится на общую сумму весов.
Но вот более удовлетворительный вариант, который просто изменяет исходный код ggplot2 stat_ecdf: https://github.com/NicolasWoloszko/stat_ecdf_weighted
ggplot2 (которая перезапишет функцию stat_ecdf). Я думаю, что пакеты будут установлены из GitHub с помощью таких команд, как library(devtools), затем install_github("NicolasWoloszko/stat_ecdf_weighted").
- person ahorn; 06.10.2020
source("https://raw.githubusercontent.com/NicolasWoloszko/stat_ecdf_weighted/master/stat_ecdf_weighted.R"). Вы можете ввести stat_ecdf (без круглых скобок), чтобы увидеть, как определяется функция до и после вызова сценария Николаса.
- person ahorn; 06.10.2020
