0

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
1
  • "data are essentially constant" と出ていますので、sd(tumor[,x] - normal[,x])0 か、ほぼ 0 の微小値になる場合がある、ということかと思います。
    – metropolis
    2022年10月5日 3:11

0

この質問に回答するには、ログインする必要があります。

のタグが付いた他の質問を参照する。