一、何为显著性检验

显著性检验的思想十分的简单，就是认为小概率事件不可能发生。虽然概率论中我们一直强调小概率事件必然发生，但显著性检验还是相信了小概率事件在我做的这一次检验中没有发生。

显著性检验即用于实验处理组与对照组或两种不同处理的效应之间是否有差异，以及这种差异是否显著的方法。

常把一个要检验的假设记作H0,称为原假设(或零假设)，与H0对立的假设记作H1，称为备择假设。

⑴在原假设为真时，决定放弃原假设，称为第一类错误，其出现的概率通常记作α；

⑵在原假设不真时，决定接受原假设，称为第二类错误，其出现的概率通常记作β。

通常只限定犯第一类错误的最大概率α，不考虑犯第二类错误的概率β。这样的假设检验又称为显著性检验，概率α称为显著性水平。

我们常用的显著性检验有t检验，卡方检验，相关性检验等，在做这一些检验时，有什么需要注意的呢？

二、正态性与P值

t检验，卡方检验，相关性检验中的pearson方法都是建立在正态样本的假设下的，所以在假设检验开始时，一般都会做正态性分析。在R中可以使用shapiro.test()。来作正态性检验。当然在norm.test包中还提供了许多其他的方法供我们选择。

P值是可以拒绝原假设的最小水平值。

三、四个重要的量

综合前面的叙述，我们知道研究显著性检验有四个十分重要的量：样本大小，显著性水平，功效，效应值。

样本大小：这个显然，样本越多，对样本的把握显然越准确，但是鉴于我们不可能拥有无限制的样本，那么多少个样本可以达到要求？今天的分享中我们可以通过R来找到答案。

显著性水平：犯第一类错误的概率，这个在做检验前我们会提前约定，最后根据P值来决定取舍。

功效：这个是在显著性检验中一般不提及但实际十分有用的量。它衡量真实事件发生的概率。也就是说功效越大，第二类错误越不可能发生。虽然显著性假设检验不提及它，但衡量假设检验的好坏的重要指标便是两类错误尽可能小。

效应值：备择假设下效应的量

四、用pwr包做功效分析

Pwr包中提供了以下函数：

下面我们来介绍以上一些函数的用法。

1、 t检验

调用格式：

pwr.t.test(n = NULL, d = NULL, sig.level =0.05, power =

NULL, type =c("two.sample", "one.sample",

"paired"),alternative = c("two.sided", "less","greater"))

参数说明：

N:样本大小

D:t检验的统计量

Sig.level：显著性水平

Power：功效水平

Type：检验类型，这里默认是两样本，且样本量相同

Alternative：统计检验是双侧还是单侧，这里默认为双侧

举例说明：已知样本量为60，单一样本t检验的统计量的值为0.2(这个可以通过t.test(data)$statistic取出来)，显著水平α=0.1，那么功效是多少呢？

R中输入命令：

[plain] view plain copy

pwr.t.test(d=0.2,n=60,sig.level=0.10,type="one.sample",alternative="two.sided")

得到结果：

One-sample t test power calculation

n = 60

d = 0.2

sig.level = 0.1

power = 0.4555818

alternative = two.sided

我们可以看到，犯第二类错误的概率在50%以上，我们应该相信这个结果吗(无论根据P值来看是拒绝还是接受)？显然不行，那么需要多少个样本才能把第二类错误降低到10%呢？

在R中输入：

[plain] view plain copy

pwr.t.test(d=0.2,power=0.9,sig.level=0.10,type="one.sample",alternative="two.sided") 得到结果：

One-sample t test power calculation

n = 215.4542

d = 0.2

sig.level = 0.1

power = 0.9