Only part of the p values were calculated for q values

[jjandshi]

Hi there,

Here is my code and error information: library(qvalue) Res_all <-read.csv("model2_com_MVPA_aC_glm_530.csv", check.names =FALSE, header = FALSE) colnames(Res_all) <- c("Taxa","beta", "se", "p")

table_median <- read.csv("/ForFDR_2710taxa_530.csv", check.names =FALSE, header = TRUE)

dat <- merge(table_median,Res_all, by.x="Taxa", by.y = "Taxa", all.x = FALSE, all.y = FALSE) phylum <- dat[dat$Rank == "phylum",] phylum[with(phylum,order(p)),] phylum$q.value <- qvalue(p =phylum$p)

Error in $<-.data.frame(*tmp*, q.value, value = list(call = qvalue(p = phylum$p), : replacement has 8 rows, data has 17

ForFDR_2710taxa_530.csv model2_com_MVPA_aC_glm_530.csv