<div class="preformatted">Hi All
I haven't used LIMMA for a while, but I remember that when I did
there was a special way to handle duplicate spots. However the
layout of the array had to be such that one replicate was printed
first for all spots and then the...I haven't used LIMMA for a while, but I remember that when I did
there was a special way to handle duplicate spots. However the
la…
Order by: Relevance
Hi
I'm trying to work with DESeq2, i already worked with edgeR and DESeq1 with relatively sucess.
They give me the same results as logFC (+ or -) but, when i work...with DESeq2, i have different results. Every gene who are positive for edgeR and DESeq1 are negative for DESeq2 and vice-versa.
DESeq2...I'll post my pipeline, thanks for your help!
bckCountTable <- read.table("DF.txt", hea…
updated 7.8 years ago • Luiz
used the following code, and got 266 genes with padj < 0.05:
table <- read.delim("test.txt")
row.names(table) <- table$Feature\_ID
count\_table <- table\[, -1\]
conds <- c("ctrl", "ctrl", "ctrl", "treatment")
cds <- newCountDataSet(count...method="blind", sharingMode="fit-only")
results <- nbinomTest(cds, "ctrl", "treatment")
In DE…
updated 9.3 years ago • wd
with readGAlignments, I usually set isDuplicate = FALSE in the scanBamParam argument to remove duplicates. I was just wondering exactly how duplicates are determined--is this based on just chromosome and the start and
updated 8.1 years ago • jfiksel
gt;The duplicated spots have random location, which means that the
number of spots between each duplicate is not the same for every gene.
This is the summary for the distances:
>
> Min. 1st Qu. Median Mean 3rd Qu. Max.
> 4.00...71.00 86.59 135.00 244.00
>
>(Distance here means number of spots between the two duplicates)
>
>The functio…
5 conditions (vehicle, tx1, tx2, tx3, tx4). When I label the groups "A", "B", "C", "D", "E" and run DESeq2 on all 5 groups at once, I get a list of 7000+ differentially expressed genes (padj < 0.05). However, when I label the groups...Shouldn't my results theoretically be the same regardless of how I named the groups, since I gave DESeq2 the same data? Am I doing something wrong/missing so…
updated 7.1 years ago • blobbyscience
vgsID , "GENE.SYMBOL")
#then I got
Error in `.rowNamesDF<-`(x, value = value) :
duplicate 'row.names' are not allowed
In addition: Warning message:
non-unique values when setting 'row.names': ........
```
I have use the...duplicate() function to check, but the console return FALSE
Does anyone have any idea how to solve this? Thanks
updated 2.1 years ago • jscl1n22
I have a RNASeq data-set which consists of 4 different conditions (Control, Treatment1, Treatment2 & Treatment3) each with duplicate samples.
I would like to identify all genes that are significantly differentially expressed across all the 4 groups...consists of 4 different conditions (Control, Treatment1, Treatment2 & Treatment3) each with duplicate samples.
I would like to ide…
updated 7.1 years ago • gupta.anuj0608
i have normalized the data using deseq2 internal normalization method. But it seems that my original input data is same as the normalized output. i don't know...on adjusted p value <0.1
Code should be placed in three backticks as shown below
```r
library(DESeq2)
library(tidyverse)
count_data <- read.csv('rawcount_data2.csv',row.names = 1)
head(count_data)
colData <- read.c…
updated 11 months ago • sajadahmad41454
the differential binding analysis of some human histone ChIP-seq samples (bowtie2 mapping to hg19, duplicate marking, MACS2 narrow peak calling, blacklisted regions filtering with bedtools intersect and Anshul Kundaje...with the Y vs H comparison:
DiffBind: Error in `.rowNamesDF<-`(x, value = value) : invalid 'row.names' length
at the point of my script that I was doing a volcano p…
updated 5.9 years ago • melnuesch
inspired by this post
\#http://seqanswers.com/forums/showthread.php?t=31036
\#install.packages("DESeq2")
\#\# try http:// if https:// URLs are not supported
\#source("https://bioconductor.org/biocLite.R")
\#biocLite("DESeq2")
\# Set working directory...head (dat)
final\_count <- read.csv("rpkm.csv\_0\_filter\_atleat\_200\_count.csv", header=TRUE, row.names=1)
library(DESeq2…
I downloaded the newest version of R as I kept getting prompted to update in RStudio. As I did, for some reason, previous work I had done following the DESeq2 vignette would not run properly giving me many errors like: 'DESeq2' package doesn't exist. I'll show my most recent errors but I can't run DESeq2 at all.
Code should be placed in three backticks as shown below
```r
library("DESeq2"…
updated 21 months ago • saguiar
Hello!
I'm working on comparing tumor and healthy samples with DESeq2. My analysis only cares about the adjusted Fold Change values that DESeq2 computes.
The input data encompasses a huge...the process up given the fact that I merely require the adjusted fold change instead of the whole DESeq2 output?
The full script that I run on the data is [here at this permalink](https://github.com/TCP…
updated 17 months ago • Luca
Hi!!
I'm having trouble to do a contrast using DESeq2. I have 12 samples, AA(2 replicates), AB(4 replicates), BA(4 replicates) and BB(2 replicates). I wanted to compare AB (or BA) with...in resultNames().
Here is my code:
<pre>
countData <- read.csv(file = dataFile, header = TRUE, row.names=1)
colData <- DataFrame(row.names=colnames(countData), condition=factor(c(rep("AA…
updated 7.5 years ago • zhuozhu132
gene
The duplicated spots have random location, which means that the number
of spots between each duplicate is not the same for every gene. This
is the summary for the distances:
Min. 1st Qu. Median Mean 3rd Qu. Max.
4.00 32.00 71.00...86.59 135.00 244.00
(Distance here means number of spots between the two duplicates)
The function duplicateCorrelation in limma can be used to…
<div class="preformatted">Hi,
I'm trying to use DESeq to know the differential expressed genes of my
datasets and i'm encountering that DESeq is not recognizing my
row.names so
i can't create my cds.
My .csv input file looks like:
transcript_id,C4,CRL_2APR10,CRL_1_15JUL11,CRL_2_15JUL11
comp1000201_c0_seq1,5.00,0.00,0.00,0.00
comp1000297_c0_seq1,7.00,0.00,0.00,0.00
comp100036_c0_seq1,0.00…
div class="preformatted">Hi ,
I am using Deseq2 for differential gene expression calculations. I
have two conditions Condition a (T) and Condition b (P),
when I run the...Thanks
in advance, code is :
> countsTable <- read.delim ("48_96_filtered.txt", header=TRUE,
row.names=1)
> pdata = data.frame(condition = factor(c( "T", "T", "T", "P", "P",
"P", "P")))
> library …
permu.size=10000, Pe=0.001)
I got this error
<pre>
Error in `.rowNamesDF<-`(x, value = value) :
duplicate 'row.names' are not allowed
In addition: Warning message:
non-unique values when setting 'row.names': 'cg00042263.ENSG00000146556
an
excellent job with limmaGUI!
My question is a simple one.
The arrays I work with tend to have duplicates on each slide. However,
they are not adjacent nor equally spaced and so when setting the
spacing of the duplicates...which appears to be compulsory)
I just make it up.
I presume this makes the correlation between duplicates meaningless
but that linear model fits are OK.
Comments?
Cheer…
Using the following code to look at differential expression in a couple of datasets. Over and over again the # of up and and down regulated genes is approximately equal. I just wonder if I'm missing something, but I can't find anything about this anywhere. It can't be an artifact of normalization. Any ideas?
library("DESeq2")
countData <- as.matrix(read.csv("gene_count_matrix.csv", r…
updated 4.5 years ago • sorrymouse
Hi Rory,
I had been previously, not marking duplicates or removing them in my ChIP-Seq datasets, and instead setting the dba.count bRemoveDuplicates to TRUE. Firstly
updated 6.1 years ago • loretta
div class="preformatted">Hi I want to compare DESeq vs DESeq2 and I am getting different number
of
DEGs which I will expect to be normal. But when I compare the 149
genes iD
that I get...with DESeq with the 869 from DESeq2 there are only ~10
genes
that are in common which I don\'92t understand (using FDR <0.05 for
both).
I want to block the Subject...using:\
\
\
DESeq\
\
library(DESeq)…
I am posting in regards to an error I've experience when running Deseq2's estimateSizeFactors. The specific error is posted below:
<pre>
> dds = DESeqDataSetFromMatrix(countData=filteredDGE...every gene contains at least one zero, cannot compute log geometric means</pre>
I've run Deseq2 several times in the past but have never experienced this error. I have trie…
updated 9.3 years ago • ssabri
<div class="preformatted">Dear List,
I had posted earlier about my problem but didn't get any response that
could
help so I am seeking your help again.
I have duplicate spots on my nimblegen array that are position
randomized so
there is no order to their placement. In Limma when we call the
duplicateCorrelation() function, it requires that the duplicate spots
have
some order to them, lik…
I am using DESeq2 for differential expression analysis. But I got big different results between the two versions while using the same...set of data as input and the same set of r code for data processing for the two DESeq2 versions. I would like to ask what is the key point changed between the two versions to make the results so different...I used DESeq2 package version 1.14.1 and pasilla version…
I am following this to make pseudobulk and perform DESeq2
https://hbctraining.github.io/scRNA-seq/lessons/pseudobulk_DESeq2_scrnaseq.html
Here they perform the DE analysis...colnames(counts) %in% rownames(cluster_metadata))])" for my 3-cell analysis because I can't have duplicate rownames for cluster_metadata. The row names of the metadata are sample names. But the samples would be the same...fo…
updated 13 months ago • puffsquash
Hi,
I am using DESeq2 for DE analysis in RNA-Seq experiment. I have 28 CLL patients raw counts and pooled male control and pooled female control...i.e. not matching with existing literature survey). e.g. hsa-mir-155 should be over-expressed but DESeq2 analysis shows under-expressed etc.
Here are my codes :
expDesign = data.frame(row.names = colnames(cts), condition = c("treated...DESeq…
updated 6.2 years ago • vivekr
Hi,
csaw allows one to remove PCR duplicates when doing read counting across windows. My question is how is this functionality actually executed.
I generally...mark duplicates using Picard Tools or Samblaster and I was of the view that csaw would remove these reads when doing read counting...wondering if csaw calls some other program (or its own method) internally for marking and remo…
lt;-read.table(file="TestFile.txt",head=TRUE,row.names=1)
expt\_design <- data.frame(row.names = colnames(total\_counts),
 ...OQ\_YQ.Sig\[complete.cases(OQ\_YQ.Sig),\]
write.table(OQ\_YQ.Sig,file="Results\_sig.txt",sep="\\t",row.names=rownames(OQ\_YQ.Sig),
&n…
Hello Michael,
I am new to DESeq2.
I have successfully used the results of function results() to do the plotMA. However, when I try
<pre>
lfcShrink...found</pre>
My codes are:
<pre>
setwd("C:/cygwin64/home/Coexpression_Nov2017")
getwd()
library(DESeq2)
library(apeglm)
library(ashr)
expected_matrix <- read.table(file = "17samples_expected_co…
updated 7.4 years ago • lychen83
software, specifically with the
GUI interface and are running into some problems. Our array has
duplicate spots in a top to bottom manner. When we try to have Limma
calculate a duplicate correlation we get the error message
TRUE)
> testFile <- tempfile()
> write.table(KarinaSusc, testFile, quote = FALSE, sep = "\t",
row.names
+ =TRUE, col.names = TRUE)
> eSet <- read.exprSet(testFile)
> eSet
Expression Set (exprSet) with
7150 genes
7 samples...sample: arbitrary numbering
2.) Second
>
KarinaSusc<-read.table("Ka_Susc.txt",header=TRUE,row.names=1,…
updated 21.4 years ago • Marcelo Luiz de Laia
I tried to follow the [DESeq2 tutorial](https://www.bioconductor.org/packages/3.7/bioc/vignettes/DESeq2/inst/doc/DESeq2.html#how-to-get-help-for...deseq2), and ended up with a error.
library("pasilla")
library("DESeq2")
pasCts <- system.file("extdata", "pasilla_gene_counts.tsv...package="pasilla", mustWork=TRUE)
cts <- as.matrix(read.csv(pasCts,sep="\t",row.names…
updated 7.2 years ago • unamourdeswann
I have missed any essential components in my script for differential gene expression analysis using DESeq2. I first filtered my data by adjusted p value (0.05) and ranked those by fold change.
The following is an example of my output...I be filtering something else? Have I missed something?
Here is the R script:
##Load DESeq2 R package
library(DESeq2)
##Load in raw re…
updated 5.7 years ago • umajs
div class="preformatted">Dear List,
I was wondering if there is a tool which can help me remove the
duplicate probes based on the location of the probes?. To be precise,
I wanted to select the probe which is closest to the 3'UTR
followed by featurecounts for raw counts matrix. when i am trying to process the raw_counts using deseq2 the results are not satisfactory to me. please suggest me if it is something wrong with my data![input normalized file...1]
this is the r-script which i am running
```r
library(DESeq2)
library(tidyverse)
count_data <- read.csv('rawcount_data2.csv',row.names = 1)
head(count_data)
co…
updated 11 months ago • sajadahmad41454
from RPKMs normalized by DESeq2. I guess this is because DESeq2 reports estimated logfc. However, sometimes sign (- or + sign) is changed. Could...res <- results(model, contrast=c("group”,"D16”,"D8"))
tmp\_df <- data.frame(res@listData)
row.names(tmp\_df) <- res@rownames
write.csv(tmp\_df, "./DESeq2\_results/Wald\_D.csv")
plotMA(res, alpha=0.…
updated 9.4 years ago • bioboy
I have repeated sampling of the same individuals I consider it would be appropriate to use Limma and Duplicate correlation (in comparisons involving paired samples) - correct?
Q2) We have discussed using Deseq2 as well, but if...I understand, the duplicate correlation function is not available there and therefore it is better to use the code above and limma in the entire...and how and where ca…
updated 3.9 years ago • sofie.haglund
calculation: featurecount assigned reads/total reads *1000000). Here is an example gene output from Deseq2:
```
Row.names baseMean log2FoldChange lfcSE pvalue padj input_rep1 input_rep2 IP_rep1 IP_rep2
WBGene00006575 70.2292642324974...input2, IP1, IP2, so this small RNA is not enriched in IP. I am not sure if the RPM or log2FC from Deseq2 give the correct result.
Here is code I used for Deseq…
updated 18 months ago • Shihui
bioconductor@stat.math.ethz.ch>
Sent: Friday, January 07, 2005 1:01 PM
Subject: random location of duplicate spots and use of limma
> >
>>
>>Hello
>>
>>There are approximately 6000 different genes on the arrays...The duplicated spots have random location, which means that the
number of
>>spots between each dup…
I could order the genes and get them next to each other. My
question is, how would I incorporate the duplicate correlation aspect
when defining the design in limma?
Thank you for your reply and help.
Yolande
[[alternative HTML
Hi
I'm trying to use DESeq2 to find differentially expressed genes between sample with replicates and sample without replicates. I have 2 condition...control and tumour. I am new to Deseq2 package, do you think my command lines are OK for my experiment?
Thanks in advance
code for find DE in tumour...nbsp;
<pre>
> countsTable <- read.delim ("/featur…
updated 8.3 years ago • Bob
I have run DeSeq2 in 2 different devices using the same count data and metadata tables, however when running the summary and results...both devices and also after round the values I get the same count table values in both,
```library("DESeq2")
# sampleTable #
csvfile <- "metadata.csv"
sampleTable <- read.csv(csvfile, row.names=1, sep=",")
sampleTable$Stage <- as.factor...da…
updated 2.2 years ago • Mafer
I have RNA-seq data from nuclear (Nuc) and cytoplasmic (cyto) fractions of both undifferentiated (UD) and differentiated (DF) cells (I have triplicates of all conditions).
Essentially I have UDNuc1, UDNuc2, UDNuc3, UDcyto1, UDcyto2, UDcyto3, DFNuc1, DFNuc2, DFNuc3, DFcyto1, DFcyto2, DFcyto3.
I have run DESeq2 on UD nuc vs UD cyto samples and DF nuc vs DF cyto samples separately but I am t…
updated 5.0 years ago • AmyNeely2023
https://support.bioconductor.org/p/79209/ and https://support.bioconductor.org/p/98476/).
I use DESeq2 in order to get normalized counts from raw TCGA-HTseq Counts.
My method is simple, put the data through DESeq2 - using design...all of my analysis (clustering, survival, gene expression heat maps) on the normalised counts from Deseq2.
However I am running into problems with huge numbe…
updated 5.5 years ago • Alex Greenshields-Watson
nbsp;
We perform our DE analysis with limma+voom packages with clean reads without PCR duplicates removal.
for testing how PCR duplicates would affect DE gene and FPKM value, we also map reads with PCR duplicates...DE analysis with the same pipeline.
It turns out that most of the DE genes are overlap in both PCR duplicates removal group and without PCR duplicates removal&…
updated 8.8 years ago • panying0713
I realize some areas may map to multiple locations in the genome, but I would like to eliminate duplicates if possible. For example, I have some DMPs that have multiple "hg19_genes_introns" entries. This is a single position...so I would like it to only have one entry. If the number of duplicates were the same for each position, this wouldn't be as much of an issue, because the proportions would …
updated 14 months ago • hayleyw
analysis performed by Anota2Seq). I am, however, struggling to understand how to accomplish this in DESeq2.
My design is as follows:
cts <- read.delim("Interaction.txt", row.names=1)
coldata <- read.delim("Sample_Sheet.txt", row.names
updated 5.6 years ago • timjoncooper
div class="preformatted">Dear all,
On my microarray I have some genes for which I have duplicate probes,
but not for all. Before I do the stats I would like to average the
duplicate probes (e.g. GB10001-RA and GB10006...ideas
what would be the most elegant way to average out expression values
based on the occurrence of duplicate entries in RG$genes$GeneName?
Limma can only average duplicates i…
updated 15.0 years ago • Tom Wenseleers
since I clearly don't know the required R code to do it, I'm here asking for help. I've generated a DESeq2 csv file containing my normalized counts organized by row (Ensembl ID's) as a R df called DSNorm. See my code below. One...be helpful? Thanks all!
Code should be placed in three backticks as shown below
```r
#load DESeq2
library("DESeq2")
#read in the above data
cts <- rea…
Hello,
I am relatively new to using DESeq2 so thank you in advanced for your help!
I have multiple samples (A, B, C, D, E) with three replicates of each sample. I would like...Here is the code I'm currently using.
```
ountData <- as.matrix((read.csv("counts.csv", row.names="gene_id")))
colData=data.frame(row.names = colnames(countData),sampleCondition=rep(c("E","E", "E", "F","F"…
updated 5.6 years ago • NG
div class="preformatted">Good morning,
I am trying to run DESeq2 with this design formula
design <- (~batch+sex+tissue)
This is what I do from a count matrix:
>library(DESeq2)
>countTable...lt;- read.table('matrix.txt', header=TRUE, row.names=1)
>ConDesign<- data.frame(row.names = colnames(countTable),
batch = factor(c("1", "2", "3"…
I am doing differential expression using DESeq2. I have 2 conditions (Veh_0h, treatment_56h) and 5 samples each.
After generating the gene count matrix from StringTie...I am doing differential expression using DESeq2. I have 2 conditions (Veh_0h, treatment_56h) and 5 samples each.
After generating the gene count matrix from StringTie output files and loading into R as `countData`, I used the fo…
updated 4.8 years ago • nattzy94
I am using DESeq2 to analysis rna-seq data with 8 biological replicates, which are paired samples. These samples are of primary cells...600px; width:630px"/>
Here is my code:
<pre>
x <- read.table("filt_counts.txt", header=T, row.names=1)
subjects=factor(c(rep(1:8, each=2)))
treat <- as.factor(rep(c("High","Low"),8))
colData <- data.frame(colnames(x),subjec…
Hi everyone,
I am using ``DESeq2`` (v1.26.0) to find DEGs between ``group1`` and ``group2``, and also between`` group2`` and ``group3``, respectly. And there are ten replicate...will be helpful for me.
Here is my code:
```
colData <- read.table('group.txt', header=TRUE, row.names=1)
readscount <- read.table('readscount.txt', header=TRUE, row.names=1)
condition <- fac…
updated 5.0 years ago • Na Chen
Hello,
I am new here.
I have certain numbers of duplicates for each Probe ID (approx. 45000) in one column and its samples data values in others 'n' columns(approx.. 230 columns
updated 3.3 years ago • biomyscience
that information on following website.[http://bioconductor.org/packages/devel/bioc/vignettes/DESeq2/inst/doc/DESeq2.html#why-un-normalized-counts][1]
But, I want to try normalizing raw count into RLE(Relative Log Expression...for your help.
```r
# setup
condition <- factor(c(c,c,c,t,t,t))
coldata <- data.frame(row.names = colnames(Counts),condition)
dds <- DESeqDataSe…
updated 2.0 years ago • SH
<div class="preformatted">Hello,
There seems to be a problem with read.ilmn() with the new version of
limma. When trying to read the control file, there is an error that
row names are not unique.
x <- read.ilmn(files="probe-profile.txt", ctrlfiles="JAlverdy-JDeF-
Dec2-13-ControlProbeData-FinalReport.txt", other.columns="Detection?)
Error in `row.names<-.data.frame`(`*tmp*`…
updated 11.2 years ago • Scott Christley
Hello all,
I tried searching for this specific problem but could not find the answer. I will happily accept a link to a previous question if such exists.
To keep this simple, I can go into more details in response, I am using DESeq2 to analyze the effect of heat stress on expression in two closely-related fish species. I want to determine if the species vary in their respons…
updated 7.2 years ago • kevin.kingsland
Hello,
I have some questions about running a likelihood ratio test in the DESeq2 package in R. My experiment has two factors: sex and population, with 8-12 replicates per condition. I wish to identify...nbsp;
countData <- read.table ("Counttable.txt", header=TRUE)
stickleDesign = data.frame (
row.names = colnames(countData),
sex = c("M", "F", "M", cont... )
pop = c("A", "A…
updated 9.2 years ago • newsomew13
Recent ...
Replies
Comment: DESeq2 - human data, weird results, troubleshooting
by
bd2000
▴
30
I didn't include all of my code the post would be too long. Here's a link to my GitHub if anyone would like to see the full code:
[Full Rst…
Answer: DESeq: Gene is not flagged as differentially expressed despite count data showin
by
James W. MacDonald
68k
The simple answer is that the within-group variability is much higher for the first gene as compared to the second. The Wald statistic is j…
Comment: DESeq: Gene is not flagged as differentially expressed despite count data showin
by
ATpoint
★
4.8k
> since there doesn't seem to be as much count variability,
No? Your CD45 samples range on normalized scale from 2 to 3255. That's 3 orde…
Comment: Ensembl site unresponsive via biomaRt
by
darylruiz6545
•
0
Thank you so much for sharing the link.
Comment: compute_fasta_digest does not match Salmon index hash
by
Alex
•
0
I use the `fastaDigest` tool from [COMBINE-lab/fastaDigest](https://github.com/COMBINE-lab/FastaDigest/), which has the function `compute_f…
Votes
Awards
• All
Popular Question to
bd2000
▴
30
Popular Question to
shepherl
4.1k
Scholar to
Michael Love
43k
Popular Question to
Rob
•
0
Popular Question to
Gordon Smyth
52k
Locations
• All
UPF, Barcelona, Spain, 3 minutes ago
United Kingdom, 50 minutes ago
United Kingdom, 1 hour ago
United States, 1 hour ago
WEHI, Melbourne, Australia, 3 hours ago
UPF, Barcelona, Spain, 3 minutes ago
United Kingdom, 50 minutes ago
United Kingdom, 1 hour ago
United States, 1 hour ago
WEHI, Melbourne, Australia, 3 hours ago
Traffic: 682 users visited in the last hour
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by the
version 2.3.6
version 2.3.6