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
The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben Martínez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...-
Channel: Articles
11-06-2024, 07:24 PM -
-
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,...-
Channel: Articles
10-18-2024, 07:11 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Today, 11:09 AM
|
0 responses
15 views
0 likes
|
Last Post
by seqadmin
Today, 11:09 AM
|
||
Started by seqadmin, Today, 06:13 AM
|
0 responses
17 views
0 likes
|
Last Post
by seqadmin
Today, 06:13 AM
|
||
Started by seqadmin, 11-01-2024, 06:09 AM
|
0 responses
30 views
0 likes
|
Last Post
by seqadmin
11-01-2024, 06:09 AM
|
||
New Model Aims to Explain Polygenic Diseases by Connecting Genomic Mutations and Regulatory Networks
by seqadmin
Started by seqadmin, 10-30-2024, 05:31 AM
|
0 responses
21 views
0 likes
|
Last Post
by seqadmin
10-30-2024, 05:31 AM
|
Leave a comment: