Question: data frame from deseq2
0
gravatar for krushnach80
8 months ago by
krushnach80280
krushnach80280 wrote:

I m plotting heatmap, based on padj value this is the code im using .

padj <- 0.05
res.index <- which(res$padj < padj)
df <- as.data.frame(
  colData(dds[res.index,])[,c("condition")])
rownames(df) <- colnames(assay(rld[res.index,]))

library(pheatmap)
pheatmap(assay(rld[res.index,]),
         cluster_rows=TRUE,
         show_rownames=FALSE, cluster_cols=TRUE,
         treeheight_col=1,
         scale="row",
         annotation_col=df)

Im getting the heatmap no issues but i want to get the value which im using for my plot , how i get that, any suggestion or help would be highly appreciated .

rna-seq R • 459 views
ADD COMMENTlink written 8 months ago by krushnach80280
1

Do you mean the rld normalized count table?

assay(rld[res.index,]
ADD REPLYlink written 8 months ago by VHahaut1.1k

yes..how do i get that?

ADD REPLYlink written 8 months ago by krushnach80280
1
gravatar for VHahaut
8 months ago by
VHahaut1.1k
Belgium
VHahaut1.1k wrote:

Ok so you mean you want to save it as a text file somewhere?

write.table(assay(rld[res.index,]), sep="\t", path, row.names=F, quote=F)

where path should be replaced by the full path where you want to save it such as:

    write.table(assay(rld[res.index,]), sep="\t", "~/Desktop/rld_normalized_count.txt", row.names=F, quote=F)
ADD COMMENTlink written 8 months ago by VHahaut1.1k

thanks a lot...i was looking to get the data so that i can use it for other analysis

ADD REPLYlink written 8 months ago by krushnach80280
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: 1554 users visited in the last hour