Any packages for normalization by read depth of the 75%tile gene?

Hi, Biostars.

I want to normalize read count data (by STAR) by read depth of the 75%tile gene like this paper.

Anyone knows useful packages (R, python ,etc) ?

(I tried by "scale" function in R, but couldn't do this.)

thanks!!

rna-seq next-gen

• 1,940 views • link

written 9.9 years ago by cfarmeri

Hi, You can try "quantile" function to calculate 75-percentile value. This link might help you: http://www.r-tutor.com/elementary-statistics/numerical-measures/percentile .

• 0 views • link

written 9.9 years ago by Persistent LABS

0 answers

No answers yet.

More posts like this

RNA seq analysis differential expression analysis
written 3.4 years ago by bioinformatics •

Hi, I'm trying to perform RNA sequencing differential expression analysis on the GEOdataset GSE111320 in the R program. https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE111320 Does anyone know how I might …
How can I be sure that raw read counts are well processed from fastq files?
written 4.8 years ago by Simon Ahn •

Hi. I'm new in bioinformatics and try to process fastq files for getting raw read count matrix. I downloaded fastq files from https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE63452 1. I …
GEO raw data download
written 6.2 years ago by kwanghoon •

Hi I'm trying to download Supplementary file(GSE109381_RAW.tar) at https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE109381 I clicked "http" link, but it will take 2 days... Is there any other way to …
Error when reading microarray idat files from illumina
written 6.8 years ago by salamandra

Hi, I tried to read an idat file with R package illuminaio: library(illuminaio) idatfile <- 'C:/Users/Tania/Documents/GSE26677/GSE26677_non-normalized.txt' idat <- readIDAT(idatfile, what='all') and also with limma package: …
How to analyze aCGH data in R?
written 7.0 years ago by Lila M

Dear all, I've tried to analyze this [data][1] that is a aCGH data in R. But I couldn't find any tutorial that can help me …
How to eliminate batch effect using removeBatchEffect function for tissue expression profiles?
written 7.7 years ago by rice.researcher

For meta-analysis work, I would like to combine two rice GEO series( [GSE27988][1] and [GSE27856][2] ) where each series analyze multiple anatomical tissue related expression …
how can I use Lumi package for GSE31699 file?
written 11.4 years ago by SCA •

I would like to use lumi package for normalize to http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE31699 I read csv file and get a matrix file. I tried to use the …
Finding The Relation Between Gene Expression And Protein Function
written 12.5 years ago by mtyler.jason •

Hi all, I have some microarray dataset like this <a href='http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE1935'>http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE1935</a>. Given this gene expression data, I want to find out which protein function is …
How Do I Find Data From A Particular Study In Dbgap?
written 15.2 years ago by Stephen

I want to follow up some GWAS hits for obesity by doing some eQTL analysis in fat and liver tissue. I found <a href="http://www.sciencedirect.com/science/article/pii/S0002929710001023">this paper</a> …
Normalize Large Number Of Cel Files
written 15.7 years ago by Mike Dewar

I'm trying to normalize a set of 508 CEL files, available on GEO: <a href="http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE15907">http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE15907</a> I've been trying to use the R package aroma affymetrix: …