Функциональная форма плотности с прибл.

Question

У меня возникла проблема при попытке приблизить функциональную форму плотности в R с approxfun().Минимальный пример следующий.

Example = seq(0.5, 1.5, 0.005)
dens = density(Example)
estim = approxfun(dens) 
plot(estim)
plot(dens) 

Проблема в том, что plot(estim) дает следующий график

, но plot(dens) дает полный график,

с осью X, охватывающей все значения, и не останавливаясь на 1, как на приблизительном графике.Следовательно, approxfun() должным образом не отражает плотность.

Что я здесь не так делаю?Я попробовал несколько конфигураций approxfun, чтобы как-то включить всю ось X, но безуспешно.

Answer 1

Ваша переменная estim - это просто функция в этой точке.Он не запоминает диапазон данных, которые использовались для его создания.При построении функции вам нужно будет указать R, где ее запускать и останавливать (по умолчанию R отображает функцию от x = 0 до x = 1).Например,

plot(estim, xlim=range(Example))

Обратите внимание, что range(Example) охватывает только наблюдаемый диапазон.Если вы хотите использовать тот же диапазон, что и для графика плотности, используйте

plot(estim, xlim=range(dens$x))