Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Creating a data.frame for DESeq2

    Hey All,

    I am a novice R user and am having a difficult time getting a data.frame set up for DESeq2. I have done the following:

    datafile <- "E:/Longitudinal Paper/DESeq/0 to 1.txt"

    > datafile

    [1] "E:/Longitudinal Paper/DESeq/0 to 1.txt"

    > countTable = read.delim(datafile, header=TRUE, row.names=1)

    > head(countTable)

    C1 C2 C3 D1 D2 D3 D4
    COX1 244969 257554 243872 235001 237881 178634 202627
    Hsp90ab1 185595 200135 228766 143571 147350 99427 116190
    CYTB 149810 151470 148100 129141 140353 86705 107290
    Eef2 142818 151153 163043 140762 134884 94678 99470
    Col4a1 126366 130972 227447 126424 122709 83403 77633
    ND1 116379 117694 111396 91746 99802 62038 76769

    > ExpDesign = data.frame(row.names=colnames(countTable),
    + condition=c("C1", "C2", "C3", "D1", "D2", "D3", "D4"),
    + libType = c("paired-end","paired-end","paired-end","paired-end","paired-end","paired-end","paired-end"))

    > ExpDesign

    condition libType
    C1 C1 paired-end
    C2 C2 paired-end
    C3 C3 paired-end
    D1 D1 paired-end
    D2 D2 paired-end
    D3 D3 paired-end
    D4 D4 paired-end

    I then try to do the following and get the subsequent errors:

    > data("countTable")
    Warning message:
    In data("countTable") : data set ‘countTable’ not found

    > countData <- counts(countTable)
    Error in (function (classes, fdef, mtable) :
    unable to find an inherited method for function ‘counts’ for signature ‘"data.frame"’

    > colData <- pData(countTable)[,c("condition", "type")]
    Error: could not find function "pData"

    What steps am I missing in between? Thanks for any help!
    Last edited by KHubbard; 10-11-2013, 12:50 PM.

  • #2
    data("countTable") doesn't mean what you think. It tries to load an R dataset called "countTable", not convert your data.frame into a format for DESeq2. In the vignette, that works, since it's just loading a specific dataset that comes with an R package. Instead, just use:

    Code:
    cds <-DESeqDataSetFromMatrix(countData = countTable, colData=ExpDesign$condition, design=~ExpDesign$condition)
    BTW, your experimental design is really just ExpDesign$condition, since they're all paired-end (so you needn't include it). Also, are your groups really unreplicated?

    Comment


    • #3
      Thanks for the reply!

      In answer to your question, all the "C's" are replicates of one another as are the "D's".

      Comment


      • #4
        Ah, good. Then you want instead:
        Code:
        condition=c("C", "C", "C", "D", "D", "D", "D")

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Exploring the Dynamics of the Tumor Microenvironment
          by seqadmin




          The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
          07-08-2024, 03:19 PM
        • seqadmin
          Exploring Human Diversity Through Large-Scale Omics
          by seqadmin


          In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
          06-25-2024, 06:43 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 07-10-2024, 07:30 AM
        0 responses
        26 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 07-03-2024, 09:45 AM
        0 responses
        201 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 07-03-2024, 08:54 AM
        0 responses
        212 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 07-02-2024, 03:00 PM
        0 responses
        193 views
        0 likes
        Last Post seqadmin  
        Working...
        X