Skip to content

Commit 50dd340

Browse files
authored
Update README.md
1 parent 3871375 commit 50dd340

File tree

1 file changed

+4
-4
lines changed

1 file changed

+4
-4
lines changed

README.md

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -28,12 +28,12 @@ Source: [inst/examples/package-lme4qtl.R](inst/examples/package-lme4qtl.R)
2828
library(lme4) # needed for `VarCorr` function
2929
library(lme4qtl)
3030
31-
### load the synthetic data:
31+
# load synthetic data set `dat40` distributed within `lme4qtl`
3232
# - table of phenotypes `dat40`
3333
# - the double kinship matrix `kin2`
3434
data(dat40)
3535
36-
### model the continiuous trait `trait1`
36+
# (1) model continiuous trait `trait1`
3737
mod <- relmatLmer(trait1 ~ AGE + SEX + (1|FAMID) + (1|ID), dat40, relmat = list(ID = kin2))
3838
3939
# get the estimation of h2
@@ -48,8 +48,8 @@ prop <- with(vf, vcov / sum(vcov))
4848
(h2 <- prop[1])
4949
#[1] `0.895419`
5050
51-
### model the binary trait `trait1bin`
52-
# - Model is nearly unidentifiable, when `(1|FAMID)` effect is added
51+
# (2) model binary trait `trait1bin`
52+
# - model is nearly unidentifiable, when `(1|FAMID)` effect is added
5353
gmod <- relmatGlmer(trait1bin ~ AGE + SEX + (1|ID), dat40, relmat = list(ID = kin2), family = binomial)
5454
```
5555

0 commit comments

Comments
 (0)