3482
правки
Изменения
м
==== Пакет <code>gamlss</code> ====
Ещё один интересный способ графического анализа представлен функцией <code>histDist</code> из пакета <code>gamlss</code>:
{{r-code|code=
<nowiki>> histDist(x, family = "NO", density = TRUE)
Family: c("NO", "Normal")
Fitting method: "nlminb"
Call: gamlssML(y = y, family = "NO", formula = x)
Mu Coefficients:
[1] -0.0462
Sigma Coefficients:
[1] 0.023
Degrees of Freedom for the fit: 2 Residual Deg. of Freedom 998
Global Deviance: 2884
AIC: 2888
SBC: 2898
</nowiki>
}}
[[Файл:Gamlss-histdist.svg|400px|центр]]
С помощью аргумента <code>family</code> можно задать семейство распределений для подгонки и сравнения<ref>Более подробную информацию о доступных семействах распределений можно получить с помощью команды <code>help("gamlss.family")</code>.</ref>.
→Графические методы
[[Файл:Car-densityPlot.svg|400px|центр]]
==== Гистограммы с наложением графика плотностей вероятнотси ====
===== Пакет <code>stats</code> =====
{{r-code|code=
<nowiki>> hist(x, freq = FALSE)
> lines(density(x))</nowiki>}}
[[Файл:Stats-hist-density.svg|400px|центр]]
Теперь наложим на наш график кривую плотностей вероятности для нормального распределения:
{{r-code|code=
<nowiki>> xfit <- seq(min(x), max(x), length = 100)
> yfit <- dnorm(xfit, mean = mean(x), sd = sd(x))
> hist(x, freq = FALSE)
> lines(density(x), col = "red")
> lines(xfit, yfit, col = "blue")
</nowiki>}}
[[Файл:Stats-density-compare.svg|400px|центр]]
===== Пакет <code>gamlss</code> =====
Более простой способ сравнение графиков плотностей вероятности представлен в функции <code>histDist</code> из пакета <code>gamlss</code>:
{{r-code|code=
<nowiki>> histDist(x, family = "NO", density = TRUE)
Family: c("NO", "Normal")
Fitting method: "nlminb"
Call: gamlssML(y = y, family = "NO", formula = x)
Mu Coefficients:
[1] -0.0462
Sigma Coefficients:
[1] 0.023
Degrees of Freedom for the fit: 2 Residual Deg. of Freedom 998
Global Deviance: 2884
AIC: 2888
SBC: 2898
</nowiki>
}}
[[Файл:Gamlss-histdist.svg|400px|центр]]
С помощью аргумента <code>family</code> можно задать семейство распределений для подгонки и сравнения<ref>Более подробную информацию о доступных семействах распределений можно получить с помощью команды <code>help("gamlss.family")</code>.</ref>.
==== Q-Q график ====
[[Файл:E1071-probplot.svg|400px|центр]]
== Многомерное нормальное распределение ==
