Rについて質問です。
初心者で申し訳ございません。
現在、腫瘍組織と正常組織を比較して、各遺伝子の発現量に差があるかを、t検定しております。
それに際してのエラーが解決できず、質問させていただきます。
func<-function(x){
result<-t.test(tumor[,x],normal[,x])
pvalue<-result$p.value
return(pvalue)
}
P_value_tt <- lapply(colnames(tumor[1:ncol(tumor)]),func)
上記Scriptを実行すると
下記のエラーが出てきます。
t.test.default(tumor[, x], normal[, x]) でエラー:
data are essentially constant
Called from: t.test.default(tumor[, x], normal[, x])
下記のようにfuncのxを仮に1に変更し実行すると、下記の通り、ある遺伝子のp値はしっかり出てきます。
funcの定義も、その後のp_value_ttの記述も問題ないと考えていますが
エラーがでる原因が分からず、難渋しております。
ご教授いただけましたら幸いです。よろしくお願いします。
result<-t.test(tumor[,1],normal3[,1])
pvalue<-result$p.value
> pvalue
[1] 0.6366935
sd(tumor[,x] - normal[,x])
が0
か、ほぼ0
の微小値になる場合がある、ということかと思います。