simulateRFLPdata: Simulate RFLP data.

Description

Simulates RFLP data for comparions of algorithms.

Usage

simulateRFLPdata(N = 10, nrBands = 3:12, bandCenters = seq(100, 800, by = 100),
                 delta = 50, refData = FALSE)

Arguments

integer: number samples which shall be simulated per number of bands.

nrBands

integer: vector of number of bands.

bandCenters

numeric: vector of band centers.

delta

numeric: uniform distribution with min = bandCenter - delta and max = bandCenter + delta is used.

refData

logical: if TRUE, additonal columns Taxonname and Accesion are generated.

Value

A data frame with N*length(nrBands) observations on the following four variables

Sample: character: sample identifier.
Band: integer: band number.
MW: integer: molecular weight.
Enzyme: character: enzyme name.

is generated. If refData = TRUE then the following two additional variables are added.

Taxonname: character: taxon name.
Accession: character: accession number.

Details

The function can be used to simulate RFLP data. For every number of band specified in nrBands a total number of N samples are generated.

First the band centers are randomly selected (with replacement) from bandCenter which form the centers of intervals of length 2*delta. From these intervals uniform random numbers are drawn leading to randomly generated RFLP data.

Examples

Run this code

# NOT RUN {
  simData <- simulateRFLPdata()
# }