Description Usage Arguments Value Examples
Run the treatment response diallel analysis.
1 2 3 |
filename |
name of data file |
savedir |
directory to save files in |
treatment |
indicate subset |
trt.colname |
indicate treatment columname |
phenotype |
indicate parameter of interest |
random |
name of random effect, or NULL |
fixed |
name of fixed effect, or NULL |
type |
specify matched pairs or other type of analysis |
BS |
whether or not to use BayesSpike |
thin |
thinning factor for chains |
lengthChains |
MCMC chain length |
numChains |
number of independent MCMC chains |
burnin |
amount of burn-in (discarded) samples from warm-up |
force |
logical indicating whether to proceed despite warnings/errors |
... |
additional arguments |
runs the analysis, generating data files.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 | ## This example may take a couple of minutes
data(FluDiData)
data(PiximusData, envir = environment())
write.csv(FluDiData, file="FluDiData.csv")
filename <- "FluDiData.csv"
args <- c("All", "D0", "NULL", "NULL", "pre", FALSE, "Trt")
treatment <- args[1]
phenotype <- args[2]
random <- args[3]
fixed <- args[4]
type <- args[5]
BS <- as.logical(args[6])
trt.colname <- args[7]
savedir <- "."
run.tr.diallel(filename=filename, savedir=savedir, treatment=treatment,
trt.colname=trt.colname, phenotype=phenotype, random=random,
fixed=fixed, type=type, BS=BS, thin=1, lengthChains=3000, numChains=5,
burnin=500, force = TRUE)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.