样本量估算及G*Power实现详细步骤(图文详解)

样本量估算及G*Power实现详细步骤(图文详解)

1.统计学基本概念

(1)假设检验、H0、H1

     假设检验(hypothesis

test)又称显著性检验(significant test),我们进行研究的目的多是为了比较总体参数差异,但实际中能够获取的是样本统计量数据,样本数据间直接比较不能准确反应总体参数情况。样本统计量比较存在差异时,可能是总体参数差异引起的,也可能是样本的抽样误差引起的。假设检验中,引入原假设(null hypothesis)H0,和与H0完全对立的备择假设(alternative hypothesis)H1

H0:总体参数无差异

H1:总体参数有差异

可以根据样本统计量的差值的数据类型、研究设计、推断目的,选择不同的检验统计量计算方法(t检验、z检验、F检验、卡方检验等),计算出得到该差值的概率P,

如果P≤α(预设的检验水准),可倾向于拒绝H0,倾向于接受H1;

如果P>α,则不足以拒绝H0,不足以接受H1。

(2)I类错误(α)、II类错误(β)

当H0成立,但因为样本的随机性,拒绝了H0,为I类错误,又称假阳性错误,检验水准α即为预先规定的允许犯I类错误的概率。

当H1成立,但因为样本的随机性,不拒绝H0,为II类错误,又称假阴性错误,用β表示(需要总体参数及样本含量才能计算出)。

正如红楼梦中所说,假作真时真亦假,其中假作真即为假阳性,I类错误;真亦假即为假阴性,II类错误。

当I类错误概率α约小,出现II类错误概率β越大,反之亦然。

增大样本量,可以同时减少α和β。

(3)检验效能

检验效能(power of a test),又称把握度,H1成立时正确拒绝H0的能力,用1-β表示。

把握度大小取决于四个方面:①两总体差异越大,1-β越大;②α越大,1-β越大;

③单侧检验比双侧检验大;④样本量越大,1-β越大。

(4)效应量

效应量(effect size),表示两总体分布的重叠程度,也就是自变量作用的大小。常见的效应量包括两总体均数的差值、两总体率的差值等。差值越大,表示两总体重叠程度越小,效应越明显。

(5)效应量与差异显著性、相关显著性的区别

差异显著性、相关显著性只是告诉我们,我们有1-α的把握认为,这差异、这相关系数是存在的,并不是完全由抽样误差造成的。但这代表差异、相关有实际的意义。尤其是大样本数据分析比较容易得到统计显著性的结果,但不意味着差异有意义。

效应量,表示的是自变量作用的大小

(6)样本量

样本量(sample size),样本量过小,假设检验效能不够,无法显现出不同处理组间的差异,样本量过大,需要增加大量人力、物力的投入。

样本量大小受以下因素影响:①α越小,样本量越大;②β越小,样本量越大;

                          ③总体间效应量越小,即两总体重叠约小,所需样本量越大

                          ④总体变异性越大,所需样本量越大

2.G*Power 软件安装及软件介绍

(1)G*Power安装

下载地址(www.gpower.hhu.de


解压缩文件夹后,双击setup,进行安装。

打开后可见如下主界面。

(2)G*Power介绍

G*Power是德国杜塞尔多夫大学开发的,专门用于计算统计功效及样本量的统计软件。该软件是完全免费的。软件可以计算t检验、F检验、χ2检验、z检验、Fisher检验等不同检验方法的样本量。

G*Power提供5种统计力分析类型:


 3.G*Power计算样本量

使用G*Power分析计算分为3步:

S1:选择合适的统计方法

S2:选择5种统计力分析类型中的1种

S3:输入分析需要的参数,点击计算即可

(1)如下图,从Test family中选择t检验、F检验、χ2检验、z检验或Fisher检验。

从Statistical test中进一步选择分类

(2)G*Power中有5种统计力分析类型可选,计算样本量选择先验分析(priori)。


(3)下框中输入分析需要的参数,点击计算即可

4.实例演示-两独立样本均数比较的样本量计算

首先看第6版医学统计学人卫版课本P159例15-4,该例子使用公式计算出了需要的样本量。


使用G*Power,选择两组独立样本t检验,输入参数:双尾、d=δ/σ=2.76/1.70=1.62352941、α=0.05、1-β=0.9、N2/N1=1,点击计算得到N1=10,N2=10,(公式算出n1=n2=8)

使用G*Power,选择两组独立样本t检验,输入参数:双尾、d=δ/σ=2.76/1.70=1.62352941、α=0.05、1-β=0.9、N2/N1=2,点击计算得到N1=7,N2=15,(公式算出n1=6,n2=12)


可发现使用G*Power得到的样本量与使用公式得到的样本近似,但比公式结果略大。可能因为计算过程中精确度选择原因导致。

©著作权归作者所有,转载或内容合作请联系作者
禁止转载,如需转载请通过简信或评论联系作者。
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 204,530评论 6 478
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 86,403评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 151,120评论 0 337
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,770评论 1 277
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,758评论 5 367
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,649评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 38,021评论 3 398
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,675评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,931评论 1 299
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,659评论 2 321
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,751评论 1 330
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,410评论 4 321
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 39,004评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,969评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,203评论 1 260
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 45,042评论 2 350
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,493评论 2 343

推荐阅读更多精彩内容