我有一系列excel文件,如下所示,其中包含一个应急表:
cns random
H3k 12 682
Not_H3k 343 30222当我使用“gdata”库使用此命令读取该文件时,它将给出下表,它不再是应急表:
random <- read.xls ("analysis.xlsx", sheet = 1, header = TRUE)
X cns random
1 UTR 12 682
2 Not_UTR 343 30222我想对这些数据进行费舍尔测试,我如何将这些data.frame数据重组成一种应变表格式?
发布于 2014-08-18 04:19:42
read.xls允许read.table的任何参数。特别是,您可以指定包含行名的列的数目。所以你可以用
random <- read.xls ("analysis.xlsx", sheet = 1, header = TRUE, row.names=1)然后继续进行测试
发布于 2014-08-18 05:46:02
您可以将fisher测试与您所拥有的dataframe:
random = structure(list(X = structure(c(2L, 1L), .Label = c("Not_UTR",
"UTR"), class = "factor"), cns = c(12L, 343L), random = c(682L,
30222L)), .Names = c("X", "cns", "random"), class = "data.frame", row.names = c(NA,
-2L))
random
X cns random
1 UTR 12 682
2 Not_UTR 343 30222
fisher.test(random[,-1])
Fisher's Exact Test for Count Data
data: ddf[, -1]
p-value = 0.1433
alternative hypothesis: true odds ratio is not equal to 1
95 percent confidence interval:
0.7892877 2.7613943
sample estimates:
odds ratio
1.550284 https://stackoverflow.com/questions/25356212
复制相似问题