Skip to main content

Table 2 Estimated variance components, heritability and correlations in Scenario 2 using the simulated data

From: Performance of the No-U-Turn sampler in multi-trait variance component estimation using genomic data

Parameter

True value

NUTS (LKJ prior)

NUTS (IW prior)

GS

REML

Mean

SE

Mean

SE

Mean

SE

Estimate

SE

Additive genetic (co)variances

 \({\sigma }_{a}^{2}(trait1)\)

0.98

1.00

0.36

1.79

0.20

1.29

0.35

1.06

0.44

 \({\sigma }_{a}^{2}(trait2)\)

4.89

5.32

1.48

5.37

1.40

5.84

1.46

5.54

1.49

 \({\sigma }_{a}(trait1, trait2)\)

0.59

0.57

0.29

0.72

0.35

0.86

0.46

0.75

0.47

Residual (co)variances

 \({\sigma }_{e}^{2}(trait1)\)

9.06

9.06

0.44

8.45

0.39

8.90

0.45

8.94

0.50

 \({\sigma }_{e}^{2}(trait2)\)

5.02

4.82

0.59

4.78

0.59

4.61

0.63

4.60

0.66

 \({\sigma }_{e}(trait1, trait2)\)

0.84

0.88

0.42

0.82

0.41

0.75

0.46

0.80

0.46

Heritabilities

 \({h}^{2}(trait1)\)

0.10

0.10

0.03

0.17

0.02

0.13

0.03

0.11

0.04

 \({h}^{2}(trait2)\)

0.50

0.52

0.08

0.52

0.08

0.55

0.08

0.54

0.09

Additive genetic correlations

 \({r}_{a}(trait1, trait2)\)

0.27

0.31

0.20

0.24

0.12

0.34

0.21

0.40

0.33

Residual correlations

 \({r}_{e}(trait1, trait2)\)

0.12

0.13

0.06

0.13

0.07

0.12

0.07

0.13

0.07

  1. NUTS No-U-Turn sampler, LKJ Lewandowski-Kurowicka-Joe, IW inverse-Wishart, GS Gibbs sampling, REML restricted maximum likelihood, SE standard error
  2. \({\sigma }_{a}^{2}\): additive genetic variance, \({\sigma }_{a}\): additive genetic covariance, \({\sigma }_{e}^{2}\): residual variance, \({\sigma }_{e}\): residual covariance, \({h}^{2}\): heritability, \({r}_{a}\): additive genetic correlation, \({r}_{e}\): residual correlation