Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Chayan: see this thread for other options in case the R solution does not work: http://seqanswers.com/forums/showthread.php?t=16505
-
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:
-
@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:
-
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:
-
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 helpsLast edited by pbseq; 01-11-2011, 06:04 AM.
Leave a comment:
-
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:
-
exporting
Ok, followup question...
Once I've done
Code:donkey <- FastqSampler(con, n=1e6, readerBlockSize=1e8, verbose=FALSE)
THanks, and apologies if this obvious.
Chris
Leave a comment:
-
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:
-
i don't think you loaded ShortRead
start R
> R
> install.packages("ShortRead")
> library("ShortRead")
> ?FastqSampler
Leave a comment:
-
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!
ChrisTags: None
Latest Articles
Collapse
-
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,...-
Channel: Articles
10-07-2024, 08:07 AM -
-
by seqadmin
Metagenomics has improved the way researchers study microorganisms across diverse environments. Historically, studying microorganisms relied on culturing them in the lab, a method that limits the investigation of many species since most are unculturable1. Metagenomics overcomes these issues by allowing the study of microorganisms regardless of their ability to be cultured or the environments they inhabit. Over time, the field has evolved, especially with the advent...-
Channel: Articles
09-23-2024, 06:35 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Today, 07:29 AM
|
0 responses
12 views
0 likes
|
Last Post
by seqadmin
Today, 07:29 AM
|
||
Genetic Barcodes and Single-Cell Sequencing Illuminate Tumor Initiation and Chemoresistance in Breast Cancer
by seqadmin
Started by seqadmin, 10-15-2024, 06:35 AM
|
0 responses
11 views
0 likes
|
Last Post
by seqadmin
10-15-2024, 06:35 AM
|
||
Started by seqadmin, 10-14-2024, 02:44 PM
|
0 responses
11 views
0 likes
|
Last Post
by seqadmin
10-14-2024, 02:44 PM
|
||
Started by seqadmin, 10-11-2024, 06:55 AM
|
0 responses
19 views
0 likes
|
Last Post
by seqadmin
10-11-2024, 06:55 AM
|
Leave a comment: