dumpData function

Dumps segregation ratio data to file for subsequent JAGS run

Given segregation ratio data provided as an object of class segRatio, data are dumped in format for use by JAGS

dumpData(seg.ratio, model, stem = "test", fix.one = TRUE,
 data.file = paste(stem, "-data.R", sep = ""))

Arguments

• seg.ratio: Object of class segRatio contains the segregation ratios for dominant markers and other information such as the number of dominant markers per individual
• model: Object of class modelSegratioMM containing mixture model information
• stem: File name stem for data file (default test )
• fix.one: Logical to fix the dosage of the observation closest to the centre of each component on the logit scale. This can greatly assist with convergence (Default: TRUE)
• data.file: Data file name which is automatically generated from stem if not specified

Returns

None.

Author(s)

Peter Baker p.baker1@uq.edu.au

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)

dumpData(sr, x)