[R] Summarizing dummy coefficients in sem package
Gus Jespersen
jesper at u.washington.edu
Tue Jun 3 23:13:51 CEST 2008
Greetings,
I am working in the sem package on a model with 3 exogenous variables (2
are nominal-categorical), and 4 endogenous, continuous variables. To
use sem with the nominal variables, I created dummy variables. Now, in
my sem output I have estimates for path coefficients for the
relationship between each level of the nominal variables and the
endogenous variables they are associated with. However, what I am
interested are the overall path coefficients from each nominal variable
to the various endogenous variables, not the relationship between each
level of this and the endogenous variables. How might I go about
achieving this?
Here is my model ('sex' and 'trt' are the nominal variables with 2 and 4
levels, respectively) followed by the relevant output:
Path Parameter
1 sexMale->nect sexnect
2 sexMale->cono sexcono
3 sexMale->herb sexherb
4 sexMale->para sexpara
5 trtControl->nect trtContnect
6 trtControl->cono trtContcono
7 trtControl->herb trtContherb
8 trtControl->para trtContpara
9 trtNoFemale->nect trtNoFemnect
10 trtNoFemale->cono trtNoFemcono
11 trtNoFemale->herb trtNoFemherb
12 trtNoFemale->para trtNoFempara
13 trtNoMale->nect trtNoMalenect
14 trtNoMale->cono trtNoMalecono
15 trtNoMale->herb trtNoMaleherb
16 trtNoMale->para trtNoMalepara
17 area->nect areanect
18 area->cono areacono
19 area->herb areaherb
20 area->para areapara
21 nect<->para nect<->para
22 nect<->herb nectherb
23 nect<->cono nectcono
24 cono<->herb conoherb
25 cono<->para conopara
26 area<->area errorarea
27 nect<->nect errornect
28 cono<->cono errorcono
29 herb<->herb errorherb
30 para<->para errorpara
31 trtControl<->trtControl errortrtControl
32 trtNoFemale<->trtNoFemale errorNoFemale
33 trtNoMale<->trtNoMale errortrtNoMale
34 sexMale<->sexMale errorsexMale
Parameter Estimates
Estimate Std Error z value Pr(>|z|)
sexnect -0.01929179 0.06376330 -0.3025532 7.6223e-01
sexcono -2.66789515 0.36084653 -7.3934344 1.4300e-13
sexherb -0.00743378 0.00566452 -1.3123412 1.8941e-01
sexpara 0.01097552 0.00980425 1.1194652 2.6294e-01
trtContnect 0.21164052 0.08736725 2.4224240 1.5417e-02
trtContcono -0.44791548 0.54528423 -0.8214349 4.1140e-01
trtContherb -0.00038585 0.00756348 -0.0510144 9.5931e-01
trtContpara -0.00015034 0.06098602 -0.0024651 9.9803e-01
trtNoFemnect 0.11658323 0.08786423 1.3268566 1.8456e-01
trtNoFemcono -0.10709039 0.55754523 -0.1920748 8.4768e-01
trtNoFemherb -0.00412156 0.00777418 -0.5301598 5.9600e-01
trtNoFempara 0.00420467 0.03406450 0.1234325 9.0176e-01
trtNoMalenect 0.23656338 0.08897775 2.6586801 7.8447e-03
trtNoMalecono 0.07633944 0.60103928 0.1270124 8.9893e-01
trtNoMaleherb -0.00156381 0.00788036 -0.1984436 8.4270e-01
trtNoMalepara 0.00910566 0.03304680 0.2755383 7.8290e-01
areanect 0.01097847 0.02250164 0.4878965 6.2562e-01
areacono 0.37431397 0.12725982 2.9413367 3.2680e-03
areaherb -0.00169721 0.00200841 -0.8450508 3.9808e-01
areapara 0.00044190 0.00365491 0.1209054 9.0377e-01
nect<->para 0.00399596 0.00425721 0.9386334 3.4792e-01
nectherb -0.00028131 0.00245373 -0.1146438 9.0873e-01
nectcono 0.52509632 0.15739193 3.3362340 8.4922e-04
conoherb 0.01380912 0.01401081 0.9856041 3.2433e-01
conopara 0.01752004 0.02413323 0.7259716 4.6786e-01
errorarea 1.90255983 0.08980542 21.1853566 0.0000e+00
errornect 0.82560204 0.03897074 21.1851759 0.0000e+00
errorcono 26.60775814 1.25597309 21.1849746 0.0000e+00
errorherb 0.00661525 0.00031226 21.1851940 0.0000e+00
errorpara 0.01967170 0.00092856 21.1852827 0.0000e+00
errortrtControl 0.19005279 0.00897095 21.1853556 0.0000e+00
errorNoFemale 0.18331925 0.00865314 21.1852797 0.0000e+00
errortrtNoMale 0.19059782 0.00899668 21.1853547 0.0000e+00
errorsexMale 0.24736468 0.01167625 21.1852846 0.0000e+00
Thanks,
Gus
More information about the R-help
mailing list