從R中已經子集的向量中找到對應的值？

Question

背景：我將數據集保存為名為D的csv文件（ 請加載到R ）：

D = read.csv("https://docs.google.com/uc?id=0B5V8AyEFBTmXWU40TUZGTURQWlU&export=download")

我使用以下函數從我的數據中獲取2個數量 （ 請提供給您的R ）：

source("https://docs.google.com/uc?id=0B5V8AyEFBTmXWTk0LWhaMkY2b3M&export=download")

這兩個數量的獲得如下：

b = BF.d.pvalue(t = D$t.value, n1 = D$n1, n2 = D$n2)
BF = b[1, ]  ;  p.value = b[2, ]

子集詳細信息：我已將大於.05 p.value及其相應的BF s進行了子集化：

 pvalue.05_1 = p.value[p.value > .05] ;
 BF.pvalue.05_1 = BF[p.value > .05]`

我進一步BF.pvalue.05_1 ，它們在1/10和1/3之間：

BF.pvalue.05_1_.1_.33 = BF.pvalue.05_1[BF.pvalue.05_1 > 1/10 & BF.pvalue.05_1 <= 1/3]

問：現在我想知道怎樣才能找到相應的p.value為BF.pvalue.05_1_.1_.33以上？

Answer 1

執行此操作的首選方法是將數據合並到dataframe ，然后使用subset命令進行過濾。

myDf = data.frame(p = p.value, BF = BF)
head(myDf)
#              p            BF
# 1 2.274873e-06  6.241835e+03
# 2 3.811612e-02  1.736017e+00
# 3 0.000000e+00 2.592434e+147
# 4 0.000000e+00 1.982820e+130
# 5 0.000000e+00  1.315152e+29
# 6 9.992007e-15  4.442134e+11 

現在，無論何時按行對數據進行子集設置 ，您都可以訪問p值和BF值。

firstSubset = subset(myDf, p > .05)
dim(firstSubset)
# [1] 175   2

secondSubset = subset(firstSubset, BF > .1 & BF < 1/3)
dim(secondSubset)
# [1] 76  2

head(secondSubset)
#            p        BF
# 28 0.8518770 0.3131790
# 34 0.9358011 0.2910234
# 35 0.9302671 0.2911639
# 52 0.6825720 0.3101911
# 88 0.7201547 0.2770751
# 96 0.6472360 0.2868055

或者，您可以同時使用兩個條件

secondSubset = subset(myDf, (BF > .1) & (BF < 1/3) & (p > .05))