第九章 对数线性模型
第一节 General过程
主要功能
实例操作
第二节 Hierarchical过程
主要功能
实例操作
第三节 Logit过程
主要功能
实例操作
对数线性模型是用于离散型数据或整理成列联表格式的计数资料的统计分析工具。在对数线性模型中,所有用作的分类的因素均为独立变量,列联表各单元中的例数为应变量。对于列联表资料,通常作χ2 检验,但χ2 检验无法系统地评价变量间的联系,也无法估计变量间相互作用的大小,而对数线性模型是处理这些问题的最佳方法。
第一节 General过程
主要功能
调用该过程可对一个或多个二维列联表资料进行非层次对数线性分析。它只能拟合全饱和模型,即分类变量各自效应及其相互间效应均包含在对数线性模型中。
返回目录
返回全书目录
实例操作
[例9-1]在住院病人中,研究其受教育程度与对保健服务满意程度的关系,资料整理成列联表后如下所示。
对保健服务满意程度
(%)
受教育程度
高
中
低
满意
不满意
65 ()
6 ()
272 ()
18 ()
41 ()
1 ()
按一般情形作χ2检验,结果显示不同受教育程度的住院病人其对保健服务满意程度无差别。但从百分比分析中可见,随受教育程度的提高,满意程度有下降的趋势;且我们还想了解受教育程度与满意程度有无交互作用和交互作用的大小。对此,必须采用对数线性模型加以分析。
数据准备
激活数据管理窗口,定义变量名:实际观察频数的变量名为freq,受教育程度和满意程度作为行、列分类变量(即独立变量),变量名分别为educ、care。输入原始数据,结果如图所示。如同第四章Crosstab过程中所述,为使列联表的频数有效,应选Data菜单的Weight Cases...项,弹出Weight Cases对话框(图),激活Weight cases by项,从变量列表中选freq点击(钮使之进入Frequency Variable框,点击OK钮即可。
图 原始数据的输入
图 频数的加权定义
统计分析
激活Statistics菜单选Loglinear中的General...项,弹出General Loglinear Analysis对话框(图)。从对话框左侧的变量列表中选care,点击(钮使之进入Factor(s)框,点击Define Range...钮,弹出General Loglinear Analysis: Define Range对话框,定义分类变量care的范围,本例为1、2,故可在Minimum处键入1,在Maximum处键入2,点击Continue钮返回General Loglinear Analysis对话框。同法将变量educ选入Factor(s)框,并定义其范围为1、3。本例要求计算各分类变量主效应和交互作用的参数估计,故点击Contrast...钮,弹出General Loglinear Analysis:Contrasts对话框,选择Display parameter estimates项,点击Continue钮返回General Loglinear Analysis对话框,最后点击OK钮即完成分析。
图 非层次对数线性模型分析对话框
结果解释
在结果输出窗口中将看到如下统计数据:
首先显示系统对403例资料进行分析,共有二个分类变量:CARE为2水平,EDUC为3水平。分析的效应有三类:满意程度(CARE)、教育程度(EDUC)和两者的交互作用(CARE BY EDUC)。系统经2次叠代后即达到相邻二次估计之差不大于规定的。
DATA Information
6 unweighted cases accepted.
0 cases rejected because of out-of-range factor values.
0 cases rejected because of missing data.
403 weighted cases will be used in the analysis.
FACTOR Information
Factor Level Label
CARE 2
EDUC 3
DESIGN Information
1 Design/Model will be processed.
Correspondence Between Effects and Columns of Design/Model 1
Starting Ending
Column Column Effect Name
1 1 CARE
2 3 EDUC
4 5 CARE BY EDUC
Note: for saturated models .500 has been added to all observed cells.
This value may be changed by using the CRITERIA = DELTA subcommand.
*** ML converged at iteration 2.
Maximum difference between successive iterations = .00000
由于本例对Model(模型)未作定义,故系统采用默认的全饱和模型,因而期望例数()与实际例数(OBS. count)相同,进而残差(Residual)、标准化残差()和校正残差()均为0。
Observed, Expected Frequencies and Residuals
Factor Code OBS. count & PCT. EXP. count & PCT. Residual Std. Resid. Adj. Resid.
CARE 1
EDUC 1 () () .0000 .0000 .0000
EDUC 2 () () .0000 .0000 .0000
EDUC 3 () () .0000 .0000 .0000
CARE 2
EDUC 1 ( ) ( ) .0000 .0000 .0000
EDUC 2 ( ) ( ) .0000 .0000 .0000
EDUC 3 ( .37) ( .37) .0000 .0000 .0000
最后输出参数估计的结果。为了唯一地估计参数,系统强行限定同一分类变量的各水平参数之和为0,故根据下列结果可推得各参数为:
λ满意 =
λ不满意 =
λ高教育程度 =
λ中教育程度 =
λ低教育程度 =
λ满意.高教育程度 =
λ满意.中高教育程度 =
λ满意.低教育程度 =
λ不满意.高教育程度 =
λ不满意.中教育程度 =
λ不满意.低教育程度 =
λ值为正,表示正效应;反之为负效应;零为无效应。分析提供的信息是:①对保健服务的满意程度高于不满意程度;②中等教育程度者的满意程度>高等教育程度者的满意程度>低等教育程度者的满意程度;③通过受教育程度与对保健服务满意程度的交互作用研究,结果表明高、中等教育未能增加人们对现有保健服务状况的满意程度。
Estimates for Parameters
CARE
Parameter Coeff. Std. Err. Z-Value Lower 95 CI Upper 95 CI
1 .15965
EDUC
Parameter Coeff. Std. Err. Z-Value Lower 95 CI Upper 95 CI
2 .19895 .29847
3 .17407 .80313
CARE BY EDUC
Parameter Coeff. Std. Err. Z-Value Lower 95 CI Upper 95 CI
4 .19895 .15834
5 .17407 .29938
返回目录
返回全书目录
第二节 Hierarchical过程
主要功能
调用该过程可对多维列联表资料进行分层对数线性分析。所谓分层即并可根据用户指定的条件,对某一或某些主效应与交互作用进行剔除,从而形成包含特定层次阶项的各种模型。
返回目录
返回全书目录
实例操作
[例9-2] 为了研究Colles骨折在不同性别中的年龄分布情况,以说明不同性别者骨折的年龄差异及其年度变化,某地收集了1978--1981年的骨折资料,数据见下表。请作对数线性模型的分析。
年龄
1978
1979
1980
1981
男
女
男
女
男
女
男
女
0—19
20--59
60--89
55
165
50
17
260
94
43
101
29
9
233
115
89
104
56
20
202
95
140
137
54
41
278
153
数据准备
激活数据管理窗口,定义变量名:实际观察频数的变量名为freq,年份、性别和年龄为分类变量,变量名分别为year、sex和age。输入原始数据,其中年份1978至1981依次为1、2、3、4,性别男为1、女为2,年龄分组依次为1、2、3。之后选Data菜单的Weight Cases...项,在Weight Cases对话框中激活Weight cases by项,从变量列表中选freq点击(钮使之进入Frequency Variable框,点击OK钮完成对频数的权重定义。
统计分析
激活Statistics菜单选Loglinear中的Hierarchical...项,弹出Hierarchical Loglinear Analysis对话框(图)。从对话框左侧的变量列表中选age,点击(钮使之进入Factor(s)框,点击Define Range...钮,弹出Hierarchical Loglinear Analysis: Define Range对话框,定义分类变量age的范围,在Minimum处键入1,在Maximum处键入9,点击Continue钮返回Hierarchical Loglinear Analysis对话框。同法将变量sex选入Factor(s)框,定义其范围为1、2;将变量year选入Factor(s)框,定义其范围为1、4。
图 层次对数线性模型分析对话框
为了更好地拟合数据,并尽可能的简单和易于解释,本例选择向后剔除法建立模型,即从所有效应均在模型中开始,然后消除那些不满足保留判据的效应。点击Model...钮,弹出Hierarchical Loglinear Analysis: Model对话框,在Model Building栏中选Use backward elimination项,点击Continue钮返回Hierarchical Loglinear Analysis对话框。
本例要求作参数估计,故点击Options...钮,弹出Hierarchical Loglinear Analysis: Options对话框,在Display for Saturated Model栏中选Parameter estimates项,点击Continue钮返回Hierarchical Loglinear Analysis对话框,之后点击OK钮即完成分析。
结果解释
在结果输出窗口中将看到如下统计数据:
首先显示,共有2540个观察例数进入分析,其中分类变量AGE为3水平,SEX为2水平,YEAR为4水平。采用全饱和模型,高阶项为年龄、性别和年份三者的交互作用。(在层次对数线性模型分析中,当指定高阶项时,即意味着包含其所属变量所有可能组合的低阶项;如本例,即包含年龄和性别的交互作用、年龄和年份的交互作用、性别和年份的交互作用、年龄的主效应、性别的主效应、年份的主效应。从最高阶到最低阶共为3阶。)
DATA Information
24 unweighted cases accepted.
0 cases rejected because of out-of-range factor values.
3 cases rejected because of missing data.
2540 weighted cases will be used in the analysis.
FACTOR Information
Factor Level Label
AGE 3
SEX 2
YEAR 4
DESIGN 1 has generating class
AGE*SEX*YEAR
Note: For saturated models .500 has been added to all observed cells.
This value may be changed by using the CRITERIA = DELTA subcommand.
The Iterative Proportional Fit algorithm converged at iteration 1.
The maximum difference between observed and fitted marginal totals is .000
and the convergence criterion is .278
系统以全饱和模型为起始,故显示各变量的实际例数、期望例数、残差和标准化残差,因期望例数与实际例数相同,进而残差、标准化残差均为0。
Observed, Expected Frequencies and Residuals.
Factor Code OBS count EXP count Residual Std Resid
AGE 1
SEX 1
YEAR 1 .00 .00
YEAR 2 .00 .00
YEAR 3 .00 .00
YEAR 4 .00 .00
SEX 2
YEAR 1 .00 .00
YEAR 2 .00 .00
YEAR 3 .00 .00
YEAR 4 .00 .00
AGE 2
SEX 1
YEAR 1 .00 .00
YEAR 2 .00 .00
YEAR 3 .00 .00
YEAR 4 .00 .00
SEX 2
YEAR 1 .00 .00
YEAR 2 .00 .00
YEAR 3 .00 .00
YEAR 4 .00 .00
AGE 3
SEX 1
YEAR 1 .00 .00
YEAR 2 .00 .00
YEAR 3 .00 .00
YEAR 4 .00 .00
SEX 2
YEAR 1 .00 .00
YEAR 2 .00 .00
YEAR 3 .00 .00
YEAR 4 .00 .00
Goodness-of-fit test statistics
Likelihood ratio chi square = .00000 DF = 0 P =
Pearson chi square = .00000 DF = 0 P =
下面,系统先显示某一阶及其更高阶交互效应为0时的似然比χ2检验概率值,因K为3时的概率值=>,故认为年龄、性别、年份三者的交互作用为0,亦即含1阶(单一变量主效应)及2阶(变量两两交互效应)的模型就能恰当地表述数据。
接着,系统又显示特定阶交互效应为0时的似然比χ2检验概率值,结果表明,单纯含1阶(单一变量主效应)或单纯含2阶(变量两两交互效应)的模型也能恰当地表述数据。
Tests that K-way and higher order effects are zero.
K DF . Chisq Prob Pearson Chisq Prob Iteration
3 6 .1964 .2007 4
2 17 .0000 .0000 2
1 23 .0000 .0000 0
Tests that K-way effects are zero.
K DF . Chisq Prob Pearson Chisq Prob Iteration
1 6 .0000 .0000 0
2 11 .0000 .0000 0
3 6 .1964 .2007 0
Note: For saturated models .500 has been added to all observed cells.
This value may be changed by using the CRITERIA = DELTA subcommand.
系统所确定的模型中各参数值如下所示,由于内容较多,各λ值如何推算及其所表示的意义,请读者参阅本章第一节。
Estimates for Parameters.
AGE*SEX*YEAR
Parameter Coeff. Std. Err. Z-Value Lower 95 CI Upper 95 CI
1 .08417 .02375
2 .1674922915 .10130 .36605
3 .07921 .13826
4 .0577506145 .05557 .16667
5 .06504 .12057
6 .05570 .02740
AGE*SEX
Parameter Coeff. Std. Err. Z-Value Lower 95 CI Upper 95 CI
1 .7059980126 .04848 .61098 .80102
2 .03276
AGE*YEAR
Parameter Coeff. Std. Err. Z-Value Lower 95 CI Upper 95 CI
1 .08417
2 .10130
3 .1339590237 .07921 .28920
4 .1990874838 .05557 .09017 .30800
5 .1982170140 .06504 .07073 .32570
6 .05570
SEX*YEAR
Parameter Coeff. Std. Err. Z-Value Lower 95 CI Upper 95 CI
1 .0471962901 .04918 .95960 .14360
2 .05818 .03615
3 .0827715134 .04734 .17556
AGE
Parameter Coeff. Std. Err. Z-Value Lower 95 CI Upper 95 CI
1 .04848
2 .7999110228 .03276 .73571 .86412
SEX
Parameter Coeff. Std. Err. Z-Value Lower 95 CI Upper 95 CI
1 .02856 .02111
YEAR
Parameter Coeff. Std. Err. Z-Value Lower 95 CI Upper 95 CI
1 .04918 .07588
2 .05818
3 .04734 .08014
系统开始对全饱和模型进行从高阶到低阶的效应项剔除。第一步,剔除3阶交互效应项(AGE*SEX*YEAR)导致χ2值为,概率为(不小于默认判据),故该效应项被剔除。
第二步,剔除2阶交互效应项,概率均小于,故2阶交互效应项不能剔除。即本例用2阶交互效应项(同时含1阶主效应项)描述模型已为最佳。
Backward Elimination (p = .050) for DESIGN 1 with generating class
AGE*SEX*YEAR
Likelihood ratio chi square = .00000 DF = 0 P =
If Deleted Simple Effect is DF . Chisq Change Prob Iter
AGE*SEX*YEAR 6 .1964 4
Step 1
The best model has generating class
AGE*SEX
AGE*YEAR
SEX*YEAR
Likelihood ratio chi square = DF = 6 P = .196
If Deleted Simple Effect is DF . Chisq Change Prob Iter
AGE*SEX 2 .0000 2
AGE*YEAR 6 .0000 2
SEX*YEAR 3 .0046 2
Step 2
The best model has generating class
AGE*SEX
AGE*YEAR
SEX*YEAR
Likelihood ratio chi square = DF = 6 P = .196
The final model has generating class
AGE*SEX
AGE*YEAR
SEX*YEAR
The Iterative Proportional Fit algorithm converged at iteration 0.
The maximum difference between observed and fitted marginal totals is .131
and the convergence criterion is .278
由于剔除了3阶交互效应项,故原全饱和模型变为层次模型,因而期望例数改变,期望例数与实际例数不同,进而残差、标准化残差均不为0。若标准化残差界于—范围内,则表示模型是恰当的。从下面的结果可知,本例的标准化残差均在—范围内,故层次模型是适合的。
Observed, Expected Frequencies and Residuals.
Factor Code OBS count EXP count Residual Std Resid
AGE 1
SEX 1
YEAR 1
YEAR 2 .62
YEAR 3 .69 .07
YEAR 4
SEX 2
YEAR 1
YEAR 2
YEAR 3
YEAR 4 .53 .08
AGE 2
SEX 1
YEAR 1 .16
YEAR 2 .31
YEAR 3
YEAR 4 .31
SEX 2
YEAR 1
YEAR 2
YEAR 3 .62
YEAR 4
AGE 3
SEX 1
YEAR 1 .30
YEAR 2
YEAR 3
YEAR 4
SEX 2
YEAR 1
YEAR 2 .67
YEAR 3
YEAR 4 .25
Goodness-of-fit test statistics
Likelihood ratio chi square = DF = 6 P = .196
Pearson chi square = DF = 6 P = .201
返回目录
返回全书目录
第三节 Logit过程
主要功能
调用此过程可完成对一个应变量与一个或多个自变量之间对数线性模型的拟合。如果分类变量未区分应变量和自变量,那么应采用本章第一、二节介绍的方法;如果应变量是二分计量,自变量是连续计量,那么应采用Logistic回归方法(详见第八章)。
返回目录
返回全书目录
实例操作
[例]在艾滋病(AIDS)相关的知识、观念、行为研究(KAB Study)中,获得了不同年龄和受教育水平的公众,对预防AIDS知识掌握程度的资料,经整理成列联表如下所示。很明显,对预防AIDS知识的掌握程度与公众的年龄和受教育水平有关,即若预防AIDS知识掌握程度为应变量,则应该受到年龄和受教育水平两个自变量的影响。下面将运用带应变量的对数线性模型进行分析。
受教育水平
年龄
预防AIDS知识掌握程度
好 一般 差
高
20-
30-
40-
50-
53
28
31
19
40
21
32
6
2
3
8
11
中
20-
30-
40-
50-
67
71
38
9
103
141
94
66
24
101
87
136
低
20-
30-
40-
50-
2
16
8
3
3
22
98
76
17
19
247
156
数据准备
激活数据管理窗口,定义变量名:实际观察频数的变量名为freq;预防AIDS知识掌握程度变量名为aids,按好、一般、差分别输入1、2、3;受教育水平变量名为educ,按高、中、低分别输入1、2、3;年龄变量名为age,20-至50-依次输入1—4。输入原始数据后选Data菜单的Weight Cases...项,在Weight Cases对话框中激活Weight cases by项,从变量列表中选freq点击(钮使之进入Frequency Variable框,点击OK钮即可。
统计分析
激活Statistics菜单选Loglinear中的Logit...项,弹出Logit Loglinear Analysis对话框(图)。从对话框左侧的变量列表中选aids,点击(钮使之进入Dependent框,点击Define Range...钮,弹出Logit Loglinear Analysis: Define Range对话框,定义应变量aids的范围,在Minimum处键入1,在Maximum处键入3,点击Continue钮返回Logit Loglinear Analysis对话框。从对话框左侧的变量列表中选age,点击(钮使之进入Factor(s)框,点击Define Range...钮,定义自变量age的范围为1、4;同法将自变量educ选入Factor(s)框,并定义其范围为1、3。本例要求计算各变量主效应和交互作用的参数估计,故点击Contrast...钮,弹出Logit Loglinear Analysis:Contrasts对话框,选择Display parameter estimates项,点击Continue钮返回Logit Loglinear Analysis对话框,最后点击OK钮即完成分析。
图 Logit对数线性模型对话框
结果解释
在结果输出窗口中将看到如下统计数据:
系统显示共有1858个观察例数进入分析,分析涉及三个变量,其中AIDS为3水平,AGE为4水平,EDUC为3水平。将产生3阶4类效应,即:预防AIDS知识掌握程度主效应(因AIDS被定义为应变量,故不再分析子变量AGE、EDUC的主效应),预防AIDS知识掌握程度分别与年龄、受教育程度的交互效应,预防AIDS知识掌握程度、年龄、受教育程度三者的交互效应。之后系统显示实际例数、期望例数、残差、标准化残差和校正残差。
DATA Information
36 unweighted cases accepted.
0 cases rejected because of out-of-range factor values.
0 cases rejected because of missing data.
1858 weighted cases will be used in the analysis.
FACTOR Information
Factor Level Label
AIDS 3
AGE 4
EDUC 3
DESIGN Information
1 Design/Model will be processed.
Correspondence Between Effects and Columns of Design/Model 1
Starting Ending
Column Column Effect Name
1 2 AIDS
3 8 AIDS BY AGE
9 12 AIDS BY EDUC
13 24 AIDS BY AGE BY EDUC
Note: for saturated models .500 has been added to all observed cells.
This value may be changed by using the CRITERIA = DELTA subcommand.
*** ML converged at iteration 2.
Maximum difference between successive iterations = .00000.
Observed, Expected Frequencies and Residuals
Factor Code OBS. count & PCT. EXP. count & PCT. Residual Std. Resid. Adj. Resid.
AIDS 1
AGE 1
EDUC 1 () () .0000 .0000 .0000
EDUC 2 () () .0000 .0000 .0000
EDUC 3 () () .0000 .0000 .0000
AGE 2
EDUC 1 () () .0000 .0000 .0000
EDUC 2 () () .0000 .0000 .0000
EDUC 3 () () .0000 .0000 .0000
AGE 3
EDUC 1 () () .0000 .0000 .0000
EDUC 2 () () .0000 .0000 .0000
EDUC 3 ( ) ( ) .0000 .0000 .0000
AGE 4
EDUC 1 () () .0000 .0000 .0000
EDUC 2 ( ) ( ) .0000 .0000 .0000
EDUC 3 ( ) ( ) .0000 .0000 .0000
AIDS 2
AGE 1
EDUC 1 () () .0000 .0000 .0000
EDUC 2 () () .0000 .0000 .0000
EDUC 3 () () .0000 .0000 .0000
AGE 2
EDUC 1 () () .0000 .0000 .0000
EDUC 2 () () .0000 .0000 .0000
EDUC 3 () () .0000 .0000 .0000
AGE 3
EDUC 1 () () .0000 .0000 .0000
EDUC 2 () () .0000 .0000 .0000
EDUC 3 () () .0000 .0000 .0000
AGE 4
EDUC 1 () () .0000 .0000 .0000
EDUC 2 () () .0000 .0000 .0000
EDUC 3 () () .0000 .0000 .0000
AIDS 3
AGE 1
EDUC 1 ( ) ( ) .0000 .0000 .0000
EDUC 2 () () .0000 .0000 .0000
EDUC 3 () () .0000 .0000 .0000
AGE 2
EDUC 1 ( ) ( ) .0000 .0000 .0000
EDUC 2 () () .0000 .0000 .0000
EDUC 3 () () .0000 .0000 .0000
AGE 3
EDUC 1 () () .0000 .0000 .0000
EDUC 2 () () .0000 .0000 .0000
EDUC 3 () () .0000 .0000 .0000
AGE 4
EDUC 1 () () .0000 .0000 .0000
EDUC 2 () () .0000 .0000 .0000
EDUC 3 () () .0000 .0000 .0000
Goodness-of-Fit test statistics
Likelihood Ratio Chi Square = .00000 DF = 0 P =
Pearson Chi Square = .00000 DF = 0 P =
下一段为拟合优度的检验。系统采用分散相似测量法(Dispersion Similarity Measure),测量值界于-1至+1之间,愈靠近| 1 |,拟合优度愈好。本例为。
Analysis of Dispersion
Dispersion
Source of Variation Entropy Concentration DF
Due to Model
Due to Residual
Total 3750
Measures of Association
Entropy = .160867
Concentration = .145879
最后,系统输出对数线性模型的各效应参数值。
由于内容较多,具体推算过程不再赘述(参阅本章第一节)。此处以AIDS主效应和AIDS与EDUC交互效应为例,演示如下:
预防AIDS知识掌握程度主效应部分,参数为
λAIDS-好 =
λAIDS-一般 =
λAIDS-差 = 0- () - =
这表明公众预防AIDS知识掌握程度一般。
预防AIDS知识掌握程度与受教育水平交互效应部分,参数为
λAIDS好-EDUC高 =
λAIDS好-EDUC中 =
λAIDS好-EDUC低 = 0 - - () =
λAIDS一般-EDUC高 =
λAIDS一般-EDUC中 =
λAIDS一般-EDUC低 = 0 - () - =
λAIDS差-EDUC高 = 0 - - () =
λAIDS差-EDUC中 = 0 - () - =
λAIDS差-EDUC低 = 0 - - =
这表明受教育水平高,预防AIDS知识掌握程度好;受教育水平低,预防AIDS知识掌握程度一般。为什么不体现受教育水平低,预防AIDS知识掌握程度差的信息呢?显然,这还需要结合年龄的因素进行分析。若用户将全部λ值都推算出来,其中会得到:
λAIDS差-AGE20-EDUC低 =
λAIDS差-AGE30-EDUC低 =
λAIDS差-AGE40-EDUC低 =
λAIDS差-AGE50-EDUC低 =
其趋势大约是年龄大的、受教育水平低的,预防AIDS知识掌握程度就较差。
Estimates for Parameters
AIDS
Parameter Coeff. Std. Err. Z-Value Lower 95 CI Upper 95 CI
1 .07013
2 .3307195684 .06115 .21087 .45057
AIDS BY AGE
Parameter Coeff. Std. Err. Z-Value Lower 95 CI Upper 95 CI
3 .5610569048 .14377 .27926 .84286
4 .4747703448 .10184 .27516 .67438
5 .10192
6 .0139480312 .13188 .10577 .27243
7 .0027846286 .09475 .02939 .18849
8 .0999432599 .08316 .26293
AIDS BY EDUC
Parameter Coeff. Std. Err. Z-Value Lower 95 CI Upper 95 CI
9 .09401 .91282
10 .08160
11 .09212 .16177
12 .0930200827 .06828 .22685
AIDS BY AGE BY EDUC
Parameter Coeff. Std. Err. Z-Value Lower 95 CI Upper 95 CI
13 .18396 .19458
14 .1990147721 .15635 .50545
15 .15180
16 .11603
17 .0248775102 .13931 .17858 .29793
18 .3092579898 .12025 .07357 .54494
19 .5096508480 .17810 .16057 .85873
20 .1850931544 .14164 .46270
21 .1964166680 .15202 .49437
22 .10424 .11606
23 .0455872550 .12925 .35269 .29893
24 .09651
1