Question: Differential expression analysis with DESeq2
0
gravatar for rhasanvandj
5 weeks ago by
rhasanvandj10
rhasanvandj10 wrote:

Hi all I am using HT-Seq count data to analyze differential expression with DESeq2 in R. I have about 60000 genes(all coding+non-coding genes) and 159 patients (two group, 0=66, 1=93). *Is this ok to use both coding and non-coding genes? or I should remove non-coding genes first? *Also, with DESeq2 I get this warning message. Is this effective on result or this is fine to have this warning? And how can I solve this?

rdata <- read.table("data.txt", header = TRUE, row.names = 1)

    library(DESeq2)

## Create metadata 
sample_org <- data.frame(row.names = colnames(rdata), c(rep("0", 66), rep("1", 93)))
colnames(sample_org) <- c("Group")

dds <- DESeqDataSetFromMatrix(countData = rdata,
                              colData = sample_org,
                              design = ~Group)


dd <- DESeq(dds)
res <- results(dd)
#Saving results
write.csv(res,"res.csv")

in this part I get warning:

dds <- DESeqDataSetFromMatrix(countData = rdata,
                         colData = sample_org,
                           design = ~Group)

The warning message is:

 Warning message:
In DESeqDataSet(se, design = design, ignoreRank) :
  some variables in design formula are characters, converting to factors
rna-seq • 168 views
ADD COMMENTlink modified 5 weeks ago by Spencer20 • written 5 weeks ago by rhasanvandj10
4
gravatar for swbarnes2
5 weeks ago by
swbarnes28.9k
United States
swbarnes28.9k wrote:

Don't worry about the warning. Unless you want to use the metadata numerically, you want them to be factors.

ADD COMMENTlink written 5 weeks ago by swbarnes28.9k
1
gravatar for Spencer
5 weeks ago by
Spencer20
Spencer20 wrote:

just ignoring it. in most cases we need to set stringsAsFactors False to prevent unexpected error

ADD COMMENTlink written 5 weeks ago by Spencer20
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2172 users visited in the last hour