我有一个包含 1893 个 p 值的数据框,我使用 R 包将其转换为 qvalues qvalue
:
pvals <- read.table("p_values.csv",header=T,sep=",",dec=".")
qvalue 手册指出,您可以通过以下方式直接绘制(使用 qplot 函数)qvalue 函数的输出:
qobj <- qvalue(pvals)
qplot(qobj)
但是,这会导致错误:
不知道如何为 qvalue 类型的对象自动选择比例。在 as.data.frame.default(x[[i]], optional = TRUE, stringsAsFactors = stringsAsFactors) 中默认为连续错误:无法将类“qvalue”强制转换为 data.frame
任何人都知道如何解决这个问题并qplot
绘制qvalue
函数的输出?