Построение BLUPS со стандартными ошибками от MCMCglmm в R - PullRequest
0 голосов
/ 28 февраля 2019

Я ищу график 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

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...