Сравнение двух групп бутстрепом
setwd('C:/R myFunction')
source('BRtwo.R')
Fdata = c(14.9, 15.7, 18.0, 15.3, 15.9)
Gdata = c(18.1, 18.5, 18.2, 19.8, 17.5)
BRt(data1 = Fdata, data2 = Gdata, 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
|