asreml软件报错:Likelihood evaluation failed with fault 333 ; trying with reduced updates

Likelihood evaluation failed with fault 333 ; trying with reduced updates

asreml软件报错:

Likelihood evaluation failed with fault 333 ; trying with reduced updates
Error in asreml(trait1 ~ Block/Rep, random = ~vm(ID, ainv), data = nested) :
ASReml failed to build MME (XXG); try increasing ‘workspace’.

十有八九:
看一下ainv有没有数值:

> head(ainv)
     Row Column Ainverse

说明构建逆矩阵出错。

狗血的解决方案:
重启R语言,在运行一下。

怀疑的bug:
ainverse间歇性神经!!!

另外的解决方式:

运行下面命令,查看ainv有无结果,如果有,就说嘛ainverse能用了。

library(asreml)
data("harvey")

ped = harvey.ped
head(ped)
str(ped)

ainv = ainverse(ped)
head(ainv)
str(ainv)

sessionInfo()

显示信息:

> library(asreml)
> data("harvey")
> 
> ped = harvey.ped
> head(ped)
  Calf   Sire Dam
1  101 Sire_1   0
2  102 Sire_1   0
3  103 Sire_1   0
4  104 Sire_1   0
5  105 Sire_1   0
6  106 Sire_1   0
> str(ped)
'data.frame':	65 obs. of  3 variables:
 $ Calf: int  101 102 103 104 105 106 107 108 109 110 ...
 $ Sire: chr  "Sire_1" "Sire_1" "Sire_1" "Sire_1" ...
 $ Dam : int  0 0 0 0 0 0 0 0 0 0 ...
> 
> ainv = ainverse(ped)
> head(ainv)
     Row Column Ainverse
[1,]   1      1 3.666667
[2,]   2      2 3.666667
[3,]   3      3 2.666667
[4,]   4      4 3.666667
[5,]   5      5 3.333333
[6,]   6      6 3.000000
> str(ainv)
 'matrix' num [1:139, 1:3] 1 2 3 4 5 6 7 8 9 10 ...
 - attr(*, "inbreeding")= Named num [1:74] 0 0 0 0 0 ...
  ..- attr(*, "names")= chr [1:74] "Sire_1" "Sire_2" "Sire_3" "Sire_4" ...
 - attr(*, "logdet")= num -18.7
 - attr(*, "dimnames")=List of 2
  ..$ : NULL
  ..$ : chr [1:3] "Row" "Column" "Ainverse"
 - attr(*, "rowNames")= chr [1:74] "Sire_1" "Sire_2" "Sire_3" "Sire_4" ...
 - attr(*, "geneticGroups")= num [1:2] 0 0
> 
> sessionInfo()
R version 4.0.3 (2020-10-10)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 10 x64 (build 18363)

Matrix products: default

locale:
[1] LC_COLLATE=Chinese (Simplified)_China.936  LC_CTYPE=Chinese (Simplified)_China.936   
[3] LC_MONETARY=Chinese (Simplified)_China.936 LC_NUMERIC=C                              
[5] LC_TIME=Chinese (Simplified)_China.936    

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] asreml_4.1.0.143 Matrix_1.2-18   

loaded via a namespace (and not attached):
 [1] lattice_0.20-41   dplyr_1.0.2       crayon_1.3.4      grid_4.0.3        R6_2.5.0         
 [6] jsonlite_1.7.1    lifecycle_0.2.0   gtable_0.3.0      magrittr_2.0.1    scales_1.1.1     
[11] pillar_1.4.7      ggplot2_3.3.2     rlang_0.4.9       rstudioapi_0.13   data.table_1.13.2
[16] generics_0.1.0    vctrs_0.3.5       ellipsis_0.3.1    tools_4.0.3       glue_1.4.2       
[21] purrr_0.3.4       munsell_0.5.0     compiler_4.0.3    pkgconfig_2.0.3   colorspace_2.0-0 
[26] tidyselect_1.1.0  tibble_3.0.4  
上一篇:最大似然估计和最大后验概率估计(MLE&MAP)


下一篇:最大似然估计