Correcting a typo in my previous post. Code should have been: FDR <- p.adjust(lrt$table$PValue, method="BH") sum(FDR < 0.05) Gordon ______________________________________________________________________ The information in this email is confidential and intend...{{dropped:4}}