MOGDx Modality Processing
Introduction
Code to process different modalities for inclusion in MOGDx pipeline
Specify Project and Dataset
project <- 'BRCA'
dataset <- 'TCGA'
trait <- 'paper_BRCA_Subtype_PAM50'
META File Generation
Meta data is (typically) located in the coldata of the mRNA gene expression experiment for TCGA datasets.
load(paste0('./../data/TCGA-',project,'/mRNA/mRNA.rda'))
# Create coldata and condition table --------------------------------------
coldata <- colData(data)
datMeta <- as.data.frame(coldata[,c('patient','race' , 'gender' , 'sample_type' , trait)])
datMeta <- datMeta[!(is.na(datMeta[[trait]])) , ]
datMeta <- datMeta[!(duplicated(datMeta[ , c('patient' , trait)])) , ]
datMeta[[trait]] <- factor(datMeta[[trait]])
rownames(datMeta) <- datMeta$patient
write.csv(datMeta , file = paste0('./../data/TCGA-',project,'/datMeta.csv'))
mRNA pre-processing
Differential Gene Expression processing pipeline for messenger RNA datasets.
# Pull in Count Matrices
load(paste0('./../data/TCGA-',project,'/mRNA/mRNA.rda'))
count_mtx <- assay(data)
colnames(count_mtx) <- substr(colnames(count_mtx) , 1, 12)
count_mtx <- count_mtx[, !(duplicated(colnames(count_mtx)))]
# Pull in Meta File
datMeta <- read.csv(paste0('./../data/TCGA-',project,'/datMeta.csv') , row.names = 1)
# Get intersection of count and meta
common_idx <- intersect(colnames(count_mtx) , rownames(datMeta))
count_mtx <- count_mtx[ , common_idx]
datMeta <- datMeta[common_idx , ]
# Perform differential expression analysis
diff_expr_res <- diff_expr(count_mtx , datMeta , trait , 500 , 'mRNA')
## [1] "Removing 2667 Genes with all 0"
## [1] "keeping 29995 genes"
## [1] "Removing 27998 Genes"
## [1] "removing outliers"
## [1] "Keeping 1047 Samples"
## [1] "Removed 36 Samples"
## [1] "performing DESeq"
## estimating size factors
## estimating dispersions
## gene-wise dispersion estimates
## mean-dispersion relationship
## final dispersion estimates
## fitting model and testing
## -- replacing outliers and refitting for 2943 genes
## -- DESeq argument 'minReplicatesForReplace' = 7
## -- original counts are preserved in counts(dds)
## estimating dispersions
## fitting model and testing

# Save differential expression results
datExpr <- diff_expr_res$datExpr
datMeta <- diff_expr_res$datMeta
dds <- diff_expr_res$dds
top_genes <- diff_expr_res$top_genes
save(datExpr, datMeta, dds, top_genes, file=paste0('./../data/',dataset,'/raw/',project,'/mRNA_processed.RData'))
miRNA preprocessing
Differential Gene Expression pipeline for micro RNA datasets
#Pull in count data
load(paste0('./../data/TCGA-',project,'/miRNA/miRNA.rda'))
# Get Count Matrices and filter for reads
read_count <- data.frame(row.names = data$miRNA_ID)
read_per_million <- data.frame(row.names = data$miRNA_ID)
for (i in 2:dim(data)[2]) {
if (i%%3 == 2) {
read_count <- cbind(read_count , data[ , i] )
}
if (i%%3 == 0) {
read_per_million <- cbind(read_per_million , data[ , i])
}
}
colname_read_count <- c()
colname_read_per_million <- c()
for (i in 2:dim(data)[2]) {
if (i%%3 == 2) {
colname_read_count <- c(colname_read_count , substr(strsplit(colnames(data)[i] , '_')[[1]][3] , 1,12))
}
if (i%%3 == 0) {
colname_read_per_million <- cbind(colname_read_per_million , substr(strsplit(colnames(data)[i] , '_')[[1]][6] , 1,12))
}
}
colnames(read_count) <- colname_read_count
colnames(read_per_million) <- colname_read_per_million
# Pull in Meta File
datMeta <- read.csv(paste0('./../data/TCGA-',project,'/datMeta.csv') , row.names = 1)
# Get Intersection of ID's
count_mtx <- read_count
count_mtx <- count_mtx[ , !(duplicated(colnames(count_mtx)))]
common_idx <- intersect(colnames(count_mtx) , rownames(datMeta))
count_mtx <- count_mtx[,common_idx]
datMeta <- datMeta[common_idx , ]
# Perform differential expression analysis
diff_expr_res <- diff_expr(count_mtx , datMeta , trait , 200 , 'miRNA')
## [1] "Removing 280 Genes with all 0"
## [1] "removing outliers"
## Warning in bicor(.): bicor: zero MAD in variable 'x'. Pearson correlation was used for individual columns with zero (or missing) MAD.
## [1] "Keeping 1015 Samples"
## [1] "Removed 50 Samples"
## [1] "performing DESeq"
## estimating size factors
## estimating dispersions
## gene-wise dispersion estimates
## mean-dispersion relationship
## final dispersion estimates
## fitting model and testing
## -- replacing outliers and refitting for 109 genes
## -- DESeq argument 'minReplicatesForReplace' = 7
## -- original counts are preserved in counts(dds)
## estimating dispersions
## fitting model and testing

