## simulate small autooctaploid data seta1 <- sim.autoMarkers(8,c(0.7,0.2,0.1),n.markers=100,n.individuals=50)##print(a1)sr <- segregationRatios(a1$markers)x <- setModel(3,8)## fit simple model in one hit## Not run:x.run <- runSegratioMM(sr, x, burn.in=200, sample=500)print(x.run$doses)## End(Not run)