R/baseContent.R

In spiderDev: Species Identity and Evolution in R--development

Defines functions baseContent

Documented in baseContent

baseContent <- function(DNAbin){
	probs <- vector()
	for (i in 1:dim(DNAbin)[1]) {
		exBf <- base.freq(DNAbin[-i, ])
		seqBf <- base.freq(DNAbin[i,], freq = TRUE)
		probs[i] <- dmultinom(seqBf, prob = exBf, log = TRUE)
	}
	probs
}

Try the spiderDev package in your browser