Hi,
Could anyone help me get some sort of a count matrix when processed RNA-seq data is extracted from GEO. I have some experience using microarray data from GEO.
For example, for GSE78220, the data has been processed by: "FASTQ files were mapped by Tophat2 Tophat BAMs were quantified and normalized using Cuffnorm (for gene analysis) and by htseq-count followed by edgeR's log CPM (for gene-set analysis) Genome_build: hg19 Supplementary_files_format_and_content: The normalized expression levels by cuffnorm"
After geting the GSE by:
g <- getGEO("GSE78220")
What would be the next step?
•
written
by
ahnje770 •
1 answer
I've been recently working on reanalyzing this data set. The data for this series was provided in the supplementary file and cannot be directly accessed by exprs(). So just download the supplementary file "GSE78220_PatientFPKM.xlsx" into your working directory.
assay <- read.csv("GSE78220_PatientFPKM.csv", row.names = 1)
For phenotypic data,
e <- getGEO("GSE78220", destdir = ".")
e <- e[[1]]
pheno <- pData(e)
Log in to answer this question.
-
In RNA-Seq data modeling process, glmfit function in R used raw RNA-Seq count data?
written by mytlsdnkr •
-
GEO RNA seq to DEG
written by unislure •
-
Low assigned alignments
written by Simon Ahn •
-
How can I be sure that raw read counts are well processed from fastq files?
written by Simon Ahn •
-
Analyzing processed RNA-seq data: is this TMM normalized?
written by Ridha •
-
How to normalize gene expression with CPM
written by concetta •
-
HTSeq-count in RNA-Seq
written by mlemusfuentes •
-
Error on HT-Seq-count while RNA quantification
written by k.kathirvel93
-
when using htseq-count, should I only count read-pairs whose two mates were both reported ?
written by Pei
-
processing Illumina HT12v4.0 expression data from GEO
written by Abhi
