内容

在这个描述我们希望为相同的模拟生成情节simulations_vignette所示。弗兰克-威廉姆斯控制限制型心肌病,但这一次。特别是,我们应用IHW-Bonferroni方法和Bonferroni方法使用相同的设置在上面的模拟。

库(ggplot2)图书馆(“网格”)图书馆(dplyr)库(cowplot)库(IHWpaper)库(“wesanderson”)
颜色< - wes_palette (“Cavalcanti1”) [c (3、5)

1空模拟

1.1弗兰克-威廉姆斯控制

null_folder < -系统。文件(“simulation_benchmarks / result_files / ihw_bonf_null”,包=“IHWpaper”) null_sim < - bind_rows(拉普兰人(file.path (null_folder, list.files (null_folder))、功能(x) readRDS (x))) % > %变异(fdr_method = ifelse (fdr_method = =“IHW-Bonferroni E3”、“IHW-Bonferroni fdr_method)) last_vals_a < - group_by (null_sim fdr_method) % > %总结(last_vals = max(罗斯福))% > %变异(last_vals = last_vals + c(+ 0.003, -0.003),标签= fdr_method颜色=颜色)panel_a < - ggplot (null_sim, aes (x =α,y =罗斯福,坳= fdr_method)) + geom_line(大小= 1.2)+ geom_abline(线型=“冲”)+ xlab(表达式(粗体(粘贴(“名义”,α))))+ ylab (fw) + scale_x_continuous(限制= c(0.01, 0.1),休息时间= seq(0.01, 0.09,长度= 5))+ # ylim(0, 0.9) +主题(情节。利润=单位(c(3 7.5。2,。2),“线”))+ scale_color_manual(值=颜色)+主题(轴。title = element_text(脸=“大胆”))panel_a < - pretty_legend panel_a (panel_a, last_vals_a, 0.102)

2效果模拟(四)

2.1弗兰克-威廉姆斯控制

ttest_folder < -系统。文件(“simulation_benchmarks / result_files / ihw_bonf_du_ttest_informative”,包=“IHWpaper”) ttest_sim < - bind_rows(拉普兰人(file.path (ttest_folder, list.files (ttest_folder))、功能(x) readRDS (x))) % > %变异(fdr_method = ifelse (fdr_method = =“IHW-Bonferroni E3”、“IHW-Bonferroni fdr_method)) last_vals_b < - group_by (ttest_sim fdr_method) % > %总结(last_vals =罗斯福(which.max (eff_size)]) % > %变异(last_vals = last_vals标签= fdr_method颜色=颜色)panel_b < - ggplot (ttest_sim, aes (x = eff_size, y =罗斯福,坳= fdr_method)) + geom_hline (yintercept = 0.1,线型=“冲”)+ geom_line(大小= 1.2)+ xlab(“效果”)+ ylab(“额外的”)+主题(情节。利润=单位(c(3 7.5。2,。2),“线”))+ scale_color_manual(值=颜色)+主题(轴。title = element_text(脸=“大胆”))panel_b < - pretty_legend panel_b (panel_b, last_vals_b, 2.52)

3权力

last_vals_c < - group_by (ttest_sim fdr_method) % > %总结(last_vals =权力[which.max (eff_size)]) % > %变异(last_vals = last_vals标签= fdr_method颜色=颜色)panel_c < - ggplot (ttest_sim, aes (x = eff_size, y =, = fdr_method)上校)+ geom_line(大小= 1.2)+ xlab(“效果”)+ ylab(“权力”)+主题(情节。利润=单位(c(3 7.5。2,。2),“线”))+ scale_color_manual(值=颜色)+主题(轴。title = element_text(脸=“大胆”))panel_c < - pretty_legend panel_c (panel_c, last_vals_c, 2.52)

4规模投资模拟

5弗兰克-威廉姆斯控制

size_investing_folder < -系统。文件(“simulation_benchmarks / result_files / ihw_bonf_wasserman_normal”,包=“IHWpaper”) size_investing_sim < - bind_rows(拉普兰人(文件。路径(size_investing_folder list.files (size_investing_folder)),函数(x) readRDS (x))) % > %变异(fdr_method = ifelse (fdr_method = =“IHW-Bonferroni E3”、“IHW-Bonferroni”, fdr_method)) last_vals_d < - group_by (size_investing_sim fdr_method) % > %总结(last_vals = log10(罗斯福(which.max (xi_max)])) % > %变异(last_vals = last_vals + c(-0.8 + 2.5) * 0.08,标签= fdr_method颜色=颜色)panel_d < - ggplot (size_investing_sim, aes (x = xi_max, y = log10(罗斯福)= fdr_method)上校)+ geom_line(大小= 1.2)+ xlab(表达式(粗体(ξ(max)))) + ylab (“log10 (fw)”) + geom_hline (yintercept = 1,线型=“冲”)+主题(情节。利润=单位(c(3 7.5。2,。2),“线”))+ scale_color_manual(值=颜色)+主题(轴。title = element_text(脸=“大胆”))panel_d < - pretty_legend panel_d (panel_d, last_vals_d, 6.02)

6权力

last_vals_e < - group_by (size_investing_sim fdr_method) % > %总结(last_vals =权力[which.max (xi_max)]) % > %变异(last_vals = last_vals + c(-0.01 + 0.01),标签= fdr_method颜色=颜色)panel_e < - ggplot (size_investing_sim, aes (x = xi_max, y =, = fdr_method)上校)+ geom_line(大小= 1.2)+ ylab(“权力”)+ xlab(表达式(粗体(ξ(max)))) +主题(情节。利润=单位(c(3 7.5。2,。2),“线”))+ scale_color_manual(值=颜色)+主题(轴。title = element_text(脸=“大胆”))panel_e < - pretty_legend panel_e (panel_e, last_vals_e, 6.02)

7复合面板:

fwer_sim_fig < - plot_grid (panel_a ggdraw (), panel_b, panel_c, panel_d, panel_e, nrow = 3,标签= c (”)”、“”、“b)”、“c)”、“d)”、“fwer_sim_fig e)”))

ggsave(情节= fwer_sim_fig文件= " fwer_simulations。pdf”,宽度= 12,高度= 16)