This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Should I remove samples after normalization of miRNA seq read counts ?

I followed https://www.bioconductor.org/help/workflows/RNAseq123/ for my rnaseq analysis. As my read counts were around 2.5 million, I had to use higher CPM. Hope this should be fine for downstream analysis.

After following command:

par(mfrow=c(1,2))                                  
lcpm <- cpm(y2, log=TRUE)                                
boxplot(lcpm, las=2, col=group$Sample, main="")            
title(main="A. Example: Unnormalised data at CPM-2",ylab="Log-cpm") 
y2 <- calcNormFactors(y2)  
y2$samples$norm.factors
lcpm <- cpm(y, log=TRUE)
boxplot(lcpm, las=2, col=group$Sample, main="")
title(main="B. Example: Normalised data at CPM-2",ylab="Log-cpm")

I got following boxplot graph. ![enter image description here][1] [1]: https://ibb.co/gWW2r6

Based on normalized data, which samples should I remove from analysis ?

rna-seq rna-seq edger mirnas
• 1,884 views • link
updated by Biostar Biostar • written by Björn • Björn

1 answer

Some of the samples look different from the others, in terms of their data distribution via the box-and-whiskers plot; however, I would reserve judgement on outliers without seeing, in addition, a PCA bi-plot and violin plot.

Kevin

• 0 views • link
written by Kevin Blighe Kevin Blighe

Log in to answer this question.

More posts like this
  • filtering lowly expressed reads
    written by Shaimaa Gamal • Shaimaa Gamal
    How do I ensure that filtration of lowly expressed counts in RNA seq data was done correctly? My goal is to find DEGs using limma …
  • WGCNA soft thresholding problem
    written by lilingjoyo • lilingjoyo
    Hi there, I'm using WGCNA to construct co-expression network. My input is z-score normalized expression from RNAseq. Here is the distribution of normalized expression cross …
  • heat map clustering
    written by Nick • Nick
    Hello everyone, During RNA-seq analysis I faced a problem, need to perform hierarchical clustering heatmap and got the following result by using following code var_genes …
  • Paired line between two observations in stripchart ?
    written by Björn • Björn
    I used following command `stripchart(normcounts["831",]~group$Groups1,vertical=TRUE,las=2,cex.axis=0.8,pch=16,col=1:5,method="jitter", ylab= "normalized counts", xlab="Patient groups", main="Main",cex.main=0.8)` and got following figure ![enter image description here][1] [1]: https://ibb.co/mt2aHS My questions: 1. how …
  • How to cluster in heatmap.2
    written by Björn • Björn
    How to cluster heatmap using different distance matrix (Manhattan or Euclidean) and split column or row-wise in following command heatmap.2(lcpm, Colv=FALSE, dendrogram="row",labRow = y2$genes$genes, Rowv …
  • How to regroup columns in R using heatmap.2 or pheatmap ?
    written by Björn • Björn
    Using following heatmap.2 command heatmap.2(highly_variable_lcpm,labRow = y2$genes$genes, Rowv = TRUE, col=rev(morecols(20)),trace="none",tracecol = "black",main="Top 20 most variable genes across samples",cex.main=0.6, ColSideColors=group.col,scale="row",margins=c(10,9)) I got the following graph …
  • boxplot in edgeR
    written by Björn • Björn
    par(mfrow=c(1,2)) lcpm <- cpm(y2, log=TRUE) boxplot(lcpm, las=2, col=group$Sample, main="") title(main="A. Example: Unnormalised data",ylab="Log-cpm") The boxplot I got was https://ibb.co/bZ8RQn Why don't I get normal boxplot …
  • How to create design matrix in edgeR
    written by Björn • Björn
    # edgeR # following are my paired samples (column Patient) with other three columns Sample Groups Patient Groups1 a baseline a B b baseline b …
  • Signature genes from RNA-seq time series (Developmental)
    written by Dataminer Dataminer
    Hi, I have a time series data in tripicates for 4 time points. My aim is to identify signature genes for each time point. What …
  • Adjusting Font.Lab R
    written by RandManP • RandManP
    <p>with the following code, xlab will be bold but ylab not. could you please help me</p> <p>pdf()</p> <p>par(font.axis=2, font.lab=2, cex.lab=1.2, cex.axis=1.2, cex.sub=1.2, mfrow=c(2,2), mar = …