AFEchidna包的简单示例4--多变量的批量分析

简单示例如下：

library(AFEchidna)

## new syntax
res22<-echidna(es0.file='fm.es0',trait=~h3+h4+h5,
                   fixed=~Trait+Trait:Rep,
                   random=~us(Trait):Mum,
                   residual=~units:us(Trait),
                   predict='Mum',
                   batch=TRUE,
                   mulT=TRUE,mulN=2)

## old syntax
trait=c('h3','h4','h5')
es0.path=r'(E:\课件\新建文件夹)'
res22<-echidna.batch(trait=trait,fixed='Trait.Rep',
                   random='us(Trait).Mum',
                   residual='units.us(Trait)',
                   predict='Mum',
                   es0.path=es0.path,mulT=T)

运行过程：

> ## new syntax
> res22<-echidna(es0.file='fm.es0',trait=trait,fixed='Trait.Rep',
+                    random='us(Trait).Fam',
+                    residual='units.us(Trait)',
+                    predict='Mum',
+                    batch=T,
+                    mulT=T,mulN=2)

Program starts running batch analysis ------ 

run 1 -- -- --:
Running Echidna for analysis:  h3 h4 

 Mon Feb 22 14:53:43 2021 
  Iteration     LogL NEDF
1         1 -5928.77 1105
2         2 -5612.69 1105
3         3 -5584.27 1105
4         4 -5548.36 1105
5         5 -5552.60 1105
6         6 -5552.97 1105
7         7 -5552.98 1105
Mon Feb 22 14:53:43 2021 LogL Converged 

run 2 -- -- --:
Running Echidna for analysis:  h3 h5 

 Mon Feb 22 14:53:44 2021 
  Iteration     LogL NEDF
1         1 -6121.71 1105
2         2 -6002.96 1105
3         3 -5875.08 1105
4         4 -5857.89 1105
5         5 -5859.26 1105
6         6 -5859.28 1105
7         7 -5859.28 1105
Mon Feb 22 14:53:44 2021 LogL Converged 

run 3 -- -- --:
Running Echidna for analysis:  h4 h5 

 Mon Feb 22 14:53