摘要
R code
d <-
expand.grid(P1 = seq(0.001, 1, length = 500),
P2 = seq(0.001, 1, length = 500))
d$`Fused P` <-
pchisq(with(d,-2 * (log(P1) + log(P2))), 2 * 2, lower.tail = F)
library(ggpubr)
library(extrafont)
loadfonts(device = "win")
cols <-
c(
'#67001f',
'#730d2c',
'#7f1938',
'#8a2545',
'#953052',
'#9f3b5f',
'#a9466c',
'#b35278',
'#bc5e85',
'#c46a91',
'#cc779d',
'#d483a9',
'#db90b4',
'#e19ebf',
'#e7acca',
'#ecb9d4',
'#f0c8de',
'#f3d6e8',
'#f6e5f1',
'#f7f4f9'
)
ggplot(d, aes(P1, P2)) +
geom_contour_filled(aes(z = `Fused P`), binwidth = 0.05, color = 8) +
coord_fixed() +
scale_fill_manual(values = cols) +
scale_x_continuous(
expand = expansion(c(0, 0)),
limits = c(0, 1),
breaks = seq(0, 1, 0.1)
) +
scale_y_continuous(
expand = expansion(c(0, 0)),
limits = c(0, 1),
breaks = seq(0, 1, 0.1)
) +
guides(fill = guide_legend(title = "Fused P")) +
theme_pubr(11, "Liberation Sans", legend = "right") +
theme(legend.text = element_text(size = 11)) +
x11(width = 7, height = 6)
许可协议
我,本作品著作权人,特此采用以下许可协议发表本作品:
|
本作品采用知识共享CC0 1.0 通用公有领域贡献许可协议授权。
|
采用本宣告发表本作品的人,已在法律允许的范围内,通过在全世界放弃其对本作品拥有的著作权法规定的所有权利(包括所有相关权利),将本作品贡献至公有领域。您可以复制、修改、传播和表演本作品,将其用于商业目的,无需要求授权。
http://creativecommons.org/publicdomain/zero/1.0/deed.enCC0Creative Commons Zero, Public Domain Dedicationfalsefalse
|
Fusing two p-values by using Fisher's method