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