Last chance! 50% off unlimited learning
Sale ends in
data(bristle3)
cross
. See read.cross
for details.There are eight phenotypes: the average and SD of the number of abdominal and sternopleural bristles in males and females for each line.
Each line is typed at 29 genetic markers on chromosome 3.
bristleX
, listeria
,
fake.bc
, fake.f2
,
fake.4way
, hyper
data(bristle3)
# Summaries
summary(bristle3)
plot(bristle3)
# genome scan for each of the average phenotypes
bristle3 <- calc.genoprob(bristle3, step=2)
out <- vector("list", 4)
for(i in 1:4) out[[i]] <- scanone(bristle3, pheno.col=2*i-1)
# Plot the results
# maximum LOD score among four phenotypes
ym <- max(sapply(out, function(a) max(a[,3])))
plot(out[[1]], ylim=c(0,ym))
for(i in 2:4)
plot(out[[i]], add=TRUE, col=c("blue","red","green")[i-1])
Run the code above in your browser using DataLab