P值的直方图

Question

我想绘制p值pv.scf。为什么我的直方图不能反映我的4个p值？

> dput(pv.scf)
c(`211124_s_at` = 0.0003175435117273, `207029_at` = 0.00882067488788698, 
`216974_at` = 0.586375149259342, `226534_at` = 0.651991817084333
)

# Plot the sorted raw P-values
png("SCF_Raw_PvaluePlot.png")
par(mfrow=c(1,2))
hist(pv.scf, cex=0.9, col=c("#c1a87d", "#7ef9ff", "#97b9ce", "#1f456e"), las=2, cex.axis=0.6, xlab="Probes", ylab="P-values", main="Raw P-values for KITLG gene\n")
names(c("211124_s_at", "207029_at", "216974_at", "226534_at"))
abline(v=.01,col=2,lwd=2)
hist(-log10(pv.scf), col=c("#c1a87d", "#7ef9ff", "#97b9ce", "#1f456e"), xlab="Probes", ylab="log10(p-values)", main="-log10(P-values) for KITLG gene between\nall groups", cex.main=0.9)
names(c("211124_s_at", "207029_at", "216974_at", "226534_at"))
abline(v= -log10(.01),col=2,lwd=2)
dev.off()

​

​

Answer 1

听起来你真的想要一个柱状图而不是柱状图。并不是所有带条形图的绘图都是一样的。尝尝这个

par(mfrow=c(1,2))
barplot(
  pv.scf, 
  cex=0.9, 
  col=c("#c1a87d", "#7ef9ff", "#97b9ce", "#1f456e"), 
  las=2, cex.axis=0.6, 
  xlab="", ylab="P-values", 
  main="Raw P-values for KITLG gene\n",
  names =c("211124_s_at", "207029_at", "216974_at", "226534_at"))
abline(v=.01,col=2,lwd=2)
barplot(
  -log10(pv.scf), 
  col=c("#c1a87d", "#7ef9ff", "#97b9ce", "#1f456e"), 
  xlab="", ylab="log10(p-values)", 
  main="-log10(P-values) for KITLG gene between\nall groups", 
  cex.main=0.9,
  las=2,
  names =c("211124_s_at", "207029_at", "216974_at", "226534_at"))
abline(v= -log10(.01),col=2,lwd=2)

这导致了

​