#Load lymphoma data
data(lymphoma)
chrom <- lymphoma[,1]
pos <- lymphoma[,2]
#pick out data for the first six samples:
cn.data <- lymphoma[,3:8]
#Create missing values in cn.data at random positions:
n <- nrow(cn.data)*ncol(cn.data)
r <- matrix(rbinom(n=n,size=1,prob=0.95),nrow=nrow(cn.data),ncol=ncol(cn.data))
cn.data[r==0] <- NA #matrix with approximately 5% missing values
mis.data <- data.frame(chrom,pos,cn.data)
#Impute missing values by constant, c=0:
imp.data <- imputeMissing(data=mis.data,method="constant")
#Impute missing values by obtained pcf-values:
pcf.est <- pcf(data=mis.data,return.est=TRUE)
imp.data <- imputeMissing(data=mis.data,method="pcf",pcf.est=pcf.est)
#Or run pcf within imputeMissing:
imp.data <- imputeMissing(data=mis.data,method="pcf")
Run the code above in your browser using DataLab