R: Построение одного ECDF поверх другого в разных цветах

Question

У меня есть пара кумулятивных эмпирических функций плотности, которые я хотел бы строить друг над другом, чтобы проиллюстрировать различия в двух кривых. Как было отмечено in a previous question функция для рисования ECDF просто plot(Ecdf()) И как я прочитал страницу тонкой ручной, я решил, что я могу построить несколько ECDFs поверх друг друга, используя что-то вроде следующего:

require(Hmisc) 
set.seed(3) 
g <- c(rep(1, 20), rep(2, 20)) 
Ecdf(c(rnorm(20), rnorm(20)), group=g) 

Однако мои кривые иногда перекрывают немного, и может быть трудно сказать, что есть что, так же, как в приведенном выше примере, который производит этот график:

Я бы очень хотел, чтобы сделать цвет этих двух CDFs разных. Однако я не могу понять, как это сделать. Какие-нибудь советы?

Answer 1

Если память служит, я сделал это в прошлом. Насколько я помню, вам нужно было обмануть его, так как Ecdf() так проклят. Я думаю, что в help(ecdf) он намекает, что это всего лишь сюжет пошаговых функций, поэтому вы можете оценить два или более ecdfs, заговорить один, а затем аннотировать через lines().

Редактировать Оказывается, это так легко, как

R> Ecdf(c(rnorm(20), rnorm(20)), group=g, col=c('blue', 'orange')) 

как страницы справки четко заявляет col= аргумент. Но я также нашел несколько сценариев, где я использовал plot.stepfun() явно.

Answer 2

Вы можете добавить каждую кривую по одному (каждый со своим стилем), например.

Ecdf(rnorm(20), lwd = 2) 
Ecdf(rnorm(20),add = TRUE, col = 'red', lty = 1)