构造并分析互动效应的线性回归模型,解释结果。
加载医院
仅保留第一列的血压,将数据存储在桌子中。
TBL =表(医院.Sex,Hospital.age,Hospital.weight,Hospital.Smoker,Hospital.blocthpressure(:,1),......'variablenames',{'性别'那'年龄'那'重量'那“吸烟者”那'血压'});
对于初始模型,请使用所有术语及其成对交互的完整模型。
mdl = stepwiselm(tbl,'互动'的)
1.去除性别:吸烟者,FSTAT = 0.050738,PVALUE = 0.8223 2.卸下体重:吸烟者,FSTAT = 0.07758,PVALUE = 0.78124 3.去除年龄:重量,FSTAT = 1.9717,PVALUE = 0.16367 4.去除性别:年龄,FSTAT= 0.32389,pvalue = 0.57067 5.去除年龄:吸烟者,Fstat = 2.4939,pvalue = 0.11768
MDL =线性回归模型:血压〜1 +年龄+吸烟*体重估计系数:估计系数估计系数______________________7.17.19.301019.0126917.19.0126重量-0.1393 0.080211 -1.7367 0.085722吸烟者_1 9.8307 1.0229 9.6102 1.2391C-15性别:重量0.2341 0.11192 2.0917 0.039162观察数:100,误差自由度:94根均匀误差:4.72 R线:0.53,调整R线:0.505 F统计与常数型号:21.2,P值= 4E-14
式形式中的最终模型是血压〜1 +年龄+吸烟+性别*重量
。该模型包括所有四种主要效应(年龄,吸烟者,性别,重量)和双向相互作用性别
和重量
。该模型对应于
在哪里
是血压
是系数
是吸烟的指标变量; 表明吸烟患者 表示非莫斯文患者
是性别的指标变量; 表示男性患者 表示女病人
是个年龄
多变的
是个重量
多变的
是错误的术语
下表显示了每个性别和吸烟组合的拟合线性模型。
从这些模型中看到,
和
当指示变量将值1相比,显示响应函数的截距变为多少变化,而与值为值0相比。
但是,当与何时需要值0时,当对性别的指示器变量取值时,显示重量变量在响应变量上的效果。您可以使用这些方法探索最终模型中的主要和交互效果linearmodel.
课程如下。
图Plotslice(MDL)
此曲线显示所有预测变量的主要效果。每个面板中的绿线显示当所有其他预测变量保持恒定时,每个面板中的变化作为预测变量的函数。例如,对于37.5岁的吸烟的雄性患者,随着患者的重量增加,预期血压会增加,鉴于其他相同。
每个面板中的虚线红色曲线显示预测响应值的95%置信度界限。
每个面板中的水平虚线显示了对应于垂直虚线的预测变量的特定值的预测响应。您可以拖动这些行以获取其他预测值值的预测响应值,如下所示。
例如,当患者是女性的,响应变量的预测值是118.3497,非镜头,40.3788磅,重139.9545磅。方括号中的值[114.621,122.079],显示了95%置信区间的下限和上限,用于估计的响应。注意,对于非墨镜的女性患者,随着重量增加,预期血压会降低,但另外一切都保持恒定。
栅格缺点(MDL)
此绘图显示主要效果。圆圈显示效果的大小,蓝线显示出主要效果的上下置信限度。例如,作为吸烟者将预期的血压增加10个单位,而不是作为一个非人的人,给予所有其他都保持不变。与女性相比,预期血压与女性相比,与女性相比增加了两个单位,给予其他预测因子保持恒定。25至50岁的年龄增加导致预期的4个单位增加,而111至202的重量变化导致预期血压的4单位减少,给予所有其他保持常数。
图PlotInteraction(MDL,'性别'那'重量'的)
该曲线显示给出另一个因素的变化的影响,另一个因素固定在一个值。
在解释互动效应时要谨慎。当所有因子组合或数据具有足够的数据时,数据高度相关时,可能难以确定改变一个因子的相互作用效果,同时保持另一个因素。在这种情况下,估计的相互作用效果是来自数据的推断。
蓝圆圈显示特定项的主要效果,如在主要效果图中。红色圆圈在一个术语中显示了一个术语的一个术语的影响。例如,在该图的下半部分,红色圆圈分别显示了女性和男性患者体重变化的影响。You can see that an increase in a female’s weight from 111 to 202 pounds causes about a 14-unit decrease in the expected blood pressure, while an increase of the same amount in the weight of a male patient causes about a 5-unit increase in the expected blood pressure, again given other predictors are held constant.
图PlotInteraction(MDL,'性别'那'重量'那'预测'的)
此图显示了随着其他预测变量变为变量的效果保持恒定。在该示例中,最后一个图显示了当在男性和女性处固定的变性时,响应变量,血压,作为重量的函数。男性和女性的线条是过境,表明体重和性别之间的强烈相互作用。当雄性患者的重量增加时,您可以看到预期的血压增加,但随着女性患者的重量增加而降低。
linearmodel.
|Fitlm.
|Stepwiselm.
|plotslice.
|丛林意见
|plotinteraction.