摘要
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)
授權條款
我,本作品的著作權持有者,決定用以下授權條款發佈本作品:
|
此檔案在創用CC 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