从R中的MCMCglmm绘制具有标准错误的BLUPS

Question

我希望从MCMCglmm对象（MCMCglmm软件包）绘制带有标准错误的BLUPS。我可以提取“ 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