||
```R
library(asreml)
```
```R
example(asreml)
```
asreml> dat <- data.frame(y=rnorm(20),x=seq(1,20))
asreml> ex.asr <- asreml(y ~ x, data=dat)
ASReml: Thu Jan 05 20:13:46 2017
LogLik S2 DF wall cpu
-11.0424 0.7404 18 20:13:46 0.0
-11.0424 0.7404 18 20:13:46 0.0
Finished on: Thu Jan 05 20:13:46 2017
LogLikelihood Converged
asreml> ## The oats data
asreml> data(oats)
asreml> oats.asr <- asreml(yield ~ Variety*Nitrogen, random = ~ Blocks/Wplots, data=oats)
ASReml: Thu Jan 05 20:13:46 2017
LogLik S2 DF wall cpu
-216.5453 306.1641 60 20:13:46 0.0
-213.5687 256.3384 60 20:13:46 0.0
-211.0520 214.6065 60 20:13:46 0.0
-209.6629 186.6835 60 20:13:46 0.0
-209.4004 179.2338 60 20:13:46 0.0
-209.3782 177.3097 60 20:13:46 0.0
-209.3779 177.0864 60 20:13:46 0.0
-209.3779 177.0833 60 20:13:46 0.0
Finished on: Thu Jan 05 20:13:46 2017
LogLikelihood Converged
很遗憾,不支持。
弃之。
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-12-27 10:49
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社