作业4.nlme:lme做分层数据分析

STATA 数据下载地址:

http://www.ats.ucla.edu/stat/stata/examples/mlm_imm/imm10.dta

数据中有很多变量,关心schid,math,homework三个。schid为学校编号。math为数学成绩。homework为花在数学作业上的时间。模型为写作业的时间影响数学成绩,同时学校是一个更高层的变量,不同学校的学生数学作业的量不同,所以学校会影响写作业的时间,从而影响数学成绩。

公式如下:(i表示与个体有关的参数,j代表与学校有关的参数)

mathij=A0i+A1i*homework+eij

#数学成绩=截距+斜率*做作业时间+测量误差(误差与学校和个人都有关系)

A0i=B00+B01*schid+u0i

#截距A0i是一个与所在学校有关的函数,=j截距B00*所在学校+误差(此误差是个体波动)

A1i=B10+B11*schid+u1i

#同上,斜率A1i=截距+斜率*学校+误差

将后两个式子代入第一个式子进行整理,把只与个人有关的项放在一起,即B00+(B10+u1i)*homework,剩下的是与学校有关的项(B01+B11*homework)*schid,以及误差项。

B00+(B01+u1i)*homework可以看成系数*1+系数*homework,对应于下文语法中的math~1+homework,,称为固定效应;与学校有关的项(B01+B11*homework)*schid对应于random=~(1+homework)|schid。



作业2.spss syntax和R分析简单效应

数据下载地址:http://core.ecu.edu/psyc/wuenschk/SPSS/ARC.sav

该数据有很多变量,此处关心三个变量Justify,Group,Gender。数据说的是不同动物被用于科研,研究者试图知道哪些变量影响我们对一个研究的伦理性的评价。其中Justify是因变量,意思是被试认为某研究符合伦理的程度;Group变量分为五个水平,即5种不同的科研用途——Cosmetic,Theory,Meat,Veterinary,Medical;Gender为被试性别,2水平。

本例关心两个简单效应,其虚无假设为:

1.Group中的Theory水平上,男性与女性的Justify评价没有差异。

2.Group中的Meat水平上,男性与女性的Justify评价没有差异。

SPSS分析过程:

UNIANOVA Justify BY Group Gender
/LMATRIX="Group2Gender12" Group*Gender 0 0 1 -1 0 0 0 0 0 0 Gender 1 -1
/LMATRIX="Group3Gender12" Group*Gender 0 0 0 0 1 -1 0 0 0 0 Gender 1 -1
/METHOD=SSTYPE(3)
/INTERCEPT=INCLUDE
/CRITERIA=ALPHA(.05)
/DESIGN= Gender Group Gender*Group

Justify BY Group Gender的格式为 因变量 BY 两个自变量

/LMATRIX=''引号内请自行给假设命名" 自变量1*自变量2  矩阵(矩阵中把需要比较的位置写成非0,此处为一个5行2列矩阵)

/DESIGN= Gender Group Gender*Group为你所定义的模型

结果1:主效应检验(主效应显著,交互不显著)

1

说明:一般来说,当交互作用显著时才关心简单效应,但是由于暂时找不到更合适的数据,为了演示做法,请假装交互作用显著。

结果2:简单效应检验(均不显著,虚无假设均成立,Group中的Theory/Meat水平上,男性与女性的Justify评价没有差异)

2 (3)

 

3 (2)

 

R分析过程:



作业1.长宽数据的转换与矩阵的提取拼接

长宽数据的转换用R完成,矩阵的提取拼接用excel完成。

宽数据来源:百度文库某PPT第39页

宽数据的源地址:http://wenku.baidu.com/link?url=G5pMx-gMQhOLbX6N2yXSY4CMsIf4yxzCjU8tJ88S9GwxTb9c1ag3tCm-1GgmFh6PkYc9jtwZfTgK97ja0veP1zV1RLXAS89jTHyjT44w9cy

宽数据及相关R文件下载地址:http://vdisk.weibo.com/s/udMx9gHGWlL7-

相关矩阵数据来源:过往研究数据

相关矩阵数据下载地址:http://vdisk.weibo.com/s/udMx9gHGWlKIN

视频观看地址:http://v.youku.com/v_show/id_XMTQ3MzY3NTM1Mg==.html

(视频转码后较模糊,建议打开高清模式或下载观看1080P版本)

视频下载地址:http://vdisk.weibo.com/s/udMx9gHGWl-SZ

所有数据及视频打包下载地址:http://vdisk.weibo.com/s/udMx9gHGWlKo-