Browse code

bug fix in check for NA's in calculateRTheta

git-svn-id: file:///home/git/hedgehog.fhcrc.org/bioconductor/trunk/madman/Rpacks/crlmm@59945 bc3139a8-67e5-0310-9ffc-ced21a209358

Rob Scharp authored on 01/11/2011 16:32:08
Showing1 changed files

... ...
@@ -2777,7 +2777,7 @@ calculateRTheta <- function(cnSet, genotype=c("AA", "AB", "BB"), batch.name){
2777 2777
 	x[x < 64] <- 64
2778 2778
 	y[y < 64] <- 64
2779 2779
 	theta[, "theta"] <- atan2(y, x)*2/pi
2780
-	if(any(is.na(theta))){
2780
+	if(any(is.na(theta[, "theta"]))){
2781 2781
 		stop("NA's in thetas.  Can not compute theta / R values for batches with fewer than 10 samples")
2782 2782
 	}
2783 2783
 	## so that 'R' is available for NP probes