just trying to launching my DE-seq2 by csv but isn't work...
1
0
Entering edit mode
3.0 years ago
Lily • 0

hello sir, im korean and currently work in dpt, bioenginearing... first of all, my major is basically experiment and im suddenly transferred into Industry-University Cooperation Foundation. its my first time coding r script. so i have to run DE-seq2 so that make my MA-plot(and have no idea any of r script knowledgement) ive tried so many time of this link : https://lashlock.github.io/compbio/R_presentation.html

and my code is right here;

library( "DESeq2" )
library(ggplot2)

countsName <- "http://bioconnector.org/workshops/data/airway_scaledcounts.csv"
download.file(countsName, destfile = "airway_scaledcounts.csv", method = "auto")

countData <- read.csv('airway_scaledcounts.csv', header = TRUE, sep = ",")
head(countData)

metaDataName <- "http://bioconnector.org/workshops/data/airway_metadata.csv"
download.file(metaDataName, destfile = "airway_metadata.csv", method = "auto")

metaData <- read.csv('airway_metadata.csv', header = TRUE, sep = ",")
metaData

dds <- DESeqDataSetFromMatrix(countData=countData, 
                              colData=metaData, 
                              design=~dex)
dds <- DESeq(dds)

?DESeq

res <- results(dds)
head(results(dds, tidy=TRUE)) #let's look at the results table

and error code emerged, but do not know what's going wrong...

Error in DESeqDataSetFromMatrix(countData = countData, colData = metaData,  : 
  ncol(countData) == nrow(colData) is not TRUE

in korea server, i cannot find how to fix this error code and i wrote this one with my lack of english skill. please help me (ㅠㅁㅠ)

DE-SEQ2 R DE-seq2 • 1.7k views
ADD COMMENT
0
Entering edit mode

I think this tutorial is deprecated as the files do not exist anymore

ADD REPLY
0
Entering edit mode

Number of samples in the data and metadata is not matching.

ADD REPLY
0
Entering edit mode

thank you so much...! have to find any other code to run TAT

ADD REPLY
0
Entering edit mode
3.0 years ago
fracarb8 ★ 1.7k

The error means that the order or the row of the metadata and the columns in your counts are not matching. DESEq2 expect the order to be the same as it doesn't do the checking and reordering for you. Order the two tables so that row 1 of the metadata table is the same as the column 1 in the count matrix, and so on.

Have a look at the current DESEq2 documentation section Count matrix input.

ADD COMMENT
0
Entering edit mode

thannk you so so much!! appreciate it!! seems better to try fixing that dds code! thank you!

ADD REPLY

Login before adding your answer.

Traffic: 1608 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6