writeJagsFile: Writes BUGS file for processing by JAGS

In polySegratioMM: Bayesian mixture models for marker dosage in autopolyploids

Description

Given the model and priors a file is written to disk for subsequent JAGS run. BUGS code contained in the model and priors objects is combined and alterered if necessary

Usage

writeJagsFile(model, priors, stem = "test")

Arguments

model	object of class modelSegratioMM specifying model parameters, ploidy etc
priors	Object of class priorsSegratioMM specifying priors
stem	File name stem for BUGS file (default “test”)

Value

None.

Author(s)

Peter Baker [e-mail]

See Also

segRatio dump

Examples

## simulate small autooctaploid data set
a1 <- sim.autoMarkers(8,c(0.7,0.2,0.1),n.markers=100,n.individuals=50)

## compute segregation ratios
sr <-  segregationRatios(a1$markers)

## set up model for 3 components of autooctoploid
x <- setModel(3,8)
x2 <- setPriors(x)

dumpData(sr, x)
inits <- setInits(x,x2)
dumpInits(inits)
##x.priors <- setPriors(x, "vague")
writeJagsFile(x, x2, stem="test")

polySegratioMM documentation built on May 2, 2019, 4:41 p.m.