BS разность средних

Анализ различий данных двух малочисленных выборок с помощью бутстрепа

# μR/h (Microroentgen per Hour) mRh

# FmRh - измерение радиционного фона

# измерение  радиционного фона стекла

FmRh = c(14.9, 15.7, 18.0, 15.3, 15.9)

GmRh = c(18.1, 18.5, 18.2, 19.8, 17.5)

# Сравнение двух независимых выборок:

t.test(FmRh, GmRh)

# p-value = 0.006981 < 0.05 нулевая гипотеза

# о равенстве средних значений двух выборок отвергнута

# радиционный фон и фон вулканического стекла различаются

# разница средних значений:

mean(FmRh) - mean(GmRh)

# -2.46 μR/h

#------------------------------------------------------

# извлекаем бустреп выборки и вычисляем их разность средних

boot_ratio = function(data1, data2){

  n1 = length(data1)

  boot1 = sample(n1, replace=TRUE)

  db1 = data1[boot1]

  n2 = length(data2)

  boot2 = sample(n2, replace=TRUE)

  db2 = data2[boot2]

 

  mean_boot = mean(db1)-mean(db2)

  return(mean_boot)

}

 

#повтор R бутстрепинг без for

#boot_ratio(data = data2)

R = 1000

out = replicate(R, boot_ratio(data1 = FmRh, data2 = GmRh))

hist(out,20,freq = F)               

lines(density(out), col = 'blue', lwd = 2)

abline(v = mean(out), lwd=2, col='red')

grid(col = 'grey70')

# разница средних значений:

mean(FmRh) - mean(GmRh)

# разность средних значений бутстреп:

mean(out)

t.test(x = out,conf.level = 0.99)$conf.int