estimate.d.RDB: Estimate d using relative distance to boundary (RDB) methods
In Stickbreaker: Fits Stickbreaking, Multiplicative and Additive Models to Data

Description Usage Arguments Details Value Examples

Estimate d using relative distance to boundary (RDB) methods

1	estimate.d.RDB(geno.matrix, fit.matrix, no.est = -100)

`geno.matrix`	Genotype matrix generated in `generate.geno.matrix` or read in.
`fit.matrix`	Fitness matrix generated in `sim.stick.data` or read in.
`no.est`	What to return when no estimate is obtained. Default is -100. Throws warning unless a number of NA is used.

The method calculates RDB for each genotype and its complement. The d.hat.RDB.all indicate the genotype pair that produces the estimate.

List:
[[1]] d.hat Is the best RDB estimate. It is the median of the positive individual RDB values.
[[2]] d.hat.RDB.all RDB estimate based on all individual estimates (not only postie values).
[[3]] d.hat.RDB.ind Vector of individual RDB values (see desalt).
[[4]] d.hat.RDB.other has alternative ways to combine the individual RDB estimates (i.e. median, mean, all values, postie values only).
d.hat.median.pos is preferred estimator (d.hat above) from the median of positive individual values.
d.hat.median is median of all values.
d.hat.mean is mean of all values.
d.hat.mean.pos is mean of positive values.

n.muts <- length(Khan.data[1,])-1
geno.matrix <- Khan.data[,seq(1, n.muts)]
fit.matrix <- as.matrix(Khan.data[,(n.muts+1)])
estimate.d.RDB(geno.matrix, fit.matrix,-100)