kbroman · August 29, 2015 14:03
diff --git a/csv2map.R b/csv2map.R
 ## simulate a cross, using genetic map from a CSV file

 # read map.csv
 maptab <- read.csv("map.csv")

 # convert to map object
 pos <- maptab[,3]
 names(pos) <- maptab[,1]
 chr <- maptab[,2]
 map <- split(pos, chr)

 # load R/qtl
 library(qtl)

 # simulate cross with that as the genetic map
 x <- sim.cross(map, n.ind=100, type="bc")
diff --git a/map.csv b/map.csv
	## simulate a cross, using genetic map from a CSV file

	# read map.csv
	maptab <- read.csv("map.csv")

	# convert to map object
	pos <- maptab[,3]
	names(pos) <- maptab[,1]
	chr <- maptab[,2]
	map <- split(pos, chr)

	# load R/qtl
	library(qtl)

	# simulate cross with that as the genetic map
	x <- sim.cross(map, n.ind=100, type="bc")
No results found