I am following this protocol for differential gene expression
https://kimlaborg.files.wordpress.com/2017/10/nprot-2016-095.pdf
In order to identify transcripts and genes with a q value <0.05 the following command in used
subset(results_transcripts,results_transcripts$qval<0.05)>subset(results_genes,results_genes$qval<0.05)
The Ballgown output is supposed to be appear in the form of table (2 tables in total, 1 for for genes and 1 for transcripts) on the screen for each of these commands as mentioned in protocol (Step16), but for my data nothing appear.
I am new to this field . I'll be really grateful if someone can guide me about how to get tables from these commands.