Генерация pdf для chisquared в R - PullRequest
0 голосов
/ 26 мая 2020

Я хочу визуализировать pdf-файл распределения chisquare в R. Я согласен с графической частью, но трудно получить фреймы данных. Я хочу PDF-файл из 100000 случайных чисел со степенями свободы 1,2,3,6.

Это то, что я пробовал.

set.seed(1)
data.frame(chisq = 0:100000) 
df1 = dchisq(x = chisq, df = 1)
df2 = dchisq(x = chisq, df = 2)
df3 = dchisq(x = chisq, df = 3)
df6 = dchisq(x = chisq, df = 6)

, а затем я получаю сообщение об ошибке:

Ошибка в dchisq (x = chisq, df = 1) : Non-numeri c аргумент математической функции

Traceback: 1. dchisq (x = chisq, df = 1)

1 Ответ

0 голосов
/ 26 мая 2020

Если вам нужно 4 разных фрейма данных, вы можете избежать использования первых двух строк кода и использовать этот код

x = 0:99999   # this way the vector has length 100000 instead of 100001
df1 = data.frame(x = x, dchisq = dchisq(x, df = 1))
df2 = data.frame(x = x, dchisq = dchisq(x, df = 2))
df3 = data.frame(x = x, dchisq = dchisq(x, df = 3))
df6 = data.frame(x = x, dchisq = dchisq(x, df = 6))

EDIT

Код после обсуждения в комментариях

ggplot(df1) +
  aes(x = x, y = dchisq) +
  geom_area(fill = "blue")
...