I just wanted to get confirmation that the following code would get me the result I wanted. So I am interested in getting the list of differential expressed genes between two groups between seizure YES and seizure NO. The expression set is the microarray data and seizure information is in the affy.clinic.sz data frame "seizure" column.
contrasts<-makeContrasts(YES-NO, levels = design.sz)
topTable(fit.sz2, adjust.method = "BH", sort.by = "logFC")
Will this list give me the list of DE genes comparing seizure YES group to NO group, where negative logFC value indicates that the gene's expression is higher in seizure NO group and positive logFC indicates that the gene's expression is higher in sz YES group?
Thank you for your help in advance!!