1

次の方法を使用して、fiser の正確なテストから p 値を取得しています。ただし、ペアごとに値「n00、n01、n10、n11」を出力する必要があります。以下のように、p値とともにマトリックスではなくテーブルとしてそれを印刷するにはどうすればよいですか?

fish <- function(y, x) {
    n00 = sum((1-x)*(1-y))
                       n01 = sum((1-x)*y)
                       n10 = sum(x*(1-y))
                       n11 = sum(x*y)
                       a = matrix(c(n00, n01, n10, n11)
                                  , nrow = 2)
                       pval = fisher.test(a)$p.value
                       return(pval)
}

chiArray <- function(x) apply(mat1, 1, fish, x)
sapply(1:nrow(mat2), function(j) chiArray(mat2[j, ]))
chisq.cna.mut.test <- sapply(1:nrow(mat2), function(j) chiArray(mat2[j, ]))

出力を次のようにしたい：

# name1_mat1 name1_mat2 n00 n01 n10 n11 pvalue
# name1_mat1 name2_mat2 n00 n01 n10 n11 pvalue

入力は (mat1) であり、この例では、mat2 が mat1 と同じであると見なします。

# ID    case1   case2   case3
# name1 0   0   0
# name2 1   0   1
# name3 0   1   1
4

1 に答える 1

0

これによって生成されたデータを使用する:

ID <- paste0("name", 1:5)
mat1 <- 1*cbind(case1=runif(5)>0.5, case2=runif(5)>0.5, case3=runif(5)>0.5)
mat2 <- 1*cbind(case1=runif(5)>0.5, case2=runif(5)>0.5, case3=runif(5)>0.5)
rownames(mat1) <- ID
rownames(mat2) <- ID

nそして、関数を変更して、p 値とともにすべての s を返すようにします。

fish <- function(y, x)
{
    n00 = sum((1-x)*(1-y))
    n01 = sum((1-x)*y)
    n10 = sum(x*(1-y))
    n11 = sum(x*y)
    a = matrix(c(n00, n01, n10, n11), nrow = 2)
    pval = fisher.test(a)$p.value
    c(a, pval)
}

あなたが望むことはこれで達成できると思います：

result <- data.frame(mat1="", mat2="", n00=0, n01=0, n10=0, n11=0, pvalue=0, stringsAsFactors=FALSE)

k <- 1

for(j in 1:nrow(mat2)) for(i in 1:nrow(mat1))
{
    result[k,1] <- paste0(rownames(mat1)[i], "_mat1")
    result[k,2] <- paste0(rownames(mat2)[j], "_mat2")
    result[k,-(1:2)] <- fish(x=mat2[j, ], y=mat1[i,])
    k <- k + 1
}

例：

> result
         mat1       mat2 n00 n01 n10 n11    pvalue
1  name1_mat1 name1_mat2   1   1   0   1 1.0000000
2  name2_mat1 name1_mat2   0   2   0   1 1.0000000
3  name3_mat1 name1_mat2   2   0   1   0 1.0000000
4  name4_mat1 name1_mat2   1   1   1   0 1.0000000
5  name5_mat1 name1_mat2   0   2   1   0 0.3333333
6  name1_mat1 name2_mat2   0   1   1   1 1.0000000
7  name2_mat1 name2_mat2   0   1   0   2 1.0000000
8  name3_mat1 name2_mat2   1   0   2   0 1.0000000
9  name4_mat1 name2_mat2   1   0   1   1 1.0000000
10 name5_mat1 name2_mat2   1   0   0   2 0.3333333
11 name1_mat1 name3_mat2   0   1   1   1 1.0000000
12 name2_mat1 name3_mat2   0   1   0   2 1.0000000
13 name3_mat1 name3_mat2   1   0   2   0 1.0000000
14 name4_mat1 name3_mat2   1   0   1   1 1.0000000
15 name5_mat1 name3_mat2   0   1   1   1 1.0000000
16 name1_mat1 name4_mat2   0   1   1   1 1.0000000
17 name2_mat1 name4_mat2   0   1   0   2 1.0000000
18 name3_mat1 name4_mat2   1   0   2   0 1.0000000
19 name4_mat1 name4_mat2   1   0   1   1 1.0000000
20 name5_mat1 name4_mat2   1   0   0   2 0.3333333
21 name1_mat1 name5_mat2   0   1   1   1 1.0000000
22 name2_mat1 name5_mat2   0   1   0   2 1.0000000
23 name3_mat1 name5_mat2   1   0   2   0 1.0000000
24 name4_mat1 name5_mat2   1   0   1   1 1.0000000
25 name5_mat1 name5_mat2   1   0   0   2 0.3333333
于 2013-06-24T15:18:28.407 に答える