R语言做wilcox分析时报错

> outTab=data.frame()
> for(i in row.names(data)){
+ 	geneName=unlist(strsplit(i,"\\|",))[1]
+ 	geneName=gsub("\\/", "_", geneName)
+ 	rt=rbind(expression=data[i,], grade=grade)
+ 	rt=as.matrix(t(rt))
+ 	wilcoxTest=wilcox.test(expression ~ grade, data=rt)
+ 	pvalue=wilcoxTest$p.value
+ 	conGeneMeans=mean(data[i,1:conNum])
+ 	treatGeneMeans=mean(data[i,(conNum+1):ncol(data)])
+ 	logFC=log2(treatGeneMeans)-log2(conGeneMeans)
+ 	conMed=median(data[i,1:conNum])
+ 	treatMed=median(data[i,(conNum+1):ncol(data)])
+ 	diffMed=treatMed-conMed
+ 	if( ((logFC>0) & (diffMed>0)) | ((logFC<0) & (diffMed<0)) ){  
+ 		  outTab=rbind(outTab,cbind(gene=i,conMean=conGeneMeans,treatMean=treatGeneMeans,logFC=logFC,pValue=pvalue))
+ 	 }
+ }
Error in wilcox.test.formula(expression ~ grade, data = rt) : 
  分组因子只能有两层
有没有人救救这个怎么解决呀