simulateRFLPdata {RFLPtools}

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

N

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.

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.

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.

Author(s)

Mohammed Aslam Imtiaz, Matthias Kohl [e-mail]

See Also

RFLPdata, RFLPref

Examples

  simData <- simulateRFLPdata()