FigShare Read Me File ARR JBH 05292020 -ARR update 06212020 -JBH update 07102020 Supplementary Figures: Figure S1. A Matrix Heatmap Heatmap of the A Matrix demonstrating relatedness between individuals. Figure S2. Correlation between A and D matrices Scatterplot of correlation between the A and D Matrices off-diagonal values Figure S3. Entry Sharing vs. Kinship between locations Comparison of entry sharing to mean kinship between Environments. Figure S4 (a-b). FA Biplot and dendrogram for Environmental clustering Biplot of first two factors of the FA model colored by environmental cluster. Dendrogram of environmental clustering using same color schematic. Figure S5. Scree Plot for Env FA clustering Scree plot for within-between variance ratio from clustering environments based on weather factors Figure S6. Stepwise Regression related to Yield 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. Figure S7. Marker Data Pipeline Workflow for processing genomic marker data from inbred parents to generate hybrid genotypes. Figure S8. Stand Boxplots by location Boxplot of stand for each environment. Outliers detected using IQR and mean-percentile criteria colored pink. Supplementary Tables: Table S1. Env -> Site year pairings Table matching each Environment name used for analysis to its site and year. Table S2. Stage 1 Covariates 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. Table S3. Stage 1 Within-Environment Selected Model Info 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. Table S4. Genetic Correlations between environments Pairwise genetic correlations between environments from Yield models fit in Echidna. 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. 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. Table S7. Loading terms for the model with best BIC from stepwise regression analysis. Supplementary Files: File S1. Supplemental Methods Description of the marker, phenotypic, and weather data cleaning and exploration done. approximately 10 pages of comprehensive data cleaning information. 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. File S3. Hybrid best linear unbiased estimators BLUEs for all traits within each environment based on selected model File S4. Echidna Code Sample Echidna code for models utilizing IDVG Variance Structures File S5. Echidna Code Sample Echidna code for models utilizing the A Matrix. File S6. Echidna Code Sample Echidna code for models utilizing the D Matrix File S7. Echidna Code Sample Echidna code for models utilizing both the A and D matrices in tandem. File S8. Environment BLUEs (averaged across all hybrids) for all traits File S9. Hybrid BLUEs (averaged across all environments) for all traits File S10. Hybrid Cluster assignments and A and D diagonal values Cluster assignments of hybrids from hierarchical clustering of marker data, along with their diagonal values from the A and D matrices. File S11. Variance Components All Traits all models 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. File S12. Three Year Loadings Factor loadings of each environment variable from the FA(10) model applied to 30-day period weather data. File S13. List of Hybrid names in order of rows and columns of the realized genomic relationship matrices File S14. Realized additive relationship matrix, TASSEL output format File S15. Realized dominance relationship matrix, TASSEL output format File S16. G2F Environmental Data Processing Functions.R This file contains functions 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. File S8. Spec-IEM Correlations This file contains correlations between data obtained from the in-field Spectrum Watchdog weather stations and data scraped from the Iowa Environmental Mesonet database. File S18. Weather variables averaged over 5-day periods for each environment. These values are in original units before scaling.