Started work on new branch of analysis analyzing expression of both colonies together
Update script to identify thermally responsive transcripts using data from both
colonies together.

Modified RxN function to take user specified formula to allow this. Dropped q-value
implementation and will add this separately.
Updated RxNseq function to (a) only return coefficients significant from anova() call, and (b) return dataframe as function
More updates to RxNseq function:
  - return each coefficient only if significant in anova() call
  - name coefficient based on model formulat in returned data.frame
Cleaning up script for gene set enrichment analysis. Separated out thermally-responsive and colony-differentiated transcripts
topGO analysis for thermal responsive genes
Add code to identify shape of response (increase, decrease, convex, concave) and further analysis of GO
Minor updates