У меня есть игрушечный образец только из четырех чисел, и я хочу построить эмпирическую функцию распределения и функцию квантиля (обратное к edf). Код выглядит следующим образом:
x <- c(-1.55, -0.67, -0.39, 0.60)
test = sort(x)
#Empirical cdf
plotout = plot(stepfun(sort(test), (0:4)/4), main="Empirical cdf")
segments(plotout$t[1], 0.4, plotout$t[3], 0.4, col="red")
segments(plotout$t[3], -0.04, plotout$t[3], 0.4, col="red")
#Empirical quantile function
plotout = plot(stepfun((1:3)/4, sort(test)), main="Empirical quantile function")
segments(0.4, 1.04*plotout$y[1], 0.4, plotout$y[2], col="red")
segments(-0.04, plotout$y[2], 0.4, plotout$y[2], col="red")
Я понял, что аргументы для построения этой функции 'stepfun' сильно отличались от графика (x, y) точек данных, например, вы не можетеиспользуйте type='l' здесь. Итак, какие аргументы доступны тогда? Если говорить точнее, то для функции квантиля было поставлено только 3 точки, как показано на рисунке результата
Как я могу поставить точку первой пропущенной точке? А что за вышеприведенные 'plotout $ t', 'plotout $ y' в коде?
Заранее большое спасибо!