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
          Choosing Between NGS and qPCR
          by seqadmin



          Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...
          10-18-2024, 07:11 AM
        • seqadmin
          Non-Coding RNA Research and Technologies
          by seqadmin




          Non-coding RNAs (ncRNAs) do not code for proteins but play important roles in numerous cellular processes including gene silencing, developmental pathways, and more. There are numerous types including microRNA (miRNA), long ncRNA (lncRNA), circular RNA (circRNA), and more. In this article, we discuss innovative ncRNA research and explore recent technological advancements that improve the study of ncRNAs.

          Nobel Prize for MicroRNA Discovery
          This week,...
          10-07-2024, 08:07 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 11-01-2024, 06:09 AM
        0 responses
        20 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 10-30-2024, 05:31 AM
        0 responses
        20 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 10-24-2024, 06:58 AM
        0 responses
        24 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 10-23-2024, 08:43 AM
        0 responses
        53 views
        0 likes
        Last Post seqadmin  
        Working...
        X