Supplemental Material for Rogers et al., 2020

Figure S1 Heatmap of the A Matrix demonstrating relatedness between individuals.<br><br>Figure S2. Scatterplot of the A and D Matrices off-diagonal values<br><br>Figure S3. Comparison of entry sharing to mean kinship between Environments.<br><br>Figure S4 (a-b). Biplot of first two factors of the FA model colored by environmental cluster. Dendrogram of environmental clustering using same color schematic.<br> <br>Figure S5. Scree plot for within-between variance ratio from clustering environments based on weather factors<br><br>Figure S6. First ten weather variables selected in stepwise regression model relating D×FA(1) environment loadings to weather variables summarized over 5-day windows. Each variable and mean environment yield are standardized.<br> <br>Figure S7. Workflow for processing genomic marker data from inbred parents to generate hybrid genotypes.<br><br>Figure S8. Boxplot of stand for each environment. Outliers detected using IQR and mean-percentile criteria colored pink.<br> <br><br><br><br>Supplementary Tables:<br>Table S1. Table matching each Environment name used for analysis to its site and year.<br> <br>Table S2. Covariates of experimental design and which year each could be fitted for/Blocking factors and residual variance models tested for each year in stage 1 analysis.<br><br>Table S3. Model chosen for each trait-environment combination including the fixed effects, random effects, and residual variance structure of the model selected using BIC. In addition, error variance from the model and both mean and plot basis heritabilities are given.<br> <br>Table S4. Pairwise genetic correlations between environments from Yield models fit in Echidna. <br><br>Table S5. ANOVA model for Yield to using weather-based environment clusters demonstrate relationship between environmental data and environment main effect and G×E variances.<br><br>Table S6. Model fit information for stepwise regression models using 5, 10, 15, and 30-day windowed environmental covariates as a predictor of parameter estimates from respective AxFA(1) and DxFA(1) models. Best fitting model for each parameter by BIC are in bold.<br><br>Table S7. Loadings terms for the model with best BIC from stepwise regression analysis.<br><br><br>File S1. Supplemental Methods. Description of the marker, phenotypic, and weather data cleaning and exploration done. approximately 10 pages of comprehensive data cleaning information.<br> <br>File S2. Stage 1 R Code. This script processes through Stage 1 modeling including multiple model types to account for experimental design in a given environment. The best model is output using BIC as a metric.<br><br>File S3. Hybrid best linear unbiased estimators BLUEs for all traits within each environment based on selected model<br><br>File S4. Sample Echidna code for models utilizing IDVG Variance Structures<br><br>File S5. Sample Echidna code for models utilizing the A Matrix.<br><br>File S6. Sample Echidna code for models utilizing the D Matrix<br><br>File S7. Sample Echidna code for models utilizing both the A and D matrices in tandem.<br> <br>File S8. Environment BLUEs (averaged across all hybrids) for all traits<br><br>File S9. Hybrid BLUEs (averaged across all environments) for all traits<br> <br>File S10. Hybrid Cluster assignments and A and D diagonal values, along with their diagonal values from the A and D matrices.<br> <br>File S11. Tables of Variance components for each trait and the model fit. For Yield, an expanded set of models were fit while all other traits fit only the IDVG models and most simple model using both the A and D realized relationship matrices.<br> <br>File S12. Factor loadings of each environment variable from the FA(10) model applied to 30-day period weather data. <br> <br>File S13. List of Hybrid names in order of rows and columns of the realized genomic relationship matrices<br><br>File S14. Realized additive relationship matrix, TASSEL output format<br><br>File S15. Realized dominance relationship matrix, TASSEL output format<br><br>File S16. R code for processing G2F weather data from cleaned files. The first function, plant, computes days since planting for each location. The second function, DailyMeans computes daily values for covariates from data taken throughout a given day at a location using Julian date format. The third function uses a while loop to allow the user to specify the sliding window size for their environmental data output. Between steps assume that the user has filled in any missing values using an appropriate source.<br><br>File S17. This files contains correlations between data obtained from the in-field Spectrum Watchdog weather stations and data scraped from the Iowa Environmental Mesonet database. <br><br><br>File S18. Weather variables averaged over 5-day periods for each environment. These values are in original units before scaling.<br><br><br>