analyze.stick.data.batch: Analyze batch data generated by 'sim.fit.stick.data.batch'
In Stickbreaker: Fits Stickbreaking, Multiplicative and Additive Models to Data

Description Usage Arguments Details Value

Analyze batch data generated by sim.fit.stick.data.batch

1 2	analyze.stick.data.batch(mut.vals, coe.vals, sig.vals, inpath, fit.methods, d.range, d.true, relative, coes)

`mut.vals`	Vector of number of mutations to simulate
`coe.vals`	Vector of stickbreaking coefficients to analyze
`sig.vals`	Vector of sigma values to analyze
`inpath`	The path to read input files from. Generally same as `outpath` from `sim.fit.stick.data.batch`.
`fit.methods`	Vector of all methods of estimating d to then fit model and output results. Accepts "MLE", "RDB", "max", "seq", "RDB.all" and "All". "All" does all methods. Default is "seq". Case sensitive. Note you should use same fit.methods (or a subset of) as used when generating data in
`d.range`	Search range for d.
`d.true`	The distance to the boundary (d) used in simulations `sim.fit.stick.data.batch`.
`relative`	`TRUE/FALSE`. Should bias and rMSE be divided by parameter value to make them relative values?
`coes`	Vector of stickbreaking coefficients for each mutation

This function assesses bias and relative root mean squared error (rMSE) from the data generated by sim.fit.stick.data.batch. It also determines the failure rate for each method. An estimate is considered invalid (a failure) when it does not fall in the range d.range. The function also returns the data frame seq.method.props that indicates for the sequential method, what proportion of the estimates come from the MLE, RDB and max methods. If "seq" is not part of fit.methods (and if fit.methods <- "All" then it is), this dataframe is empty but still returned.

List:
[[1]] error.table A data frame with the parametric conditions and the failure rate, bias and rMSE for the parameters d, the stickbreaking coefficients, and sigma.
[[2]] seq.method.props A data frame with the proportion of sequential estimates coming from the MLE RDB and max methods for each parametric condition. If "seq" is not among the fit.methods, the dataframe is returned but empty.

Stickbreaker documentation built on May 29, 2017, 9:01 a.m.