dumpData {polySegratioMM}

Dumps segregation ratio data to file for subsequent JAGS run

Description

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

Usage

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

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)

dumpData(sr, x)