R语言统计分析——功效分析3（相关、线性模型）

参考资料：R语言实战【第2版】

1、相关性

        pwr.r.test()函数可以对相关性分析进行功效分析。格式如下：

        pwr.r.test(n=, r=, sig.level=, power=, alternative=)

其中，n是观测数目，r是效应值（通过线性相关系数衡量），sig.level是显著性水平，power是功效水平，alternative指定显著性检验是双边检验（"two.sided"）还是单边检验（"less"或"greater"）。

        假设正在研究抑郁与孤独的关系，我们的零假设和备择假设分别是：

H0：

HA：

其中，是两个心理变量的总体相关性大小。我们设定显著性水平为0.05，而如果H0是错误的，我们想有90%的信心拒绝H0，那么研究需要多少观测呢？代码如下：

# 加载pwr包
library(pwr)
# 相关性功效分析
pwr.r.test(
  r=0.25,
  sig.level = 0.05,
  power=0.9,
  alternative = "greater"
)

        结果表明：我们需要134个受试者来评价抑郁与孤独的关系，以便在零假设为假的情况下有90%的信息拒绝它。

2、线性模型

        多于线性模型（比如多元回归），pwr.f2.test()函数可以完成相应的功效分析，格式为：

        pwr.f2.test(u=, v=, f2=, sig.level=, power=)

其中，u和v分别是分子自由度和分母自由度，f2是效应值。

 其中，R^2=多重相关性的总体平方值

 其中，=集合A中变量对总体方差的解释率；=集合A和B中变量对总体方差的解释率。

        当要评价一组预测变量对结果的影响程度，适宜采用第一个公式来计算f2；当药评价一组预测变量对结果的影响超过第二组变量（协变量）多少时，适宜采用第二个公式。

        现在假设我们想研究老板的领导风格对员工满意度的影响，是否超越薪水和工作小费对员工满意度的影响。领导风格可用四个变量来评估，薪水和小费与三个变量有关。过去的经验表明，薪水和消费能够解释约30%的员工满意度方差。从而现实出发，领导风格至少能解释35%的方差。假定显著性水平为0.05，那么在90%的置信度情况下，我们需要多少受试者才能得到这样的方差贡献率呢？

        此处，sig.level=0.05,power=0.9,u=3（总预测变量数减去集合B中的预测变量数），效应值为f2=(0.35-0.30)/(1-0.35)=0.0769。

# 线性模型功效分析
pwr.f2.test(
  u=3,
  f2=0.0769,
  sig.level=0.05,
  power=0.9
)

        在多元回归分析中，分母的自由度等于N-k-1，N是总观测数，k是观测变量数。本例中，N-7-1=185，即需要的样本大小为N=185+7+1=193。