Анализ различий в двух группах бутстреп подходом
Решим задачу. В нашем распоряжении две выборки массы грецких орехов у одного и того же дерева урожая этого и прошлого года. Объемы выборок отличаются. Необходимо выяснить различаются ли группы между собой (или более точно имеют ли группы статистически значимые отличия между собой). Решить такую задачу просто, если применить написанную мной функцию BRtwo.R как всегда не забудем поместить эту функцию в корневой каталог и указать путь к нему.
setwd('C:/R myFunction')
source('BRtwo.R')
m3 = c(11.95, 10.95, 12.56, 9.16, 10.08, 10.72, 9.29, 6.15,
5.12, 11.16, 9.51, 10.73, 10.44, 10.10, 10.10, 11.35,
13.05, 12.09, 11.09, 7.86, 5.12, 6.41, 8.76,
10.40, 11.36, 6.07, 11.11, 9.29, 10.25, 11.72,
9.82, 8.07, 9.20, 9.65, 9.87, 12.39, 10.72, 11.84,
9.81, 10.10, 11.04, 10.08, 12.91, 10.30, 12.80,
9.40, 8.65, 10.21, 9.50, 10.97, 11.65, 8.27, 9.90,
9.18, 9.83, 10.41, 10.24, 11.26, 12.37, 9.73, 6.10,
5.84, 12.00, 8.70, 10.15, 10.62, 5.21, 11.34, 9.39,
9.96, 4.80, 10.97, 11.51, 12.44, 11.17, 12.43, 9.41,
10.78, 11.05, 10.85, 10.68, 9.02, 12.37, 10.43,
11.96, 10.93, 10.59, 10.76, 8.27, 9.09, 9.32, 9.36,
10.51, 11.70, 11.76, 10.72, 7.23, 9.34, 9.39, 11.53,
9.85, 10.68, 11.19, 11.29, 10.14, 12.88, 8.85, 10.88,
11.17, 13.30, 10.45, 10.13, 12.07, 10.57, 10.20, 9.74,
8.09, 10.62, 11.12, 13.98, 11.87, 13.48, 11.00, 5.86,
11.24, 10.14, 10.58, 11.04, 9.03, 9.67, 10.73, 10.21,
9.99, 10.98, 5.49, 10.28, 10.70, 13.00, 11.98, 10.04,
10.79, 9.65, 10.05, 13.12, 9.41, 4.32, 10.77, 11.40)
m2 = c(11, 10.5, 10.7, 10.9, 7.1, 10.4, 11, 11.9, 10, 9.4, 10.1, 11.7,
10.5, 10.3, 10.5, 10.3, 10.9, 10.5, 10, 11.3, 12.5, 10.4, 9.9,
11.8, 6.2, 9.4, 10.3, 5.9, 9.9, 10.2, 10.6, 8.9, 9.9, 8.2, 11.1,
11.1, 6.2, 10.6, 9.7, 9.7, 10.3, 11.8, 10.3, 12, 11.1, 10.9, 9.5,
11.1, 11.2, 10, 8.6, 10.8, 10.5, 5.9, 10.6, 10.1, 9.7, 11.2, 11.2,
11.6, 11.2, 10.5, 5.6, 12.3, 10.2, 11.4, 10.8, 10, 8.9, 10.9,
9.2, 10.6, 11.9, 10.8, 10.6, 10.2, 11.3, 15.8, 10.5, 12.2, 10.8,
11.7, 6.7, 12.3, 9.4, 8, 10.1, 9.3, 9.8, 11.6, 9.8, 10.2, 8.9,
10.4, 10.4, 8.5, 10.4, 10.5, 10.1, 9.7, 9.5, 10.7, 12.1, 8.9,
10.3, 4.1, 10.2, 9.8, 9.8, 10.9, 9.4, 8, 11.2, 7.6, 11.2,
10.6, 11.2, 11.1, 9.7, 10.9)
BRt(data1 = m2, data2 = m3, R = 1000,CLevel = 0.95)
> BRt(data1 = m2, data2 = m3, R = 1000,CLevel = 0.95) описательная статистика исходных данных средние значения, стандартные отклонения --------------------------------------- | data1 | data2 | Mean1 SD1 | Mean2 SD2 10.172 1.556 == 10.188 1.836 различия в данных статистически незначимы 0.442624714056 |
Как видим все не просто а очень просто! Получили на выходе небольшой статистический отчет: различия в данных статистически незначимы. Что собственно и ожидалось.