diff options
Diffstat (limited to 'scripts')
| -rw-r--r-- | scripts/wgcna_analysis.R | 20 |
1 files changed, 19 insertions, 1 deletions
diff --git a/scripts/wgcna_analysis.R b/scripts/wgcna_analysis.R index a8170b6..8e90d7d 100644 --- a/scripts/wgcna_analysis.R +++ b/scripts/wgcna_analysis.R @@ -10,10 +10,28 @@ wgcnaRawData <- read.csv(file = "wgcna_data.csv") # transform expressionData -datExpr <- as.data.frame(t(wgcnaRawData)); +dataExpr <- as.data.frame(t(wgcnaRawData)); +# data cleaning + +# adopted from docs +gsg = goodSamplesGenes(dataExpr, verbose = 3); + + + +if (!gsg$allOK) +{ +# Optionally, print the gene and sample names that were removed: +if (sum(!gsg$goodGenes)>0) +printFlush(paste("Removing genes:", paste(names(datExpr0)[!gsg$goodGenes], collapse = ", "))); +if (sum(!gsg$goodSamples)>0) +printFlush(paste("Removing samples:", paste(rownames(datExpr0)[!gsg$goodSamples], collapse = ", "))); +# Remove the offending genes and samples from the data: +dataExpr <- dataExpr[gsg$goodSamples, gsg$goodGenes] +} + |