我正在尝试一些模拟的p值分布。
使用hist()
绘制数据时,结果与预期一致。 p值均匀分布,峰值接近于零('信号')。
data = read_tsv("./data_vector.tsv")
hist(data$p_values, breaks = seq(0, 1, by=1/30))
但是,当使用ggplot / qplot执行相同操作时,左侧的峰值将丢失:
qplot(data$p_values, geom="histogram", bins=30)
我出错了什么?我本以为这两个命令是等价的。
数据&代码: 我的输入数据和Rmarkdown报告可从此gist
获得