全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 Stata专版
2370 2
2024-01-09


方差分析(analysis of variance,ANOVA)是在多个总体或者组(group)之间进行比较的方法。在方差分析中,从每个组中选取样本来检验两个或更多组中差异的影响。分组的条件被称为因素(factor),有时也称兴趣因素。构成因素的不同水平(level)类似于属性变量下包括的类别。


在最简单的单因素方差分析(one-way ANOVA,又称为完全随机设计)中,只检验一种因素,构成因素的不同水平称为变量分组的基础。单因素方差分析是一个包括两部分的过程。首先确定各组平均值之间是否有显著性差异。如果拒绝了零假设代表的各组平均值之间没有差别,那么再以第二种方法继续寻找哪些组的平均值明显的不同于其他组的平均值。


oneway x1 x2 //x2因素作用下的不同水平,x1变量在每组中有相对应的数值,执行后可以得到ANOVA的分析结果


/*df是自由度   degree of freedom

SS:离均差平方和 也就是变量中每个数据点与变量均值差的平方和

MS: mean square =SS/DF

F 值  = 组间差异 /组内差异  处理的MS/误差的MS MSA/MSW

造成波动的原因可分成两类,,一是不可控的随机因素:误差的MS-MSA,另一是研究中施加的对结果形成影响的可控因素处理的MS-MSW。

看F的p值  p远小于<0.05说明至少有两个组是有显著差别的

最后一行是方差齐性检验(test for equal variances),p值>0.05说明满足方差齐性*/


oneway x1 x2, bonferroni//如果至少两个组显著有差别,用bonferroni检验两两比较是哪两个组

//在结果中两两比较组之间的p如果<0.05,则说明有显著差异,正负号说明是显著大还是小于另一个


oneway x1 x2, bonferroni tabulate // 看每组详细的频数统计表格

codebook x1 x2//查看每个变量描述性统计信息内容


要使用单因素方差分析F检验,你必须对所要考察的数据作出三个假设:


样本选择的随机性与独立性(randomness and independence)

样本来源的c个组的正态性(normality)

同方差(homegeneity of variance)(c个组的方差相同)


如果每组中都有相同大小的样本容量,在F分布基础上推断就不会被异方差问题严重影响。可以使用同方差Levene检验来检验c个组的方差是否相等。

有许多方法可以用来确定c个均值的哪些是显著不同的,一个最常用的方法是单因素方差分析Tukey-Kramer多重比较过程(Tukey-Kramer multiple comparisons procedure for one-way ANOVA)


二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2024-1-10 13:03:50
太棒啦
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2024-1-11 01:39:37
这里句一个利用面板数据做ANOVA检验的STATA例子

/*导入数据*/
import excel "/Users/data.xlsx", sheet("Sheet1") firstrow
//这里要注意,要做ANOVA的样本数据一般都是宽数据,要想在STATA里做处理,需要将其转化成长数据,所以下面对面板数据做下处理,假设作用因素变量有五个值,每个值作用10个样本(及另一个变量)

/*假设变量名为字符串,对变量名做预处理*/
rename Kidney EffectKidney
rename Shrimp EffectShrimp
rename ChickenLiver EffectChickenLiver
rename Salmon EffectSalmon
rename Beef EffectBeef
/*生成一个新的变量ID,方便做数据转换*/
gen ID= ., before (EffectKidney)
replace ID = _n if _n <= 10 //对变量赋值1-10

reshape long Effect, i(ID) j(Type) string //数据转换,从wide变long数据,这里的string是因为转换后的Type变量下的数值为字符串类型,否则会报错
//两篇reshpe命令的解释讲解文章https://mp.weixin.qq.com/s/RLaU1vkAlmUxGIQL7IqUrA
//https://mp.weixin.qq.com/s/lZ_psi8eRYkBfRSLUCdUOQ

/*对Type变量值做调*/
replace Type = "1" if Type == "Beef"
replace Type = "2" if Type == "ChickenLiver"
replace Type = "3" if Type == "Kidney"
replace Type = "4" if Type == "Salmon"
replace Type = "5" if Type == "Shrimp"
destring Type,replace //去字符类型

/*最后进行ANOVA检验*/
oneway Effect Type



二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群