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 (ㅠㅁㅠ)