2017-02-28 38 views
0

我目前正在試圖運行既與X數據集和Y數據集Spearman秩相關檢驗,均含有3行(個人)。我可以用cor()運行spearman並獲取值,它們全部是:-1,-.5,-5或1.這對我來說看起來不太合適。我沒有任何0數據集。然而,當我使用rcorr(),它給了我一個錯誤:斯皮爾曼等級相關與少於4觀察?

Error in rcorr(BPT2, y = FunT2, type = "spearman") : 
must have >4 observations 

我比較細菌門在腸道(我X)代謝讀出(我Y

所以我的問題:

  1. 是我得到準確的結果?

  2. 我應該使用在首位Spearman等級相關?

謝謝!

+0

如果你想幫助排除故障代碼,您應該提供問題的[重複的例子(http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example)。如果你問的是一個正確的統計分析,那麼這類問題屬於[stats.se],而不是Stack Overflow(如果你曾經使用統計標籤,這是一個好兆頭,你可能在錯誤的網站作爲標籤描述本身如是說。) – MrFlick

回答

0

你是正確的,只有四種可能的輸出。這是因爲對於3次觀測,秩相關只有6種情況是可能的。

如果我們認爲我們的X常數1:3,有y的6個可能的等級值(逐行讀取明智):

Var1 Var2 Var3 
1  3 2 1 
2  2 3 1 
3  3 1 2 
4  1 3 2 
5  2 1 3 
6  1 2 3 

當你計算X和Y之間的每個相關,只有以下可能的回報:

apply(df, 1, function(x){cor(1:3, x)}) 

    1 2 3 4 5 6 
-1.0 -0.5 -0.5 0.5 0.5 1.0 

由於兩個重複,您會得到4個值。

所以 - 這是數學上可以計算,但它不是在描述的分佈是非常有用的。