seminar CI

### Доверительные интервалы ###

# загружаем данные 
dat <- read.csv("https://allatambov.github.io/psms/chip-n-dale.csv")

# смотрим на их описание
str(dat)
summary(dat)

# удаляем пропущенные значения
dat <- na.omit(dat)

# строим гистограммы
hist(dat$cheer, 
     col = "red", 
     border = "yellow",
     main = "Уровень веселости",
     xlab = "Веселость (в % от Дейла)",
     ylab = "Частота")

hist(dat$grump, 
     col = "burlywood1", 
     border = "firebrick",
     main = "Уровень ворчливости",
     xlab = "Ворчливость (в % от Чипа)",
     ylab = "Частота")

# строим ящики с усами
boxplot(dat$cheeze, 
        col = "yellow", 
        border = "sienna2")

boxplot(dat$tech, 
        col = "plum1")

# выбираем нужные строки
chip <- subset(dat, dat$assoc == "Chip")

# доверительный интервал для доли
Nchip <- nrow(chip)
N <- nrow(dat)

install.packages("DescTools")  # установить библиотеку
library(DescTools) # обратиться к ней

# 90%-ный интервал
ci90 <- BinomCI(Nchip, N, conf.level = 0.90)
ci90 * 100

# 95%-ный интервал
ci95 <- BinomCI(Nchip, N, conf.level = 0.95)
ci95

# доверительный интервал для среднего 
mci95 <- MeanCI(dat$reaction)
mci95

mci99 <- MeanCI(dat$reaction, conf.level = 0.99)
mci99

# сравним длины интервалов
mci95[3] - mci95[2]
mci99[3] - mci99[2]