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 (ㅠㅁㅠ)
I think this tutorial is deprecated as the files do not exist anymore
Number of samples in the data and metadata is not matching.
thank you so much...! have to find any other code to run TAT