Learn R Programming

seqTools (version 1.6.0)

sim_fq: sim_fq: Performs an experimental series of separation capabilities of hierarchical clustering (HC) based on DNA k-mers in FASTQ files using simulated DNA content.

Description

Writes compressed FASTQ files where sequence sections contain concatenated k-mers which are uniformly distributed in the range of k-mers for given k. The function first writes a batch of randomly FASTQ files containing randomly simulated DNA sequence. In a second step the function repeatedly writes FASTQ files with random DNA sequence where a fraction of the reads is 'contaminated' with given DNA k-mers. In a third step, for each set of simulated and contaminated files, a hierarchical cluster (HC) tree based on DNA k-mers is calculated. For each set of files, the size of the smaller fraction in the first half of the tree is counted (perc). The value can be used as measure for separation capability of the HC algorithm.

Usage

sim_fq(nRep=2, nContamVec=c(100, 1000), grSize=20, nSeq=1e4, k=6, kIndex=1365, pos=20)

Arguments

nRep
numeric. Number of replicates for each combination of each nContamVec value
nContamVec
numeric. Vector with nContam (absolute number of contaminated reads) values.
grSize
numeric. Number FASTQ files in control and contamination group.
nSeq
numeric. Number of reads per FASTQ file.
k
numeric. k value used in fastqq function.
kIndex
numeric. k-mer index of inserted k-mer(s). The k-mer index can be retreaved for a given k-mer with 'kMerIndex'. Default value is 1365 (="CCCCCC").
pos
numeric. Determines at which position in sequence the k-mer is inserted. 1-based (1=first position).

Value

containing results of the counted perc values for each repetition of the simulation.

Details

The function is intended to be used as explorative tool (not for routine quality assessment). There are some files written and there will be a lot of output on the terminal. It is therefore recommended to switch to a separate working directory and to run this function on a separate terminal. The function is not exported.

Examples

Run this code
kMerIndex("CCCCCC")
## Not run: res <- seqTools:::sim_fq(nRep=2, nContamVec=c(10, 100),
#                     grSize=4, nSeq=1e2)## End(Not run)

Run the code above in your browser using DataLab