R/aa.conservation_first.R

In VLF: Frequency Matrix Approach for Assessing Very Low Frequency Variants in Sequence Records

Defines functions aa.conservation_first

Documented in aa.conservation_first

aa.conservation_first <-
function(modal, p, seqlength){
	conserve <- 0
	for(i in 1:seqlength){
		if(modal[i] >= p){
			conserve <- conserve + 1
		}
	}
	return(conserve)
}

Try the VLF package in your browser