Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • R script help

    Hey everyone,

    I am working on an R script that will pull contig ID and ranges of differentially expressed transcripts out of the genomic reference I used as input for Cuffdiff. I have the script working, but when it runs into transcripts with the same contig ID (Different ranges) it will only count the first one it encounters and skips the redunancy.

    The format of my contig name and range is formatted like this:

    NAME START END
    gi340012932gbAFJA01012285.1 0 524
    gi340012938gbAFJA01012279.1 10334 12405
    gi340012938gbAFJA01012279.1 2260 7496
    gi340012938gbAFJA01012279.1 69 2164
    gi340012938gbAFJA01012279.1 9069 10084
    gi340013005gbAFJA01012212.1 5696 9927
    gi340013006gbAFJA01012211.1 4079 11139

    As you can see, there is redundancy in the contig ID, but the ranges are different.



    I am getting warning messages looking something like this:

    Warning messages:
    1: In St:Ed : numerical expression has 2 elements: only the first used
    2: In St:Ed : numerical expression has 2 elements: only the first used
    3: In St:Ed : numerical expression has 2 elements: only the first used
    4: In St:Ed : numerical expression has 2 elements: only the first used

    I have copied my script here:

    setwd("C:/Users/Insect1/Documents/Alex/fasta/Treatment_Specific_Analysis")
    files <- list.files(".") # Will list all your files
    library("seqinr") # load the seqinr package
    dnafile <- ("Mrotcontig.fsa") # change the fasta file into an R object
    data <- read.csv("R_Index_N6NF_Apr_9_2013.csv") #reads the csv file "Index" into R
    myseq <- read.fasta(file = dnafile) # creates a fasta file
    for(i in 1:length(data$NAME)){seqname <-data$NAME[i]
    #seqnameB <- noquote(seqname)
    seqB <- getSequence(eval(parse(text=paste("myseq$", seqname, sep = "")))) # the eval and paste functions are needed to construct the search oommand "myseq$MORT00010950" without spaces or quotes mark which make the command nonfunctional
    St <- data[which(data$NAME == seqname), 2]
    Ed <- data[which(data$NAME == seqname), 3]
    trimmedseq <- seqB[St:Ed]
    seqnameC <- paste(seqname,St,"-",Ed)
    write.fasta(sequences = trimmedseq, names = seqnameC, file.out = "Blast.txt", open = "a")
    }
    Blast <-read.fasta("blast.txt", as.string = TRUE)

    Is there an easy fix for making the script address the redundancy in the contig IDs?

    Thanks,
    Alex

  • #2
    Hi Bioinformaticsnewb,
    My experience in R is below standard. Intuitively, what I'll do which you can try while awaiting a better suggestion would be to rename the ID like so:
    gi340012932gbAFJA01012285.1_0-524
    gi340012938gbAFJA01012279.1-10334-12405
    gi340012938gbAFJA01012279.1_2260-7496

    This will mitigate the effect of redundancy.

    HTH

    Comment


    • #3
      R script help

      Hi,

      I'm not sure that I follow what you're trying to do, but

      Is your 'data' object a data frame, and do columns 2 and 3 have names?

      Does the 'R_Index_ ..' file contain the data that you showed an example of with
      columns for NAME START END ?

      One of the components of 'data' appears to be called 'NAME', because
      you refer to data$NAME[i] at the beginning of the for loop.

      I would try changing

      Code:
      St <- data[which(data$NAME == seqname), 2]
      Ed <- data[which(data$NAME == seqname), 3]
      to

      Code:
      St <-data$START[i]
      Ed <- data$END[i]
      but that depends on what is in your 'index' file.

      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, Yesterday, 06:37 PM
      0 responses
      10 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, Yesterday, 06:07 PM
      0 responses
      9 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-22-2024, 10:03 AM
      0 responses
      49 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-21-2024, 07:32 AM
      0 responses
      67 views
      0 likes
      Last Post seqadmin  
      Working...
      X