2016-11-09 6 views
1

Я хотел бы сделать график функции квантили стандартного нормального распределения с помощью R. Я использую командуPlotting нормальный квантиль функции с R

plot(qnorm,ylim=c(-4,4)) 

и я получаю этот график

enter image description here

Однако график, который я получаю, усечен примерно в $ -2 $ и около $ 2 $. Я хотел бы увидеть что-то вроде этой картины

enter image description here

Последняя картина была сгенерирована с помощью команды

plot(1:99999/10^5,qnorm(1:99999/10^5),type="l") 

, но это неэффективно и медленно.

Почему R делает усеченный граф? Какую команду следует использовать для создания графика без усечения?

Любая помощь очень ценится!

+2

Смотрите страницу справки для 'curve' и обратить особое внимание на' n' аргумента. Для эффективного построения графика используйте последний метод, но выберите значения x, расположенные на некотором расстоянии (приблизительно) обратно пропорционально производной функции, которую вы рисуете. – whuber

ответ

1

Этого должно быть достаточно:

plot(qnorm,from=pnorm(-4),to=pnorm(4)) 

Если вы хотите его немного более гладкой, где радиус кривизны мал, я предложил бы увеличить параметр n немного, например:

plot(qnorm,from=pnorm(-4),to=pnorm(4),n=201) 

см ?plot.function и ?curve