dumpData: Dumps segregation ratio data to file for subsequent JAGS run

In polySegratioMM: Bayesian mixture models for marker dosage in autopolyploids

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)

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