Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • GenoMax
    replied
    Chayan: see this thread for other options in case the R solution does not work: http://seqanswers.com/forums/showthread.php?t=16505

    Leave a comment:


  • chayan
    replied
    i need to subsampling my fastq file of Iontorrent shotgun reads. I install R-2.15 and going by the previous suggesitios in this thread i tried to install ShortRead by > install.packages("ShortRead") but get the folowing error
    Warning message:
    In getDependencies(pkgs, dependencies, available, lib) :
    package ‘Shortread’ is not available

    ne help regarding this and if any other ways to random sampling of my fastq file. Thanx in advannce

    Regards
    Chayan

    Leave a comment:


  • westerman
    replied
    @shanebrubaker:

    I suspect that you will get a better response if you ask your questing in a different thread with a title that contains "Digital Normalization" instead of burying it inside this thread. As for the answer to your question, I do not know. I am just starting to explore the program and may find out the answer later today.

    Leave a comment:


  • shanebrubaker
    replied
    On a related note, I was trying to reduce my data size by using khmer, from a paper on Digitial Normalization. Does anyone have experience with this tool? I am getting an error when using it where it says I have no paired reads ... but I do. Are there simliar tools?

    Leave a comment:


  • pbseq
    replied
    Hi chrisbala,
    in ShortRead package ( manual :
    http://www.bioconductor.org/packages.../ShortRead.pdf
    )

    there is writeFastq method that looks like what you need:
    from manual:

    writeFastq signature(object = "ShortReadQ", file = "character", mode="character",
    ...): Write object to file in fastq format. mode defaults to ‘w’. This creates a new
    file, or fails if file already exists. Use mode="a" to append to an existing file. file is
    expanded using path.expand.


    hope it helps
    Last edited by pbseq; 01-11-2011, 06:04 AM.

    Leave a comment:


  • ewilbanks
    replied
    You can also use the mothur package which has a similar subsampling feature. The examples shows fasta only but I think it can handle paired fasta/qual files.

    Leave a comment:


  • chrisbala
    replied
    exporting

    Ok, followup question...

    Once I've done

    Code:
    donkey <- FastqSampler(con, n=1e6, readerBlockSize=1e8, verbose=FALSE)
    How do I write the sampled sequences to a fastq file?

    THanks, and apologies if this obvious.

    Chris

    Leave a comment:


  • chrisbala
    replied
    haha, I'm a bit of a newbie but not quite THAT bad...

    I was testing it out on my Mac/GUI version of R. That does not seem to work.

    Linux version seems fine. Maybe I'm still missing something ... but thats my guess as to what was going on.

    Thanks though!

    Leave a comment:


  • Zigster
    replied
    i don't think you loaded ShortRead
    start R
    > R
    in R:
    > install.packages("ShortRead")
    > library("ShortRead")
    > ?FastqSampler

    Leave a comment:


  • chrisbala
    started a topic Subsampling Fastq

    Subsampling Fastq

    Anyone know an easy way to randomly subsample from a fastq (Illumina data)?

    ShortReads seems to have a function FastqSampler, but I can't seem to make it work (?FastqSampler gives No documentation for 'FastqSampler' in specified packages and libraries:
    you could try '??FastqSampler')

    Thanks!

    Chris

Latest Articles

Collapse

  • seqadmin
    Recent Advances in Sequencing Analysis Tools
    by seqadmin


    The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
    05-06-2024, 07:48 AM
  • seqadmin
    Essential Discoveries and Tools in Epitranscriptomics
    by seqadmin




    The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...
    04-22-2024, 07:01 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 05-14-2024, 07:03 AM
0 responses
19 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-10-2024, 06:35 AM
0 responses
42 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-09-2024, 02:46 PM
0 responses
53 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-07-2024, 06:57 AM
0 responses
42 views
0 likes
Last Post seqadmin  
Working...
X