... | ... | @@ -5,7 +5,7 @@ Fitting a Cholesky model (10 genetic and 10 residual parameters) to the data usi |
|
|
#Please downloaded externally
|
|
|
#load("G.large.RData")
|
|
|
#load("ph.large.RData")
|
|
|
#fit <- gsem.fit(ph.large, G.large, LogL = TRUE, estSE = TRUE)
|
|
|
#fit <- grmsem.fit(ph.large, G.large, LogL = TRUE, estSE = TRUE)
|
|
|
|
|
|
> load("fit.large.RData")
|
|
|
> fit.large
|
... | ... | @@ -97,7 +97,7 @@ attr(,"class") |
|
|
The variance can be obtained as
|
|
|
|
|
|
```{r eval = FALSE}
|
|
|
> var.fit <- gsem.var(fit)
|
|
|
> var.fit <- grmsem.var(fit)
|
|
|
> print(var.fit)
|
|
|
$VA
|
|
|
1 2 3 4
|
... | ... | |