AAfun4包的新函数

主要针对原有的asreml.batch()，增加了两个函数：repin和update，用于计算新的遗传参数，或者重新运行新的批量分析。最新版的asreml.batch()，基本上能涵盖常见所有的可用模型，比如家系模型，个体模型，空间模型，基因组选择模型，及其多变量模型。

repin用于计算批量分析结果的新参数，无需重新调用Asrem程序。

update用于重新运行批量分析，只需输入相应的新模型项，比如固定效应、随机效应和误差效应等，无需重复输入原有的其它参数。同时，运行的新结果，可直接传递给repin计算遗传参数。

1 repin

用法：

repin(object, mulp, asrV = 4)
     -- object：批量分析运行的结果；
     --mulp：新的参数计算公式；
     --asrV：asreml版本，默认为最新版4.1。

简单示范如下：

// example
pkgs <- c('asreml','AAfun4')
lapply(pkgs,library,character.only=TRUE)
data(PrSpa)
df1<-subset(PrSpa,Spacing==3)

mulp1<-c(h2 ~ 4 * V1/(V1+V2))
run1<-asreml.batch(data=df1,factorN=1:5,traitN=c(9:13),
             FMod=y~1+Rep+Plot,
             RMod=~Fam,
             asrV=4,all.result=T,
             #family=asr_gaussian(),
             mulp=mulp1)
repin(run1,mulp=c(h2 ~ 4 * V1/(V1+V2),
                             H2 ~  V1/(V1+V2/4)))

结果如下：

// result
> run1<-asreml.batch(data=df1