假设检验-两服从正态分布的独立总体均值检验

19 篇文章 9 订阅

传送:随机变量概率分布函数汇总-离散型分布+连续型分布

关注的结果变量为连续型组间比较(两组数据必须是独立的),并假设其呈现正态分布。首先判断是否为正态分布qqnorm(x1);qqline(x1)

传送:假设检验-KS检验

传送:假设检验-W检验

传送:假设检验-单样本t检验

假设条件:X,Y是两个独立的正态总体,X-N(\mu _{1},\sigma _{1}^{2}),Y-N(\mu _{2},\sigma _{2}^{2}),X1,X2...Xn是来自X的样本,Y1,Y2...Yn是来自Y的样本。样本的均值分别是\bar{X},\bar{Y},方差分别为S_{1}^{2},S_{2}^{2}

一、两正态总体方差均已知

当两个正态总体方差均已知时,在原假设\mu _{1}=\mu _{2}条件下,构造服从正态分布的检验统计量

双侧检验的拒绝域为|Z|\geq Z_{\alpha /2},单侧检验的拒绝域为Z\geq Z_{\alpha }Z\leq- Z_{\alpha}

	z.test2=function(x,y,sigma1,sigma2,alternative="two.sided"){
	n1=length(x);n2=length(y)
	result=list()
	mean=mean(x)-mean(y)
	z=mean/sqrt(sigma1^2/n1+sigma2^2/n2)	#构造z统计量
	options(digits=4)
	result$mean=mean;result$z=z
	result$P=2*pnorm(abs(z),lower.tail=FALSE)	#计算落入拒绝域的概率
	#单侧检验-重新计算P值
	if (alternative="greater")	#H0:µ1≤µ2,H1:µ1>µ2
    result$P=pnorm(z)    #参考[两正态总体方差未知但相等]的统计量表达式理解
	else if (alternative="less") result$P=pnorm(z,lower.tail=FALSE)
	result
	}

二、两正态总体方差未知但相等

当原假设为真时,构造服从t分布的检验统计量

双侧检验的拒绝域是|T|\geq t_{\alpha /2}(n_{1}+n_{2}-2),单侧检验的拒绝域是T\geq t_{\alpha}(n_{1}+n_{2}-2)或者T\leq - t_{\alpha}(n_{1}+n_{2}-2)

t.test(x1,x2,var.equal=T,conf.level=0.95)	#默认条件下是方差不相等
t.test(y~x,data)	#x为一个二分变量

三、两正态总体方差未知且不等

根据样本方差,构造服从t分布的检验统计量

双侧检验的拒绝域是|T|\geq t_{\alpha /2}(\nu ),单侧检验的拒绝域是T\geq t_{\alpha}(\nu )或者T\leq - t_{\alpha}(\nu )

t.test(x1,x2,var.equal=F,conf.level=0.95)

四、配对样本t检验

X1,X2..Xn是X的样本,Y1,Y2..Yn是Y的样本,令Z_{i}=X_{i}-Y_{i} (i=1,2,3,...),记\mu =\mu _{1}-\mu _{2}\sigma ^{2}=\sigma _{1}^{2}+\sigma _{2}^{2},则Z1,Z2,...Zn服从正态分布总体Z-N(\mu ,\sigma ^{2}),构造统计量

t.test(before,after,paired=TRUE) #非独立样本的t检验组间的差异呈现正态分布

 

  • 2
    点赞
  • 29
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
进行单正态总体均值假设检验可以使用 t 检验。假设总体服从正态分布,且方差未知样本量较小(小于30),可以使用 t 分布进行假设检验。 以下是一个单正态总体均值假设检验的例子: 假设现有一组数据 x,样本量为 n,要检验均值是否等于 mu0。则假设检验的步骤如下: 1. 建立假设。原假设 H0:mu = mu0,备择假设 Ha:mu != mu0。 2. 确定显著性水平 alpha。通常取 0.05 或 0.01。 3. 计算 t 统计量: ``` t <- (mean(x) - mu0) / (sd(x) / sqrt(n)) ``` 其中,mean(x) 是样本均值,sd(x) 是样本标准差,n 是样本量。 4. 计算自由度 df: ``` df <- n - 1 ``` 5. 确定临界值: ``` t.crit <- qt(1 - alpha/2, df) ``` 其中,qt 函数是 t 分布的分位数函数,1 - alpha/2 表示双侧检验的置信水平。 6. 判断决策。如果 |t| > t.crit,则拒绝原假设,接受备择假设;否则,无法拒绝原假设。 以下是一个完整的 R 代码示例: ``` # 建立假设 mu0 <- 100 alpha <- 0.05 x <- c(98, 102, 99, 101, 97, 103) # 计算 t 统计量 n <- length(x) t <- (mean(x) - mu0) / (sd(x) / sqrt(n)) # 计算自由度 df <- n - 1 # 确定临界值 t.crit <- qt(1 - alpha/2, df) # 判断决策 if (abs(t) > t.crit) { cat("拒绝原假设,接受备择假设。\n") } else { cat("无法拒绝原假设。\n") } ``` 在这个例子中,我们假设总体均值为 100,样本量为 6,样本数据为 x = c(98, 102, 99, 101, 97, 103)。根据计算,t 统计量为 -0.4714,自由度为 5,临界值为 2.571。由于 |t| < t.crit,因此无法拒绝原假设,即我们不能得出总体均值不等于 100 的结论。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值