These are the data and code to go with "Rapid evolution of thermal tolerance and phenotypic plasticity in variable environments" Figure 01 takes the following data/scripts: 20210610_Thally02_Figure01_plot_and_stats.R with track_keeper.csv and corr_Response growth .csv . These files contain growth rates per transfers for all selection environments throughout the experiment and growth rates in correlated environments during reciprocal transplants, respectively. Figure 02 takes the following data/scripts: 20210610_Thally02_Figure02_plot_and_stats.R with 20161120_res_logis_t000.csv and 20161123_resloglint300.csv . These files contain the output of the shapes of the growth curves (i.e. information on lag time , growth at µmax, K etc) for all samples in all selection environments at t0 and t300, respectively The remaining figures - position not clear at time of submission - take the following data/script. For plasticity in FRRF data, the script 20181204_FRRF_plasticity.R takes the FRRF raw data contained in allfvfmdata_thally_t300_t000.csv . Extracted parameters are in files CvaluesThally02.csv, psiPSI_slope_intercept.csv, and rP_extracted_values.csv and can be analysed using the R script extracted parameter plots.R . R script FRRF visualisation only .R is for visualisation only, as the title suggests. For comparing plasticity/growth , the data are in 20170327_giantbigtable.csv , and can be visualised/analysed in plast vs growth.R In order to recreate the AMOVAS based on SNVs, use all_variants_fixed-only_using_5x_depth_threshold.cvs with amova thally02.R For additional information, please contact elisa.schaum@uni-hamburg.de