#Load the lymphoma data set:
data(lymphoma)
#Take out a smaller subset of 3 samples (using subsetData):
sub.lymphoma <- subsetData(lymphoma,sample=1:3)
#Run pcf:
seg <- pcf(data=sub.lymphoma,gamma=12)
#Make a matrix with two positions and chromosomes for which we want to
#interpolate the pcf-estimate:
pos <- c(2000000,50000000)
chr <- c(1,2)
x <- cbind(chr,pos)
#Interpolate
int.pcf <- interpolate.pcf(seg,x)
Run the code above in your browser using DataLab