# NOT RUN {
# this example is from a haplotype mapping problem submitted by Mao Jianfeng
# }
# NOT RUN {
hapIntList<-
getIntersectList(3,xnames=c("hap.Pd","hap.Pt","hap.Py"))
# enter the data as follows:
# Number of elements in hap.Pd - 1: 27.586
# Number of elements in hap.Pt - 1: 20.689
# Number of elements in hap.Py - 1: 31.035
# Number of elements in hap.Pd-hap.Pt - 1: 10.345
# Number of elements in hap.Pd-hap.Py - 1: 10.345
# Number of elements in hap.Pt-hap.Py - 1: 0
# Number of elements in hap.Pd-hap.Pt-hap.Py - 1: 0
# Total number of elements - 1: 100
# }
# NOT RUN {
hapIntList<-structure(list(structure(c(27.586, 20.689, 31.035),
.Names = c("hap.Pd","hap.Pt","hap.Py")),
structure(c(10.345, 10.345, 0),
.Names = c("hap.Pd-hap.Pt","hap.Pd-hap.Py","hap.Pt-hap.Py")),
structure(0, .Names = "hap.Pd-hap.Pt-hap.Py"),100),
class = "intersectList")
intersectDiagram(hapIntList)
# }
Run the code above in your browser using DataLab