|
|
The factorial co-heritabilities can be estimated as shown below (shortened output). Note that for the first genetic factor in a Cholesky model, the factorial co-heritability is always **one**, as per definition, this factor explains the entire genetic variation of the trait. Note that `Vi` is the estimated explained trait variance by a factor and corresponds to the squared path coefficient.
|
|
|
|
|
|
```{r eval=FALSE}
|
|
|
> gsem.fcoher(fit)
|
|
|
$fcoher.model.out
|
|
|
|
|
|
#fcoher.model.out[,c(1,7:10)]
|
|
|
label Vi Vi.se FCOHER FCOHER.se
|
|
|
1 a11 3.484351e-01 0.0223563431 1.0000000000 0.000000000
|
|
|
2 a21 1.650216e-01 0.0181433759 0.3059306238 0.027526691
|
|
|
3 a31 1.699770e-01 0.0179853613 0.3167175772 0.025135935
|
|
|
4 a41 6.309215e-02 0.0127378513 0.0957606594 0.018914180
|
|
|
5 a22 3.743871e-01 0.0160084532 0.6940693762 0.027526691
|
|
|
6 a32 1.701516e-01 0.0125017124 0.3170428312 0.024459178
|
|
|
7 a42 1.652393e-01 0.0164616136 0.2507985687 0.023183443
|
|
|
8 a33 1.965547e-01 0.0120137852 0.3662395917 0.021786386
|
|
|
9 a43 2.455148e-04 0.0007739957 0.0003726401 0.001174465
|
|
|
10 a44 4.302756e-01 0.0194021010 0.6530681318 0.023536859
|
|
|
11 e11 6.308243e-01 0.0198780916 NA NA
|
|
|
12 e21 8.424458e-02 0.0075712747 NA NA
|
|
|
13 e31 1.982422e-01 0.0108236256 NA NA
|
|
|
14 e41 4.515762e-05 0.0001572113 NA NA
|
|
|
15 e22 2.715217e-01 0.0092143864 NA NA
|
|
|
16 e32 2.982679e-03 0.0010945112 NA NA
|
|
|
17 e42 3.323062e-04 0.0004363328 NA NA
|
|
|
18 e33 1.645033e-01 0.0060606482 NA NA
|
|
|
19 e43 5.184572e-04 0.0005643313 NA NA
|
|
|
20 e44 2.182086e-01 0.0084150015 NA NA
|
|
|
|
|
|
```
|
|
|
|
|
|
Similarly, the factorial co-environmentalities is shown below (shortened output). Again, note that for the first residual factor in a Cholesky model, the factorial co-environmentality is always **one**, as per definition this factor explains the entire residual variation of the trait.
|
|
|
|
|
|
```{r eval=FALSE}
|
|
|
> gsem.fcoher(fit)
|
|
|
$fcoher.model.out
|
|
|
|
|
|
#fcoher.model.out[,c(1,7:8,13:14)]
|
|
|
label Vi Vi.se FCOENV FCOENV.se
|
|
|
1 a11 3.484351e-01 0.0223563431 NA NA
|
|
|
2 a21 1.650216e-01 0.0181433759 NA NA
|
|
|
3 a31 1.699770e-01 0.0179853613 NA NA
|
|
|
4 a41 6.309215e-02 0.0127378513 NA NA
|
|
|
5 a22 3.743871e-01 0.0160084532 NA NA
|
|
|
6 a32 1.701516e-01 0.0125017124 NA NA
|
|
|
7 a42 1.652393e-01 0.0164616136 NA NA
|
|
|
8 a33 1.965547e-01 0.0120137852 NA NA
|
|
|
9 a43 2.455148e-04 0.0007739957 NA NA
|
|
|
10 a44 4.302756e-01 0.0194021010 NA NA
|
|
|
11 e11 6.308243e-01 0.0198780916 1.0000000000 0.0000000000
|
|
|
12 e21 8.424458e-02 0.0075712747 0.2367975722 0.0172339750
|
|
|
13 e31 1.982422e-01 0.0108236256 0.5420479497 0.0161908960
|
|
|
14 e41 4.515762e-05 0.0001572113 0.0002061008 0.0007175779
|
|
|
15 e22 2.715217e-01 0.0092143864 0.7632024278 0.0172339750
|
|
|
16 e32 2.982679e-03 0.0010945112 0.0081554556 0.0029995986
|
|
|
17 e42 3.323062e-04 0.0004363328 0.0015166558 0.0019944184
|
|
|
18 e33 1.645033e-01 0.0060606482 0.4497965946 0.0162294819
|
|
|
19 e43 5.184572e-04 0.0005643313 0.0023662548 0.0025782116
|
|
|
20 e44 2.182086e-01 0.0084150015 0.9959109886 0.0034506100
|
|
|
``` |
|
|
\ No newline at end of file |