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.