Applying Spatial Statistical Analysis to Ordinal Data for Soybean Iron Deficiency Chlorosis

Abstract

Accounting for field variation patterns plays a crucial role in interpreting phenotype data and, thus, in plant breeding. Several spatial models have been developed to account for field variation. Spatial analyses show that spatial models can successfully increase the quality of phenotype measurements and subsequent selection accuracy for continuous data types such as grain yield and plant height. The phenotypic data for stress traits are usually recorded in ordinal data scores but are traditionally treated as numerical values with normal distribution, such as iron deficiency chlorosis (IDC). The effectiveness of spatial adjustment for ordinal data has not been systematically compared. The research objective described here is to evaluate methods for spatial adjustment of ordinal data, using soybean IDC as an example. Comparisons of adjustment effectiveness for spatial autocorrelation were conducted among eight different models. The models were divided into three groups: Group I, moving average grid adjustment; group II, geospatial autoregressive regression (SAR) models; and Group III, tensor product penalized P-splines. Results from the model comparison show that the effectiveness of the models depends on the severity of field variation, the irregularity of the variation pattern, and the model used. The geospatial SAR models outperform the other models for ordinal IDC data. Prediction accuracy for the lines planted in the IDC high-pressure area is 11.9% higher than those planted in low-IDC-pressure regions. The relative efficiency of the mixed SAR model is 175%, relative to the baseline ordinary least squares model. Even though the geospatial SAR model is the best among all the compared models, the efficiency is not as good for ordinal data types as for numeric data.