Skip to content

Commit

Permalink
Rcode
Browse files Browse the repository at this point in the history
  • Loading branch information
Mahsa-Ehsanifard committed Jul 18, 2024
1 parent e23d215 commit 05edd78
Showing 1 changed file with 21 additions and 0 deletions.
21 changes: 21 additions & 0 deletions geoquery.R
Original file line number Diff line number Diff line change
@@ -0,0 +1,21 @@
library(GEOquery)
library(limma)
library(Biobase)
library(dplyr)

data <- getGEO(GEO = "GSE58831", GSEMatrix = T, getGPL = T, AnnotGPL = T)
assay <- data$GSE58831_series_matrix.txt.gz@assayData$exprs
pheno <- data$GSE58831_series_matrix.txt.gz@phenoData@data
feature <- data$GSE58831_series_matrix.txt.gz@featureData@data

#for exchanging prob ids into gene symbols of feature file of GEO
#removing repetitive genes in data
exdat <- assay[!duplicated(feature$`Gene symbol`[match(rownames(assay), rownames(feature))]),]
#exchanging the unique rownames into gene symbols
rownames(exdat) <- feature$`Gene symbol`[!duplicated(feature$`Gene symbol`)]






0 comments on commit 05edd78

Please sign in to comment.