- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我正在计算 R 数据框中每行的超几何测试,希望是正确的。
其中第 1 列是基因 (microRNA) 的名称,“Total_mRNA”列是基因组中总共存在多少 mRNA,因此不会改变。 “Total_targets_targets”列是如果所有 mRNA 都存在的话,每个 microRNA 可以靶向的 mRNA 数量。然而,对于这个例子,仅存在“subset_mRNA”(该数字也始终相同),并且其中我知道每个 microRNA 可以针对“subset_targets”有多少个 mRNA。
为了确定每个 microRNA 的目标与背景(总 mRNA 和针对它们的总 microRNA)相比是否富集,我对每行执行超几何测试,如下所示:
phyper(targets-in-subset, targets-in-bkgd, failure-in-bkgd, sample-size-subset, lower.tail= FALSE)
dput(df1)
structure(list(Genes_names = c("microRNA-1", "microRNA-2", "microRNA-3",
"microRNA-4", "microRNA-5", "microRNA-6", "microRNA-7", "microRNA-8",
"microRNA-9", "microRNA-10"), Total_mRNAs = c(61064L, 61064L,
61064L, 61064L, 61064L, 61064L, 61064L, 61064L, 61064L, 61064L
), Total_targets_targets = c(1918L, 7807L, 3969L, 771L, 2850L,
1355L, 1560L, 2478L, 1560L, 2478L), subset_mRNAs = c(17571L,
17571L, 17571L, 17571L, 17571L, 17571L, 17571L, 17571L, 17571L,
17571L), subset_targets = c(544L, 2109L, 1137L, 213L, 793L, 394L,
430L, 686L, 430L, 686L)), class = "data.frame", row.names = c(NA,
-10L))
df1$pvalue <- phyper(df1$subset_targets, df1$Total_targets_targets, df1$Total_mRNAs-df1$Total_targets_targets, df1$subset_mRNAs, lower.tail= FALSE)
现在的问题是 Bonferroni 如何纠正这个值?这个计算理论上正确吗?
最佳答案
如果您没有大量样本,请避免这种痛苦,只需使用 Fisher 测试并使用 p.adjust
进行 bonferroni:
library(broom)
result = lapply(1:nrow(df1),function(i){
not_target_subset = df1$Total_targets_targets[i] - df1$subset_targets[i]
not_subset = df1$Total_mRNAs[i] - df1$subset_mRNAs[i] - not_target_subset
M = cbind(c(df1$subset_targets[i],df1$subset_mRNAs[i]-df1$subset_targets[i]),
c(not_target_subset,not_subset))
res = data.frame(Genes_names=df1$Genes_names[i],
tidy(fisher.test(M,alternative="greater")))
return(res)
})
result= do.call(rbind,result)
result$padj = p.adjust(result$p.value,"bonferroni")
您的超几何代码略有偏差。请注意,您正在进行单侧超几何测试。
您可以查看这个post for how to put the tables into phyper和 this for why you need the -1 。因此我们计算超几何 p 值:
result$hyper_p = with(df1,
phyper(subset_targets-1,subset_mRNAs,Total_mRNAs-subset_mRNAs, Total_targets_targets, lower.tail= FALSE)
)
你可以看到它相符:
Genes_names estimate p.value conf.low conf.high
1 microRNA-1 0.9793710 0.6655527 0.8984025 Inf
2 microRNA-2 0.9047305 0.9998968 0.8647701 Inf
3 microRNA-3 0.9933480 0.5791759 0.9350214 Inf
4 microRNA-4 0.9441864 0.7722712 0.8229140 Inf
5 microRNA-5 0.9520878 0.8789562 0.8863785 Inf
6 microRNA-6 1.0151760 0.4119600 0.9168998 Inf
7 microRNA-7 0.9404619 0.8641420 0.8539585 Inf
8 microRNA-8 0.9454359 0.8942082 0.8756678 Inf
9 microRNA-9 0.9404619 0.8641420 0.8539585 Inf
10 microRNA-10 0.9454359 0.8942082 0.8756678 Inf
method alternative hyper_p
1 Fisher's Exact Test for Count Data greater 0.6655527
2 Fisher's Exact Test for Count Data greater 0.9998968
3 Fisher's Exact Test for Count Data greater 0.5791759
4 Fisher's Exact Test for Count Data greater 0.7722712
5 Fisher's Exact Test for Count Data greater 0.8789562
6 Fisher's Exact Test for Count Data greater 0.4119600
7 Fisher's Exact Test for Count Data greater 0.8641420
8 Fisher's Exact Test for Count Data greater 0.8942082
9 Fisher's Exact Test for Count Data greater 0.8641420
10 Fisher's Exact Test for Count Data greater 0.8942082
关于r - 如何使用 Bonferroni 校正计算数据框中每一行的超几何检验,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/64969842/
SQL 和一般开发的新手,我有一个表(COUNTRIES),其中包含字段(INDEX、NAME、POPULATION、AREA) 通常我添加一个客户端(Delphi)计算字段(DENSITY)和 On
我想使用 calc(100%-100px),但在我的 demo 中不起作用由于高度只接受像素,因此如何将此百分比值转换为像素。 最佳答案 以下将为您提供高度: $(window).height();
我正在尝试在 MySQL 中添加列并动态填充其他列。 例如我有一张表“数字”并具有第 1 列、第 2 列、第 3 列,这些总数应填充在第 4 列中 最佳答案 除非我误解了你的问题,否则你不只是在寻找:
我想返回简单计算的结果,但我不确定如何执行此操作。我的表格如下: SELECT COUNT(fb.engineer_id) AS `total_feedback`, SUM(fb.ra
我一直在尝试做这个程序,但我被卡住了,我仍然是一个初学者,任何帮助将不胜感激。我需要程序来做 打印一个 10 X 10 的表格,其中表格中的每个条目都是行号和列号的总和 包含一个累加器,用于计算所有表
这个计算背后一定有一些逻辑。但我无法得到它。普通数学不会导致这种行为。谁能帮我解释一下原因 printf ("float %f\n", 2/7 * 100.0); 结果打印 1.000000 为什么会
我想计算从 0 到 (n)^{1/2} - 1 的数字的 AND每个数字从 0 到 (n)^{1/2} - 1 .我想在 O(n) 中执行此操作时间,不能使用 XOR、OR、AND 运算。 具体来说,
如何在 Excel 中将公式放入自定义数字格式?例如(出于说明目的随机示例), 假设我有以下数据: 输入 输出 在不编辑单元格中的实际数据的情况下,我想显示单元格中的值除以 2,并保留两位小数: 有没
每次我在 Flutter 应用程序中调用计算()时,我都会看到内存泄漏,据我所知,这基本上只是一种生成隔离的便捷方法。我的应用程序内存占用增加并且在 GC 之后永远不会减少。 我已将我的代码简化为仅调
我有数字特征观察 V1通过 V12用于目标变量 Wavelength .我想计算 Vx 之间的 RMSE列。数据格式如下。 每个变量“Vx”以 5 分钟的间隔进行测量。我想计算所有 Vx 变量的观测值
我正在寻找一种使用 C 语言计算文件中未知字符数的简单方法。谢谢你的帮助 最佳答案 POSIX 方式(可能是您想要的方式): off_t get_file_length( FILE *file ) {
我正在使用 Postgres,并且我正试图围绕如何在连续日期跨度中得出第一个开始日期的问题进行思考。例如 :- ID | Start Date | End Date =================
我有一个订单表格,我在其中使用 jQuery 计算插件来汇总总数。 此求和工作正常,但生成的“总和”存在问题。总之,我希望用逗号替换任何点。 代码的基础是; function ($this) {
我在使用 double 变量计算简单算术方程时遇到问题。 我有一个具有 double 属性 Value 的组件,我将此属性设置为 100。 然后我做一个简单的减法来检查这个值是否真的是 100: va
我在这里看到了一些关于 CRC 32 计算的其他问题。但没有一个让我满意,因此是这样。 openssl 库是否有任何用于计算 CRC32 的 api 支持?我已经在为 SHA1 使用 openssl,
当我在PHP日期计算中遇到问题时,我感到惊讶。 $add = '- 30 days'; echo date('Y-m-01', strtotime($add)); // result is 2017-
我正在使用 javascript 进行练习,我编写了这个脚本来计算 2 个变量的总和,然后在第三个方程中使用这个总和!关于如何完成这项工作的任何想法都将非常有用! First Number:
我有一个来自EAC的提示单和一个包含完整专辑的FLAC文件。 我正在尝试制作一些python脚本来播放文件,因为我需要能够设置在flac文件中开始的位置。 如何从CueSheet格式MM:SS:FF转
这个问题已经有答案了: Adding two numbers concatenates them instead of calculating the sum (24 个回答) 已关闭去年。 我有一个
4000 我需要上面字段 name="quantity" 和 id="price" 中的值,并使用 javascript 函数进行计算,并将其显示在字段 id= 中仅当我单击计算按钮时才显示“总
我是一名优秀的程序员,十分优秀!