自编函数做t检验

# 完整的检验作业

# 题目要求：
# 比较两个总体均值差异显著性检验的函数
# 要求函数输入参数只有两个，1.data 2.alph=0.05
# 输出结果：样本均值，样本方差，合并方差，检验统计量，t分布临界值，结论

# 检验函数编写

# 法1 自写函数t.test检验函数
# 这里传入数据的列名为a,b
t_test_fun = function(data,alph=0.05){
  # 列数据作为计算对象
  a = data$a
  b = data$b
  # 1.计算基本量
  # 1.1向量a,b长度
  n1 = length(a)
  n2 = length(b)
  # 1.2计算a,b均值
  a_bar = mean(a)
  b_bar = mean(b)
  message('a组均值：',a_bar,'\n')
  message('b组均值：',b_bar,'\n')
  # 计算a,b方差
  var1 = var(a)
  var2 = var(b)
  message('a组方差：',var1,'\n')
  message('b组方差：',var2,'\n')
  # 求合并方差(方差齐性)
  concat_vars = (var1*(n1-1)+var2*(n2-1))/(n1+n2-2)
  message('联合方差为：',concat_vars,'\n')  
  # 计算t统计量
  t_statistic = (a_bar-b_bar)/sqrt(concat_vars*(1/n1+1/n2))
  message('t统计量为：',t_statistic,'\n')
  # 计算自由度
  v = (n1+n2-2)
  message('自由度为：',v,'\n')
  #t分布临界值
  t_value = qt(alph,v)
  if (abs(t_statistic) > abs(t_value)){
    print('拒绝原假设，两组数据之间存在显著性差异')
  }else{
    print('接受原假设，两组数据之间不存在显著性差异')
  }
}

# 产生随机数
set.seed(0) 
A = round(runif(100,20,40))
B = round(runif(100,20,40))
data1 = data.frame(a = A,b = B)
t_test_fun(data1,alph = 0.05)

# 法2 验证
t.test(data1$a,data1$b,alternative = 'two.sided',
       conf.level = 0.95)