checkpeaks {KinMixLite}

Check and modify database to have positive frequencies for all observed peaks/alleles

Description

Check whether database has positive frequencies for all peaks/alleles observed in epg and genotype profiles, and optionally modify db by addition of small positive frequencies so that it does, followed by renormalisation of frequencies for each allele to sum to 1.

Usage

checkpeaks(x,db,fix=0.003)

Arguments

x

data frame, the epg or genotype profile; see formats.

db

data frame, the db; see formats.

fix

numeric: if positive, increment to db frequency for each identified discrepant peak

Value

(possibly modified) db

Author(s)

Peter Green ([e-mail])

Examples

data(test2data)

db<-checkpeaks(epg,db)
db<-checkpeaks(Cgt,db)

Xgt<-data.frame(marker=c('D10','D12'),allele1=c(8,13),allele2=c(13,10))
db<-checkpeaks(Xgt,db)
db