R/polymorphic.R

In msap: Statistical analysis for Methylation-sensitive Amplification Polymorphism data

#polymorphic.R
# Andrés Pérez-Figueroa

# From a binary matrix returns a list with indices of those polymorphic columns (i.e. at least two ocurrences for every state)
polymorphic <- function(col){
	tcol <- table(col)
	if(length(tcol)>2){
		cat("ERROR! data matrix in 'polymorphic' is not binary\n")
	} 
	if(length(tcol)==1) return(FALSE)
	if(length(tcol)==2){
		if(min(tcol)>0) return (TRUE)
		else return (FALSE)
	}
}