Closed famuvie closed 8 years ago
As a workaround, we can use the following hack.
For example, say that we wanted to use the OPTION maxrounds 0
for AIREML,
library(breedR)
remlf90 <- eval(parse(text = gsub('"sol se"', '"sol se", "maxrounds 0"', deparse(remlf90))))
environment(remlf90) <- environment(breedR::remlf90)
afterwards, rm(remlf90)
recovers the regular operation.
Include a new argument
options
(or perhapscontrol.reml
) of the functionremlf90()
that passes directly the content to theOPTIONS
of the (AI)REML programs. This allows to manually take advantage of some options such asmaxrounds
.