Code for parallelizing R/qtlhot. See installed parallel directory for proper use. There is apparently an S3 parallel method, so doc has to be as shown below, even though it is called as parallel.qtlhot.
# S3 method for qtlhot
parallel(x, data = 1, …, dirpath = ".")
qtlhot.phase0(dirpath, init.seed = 92387475, len = rep(400,16), n.mar = 185, n.ind = 112,
n.phe = 100, latent.eff = 0, res.var = 1, lod.thrs, …)
big.phase0(dirpath, cross, trait.file, trait.matrix, droptrait.names = NULL,
keeptrait.names = NULL, lod.thrs, sex = "Sex", trait.index,
batch.effect = NULL, size.set = 250, offset = 0, subset.sex = NULL, verbose = TRUE)
phase of parallel processing (1,2,3)
index for parallel processing (1,2,...)
additional arguments passed along
directory path as character string
initial seed for pseudorandom number generation
vector of chromosome lengths for simulated map
number of markers for simulated map
number of individuals for simulated cross
number of phenotypes for simulated phenotypes
size of latent effect
magnitude of residual variance
vector of LOD thresholds to examine
object of class cross
character string name of trait file
character string name of trait matrix
vector of character strings for traits to drop
(none if NULL
)
vector of character strings for traits to keep
(keep all if NULL
)
character string name of phenotype for sex
vector of character strings for trait names
character string for batch effect (none if NULL
)
maximum size of set of traits to scan at one time
offset for name of trait RData files
string of sex to subset on (both sexes if NULL
)
verbose output if TRUE