# load the aberration profiles for three data types
example.molecular.data.dir <- paste0(path.package('iSubGen'),'/exdata/');
molecular.data <- list();
for (i in c('cna','snv','methy')) {
molecular.data[[i]] <- load.molecular.aberration.data(
paste0(example.molecular.data.dir,i,'_profiles.txt'),
patients = c(paste0('EP00',1:9), paste0('EP0',10:30))
);
}
# calculate scaling factors for all three data types
scaling.factors <- calculate.scaling(molecular.data);
# save the scaling factors to file
write.scaling.factors(scaling.factors, example.molecular.data.dir);
Run the code above in your browser using DataLab