分析思路 - EmpowerStats.PDF

重复测量数据分析 2.1 例2.1 对5 个国家过去若干年疟疾发病率与肿瘤死亡率关系分析，数 据如下表所示: Country Year Malaria Male Female Income El Salvador 1958 3.8513 22.53 38.97 low El Salvador 1959 7.0503 23.76 41.64 low …… El Salvador 2008 0.0055 57.1 63.51 low Greece 1961 0.0177 116.45 74.94 high Greece 1962 0.0156 118.97 74.59 high …… Greece 1981 0.0066 142.7 79.93 high Mexico 1969 0.9556 51.76 64.21 high Mexico 1970 1.0863 54.54 66.52 High …… Mexico 2008 0.0214 70.51 64.4 high Nicaragua 1973 1.8888 17.15 21.9 low Nicaragua 1974 5.2444 16.44 28.32 low …… Nicaragua 2008 0.1395 56.57 58.19 low Panama 1955 2.0722 57.86 60.1 low Panama 1956 3.2719 59.22 64.18 low …… Panama 2008 0.2248 89.59 68.02 low 点击这里下载 malaria_cancer.xls (注: 该数据非真实数据，仅供练习用，不 能用来进行科学推断) 数据特点 （1）每个国家有资料的年份(t)不完全相同；（2）X 变量（疟疾发病率） 与Y 变量（男、女肿瘤死亡率）都随时间在变化。 分析思路 混合线性模型原理参考例1.1 1、如何拟合X、t 与Y 的关系 这里每个国家为一个观察对象，国家之间基础的Y 水平不同，可以 用混合效应模型，设置截距为随机效应。同前所述，我们可以建一个简 单模型： Y = β + β*YEAR + β*X +e i 0i 1 2 i 这个模型中，我们假定了Y 随t （YEAR）是直线性的变化，如果不做 这样的假定，用非参数的曲线来拟合t,就把模型简单地变为： Y= β +s(YEAR) + β*X +e 0i 2 这里s(YEAR)即表示同spline smoothing 拟合Y 与YEAR 的关系。当 然X 与Y 的关系也可以用曲线来拟合。 Y= β +s(YEAR) +s(X) +e 0i 2、两分类变量INCOME 与X 的交互作用 这里有个INCOME 变量，是否INCOME 高与低不同，X 对Y 的作用也不 同呢？可以用模型： Y= β +s(YEAR) + β*