sim_fq: sim_fq: Performs an experimental series of separation...
In seqTools: Analysis of nucleotide, sequence and quality content on fastq files.

Description Usage Arguments Details Value Author(s) Examples

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.

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

`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 'getKmerIndex'. Default value is 1365 (="CCCCCC").
`pos`	`numeric`. Determines at which position in sequence the k-mer is inserted. 1-based (1=first position).

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.

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

Wolfgang Kaisers