Update 6 files

- /Data Annotation/GO cellular component.png
- /Data Annotation/GO biological function.png
- /Data Annotation/GO molecular function.png
- /Data Annotation/biomart_for_annotation.R
- /Data Annotation/GO_enrichment.R
- /Data Annotation/mrna.csv

Data Annotation/GO_enrichment.R

+library(clusterProfiler)
+library(org.Hs.eg.db)
+library(AnnotationDbi)
+
+data <- read.csv("C:/Users/Emre/Desktop/go.csv", sep = ";", header = TRUE)
+
+genes_to_test <- data$ensembl_gene_id
+
+GO_results <- enrichGO(gene = genes_to_test, OrgDb = "org.Hs.eg.db", keyType = "ENSEMBL", ont = "BP")
+as.data.frame(GO_results)
+
+fit <- plot(barplot(GO_results, showCategory = 15))
+fit
+
+
+GO_results2 <- enrichGO(gene = genes_to_test, OrgDb = "org.Hs.eg.db", keyType = "ENSEMBL", ont = "MF")
+as.data.frame(GO_results2)
+
+fit2 <- plot(barplot(GO_results2, showCategory = 15))
+fit2
+
+
+GO_results3 <- enrichGO(gene = genes_to_test, OrgDb = "org.Hs.eg.db", keyType = "ENSEMBL", ont = "CC")
+as.data.frame(GO_results3)
+
+fit3 <- plot(barplot(GO_results3, showCategory = 15))
+fit3
+

Data Annotation/biomart_for_annotation.R

+library(biomaRt)
+ensembl <- useMart("ensembl", dataset = "hsapiens_gene_ensembl")
+mrna <- read.csv("C:/Users/Emre/Desktop/mrna.csv")
+gene_list <- mrna$Gene_symbol
+gene_list
+
+ensembl_ids <- getBM(attributes = c("ensembl_gene_id", "external_gene_name"),
+                     filters = "external_gene_name",
+                     values = gene_list,
+                     mart = ensembl)
+
+print(ensembl_ids)
+
+library(writexl)
+write_xlsx(ensembl_ids, path = "C:/Users/Emre/Desktop/go.csv")
+
+