Ответы, отчеты, курсовые, лекции 1-3 курс (1, 2, 3, 4, 5 семестр) ЧГУ КТ-43-15 / 3 семестр / R Language лабы (Сковорцов) / 6

#6.1
Задание 6.1. Воспользуемся еще одним встроенным набором данных в R -
ToothGrowth. Данные позволяют исследовать рост зубов у морских свинок
в зависимости от дозировки витамина C и типа потребляемых продуктов.
Сравните среднее значение длины зубов свинок, которые потребляли
апельсиновый сок (OJ) с дозировкой 0.5 миллиграмм, со средним значением
длины зубов свинок, которые потребляли аскорбиновую кислоту (VC) с
дозировкой 2 миллиграмма.
Значение t - критерия сохраните в переменную t_stat.

df <- read.csv("ToothGrowth.csv")
str(df)
lenght_VC<-ToothGrowth$len[ToothGrowth$supp=="VC"&ToothGrowth$dose==0.5]
lenght_OJ<-ToothGrowth$len[ToothGrowth$supp=="OJ"&ToothGrowth$dose==2]
t_stat -> t.test(lenght_VC,lenght_OJ)$statistic
t_stat

#6.2
Задание 6.2. В этом задании нужно проверить гипотезу о равенстве средних двух
выборок. Загрузите один из доступных наборов данных по ссылке:
https://goo.gl/VSIZ9c. Сначала с помощью теста Бартлетта проверьте гомогенность
дисперсий двух выборок. В случае, если дисперсии значимо не отличаются (с уровнем
0.05), примените тест Стьюдента, иначе - непараметрический тест (Манна-Уитни).
Обратите внимание, что по умолчанию в t.test стоит var.equal = FALSE, так как мы
будем применять его только в случае гомогенности (однородности) дисперсий,
измените значение этого параметра на var.equal = TRUE. Сохраните в переменную
main_stat получившийся p-value, с точностью четыре знака после запятой.

t <- read.table("dataset_1.txt")
bt <- bartlett.test(V1 ~ V2, t)$p.value
bt
if (bt == 0.05)
{
main_stat <- round(t.test(V1 ~ V2, t, var.equal = TRUE)$p.value, 4)
}else
{
main_stat <- round(wilcox.test(V1 ~ V2, t)$p.value, 4)
}
main_stat