Обработка данных эксперимента

x = c(7, 11, 8, 11, 10, 8, 9, 11, 7, 8, 11, 10)
CI = 98
edp = function(data, CI){
  
# Обработка данных эксперимента
# experience data processing edp
n = length(data)
Mean = mean(data)
SD = sd(data)
mx = SD/sqrt(n)# ошибка средн арифм
cat('среднее значение:', Mean,'\n')
cat('\n')
cat('ошибка средней:', mx,'\n')
cat('\n')
cat(Mean, "-", mx, "  ",Mean, "+", mx,'\n')
cat(Mean-mx, '...', Mean+mx,'\n')
cat('\n')
# Calculate the confidence interval
ans <- t.test(data,conf.level = CI/100)
# Extract the confidence interval
ci <- ans$conf.int
# Print the confidence interval
cat('Генеральная средняя с вероятностью', CI, '%\n')
cat('находиться в интервале\n')
print(ci[1:2])
}
x=c(3,4,5,4, 5, 2, 5, 5)
edp(data = x, CI = 95)
 

среднее значение: 4.125 

ошибка средней: 0.3980982 

4.125 - 0.3980982 4.125 + 0.3980982 
3.726902 ... 4.523098 

Генеральная средняя с вероятностью 95 %
находиться в интервале
[1] 3.183647 5.066353