if (FALSE) {
# Load three data types and create an autoencder for each
example.molecular.data.dir <- paste0(path.package('iSubGen'),'/exdata/');
molecular.data <- list();
ae.result <- 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))
);
ae.result[[i]] <- create.autoencoder(
data.type = i,
data.matrix = molecular.data[[i]],
encoder.layers.node.nums = c(10,2)
)$autoencoder;
}
# Create a matrix of the bottleneck layers
irf.matrix <- create.autoencoder.irf.matrix(
data.types = names(molecular.data),
data.matrices = molecular.data,
autoencoders = ae.result
);
}
Run the code above in your browser using DataLab