I ran the following loop on my SNP data and got an error message as
indicated
for (i in genenames){
+   current <- fst1[which(fst1$Gene == i),]
+   num <- nrow(current)
+   fst <- max(current$fst)
+   position <- mean(current$pos)
+   nposition <- mean(current$newpos)
+   numhigh <- nrow(current[which(current$fst > threshold),])
+   colors <- mean(current$colors)
+   output <- matrix(NA,nrow=1,ncol=8)
+   numthresh <- paste("# SNPs > Fst = ", threshold, sep="")
+   colnames(output) <- c("gene", "gene_old", "pos", "newpos", "# Snps",
numthresh, "Max.Fst", "colors")
+   output[1,1] <- i
+   output[1,2] <- as.character(current[1, "gene_old"])
+   output[1,3] <- position
+   output[1,4] <- nposition
+   output[1,5] <- num
+   output[1,6] <- numhigh
+   output[1,7] <- fst
+   output[1,8] <- colors
+   maxfstgene <- rbind(maxfstgene, output)
+ }
Error in output[1, 2] <- as.character(current[1, "gene_old"]) :
  replacement has length zero
In addition: Warning message:
In mean.default(current$pos) :
  argument is not numeric or logical: returning NA
--

-- 
* "The information contained in this email and any attachments may be 
legally privileged and confidential. If you are not an intended recipient, 
you are hereby notified that any dissemination, distribution, or copying of 
this e-mail is strictly prohibited. If you have received this e-mail in 
error, please notify the sender and permanently delete the e-mail and any 
attachments immediately. You should not retain, copy or use this e-mail or 
any attachments for any purpose, nor disclose all or any part of the 
contents to any other person."*

	[[alternative HTML version deleted]]


