diff --git a/.~lock.Schrick-Noah_CS-6643_Lab-4.docx# b/.~lock.Schrick-Noah_CS-6643_Lab-4.docx# new file mode 100644 index 0000000..111e6d0 --- /dev/null +++ b/.~lock.Schrick-Noah_CS-6643_Lab-4.docx# @@ -0,0 +1 @@ +,noah,NovaArchSys,29.09.2022 14:42,file:///home/noah/.config/libreoffice/4; \ No newline at end of file diff --git a/Schrick-Noah_CS-6643_Lab-4.R b/Schrick-Noah_CS-6643_Lab-4.R new file mode 100644 index 0000000..bda527f --- /dev/null +++ b/Schrick-Noah_CS-6643_Lab-4.R @@ -0,0 +1,32 @@ +# Lab 4 for the University of Tulsa's CS-6643 Bioinformatics Course +# Differential Expression +# Professor: Dr. McKinney, Fall 2022 +# Noah L. Schrick - 1492657 + +## Set Working Directory to file directory - RStudio approach +setwd(dirname(rstudioapi::getActiveDocumentContext()$path)) + +#### Part A: Preparing Data +load("sense.filtered.cpm.Rdata") +# load phenotype (mdd/hc) data +subject.attrs <- read.csv("Demographic_symptom.csv", + stringsAsFactors = FALSE) + +if (!require("dplyr")) install.packages("dplyr") +library(dplyr) +# grab intersecting X (subject ids) and Diag (Diagnosis) from columns +phenos.df <- subject.attrs %>% + filter(X %in% colnames(sense.filtered.cpm)) %>% + dplyr::select(X, Diag) +mddPheno <- as.factor(phenos.df$Diag) + +# Normalized and transform +library(preprocessCore) +mddExprData_quantile <- normalize.quantiles(sense.filtered.cpm) +mddExprData_quantileLog2 <- log2(mddExprData_quantile) +# attach phenotype names and gene names to data +colnames(mddExprData_quantileLog2) <- mddPheno +rownames(mddExprData_quantileLog2) <- rownames(sense.filtered.cpm) + +length(rownames(sense.filtered.cpm)) +#### Part B: Filter noise genes \ No newline at end of file diff --git a/Schrick-Noah_CS-6643_Lab-4.docx b/Schrick-Noah_CS-6643_Lab-4.docx new file mode 100644 index 0000000..3d88b9c Binary files /dev/null and b/Schrick-Noah_CS-6643_Lab-4.docx differ