|
|
Fitting a Cholesky model (10 genetic and 10 residual parameters) to the data using 4 cores, sharing memory across 24 cores, required 34h using R MKL 3.6.3. Therefore, the model was pre-fitted and saved (`fit.large.RData`) to illustrate `grmsem` follow-up functions. Note that some output has been omitted for clarity.
|
|
|
|
|
|
```{r eval = FALSE}
|
|
|
#Do not run:
|
|
|
#load("G.large")
|
|
|
#load("ph.large")
|
|
|
#Do not run!
|
|
|
#Please downloaded externally
|
|
|
#load("G.large.RData")
|
|
|
#load("ph.large.RData")
|
|
|
#fit <- gsem.fit(ph.large, G.large, LogL = TRUE, estSE = TRUE)
|
|
|
|
|
|
> load("fit.large")
|
|
|
> fit
|
|
|
> load("fit.large.RData")
|
|
|
> fit.large
|
|
|
$model.in
|
|
|
part label value freepar
|
|
|
1 a a11 0.01686961 1
|
... | ... | |