Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • error in DESeq analysis

    Hello,

    I have RNA-seq data, which I am trying to analyze with DESeq. My file (.csv) appears to be correct

    >head(myfile)
    "
    gene_id VZ_w13 VZ_w14a VZ_w14b VZ_w15a VZ_w15b VZ_w16a
    1 ENSG00000253101 0 0 0 0 0 0
    ..."

    However, when I try

    >cds <- newCountDataSet(myfile,conds)

    I get the following error message:

    "Error in newCountDataSet(myfile, conds) : The countData is not integer."

    The problem, as far as I can tell, is that my data are numerical, not integer, because when I run

    >str(myfile)
    "'data.frame': 53507 obs. of 14 variables:
    $ VZ_w13 : num 0 0 0 0 8 0 0 0 0 0 ..."

    Does anyone have a way to make my data integer, not numerical? As you can see, the data are in fact integers.

    Thanks,
    Stephan

  • #2
    What about as.integer

    If that is really the problem then applying as.integer to your columns should work (note that: as.integer(12.9) == 12)

    Cheers,
    Paul

    "You are only young once, but you can stay immature indefinitely."

    Comment


    • #3
      Excuse my ignorance. How do I apply as.integer to my columns?

      Comment


      • #4
        Example

        Hi,

        here a toy example, which you can just transfer to your data-set:

        Code:
        > data.frame(SomeValue=seq(0,2,0.4))
          SomeValue
        1       0.0
        2       0.4
        3       0.8
        4       1.2
        5       1.6
        6       2.0
        > the.data <- data.frame(SomeValue=seq(0,2,0.4))
        > the.data$SomeValue
        [1] 0.0 0.4 0.8 1.2 1.6 2.0
        > as.integer(the.data$SomeValue)
        [1] 0 0 0 1 1 2
        > the.data$SomeValue <- as.integer(the.data$SomeValue)
        > the.data$SomeValue
        [1] 0 0 0 1 1 2
        Basically you just apply the function to the column and then replace the column with the result. If you have many columns you might want to use the apply function, which you would do like this:
        Code:
        > the.data <- data.frame(SomeValue=seq(0,2,0.4), AnotherValue=seq(4,6,0.4))
        > the.data
          SomeValue AnotherValue
        1       0.0          4.0
        2       0.4          4.4
        3       0.8          4.8
        4       1.2          5.2
        5       1.6          5.6
        6       2.0          6.0
        > apply( the.data, 2, as.integer )
             SomeValue AnotherValue
        [1,]         0            4
        [2,]         0            4
        [3,]         0            4
        [4,]         1            5
        [5,]         1            5
        [6,]         2            6
        When in doubt, use the R-help system to find out about things.

        Cheers,
        Paul

        "You are only young once, but you can stay immature indefinitely."

        Comment


        • #5
          Hi Stephen, looks like the first column of your data are the geneIDs (not integers). Just remove this column and it should be ok. You could use the row.names option when reading the file.

          Comment


          • #6
            Got it working now. Thanks for your help.

            Comment

            Latest Articles

            Collapse

            • seqadmin
              Strategies for Sequencing Challenging Samples
              by seqadmin


              Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
              03-22-2024, 06:39 AM
            • seqadmin
              Techniques and Challenges in Conservation Genomics
              by seqadmin



              The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

              Avian Conservation
              Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
              03-08-2024, 10:41 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, 03-27-2024, 06:37 PM
            0 responses
            13 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-27-2024, 06:07 PM
            0 responses
            11 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-22-2024, 10:03 AM
            0 responses
            53 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-21-2024, 07:32 AM
            0 responses
            69 views
            0 likes
            Last Post seqadmin  
            Working...
            X