This work is licensed under a Creative Commons Attribution 4.0 International License.
Completely re-working analysis, splitting thermally-responsive transcripts to examine differences among transcripts with colony x temperature interaction.
Nice tip on how to identify inflection point.
Found that my “Intermediate” expressed transcripts included those with the correct curvature, but that had maximum expression below (<10) or above (>31) the ‘intermediate’ interval. Felt rather foolish - this is exactly the issue with identifying stabilizing selection highlighted by Mitchell-Olds & Shaw (1987). At least I’m not the first (or last) to make this mistake.
Modified function to only call a transcript “Intermediate” if curvature is correct and the maximum is between 10 and 31C.
Interestly, also found that there’s an R function named after Mitchell-Olds & Shaw (1987) to check if maximum or minimum of quadratic function is within studied interval.
This work is licensed under a Creative Commons Attribution 4.0 International License.