对数线性模型在高维列联表资料分析中的应用思考与练习参考答案.doc

对数线性模型在高维列联表资料分析中的应用思考与练习参考答案一、最佳选择题1. 对数线性模型优越性在于（ E ）。A. 可有效地分析高维列联表资料B. 可进行多个分类变量间关系的探索性分析C. 既可分析反应变量与解释变量间的关系，也可分析解释变量间的关系D. 可有效避免混杂因素的影响，提高分析效率E. 以上答案都对 2. 对数线性模型中对主效应解释正确的是（ D ）。A. 主效应大小反映了自变量对因变量的影响B. 主效应反映了因变量对期望频数的贡献C. 主效应反映了自变量对期望频数的贡献D. 对主效应的分析无实际意义 E. 主效应无统计学意义时，不必再拟合其他谱系模型3. 对数线性模型中对交互效应解释，描述不正确的是（ B ）。A. 交互效应指某几个自变量对因变量的联合作用B. 某两维交互效应为零，则包含该效应的三维交互效应必为零 C. 某三维交互效应不为零，并不表明该效应包含的所有两维交互效应也不为零D. 高维列联表中，某单元格的频数越大，交互效应越大E. 交互效应的大小受主效应大小的影响4. 列联表资料中存在有序分类变量时，正确的处理方法是（ C ）。A. 一律按无序分类变量引入模型B. 一律按有序分类变量引入模型C. 按有序还是无序变量引入，须视拟合效果和专业解释而定 D. 此种情形，不宜拟合对数线性模型E. 以上答案都不对5. 不能反映模型拟合效果的统计量或方法是（ D ）。A. Pearsons B. 残差分析 C. 似然比统计量D. 参数检验 E. AIC二、思考题1. 简述对数线性模型的分析思路。答： 判断适用条件，包括设计类型、样本含量。 拟合饱和模型，通过参数初步了解可能有统计学意义的效应项。 采用逐步法，通常采用后退法筛选最优模型。获得最优模型的参数估计结果。 结合专业解释结果。2. 对于高维列联表，采用对数线性模型比分层分析有什么优点？答： 提高分析效率。 可同时考察更多的变量的关系。 可检验多个分类变量间是否存在高维交互效应。3. 对比对数线性模型与Logistic模型的应用场合。答： 对数线性模型不属于因果关系模型，通常用于高维列联表资料的探索性分析。对数线性模型通常要求变量都为定性变量，如果为定量变量，需离散化。 Logistic模型描述一个分类变量与多个变量之间的关系，有明确的反应变量与解释变量，用于因果关系的描述，所以常用于验证性的分析。Logistic模型对自变量较宽容，可以是定量的，也可以是定性的，分布方面也没有特别的要求。三、计算题1. 引用第9章例9-5的数据，有研究表明不同国籍人的血型是不同的。现有2 500例不同国籍人的血型分布资料教材表20-21，请问国籍与血型是否有关？教材表20-21 2 500例不同国籍人的血型分布国籍血 型合 计OABAB美国人450410100401 000中国人3002503501001 000挪威人190250 40 20 500合 计9409104901602 500请用对数线性模型分别拟合主效应模型与饱和模型，并与第9章关联性分析结果对比，看看彼此之间有何联系？解：SPSS数据输入格式如下（练习图20-1）：练习图20-1 SPSS数据输入格式主要分析结果如下：（1）主效应分析结果Goodness-of-Fit Tests(a,b) ValuedfSig.Likelihood Ratio331.1476.000Pearson Chi-Square332.9676.000a Model: Poissonb Design: Constant + 国籍 + 血型（2）饱和模型分析结果Parameter Estimates(b,c)ParameterEstimateStd. ErrorZSig.Constant3.020.22113.676.000血型 = 12.229.2329.591.000血型 = 22.503.23010.896.000血型 = 3.681.2712.512.012国籍 = 1.681.2712.512.012国籍 = 21.590.2426.560.000血型 = 1 * 国籍 = 1.180.284.632.527血型 = 1 * 国籍 = 2-1.134.259-4.371.000血型 = 2 * 国籍 = 1-.187.283-.661.508血型 = 2 * 国籍 = 2-1.590.258-6.155.000血型 = 3 * 国籍 = 1.228.329.693.488血型 = 3 * 国籍 = 2.568.2941.935.053a This parameter is set to zero because it is redundant.b Model: Poissonc Design: Constant + 血型 + 国籍 + 血型 * 国籍（3）结果分析主效应分析结果提示，模型拟合效果不好，说明血型与国籍两个变量并不独立，拟合优度Pearson =332.967，与第9章所求统计量完全相等，读者可自行考虑为什么。饱和模型的参数估计结果提示，两维交互效应的参数估计值有统计学意义，表示两变量具有关联性。2. Sewell与Shan对Wisconsin州的10 317名高中生是否计划上大学的一项调查，调查变量包括家庭的社会经济地位(socioeconomic status,低/中低/中高/高)、智力水平(intelligence,低/中下/中上/高)、是否计划入学(plan college，是/否)、性别(sex,男/女)和父母重视鼓励程度(parental encouragement,低/高)。5个变量分别以T、I、P、S、E表示，构成44222的高维表，见教材表20-22。试分析与“是否计划入学”有关的因素。教材表20-22 对Wisconsin州的10 317名高中生是否计划上大学的调查资料计划入学(P)是否性别(S)男女男女父母重视(E)高低高低高低高低地位(T)智力(I)低低 4 13 5 9349 64455 44中低 9 3315 14207 72312 47中高12 38 8 20126 54216 35高10 8913 28 67 43 96 24中低低 2 2711 29232 84285 61中低 7 6419 47201 95236 88中高12 9312 62115 92164 85高1711915 72 79 59113 50中高低 8 47 7 36164 91163 72中低 4 7413 75120110193 90中高1714812 91 92100174100高 619820142 42 73 81 77高低 4 391636 48 57 50 58中低 5123 5110 47 90 70 76中高 922413230 41 65 47 81高 841413360 17 54 49 98解：主要分析结果：用后退法获得最优模型PETI，SET，STI，PSI，偏关联似然比检验如下：If Deleted Simple Effect is DF L.R. Chisq Change ProbP*E*T*I 9 23.581 .0050S*E*T 3 17.468 .0006S*T*I 9 22.788 .0067P*S*I 3 32.371 .0000分析结果提示：各变量对入学态度都有不同的影响，其中智力、社会地位、性别与入学态度具有联合作用，父母鼓励重视程度对学生智力及是否入学具有较强的影响。模型参数估计结果（剔除了冗余项）如下表：Parameter EstimatesParameterEstimateStd. ErrorZSig.Constant4.523(a) E = 1-.716.152-4.700.000I = 1-.537.145-3.691.000I = 2-.160.130-1.231.218I = 3-.108.129-.839.401P = 11.394.09714.392.000S = 1-.429.097-4.411.000T = 1-1.317.168-7.848.000T = 2-.556.137-4.044.000T = 3-.125.124-1.006.314E = 1 * I = 1.719.2033.548.000E = 1 * I = 2.528.1912.756.006E = 1 * I = 3.341.2001.702.089P = 1 * E = 1-2.732.266-10.262.000S = 1 * E = 1-.333.114-2.933.003E = 1 * T = 12.068.2189.493.000E = 1 * T = 21.481.1977.509.000E = 1 * T = 3.739.1973.750.000P = 1 * I = 1-1.693.188-8.988.000P = 1 * I = 2-1.040.149-6.972.000P = 1 * I = 3-.419.142-2.952.003S = 1 * I = 1.558.1573.555.000S = 1 * I = 2.538.1383.907.000S = 1 * I = 3.164.1361.206.228T = 1 * I = 1.978.2274.305.000T = 1 * I = 2.765.2153.549.000T = 1 * I = 3.488.2232.192.028T = 2 * I = 1.736.2013.662.000T = 2 * I = 2.600.1853.244.001T = 2 * I = 3.586.1853.172.002T = 3 * I = 1.268.1931.386.166T = 3 * I = 2.319.1741.827.068T = 3 * I = 3.290.1731.681.093P = 1 * S = 1.512.0915.633.000P = 1 * T = 1-1.188.180-6.620.000P = 1 * T = 2-1.129.151-7.491.000P = 1 * T = 3-.843.133-6.329.000S = 1 * T = 1.969.1436.776.000S = 1 * T = 2.492.1204.108.000S = 1 * T = 3.262.1062.460.014P = 1 * E = 1 * I = 11.549.3923.953.000P = 1 * E = 1 * I = 2-.069.436-.159.874P = 1 * E = 1 * I = 3.238.370.644.519E = 1 * T = 1 * I = 1.411.2711.518.129E = 1 * T = 1 * I = 2.045.262.173.863E = 1 * T = 1 * I = 3.082.276.298.765E = 1 * T = 2 * I = 1.003.254.013.990E = 1 * T = 2 * I = 2-.208.243-.856.392E = 1 * T = 2 * I = 3-.450.253-1.783.075E = 1 * T = 3 * I = 1.267.2561.041.298E = 1 * T = 3 * I = 2.155.245.632.528E = 1 * T = 3 * I = 3.173.253.683.494P = 1 * E = 1 * T = 1.324.380.852.394P = 1 * E = 1 * T = 2.433.3491.242.214P = 1 * E = 1 * T = 3.431.3571.207.227S = 1 * E = 1 * T = 1-.552.149-3.701.000S = 1 * E = 1 * T = 2-.098.138-.709.479S = 1 * E = 1 * T = 3-.236.139-1.700.089P = 1 * S = 1 * I = 1-.768.158-4.867.000P = 1 * S = 1 * I = 2-.540.132-4.106.000P = 1 * S = 1 * I = 3-.177.126-1.403.161P = 1 * T = 1 * I = 1.054.331.164.870P = 1 * T = 1 * I = 2-.077.269-.286.775P = 1 * T = 1 * I = 3-.427.265-1.613.107P = 1 * T = 2 * I = 1.611.2642.316.021P = 1 * T = 2 * I = 2.291.2181.333.182P = 1 * T = 2 * I = 3-.167.210-.795.427P = 1 * T = 3 * I = 1.618.2412.561.010P = 1 * T = 3 * I = 2.209.1991.049.294P = 1 * T = 3 * I = 3-.140.190-.737.461S = 1 * T = 1 * I = 1-.501.196-2.552.011S = 1 * T = 1 * I = 2-.587.183-3.217.001S = 1 * T = 1 * I = 3-.320.185-1.727.084S = 1 * T = 2 * I = 1-.396.185-2.14