Построение обратного ecdf с использованием ggplot2

Question

Я пытаюсь построить обратный ecdf (что = "1-F" эквивалент функции Ecdf в Hmisc) с помощью ggplot.

Вот данные:

dat = data.frame(ID = rep(c(1001:1015),each=3), 
      Month = rep(c("April", "July", "October")), 
      Value = runif(45,2.0,14.6)) 

Я мог бы построить ecdf для каждого месяца с помощью ggplot2:

ggplot(dat, aes(Value, colour = Month)) + stat_ecdf() 

Однако, мне нужно построить обратную ecdf для каждого месяца. Я пробовал несколько способов, но не мог это решить. Любая помощь ценится!

Бест,

дб

Answer 1

Лучше, если вы можете обмениваться данными, но попробовать:

series <- rnorm(1000, 200, 50) 
n<-ceiling(max(series)) 
qplot(1:n,1-(ecdf(series)(1:n)))+geom_step() 

Answer 2

Еще два способа сделать это, как я недавно работал над этим вакцинных клинических испытаний :

1. Использование Hmisc Ecdf. Это прямолинейно, и он показывает это, хотя трудно понять детали изменения разных элементов графика.

2. Рассчитать кумулятивное распределение, а затем 1-кумулятивное восстановление в обратном порядке. Постройте обратное, используя ggplot2, используя geom_step, если вам нравится функция шага на графике. Ниже функция будет использовать ecdf от базового г, чтобы дать вам кумулятивное распределение, а затем 1-накопленный:

   rcdf <- function (x) { 
   cdf <- ecdf(x) 
   y <- cdf(x) 
   xrcdf <- 1-y 
   } 
   

в приведенном выше rcdf является определенным пользователем функцией, определенной с помощью ecdf.