# Save differential expression results
datExpr <- diff_expr_res$datExpr
datMeta <- diff_expr_res$datMeta
dds <- diff_expr_res$dds
top_genes <- diff_expr_res$top_genes
save(datExpr, datMeta, dds, top_genes, file=paste0('./../data/',dataset,'/raw/',project,'/miRNA_processed.RData'))
DNAm preprocessing
Elastic Net Regression and Processing pipeline for DNA methylation datasets
# Load CpG Counts ---------------------------------------------------------
load(paste0('./../data/TCGA-',project,'/DNAm/DNAm.rda'))
count_mtx <- assay(data)
to_keep = complete.cases(count_mtx) #removed 191928 cpg sites
length(to_keep) - sum(to_keep)
## [1] 191928
count_mtx <- t(count_mtx[to_keep,])
rownames(count_mtx) <- substr(rownames(count_mtx) , 1,12)
# Compute the variance across CpG sites
cpg_variances <- colVars(count_mtx)
# Sort the variances in descending order and get the indices
sorted_indices <- order(cpg_variances, decreasing = TRUE)
# Select the top 200000 most variable CpG sites
num_top_cpg <- 200000
top_cpg_indices <- sorted_indices[1:num_top_cpg]
count_mtx <- count_mtx[ , top_cpg_indices]
# Pull in Meta File
datMeta <- read.csv(paste0('./../data/TCGA-',project,'/datMeta.csv') , row.names = 1)
# Get Intersection of ID's
count_mtx <- count_mtx[!(duplicated(rownames(count_mtx))) , ]
common_idx <- intersect(rownames(count_mtx) , rownames(datMeta))
count_mtx <- count_mtx[common_idx , ]
datMeta <- datMeta[common_idx , ]
# Run glmnet to get CpG's associated with phenotypes of interest ------------
phenotypes <- datMeta[,c('patient' , trait )]
colnames(phenotypes)
## [1] "patient" "paper_BRCA_Subtype_PAM50"
traits <- c(trait)
traitResults <- lapply(traits, function(trait) {
cvTrait(count_mtx, phenotypes, traits, nFolds = 10)
})
## [1] "Removing rows with missing paper_BRCA_Subtype_PAM50 from training data."
## [1] "Fitting lasso model"
## Training
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
##
## Call: cv.glmnet(x = trainMethyl, y = as.factor(trainPhenotypes[[trait]]), type.measure = "class", parallel = TRUE, trace.it = 1, seed = 42, family = "multinomial", alpha = 1, nFolds = nFolds)
##
## Measure: Misclassification Error
##
## Lambda Index Measure SE Nonzero
## min 0.03405 48 0.2141 0.01418 17
## 1se 0.05421 38 0.2282 0.01250 10
cpg_sites <- c()
for (res in traitResults) {
trait_coefs <- coef(res$model , s = "lambda.min")
cpg_sites_tmp <- c()
for (coefs in trait_coefs) {
class_coefs <- rownames(coefs)[which(coefs != 0)]
class_coefs <- class_coefs[2:length(class_coefs)]
cpg_sites_tmp <- unique(c(cpg_sites_tmp ,class_coefs ))
}
cpg_sites[[res$trait]] <- cpg_sites_tmp
}
# Save CpG sites and expression data
datExpr <- count_mtx
save(cpg_sites , datExpr , datMeta , file = paste0('./../data/',dataset,'/raw/',project,'/DNAm_processed.RData'))
Protein (RPPA) pre-processing
Elastic Net Regression and Processing pipeline for proteomic datasets
load(paste0('./../data/TCGA-',project,'/RPPA/RPPA.rda'))
count_mtx <- t(data[ , 6:ncol(data) ])
colnames(count_mtx) <- data$peptide_target
rownames(count_mtx) <- substr(rownames(count_mtx) , 1,12)
# Load Meta Data -----------------------------------------------------------
datMeta <- read.csv(paste0('./../data/TCGA-',project,'/datMeta.csv') , row.names = 1)
# Subset to match ID's ----------------------------------------------------
count_mtx <- count_mtx[!(duplicated(rownames(count_mtx))) , ]
common_idx <- intersect(rownames(count_mtx) , rownames(datMeta))
count_mtx <- count_mtx[common_idx , ]
datMeta <- datMeta[common_idx , ]
# Run glmnet R2 regression to identify proteins of interest -------------------------------
count_mtx <- count_mtx[,colSums(is.na(count_mtx))<0.5*nrow(count_mtx)] #remove columns with more than 50% NA
for(i in 1:ncol(count_mtx)){
count_mtx[is.na(count_mtx[,i]), i] <- mean(count_mtx[,i], na.rm = TRUE)
}
phenotypes <- datMeta[,c('patient' , trait , 'race' , 'gender')]
colnames(phenotypes)
## [1] "patient" "paper_BRCA_Subtype_PAM50" "race" "gender"
traits <- c(trait)
traitResults <- lapply(traits, function(trait) {
cvTrait(count_mtx, phenotypes, trait, nFolds = 10)
})
## [1] "Removing rows with missing paper_BRCA_Subtype_PAM50 from training data."
## [1] "Fitting lasso model"
## Training
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
##
## Call: cv.glmnet(x = trainMethyl, y = as.factor(trainPhenotypes[[trait]]), type.measure = "class", parallel = TRUE, trace.it = 1, seed = 42, family = "multinomial", alpha = 1, nFolds = nFolds)
##
## Measure: Misclassification Error
##
## Lambda Index Measure SE Nonzero
## min 0.01034 36 0.1968 0.01027 24
## 1se 0.02877 25 0.2025 0.01197 13
protein_sites <- c()
for (res in traitResults) {
trait_coefs <- coef(res$model , s = "lambda.min")
protein_sites_tmp <- c()
for (coefs in trait_coefs) {
class_coefs <- rownames(coefs)[which(coefs != 0)]
class_coefs <- class_coefs[2:length(class_coefs)]
protein_sites_tmp <- unique(c(protein_sites_tmp ,class_coefs ))
}
protein_sites[[res$trait]] <- protein_sites_tmp
}
# Save proteins and expression data
datExpr <- count_mtx
save(protein_sites , datExpr , datMeta , file = paste0('./../data/',dataset,'/raw/',project,'/RPPA_processed.RData'))
CNV pre-processing
Elastic Net Regression and Processing pipeline for Copy Number Variation datasets
# Read in Count Matrix
load(paste0('./../data/TCGA-',project,'/CNV/CNV.rda'))
count_mtx <- t(assay(data))
rownames_mtx <- c()
for (name in strsplit(rownames(count_mtx) , ',')) {
rownames_mtx <- c(rownames_mtx , substr(name[1] ,1, 12))
}
rownames(count_mtx) <- rownames_mtx
count_mtx <- count_mtx[,colSums(is.na(count_mtx))<0.5*nrow(count_mtx)] #remove columns with more than 50% NA
# Read in Meta Data
datMeta <- read.csv(paste0('./../data/TCGA-',project,'/datMeta.csv') , row.names = 1)
#Intersect to common IDs
count_mtx <- count_mtx[!(duplicated(rownames(count_mtx))) , ]
common_idx <- intersect(rownames(count_mtx) , rownames(datMeta))
count_mtx <- count_mtx[common_idx , ]
datMeta <- datMeta[common_idx , ]
# Replace NA's with 0's
count_mtx[is.na(count_mtx)] <- 0
# Perform log transform on count matrix to give normal distribution resemblance
count_mtx_log <- log(count_mtx)
# Run glmnet R2 regression to get CNV's of interest
phenotypes <- datMeta[,c('patient' , trait , 'race' , 'gender')]
colnames(phenotypes)
## [1] "patient" "paper_BRCA_Subtype_PAM50" "race" "gender"
traits <- c(trait)
traitResults <- lapply(traits, function(trait) {
cvTrait(count_mtx_log, phenotypes, traits, nFolds = 10)
})
## [1] "Removing rows with missing paper_BRCA_Subtype_PAM50 from training data."
## [1] "Fitting lasso model"
## Training
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
## | | | 0% | |= | 1% | |=== | 2% | |==== | 3% | |===== | 4% | |====== | 5% | |======== | 6% | |========= | 7% | |========== | 8% | |============ | 9% | |============= | 10% | |============== | 11% | |================ | 12% | |================= | 13% | |================== | 14% | |==================== | 15% | |===================== | 16% | |====================== | 17% | |======================= | 18% | |========================= | 19% | |========================== | 20% | |=========================== | 21% | |============================= | 22% | |============================== | 23% | |=============================== | 24% | |================================ | 25% | |================================== | 26% | |=================================== | 27% | |==================================== | 28% | |====================================== | 29% | |======================================= | 30% | |======================================== | 31% | |========================================== | 32% | |=========================================== | 33% | |============================================ | 34% | |============================================== | 35% | |=============================================== | 36% | |================================================ | 37% | |================================================= | 38% | |=================================================== | 39% | |==================================================== | 40% | |===================================================== | 41% | |======================================================= | 42% | |======================================================== | 43% | |========================================================= | 44% | |========================================================== | 45% | |============================================================ | 46% | |============================================================= | 47% | |============================================================== | 48% | |================================================================ | 49% | |================================================================= | 50% | |================================================================== | 51% | |==================================================================== | 52% | |===================================================================== | 53% | |====================================================================== | 54% | |======================================================================== | 55% | |========================================================================= | 56% | |========================================================================== | 57% | |=========================================================================== | 58% | |============================================================================= | 59% | |============================================================================== | 60% | |=============================================================================== | 61% | |================================================================================= | 62% | |================================================================================== | 63% | |=================================================================================== | 64% | |==================================================================================== | 65% | |====================================================================================== | 66% | |======================================================================================= | 67% | |======================================================================================== | 68% | |========================================================================================== | 69% | |=========================================================================================== | 70% | |============================================================================================ | 71% | |============================================================================================== | 72% | |=============================================================================================== | 73% | |================================================================================================ | 74% | |================================================================================================== | 75% | |=================================================================================================== | 76% | |==================================================================================================== | 77% | |===================================================================================================== | 78% | |======================================================================================================= | 79% | |======================================================================================================== | 80% | |========================================================================================================= | 81% | |=========================================================================================================== | 82% | |============================================================================================================ | 83% | |============================================================================================================= | 84% | |============================================================================================================== | 85% | |================================================================================================================ | 86% | |================================================================================================================= | 87% | |================================================================================================================== | 88% | |==================================================================================================================== | 89% | |===================================================================================================================== | 90% | |====================================================================================================================== | 91% | |======================================================================================================================== | 92% | |========================================================================================================================= | 93% | |========================================================================================================================== | 94% | |============================================================================================================================ | 95% | |============================================================================================================================= | 96% | |============================================================================================================================== | 97% | |=============================================================================================================================== | 98% | |================================================================================================================================= | 99% | |==================================================================================================================================| 100%
##
## Call: cv.glmnet(x = trainMethyl, y = as.factor(trainPhenotypes[[trait]]), type.measure = "class", parallel = TRUE, trace.it = 1, seed = 42, family = "multinomial", alpha = 1, nFolds = nFolds)
##
## Measure: Misclassification Error
##
## Lambda Index Measure SE Nonzero
## min 0.01129 61 0.2571 0.009665 119
## 1se 0.01798 51 0.2665 0.015134 74
cnv_sites <- c()
for (res in traitResults) {
trait_coefs <- coef(res$model , s = "lambda.min")
cnv_sites_tmp <- c()
for (coefs in trait_coefs) {
class_coefs <- rownames(coefs)[which(coefs != 0)]
class_coefs <- class_coefs[2:length(class_coefs)]
cnv_sites_tmp <- unique(c(cnv_sites_tmp ,class_coefs ))
}
cnv_sites[[res$trait]] <- cnv_sites_tmp
}
rm(count_mtx_log)
# Save CNV's and expression
datExpr <- count_mtx
save(cnv_sites , datExpr , datMeta , file = paste0('./../data/',dataset,'/raw/',project,'/CNV_processed.RData'))