Я ищу график BLUPS со стандартными ошибками от объекта MCMCglmm (пакет MCMCglmm).Я могу извлечь "BLUPS", но не уверен, как извлечь стандартные ошибки.Ниже приведен код, используемый для извлечения BLUPS.Любой совет будет высоко ценится.
Value = colMeans(mcmc_1$Sol)) %>%
separate(Trait, c("Trait","Type","ID" ), sep = "\\.", fill = "right") %>%
filter(Type == "ID") %>%
filter(Trait %in% c("trait1", "trait2", "trait3")) %>%
select(-Type) %>%
spread(Trait, Value)
# A tibble: 72 x 4
ID trait1 trait2 trait3
<chr> <dbl> <dbl> <dbl>
1 1 -0.353 0.296 -0.122
2 10 0.0794 -0.0537 0.384
3 11 -0.289 0.186 -0.120
4 12 -0.0657 -0.175 0.558
5 13 -0.442 -0.133 0.0189
6 14 0.384 0.0833 0.0862
7 15 -0.793 0.0599 -0.265
8 16 -1.32 -0.313 -0.442
9 17 -0.329 -0.121 -0.604
10 18 1.20 0.356 -0.574
# ... with 62 more rows