R/filterXY.R

In wateRmelon: Illumina 450 methylation array normalization and metrics

filterXY <-
function(methLumi_data){

	probeXY <- fData(methLumi_data)$TargetID[ which(fData(methLumi_data)$CHR=="Y" | fData(methLumi_data)$CHR=="X") ]
	indexXY <- which(is.element(featureNames(methLumi_data), probeXY))
	rm(probeXY)
	if(length(indexXY) > 0) methLumi_data <- methLumi_data[-indexXY,]
	
	print(paste("Nb probes on X & Y chr to remove: ", length(indexXY), sep=""))
	
	return(methLumi_data)
}