Построение гистограммы и оценки плотности распределения

 hist(x, breaks = "Sturges", freq = NULL, probability = !freq, include.lowest = TRUE, right = TRUE, density = NULL, angle = 45, col = NULL, border = NULL, main = paste("Histogram of" , xname), xlim = range(breaks), ylim = NULL, xlab = xname, ylab, axes = TRUE, plot = TRUE, labels = FALSE, nclass = NULL, warn.unused = TRUE, ...)

x – вектор значений, для которых строится гистограмма

breaks – границы интервалов разбиения, вычисляются автоматически или задаются вектором пользовательских значений

freq – по умолчанию равен TRUE, что соответствует расчету абсолютных частот. Если FALSE, то рассчитываются оценки вероятностей

probability – противоположно freq right – если TRUE, то строится гистограмма с закрытой правой границей

plot – по умолчанию TRUE, если FALSE, то не выводит график (используется, например, если нужно получить только частоты)

Обращение к различным значениям функции осуществляется с помощью символа “$”, например: hist(…)$counts – полученные значения абсолютных частот hist(…)$breaks – построенные границы интервалов разбиения

hist(…)mids – середины интервалов разбиения density(x) – оценка плотности распределения

plot(density(x)) – график оценки плотности

ecdf(x) – эмпирическая функция распределения

plot(ecdf(x)) – график эмпирической функции распределения

Пример:

x=rnorm(200,3,2) # генерируем выборку

hist(x,freq=FALSE,density=6,angle=45, col="blue") # строим гистограмму относительных частот. Здесь density, angle и col являются графическими параметрами, определяющими плотность и угол наклона штриховки и цвет графика.

x1=seq(-3,9,0.01) # задаем диапазон изменения х для построения графика плотности

y1=dnorm(x1,3,2) # вычисляем значения плотности при известных параметрах

lines(x1,y1,col='red') # строим график плотности

hist(x,plot='FALSE')$breaks # границы интервалов разбиения [1] -2 -1 0 1 2 3 4 5 6 7 8 60

hist(x,plot='FALSE')$counts # абсолютные частоты по интервалам