The effect of spatial lag on modeling geomatic covariates using analysis of variance

The effect of spatial lag on modeling geomatic covariates using analysis of variance | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article The effect of spatial lag on modeling geomatic covariates using analysis of variance Aquiles Enrique Darghan Contreras, Darlley Stiven Taborda Lozada, and 3 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-3243407/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 22 Jul, 2024 Read the published version in Applied Geomatics → Version 1 posted 4 You are reading this latest preprint version Abstract In recent years, statistical methods involving spatial considerations have been developed, for example, those incorporating data with some type of georeferencing. The descriptive part of geographic information systems currently provides many visualization and analytic tools; however, the latter is still quite limited. In this sense, research of a spatial nature is seen as combining non-spatial statistical methods for inferential treatment that can certainly invalidate the excellent capture work with advanced tools such as those observed every day in the geomatic context. This prompted the current document, drawing attention to how geomatic information analyzed with statistical methods that imply independence in modeled observations can be invalid. The Moran index is compared with a proposal for a spatial lag coefficient in the context of experimental design so that users of variance analysis do not apply this well-known procedure in a ritualistic way, perhaps revising some assumptions and perhaps ignoring more important ones. The distortion of the p value generated from the analysis of variance is clear in the presence of spatial dependence. In this case it is associated with the lag or spatial overlap. The methodology is simple to adopt in other experimental designs with the simple consideration of the design matrix and its reparameterization and the choice of the appropriate weight matrix. This will allow users to reconsider the traditional method of analysis and incorporate some methodology to support spatial dependency structures. spatial lag spatial overlapping spectral index sensing Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Introduction Analysis of variance (ANOVA) is a very useful tool in the biological and industrial context, and it has undoubtedly gained ground in modern fields such as geomatics, where data captured with remote or proximal sensors predominate. ANOVA has been one of the standard tools used to analyze experimental data for decades; however, currently there are so many application modalities that it is not enough to appear in the methodology that applied an ANOVA (AOV, ANAVAR, ANDEVA) to have complete clarity of what is behind the method and its proper use, since aspects such as balancing, blocking, the presence of covariates, temporal measurements, geo-referenced information, discrete or categorical response(s) of univariate or multivariate nature generate their own analytic strategies for estimating parameters of the associated model (Christensen, 2019 ). Currently, remote sensing data capture technologies, and drone images or special cameras allow the creation of spectral indexes (Adak et al., 2021 ) that are associated with data that usually have very different structure to the experimental design, obligating users to do some kind of aggregation to experimental plots or other units in order to be able to treat this information as a source of additional variation in the ANOVA or simply as the same response. These data that involve some form of georeferencing, may also incorporate some spatial dependency structure (and many other temporal and space-times). Regarding different aspects such as analysis of vegetation cover, climatic and meteorological variables, among others, it is important to highlight that these vary due to many factors, including the time scale, therefore, it is necessary to develop evaluations with updated information from two or more periods, which allow comparisons to identify temporary variations of the object of study, before the attribution of the cause of said change. This process is known as multitemporal evaluation, which consists of comparing the object in the same geographical area in different periods, in order to know its variation, for this there are several evaluation methods (Gao et al., 2012 ). If these aspects are overlooked in the inferential process, it could invalidate its application, especially if they are analyzed with the methods usually studied in basic courses of experimental design (Kenny, 1995 ; Kenny & Judd, 1986 ; Zimmerman et al., 2007 ; Zimmerman & Zumbo, 2010 ). It is surprising that many users of the technique incompletely report aspects for consideration prior to the ANOVA (experimental design) in the ANOVA’s own results; and they report insufficient details of analysis after ANOVA (interactions or multiple comparisons) to be certain of the scope and interpretation of the results (Acutis et al., 2012 ). This allows its statistical software package to set the tone for substantive data analysis decision-making usurping in some way the preponderant role of the users, perhaps assuming that by the fame of the software or its recognition it can be taken for granted that everything developed there is directly applicable or interpretable. Several documents highlight the misuse of ANOVA in two-way designs due to heterogeneous variance (Rowell & Walters, 2006 ) highlighting defective data (Pearce, 2006 ). Defective data includes the non-specification of the type of ANOVA (Abubakar et al., 2022 ), confusion in the choice of the type of effects (fixed or random) (Bennington & Thayne, 1994 ), assumption too quickly of the weak argument that not rejecting the null hypothesis implies a strong conclusion (Rong, 2000 ), inappropriate use of the terms associated with the error in the construction of the statistical F or the type of sums of squares (Li & Lomax, 2011 ). Such defective data even includes the belief, supported in recognized literature that only a couple of assumptions are necessary to make use of the technique (Acutis et al., 2012 ). This ignores the supposed key to independence (Lindman, 1992 ) in which it is closely linked to the reality of the data captured with sophisticated sensors or equipment related to agriculture, such as performance monitors, spectrometric or other techniques for capturing data from radar or satellite images, or images from spectral or hyperspectral cameras, etc. All this information is usually converted into indices that usually relate to the physical-chemical and microbiological properties of soils or water as well as the vegetation, specifically in the development of the crop. In turn, this information is quite often incorporated into the structure of the experimental design in order to associate the data with each observation or experimental unit and thus, using ANOVA, evaluate the effect of the treatments (Atik & Akdemir, 2022 ; Barbosa et al., 2020 ; Ding et al., 2020 ; Firozjaei et al., 2020 ; Stoy et al., 2022 ; Volcani et al., 2005 ; Ya’acob et al., 2014 ). In the documents cited above, how SIMPLE (Simple, Informative, Meaningful, Powerful, Logical, Effective,) (Mcintosh, 2015 ) were these analyses in general? In the case of data which in many cases have been shown to have spatial or temporal dependency, shouldn´t an adequate description be made of the essential aspects of the ANOVA to involve the dependency structure present in the data? Knowing the weaknesses of ANOVA in certain contexts and of certain types of data can very surely enhance its benefits. A form of spatial dependence can arise with remote sensor data or similar technologies, such as spatial lag or overlap (Shukl & Subrahmanyam, 1999 ), in order to find data associated with experiments where some design model was adjusted for the comparison of treatments, using as response(s) or some performance indicator evaluated in field or some estimate generated by some type of sensor (Wang et al., 2016 ). Some of these indices are added to each polygon that make up the experimental plot to finally incorporate them into the model, sometimes as a response and sometimes as covariates (if treatments do not affect this covariate, or a multiple-slope model has been specified). The literature points out how these indices show some kind of spatial dependence (dos Santos et al., 2021 ; H. Zhang et al., 2011 ); however, the rituality of traditional analysis mainly by knowledge of the methods that usually do not treat spatial dependence allows many excellent data captures, invalidating the work of field or laboratory tests because of the analytic technique (Gotway & Cressie, 1990 ). With all the reasons and documentation presented, it is important to understand how ANOVA can be invalidated, at least in the context of spatial lag; but the same results extend to other situations that may have the same effect, such as edge effects, competition, poor specification of a model, etc. (Christensen & Bedrick,1997). Using conditional simulations, a Gaussian variogram was generated with fully defined parameters, and 1000 simulations were obtained to generate spatial dependence. The final goal was to show the invalidity of the ANOVA when spatial lag was presented using the Moran index and the spatial lag coefficient, although the other assumptions usually found as the most important in the context of the ANOVA are fulfilled. Methodology In the linear models treated in experiment design, Gauss-Markov dominates, written as \(\mathbf{y} = \mathbf{X}\varvec{\beta } + \varvec{\epsilon }\) ; \(E\left(\varvec{\epsilon }\right)=0\) and \(Cov\left(\varvec{\epsilon }\right)={}^{2}\varvec{V}\) , where \(\varvec{V}\) is a known matrix representing some correlation structure between residuals. The different models depend on the assumptions made about \(\varvec{V}\) , and it is precisely the simplest one dealt with in the design where it is assumed that V = I (an identity matrix), that is, it assumes independence and identical distribution, something that cannot be assumed so quickly when data is associated with information from remote sensing technologies. To deal with this limitation of the model, an extended model by spatial lag is used written as: $$\mathbf{y} = \mathbf{X}\varvec{\beta } + {{\alpha }}_{\text{s}}\mathbf{W} \mathbf{X}\varvec{\beta } + \varvec{\epsilon }$$ 1 where \(\mathbf{y}\) is a random vector of length n denoting the response of n experimental units where the response associated with the spectral index, or any other index could usually be involved. \(\mathbf{X}\) is the known design matrix of dimension n × p and incomplete range (as a design model), \(\varvec{\beta }\) is a vector of unknown parameters of length p , which, in the case that develops consists of the effect of blocks and treatments, \({{\alpha }}_{\text{s}}\) is the coefficient associated with the effect of spatial lag (overlap) and \(\mathbf{W}\) = ( w ij ) is the matrix of spatial weights of dimension n × n , where w ij represents the effect of unit j on unit i . A conditional simulation was generated from a Gaussian variogram (Gräler et al, 2019) with parameters as described in the following R script: The 1000 simulations were assigned to a 10m×6m grid for 60 polygon-associated cells representing the experimental units. With this matrix and their respective latitude and longitude coordinates artificially created for the grid (per row[y] and column[x]), three treatments were assigned in two blocks to associate the ANOVA to the simple factorial design in generalized and random complete blocks without interaction. The random generation of this data set ensured that in the 1000 runs of the ANOVA there was no effect of the treatments nor that the blocking was efficient in not distorting the distribution of the test statistic (F) as regards the no centrality parameter. The image generated for the grid with treatment-block-repeat assignment (T.B.R) is shown in Fig. 1 . Each cell associated with each polygon was labeled with the respective treatment (T:{1,2,3}), blocking ratio (B:{1,2}), and repetition (R:{1,2,...,10}). To generate the spatial lag, the cells were spatially varied (exchanging) by two in two in the combinations of T.B.R within the grid to modify the distance between these units and change the pattern of the neighborhood and, therefore, of the spatial dependence. In Fig. 2 , each image with different ordering (ord0[data from the original variogram with high spatial dependence] up to ord60 [data with all grid elements exchanged and thus with low spatial dependence]) represents the same T.B.R response. For all systems, the Moran index was calculated (J. Zhang et al., 2014 ) and applied to the residual design model once the ANOVA ran in all simulations, thus creating different scenarios of spatial dependence. This index was calculated with the ape library of R (Paradis & Schliep, 2019 ). In each simulation, the ANOVA F statistic was extracted, normality and homoscedasticity of residues were tested for treatments. The same O T lag statistic was calculated based on the Chi-square distribution (with a degree of freedom) that appears in Darghan et al. ( 2012 ) and that is annexed in this document in a script of R for use, as well as the coefficient associated with the spatial lag ( \({{\alpha }}_{\text{s}}\) ) which was obtained from: $${ \widehat{{\alpha }}}_{\text{s}}=\frac{\mathbf{y}\mathbf{{\prime }}{\mathbf{M}}_{0}\varvec{W}(\mathbf{I}-{\mathbf{M}}_{0})\mathbf{y}}{\varvec{y}\varvec{{\prime }}{\mathbf{M}}_{0}\varvec{W}\varvec{W}{\mathbf{M}}_{0}\varvec{y}} \left(2\right)$$ where \({\mathbf{M}}_{0}\) is associated with the design matrix having already repaired the model in the Christensen ( 2011 ) style, with \({\mathbf{M}}_{0}={\varvec{X}}_{0}{\left({\varvec{X}}_{0}\varvec{{\prime }}{\varvec{X}}_{0}\right)}^{-1}{\varvec{X}}_{0}^{\varvec{{\prime }}}\) where \({\varvec{X}}_{0}=\varvec{X}\varvec{Z}\) being \(\varvec{Z}\) generated with the " mz " function of R available in (GitHub1: supplementary material). This matrix depends on the selected design, and without loss of generality, it will be enough to modify the design matrix (in balanced cases) and thus obtain this special matrix for the new design. The equation in (2) also depends on the weight matrix and in this case the inverses of the Euclidean distances were used with type C standardization according to R spdep library (Bivand, 2022 ). For all spatial variations all these measurements were calculated in addition to the ANOVA (AOV function of R of the library stats (R Core Team, 2022 )), where the residuals of the model were extracted to test normality (function Shapiro.test of the library stats of R) and homoscedasticity test in the treatments (Bartlett.test function of R stats ). All the information formed a dimension matrix (31000×7) that appears available in github2 (supplementary material- 31 ordering scenarios for 1000 simulations and seven responses). With all these data diagrams that were made that related the p value of the Moran index and p value of the ANOVA statistic F with the statistic O T , the spatial lag coefficient in each treatment was arranged to demonstrate how this order was related to the Moran index and the lag coefficient, thus highlighting the little or lack of usefulness of the ANOVA in the presence of a dependence or spatial lag. Simulation Results Figure 3 is one of the most important in this document because it shows the effect of the spatial lag in the distributional behavior of the F statistic. The monotonous continuous line represents the theoretical distribution of the F statistic (central) with 2 and 56 degrees of freedom (for the established design). The dotted vertical line represents the quantile value for a 95% confidence level and the x-axis represents the values of the F obtained in the 1000 simulations. The histogram describes the observed distributional behavior of the F statistic in the presence of spatial lag. It is evident how the theoretical curve F cannot be used to obtain the p value of the ANOVA even though practically all the generated F values do not reject the null hypothesis of the null effect of treatments. Even in the presence of spatial lag, both scenarios lead to no rejection, but the p value, usually generated from the ANOVA is not the measure that represents what is occurring. This would also apply if the rejection of the null hypothesis was even due to a difference in treatment means. Figure 4 shows the probabilities associated with theoretical F and simulated F in the ord0, ord30 and ord60 systems to indicate the spatial dependence according to the Moran index. Already in the last order the spatial lag has been lost and therefore the p value of Moran is greater than 5% in almost all simulations. The dotted line represents the points where both probabilities coincide and happens when the statistic F in both cases tends to zero. It is also evident that it is more likely not to reject the null hypothesis of equality of means in spatial lag scenarios, making it easier to obtain false negatives in the presence of this effect. Figure 5 is like Fig. 4 except that the coefficient associated with the spatial lag for the same orders (ord0, ord30 and ord60) is highlighted here. The estimated coefficient of lag is near zero in spatial variations of the greatest occurrence of re-orders. Figure 6 shows the O T statistic for spatial lag for the null hypothesis Hoα s = 0. This hypothesis is tested with one degree of freedom. Also shown is a vertical line (for the 95% confidence level) in different systems to delimit in the x-axis the quantum of the theoretical F. The crossing of these lines yields four regions: A) [spatial lag and no treatment effect], B) [spatial lag and no treatment effect], C) [spatial lag and no treatment effect] and D) [spatial lag and no treatment effect]. From orders higher than or equal to 36 (ord36) practically all the data fall in region B of absence of spatial lag and null effect of the treatments. In the ord0, where the data with Moran index indicate spatial dependence, for a fraction of these data there is a spatial lag (B) and in the other there is none (A). To better understand these results, Table 1 describes the proportion of data in each region. Table 1 Distribution of the proportion of data by region in Fig. 6 . Region Order A B C D Ord0 49.7 50.2 0.0 0.1 Ord12 37.9 62 0.0 0.1 Ord24 30.1 69.8 0.1 0.0 Ord36 2.7 97.2 0.0 0.1 Ord48 2.8 97.1 0.0 0.1 Ord60 1.3 98.6 0.0 0.1 From Table 1 it is evident that from the Ord36 most of the data falls into B and at most 10% between regions C and D. This clearly shows as the spatial lag is lost; the spatial dependence is lost even though the ANOVA shows no evidence in the data to reject the null hypothesis of null effect. Finally, Fig. 7 separates by ordering the behavior of the estimated spatial lag coefficient \({\widehat{\alpha }}_{s}\) where an average at zero and greater variability is clearly perceived when there is greater spatial dependence and thus spatial lag. The clearest points indicate the cases with no spatial dependence when the estimated values of \({\widehat{\alpha }}_{s}\) oscillate approximately between − 0.01 and 0.01, maintaining a similar variability since the ord36, when the Moran index already suggests spatial independence. Although ANOVA yields in practice results associated with not rejecting the hypothesis of the null effect of treatments, the p value of the test is not appropriate as it does not reflect the probability of obtaining test results for the F statistic, at least as extreme as the result observed, on the assumption that the null hypothesis is correct. Spatial dependence has a direct impact on ANOVA validity. This violation of the assumption of independence by spatial lag (although it can be violated in many ways (Christensen, 1997 ) may have several important implications in the ANOVA, highlighting here the underestimation of the p-value of Exhibit F that affects the validity of the associated hypothesis tests and thus the interpretation of results, since the spatial dependence introduces a structure in the data that is not taken into account in the conventional ANOVA that is usually taught in basic courses of experimental design. Although there are different methods for detecting spatial dependency, surely the Moran index or some of its modifications are the most used and are already programmed in several statistical packages. In the case of the lag coefficient, it could also be a convenient measure adapted to each experimental design using the design matrix. As this represents the ratio of two independent quadratic forms, it is possible not only to obtain their moments but even their distributional behavior, making this a measure of alternative interest to judge spatial lag in the context of experimental design, provided that a weight matrix is defined ensuring symmetry in the matrices associated with the quadratic forms. Conclusions Spatial dependence represents a significant violation of the assumption of independence of the residues in the ANOVA. This assumption assumes that observations are independent of each other, that is, that the value of an observation is not related to the value of another nearby observation in space. When spatial dependence exists, close observations tend to be more like each other than more distant observations, meaning that the residues of neighboring observations are correlated and do not meet the assumption of independence. Since the observations taken by remote or proximal sensors of physic-chemical or microbiological properties of soils and vegetation can be grouped in space by evolutionary processes, degradation, crop or soil management or other agents can generate similar predictor variables. This similarity generates a structure in the correlation of observations that manifests as spatial dependence. If this dependence is not considered and in this case is caused by spatial lag in treatments, there is a risk of obtaining incorrect estimates of the p value associated with the ANOVA F test statistic; therefore, this leads to inappropriate inferences and, therefore, to interpretations of little or no value. Spatial autocorrelation can lead to overestimation or underestimation of the significance of treatment effects that can lead to erroneous acceptance or incorrect rejection of associated hypotheses. Addressing spatial dependence in the analysis of variance by Moran's index or by the lag coefficient for a specified experimental design can help to identify the problem and to take the necessary measures using the appropriate procedures or models for an estimate of some response or an appropriate comparison of treatments. Spatial dependence represents a critical violation of the independence assumption in the ANOVA and should not be omitted in analyses where studies have proven some form of dependence, as is the case in many of the phenomena that are modeled in the context of geomatics. Declarations Conflict of Interest The authors declare that there is no conflict of interest. Data Availability The data as well as the R scripts to generate the simulations can be found at: GitHub 1 : https://github.com/CarlosRivera1212/rezago_aov GitHub 2 : https://github.com/darklley/Rezago References Abubakar, B., Uthman, Y. A., Jatau, A. I., Danbatta, A., Nuhu, H. N., & Mustapha, M. (2022). Misuse of analysis of variance in African biomedical journals: a call for more vigilance. Bulletin of the National Research Centre 2022 46:1 , 46 (1), 1–8. https://doi.org/10.1186/S42269-022-00924-8 Acutis, M., Scaglia, B., & Confalonieri, R. (2012). Perfunctory analysis of variance in agronomy, and its consequences in experimental results interpretation. European Journal of Agronomy , 43 , 129–135. https://doi.org/10.1016/J.EJA.2012.06.006 Adak, S., Bandyopadhyay, K., Sahoo, R., Mridha, N., Shrivastava, M., & Purakayastha, T. (2021). Prediction of Wheat Yield Using Spectral Reflectance Indices Under Different Tillage, Residue and Nitrogen Management Practices. Current Science , 121 , 402. https://doi.org/10.18520/cs/v121/i3/402-413 Atik, M., & Akdemir, B. (2022). Spatial Variability of NDVI and SPAD for Variable- and Fixed-Rate Fertilization in Olive Orchards. Erwerbs-Obstbau , 64 (4), 647–655. https://doi.org/10.1007/S10341-022-00748-W/METRICS Barbosa, J. de A., de Faria, R. T., Coelho, A. P., Dalri, A. B., & Palaretti, L. F. (2020). Nitrogen fertilization management in white oat using spectral indices. Pesquisa Agropecuária Tropical , 50 , 1–9. https://doi.org/10.1590/1983-40632020V5064924 Bennington, C. C., & Thayne, W. V. (1994). Use and misuse of mixed model analysis of variance in ecological studies. Ecology , 75 (3), 717–722. https://doi.org/10.2307/1941729 Bivand, R. (2022). R Packages for Analyzing Spatial Data: A Comparative Case Study with Areal Data. Geographical Analysis , 54 (3), 488–518. https://doi.org/10.1111/GEAN.12319 Christensen, R.,Bedrick,E.J. (1997). Testing the independence assumption in linear models. Journal of the American Statistical Association, 92, 1006–1016. Christensen, R. (2011). Plane Answers to Complex Questions . https://doi.org/10.1007/978-1-4419-9816-3 Christensen, R. (2019). Advanced Linear Modeling . https://doi.org/10.1007/978-3-030-29164-8 Darghan, A., Surendra, P., & Goitía, A. (2012). An application of the overlapping test in a two-way classification reparameterized model. Revista de La Facultad de Agronomia , 29 , 124–137. Ding, Y., Jiang, Y., Yu, H., Yang, C., Wu, X., Sun, G., Fu, X., & Dou, X. (2020). Measurement Method for Height-Independent Vegetation Indices Based on an Active Light Source. Sensors 2020, Vol. 20, Page 1830 , 20 (7), 1830. https://doi.org/10.3390/S20071830 dos Santos, R. A., Filgueiras, R., Mantovani, E. C., Fernandes-Filho, E. I., Almeida, T. S., Venancio, L. P., & da Silva, A. C. B. (2021). Surface reflectance calculation and predictive models of biophysical parameters of maize crop from RG-NIR sensor on board a UAV. Precision Agriculture , 22 (5), 1535–1558. https://doi.org/10.1007/S11119-021-09795-X/METRICS Firozjaei, M. K., Fathololoumi, S., Mijani, N., Kiavarz, M., Qureshi, S., Homaee, M., & Alavipanah, S. K. (2020). Evaluating the Spectral Indices Efficiency to Quantify Daytime Surface Anthropogenic Heat Island Intensity: An Intercontinental Methodology. Remote Sensing 2020, Vol. 12, Page 2854 , 12 (17), 2854. https://doi.org/10.3390/RS12172854 Gao, J., Li, X., and Brierley, G. (2012). Topographic influence on wetland distribution and change in Maduo County, Qinghai-Tibet Plateau, China. Journal of Mountain Science, 9(3):362{371. Gotway, C. A., & Cressie, N. A. C. (1990). A Spatial Analysis of Variance Applied to Soil-Water Infiltration. Water Resources Research , 26 (11). https://doi.org/10.1029/WR026i011p02695 Gräler,B.,Pebesma,E.,Heuvelink,G.(2016). Spatio-Temporal interpolation using gstat. The R Journal, 8(1), pp. 204-218. Kenny, D. A. (1995). The effect of nonindependence on significance testing in dyadic research. Personal Relationships , 2 (1), 67–75. https://doi.org/10.1111/J.1475-6811.1995.TB00078.X Kenny, D. A., & Judd, C. M. (1986). Consequences of Violating the Independence Assumption in Analysis of Variance. Psychological Bulletin , 99 (3), 422–431. https://doi.org/10.1037/0033-2909.99.3.422 Li, J., & Lomax, R. G. (2011). Analysis of Variance: What Is Your Statistical Software Actually Doing? Http://Dx.Doi.Org/10.1080/00220973.2010.481684 , 79 (3), 279–294. https://doi.org/10.1080/00220973.2010.481684 Lindman, H. R. (1992). Analysis of Variance in Experimental Design . https://doi.org/10.1007/978-1-4613-9722-9 Mcintosh, M. S. (2015). Can Analysis of Variance Be More Significant? Agronomy Journal , 107 (2), 706–717. https://doi.org/10.2134/AGRONJ14.0177 Paradis, E., & Schliep, K. (2019). ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics , 35 (3), 526–528. https://doi.org/10.1093/BIOINFORMATICS/BTY633 Pearce, S. C. (2006). Defective data in the analysis of variance. Http://Dx.Doi.Org/10.1080/02664768600000022 , 13 (2), 139–147. https://doi.org/10.1080/02664768600000022 R Core Team. (2022). R: A Language and Environment for Statistical Computing . https://www.R-project.org/ Rong, Y. (2000). Statistical Methods and Pitfalls in Environmental Data Analysis. Environmental Forensics , 1 (4), 213–220. https://doi.org/10.1006/ENFO.2000.0022 Rowell, J. G., & Walters, D. E. (2006). A critical appraisal of the analysis of variance for two-way classifications. Http://Dx.Doi.Org/10.1080/02664768500000021 , 12 (2), 157–169. https://doi.org/10.1080/02664768500000021 Shukl, G. K., & Subrahmanyam, G. S. V. (1999). A note on an exact test and confidence interval for competition and overlap effects. Biometrics , 55 (1), 273–276. https://doi.org/10.1111/J.0006-341X.1999.00273.X Stoy, P. C., Khan, A. M., Wipf, A., Silverman, N., & Powell, S. L. (2022). The spatial variability of NDVI within a wheat field: Information content and implications for yield and grain protein monitoring. PLOS ONE , 17 (3), e0265243. https://doi.org/10.1371/JOURNAL.PONE.0265243 Volcani, A., Karnieli, A., & Svoray, T. (2005). The use of remote sensing and GIS for spatio-temporal analysis of the physiological state of a semi-arid forest with respect to drought years. Forest Ecology and Management , 215 (1–3), 239–250. https://doi.org/10.1016/J.FORECO.2005.05.063 Wang, C., Nie, S., Xi, X., Luo, S., Sun, X., Zhang, J., Lin, X., Baghdadi, N., Gloaguen, R., & Thenkabail, P. S. (2016). Estimating the Biomass of Maize with Hyperspectral and LiDAR Data. Remote Sensing 2017, Vol. 9, Page 11 , 9 (1), 11. https://doi.org/10.3390/RS9010011 Ya’acob, N., Ismail, N. S., Mustafa, N., & Yusof, A. L. (2014). Investigation of flood pattern using ANOVA statistic and remote sensing in Malaysia. IOP Conference Series: Earth and Environmental Science , 20 (1), 012030. https://doi.org/10.1088/1755-1315/20/1/012030 Zhang, H., Lan, Y., Lacey, R., Hoffmann, W. C., & Westbrook, J. K. (2011). Spatial Analysis of NDVI Readings with Different Sampling Densities. Transactions of the ASABE , 54 (1), 349–354. https://doi.org/10.13031/2013.36244 Zhang, J., Atkinson, P., & Goodchild, M. F. (2014). Scale in spatial information and analysis. Scale in Spatial Information and Analysis , 1–347. Zimmerman, D. W., Williams, R. H., & Zurabo, B. D. (2007). Effect of nonindependence of sample observations on some parametric and nonparametric statistical tests. Http://Dx.Doi.Org/10.1080/03610919308813123 , 22 (3), 779–789. https://doi.org/10.1080/03610919308813123 Zimmerman, D. W., & Zumbo, B. D. (2010). Correction for Nonindependence of Sample Observations in ANOVA F Tests. Http://Dx.Doi.Org/10.1080/00220973.1992.9943874 , 60 (4), 367–381. https://doi.org/10.1080/00220973.1992.9943874 Cite Share Download PDF Status: Published Journal Publication published 22 Jul, 2024 Read the published version in Applied Geomatics → Version 1 posted Reviewers agreed at journal 01 Sep, 2023 Reviewers invited by journal 25 Aug, 2023 Editor assigned by journal 11 Aug, 2023 First submitted to journal 09 Aug, 2023 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-3243407","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":228833104,"identity":"f21d9d3f-ab07-45d0-8d51-6cd21b27b412","order_by":0,"name":"Aquiles Enrique Darghan Contreras","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA/klEQVRIiWNgGAWjYJACCSBOALM+ADEfVIQ4LYwzgAQbSVqYeYjRwt+efPDGzz21efyzmx8+tm27I8/GwHzwNg+DnRxOG848S7bseXa8WOLOMWPj3LZnhm0MbMnWPAzJxjituZFjJsFz4Fhiw40cNunctsOMbQw8ZtI8DAcSG3DokL+R/03yD1DLfJAWy7bD9m0M/N/wajEAqeQ5UJO4AcRgbDucCLSFDa8WwzPPjK1lDhxI3Hgjzdiw59zh5DZmNmPLOQa4/SJ3PPnhzTcH6hLn3Uh++OBH2WHbfvbmhzfeVOAOMWiMHEYSYAY7GLcGqJY6fCpGwSgYBaNgpAMABRRXsQmn0OwAAAAASUVORK5CYII=","orcid":"https://orcid.org/0000-0001-5790-1684","institution":"Universidad Nacional de Colombia","correspondingAuthor":true,"prefix":"","firstName":"Aquiles","middleName":"Enrique Darghan","lastName":"Contreras","suffix":""},{"id":228833105,"identity":"73cde751-a632-4165-901a-4d9179bf8fbf","order_by":1,"name":"Darlley Stiven Taborda Lozada","email":"","orcid":"","institution":"Universidad Nacional de Colombia - Sede Bogotá: Universidad Nacional de Colombia","correspondingAuthor":false,"prefix":"","firstName":"Darlley","middleName":"Stiven Taborda","lastName":"Lozada","suffix":""},{"id":228833106,"identity":"65733a59-58fa-4df9-aff8-04041cb766a6","order_by":2,"name":"Nair Jose Gonzalez Sotomayor","email":"","orcid":"","institution":"Universidad Nacional de Colombia - Sede Bogotá: Universidad Nacional de Colombia","correspondingAuthor":false,"prefix":"","firstName":"Nair","middleName":"Jose Gonzalez","lastName":"Sotomayor","suffix":""},{"id":228833107,"identity":"be52725a-7ef6-42cf-8160-163ac91ce3e6","order_by":3,"name":"Carlos Armando Rivera Moreno","email":"","orcid":"","institution":"Universidad Nacional de Colombia - Sede Bogotá: Universidad Nacional de Colombia","correspondingAuthor":false,"prefix":"","firstName":"Carlos","middleName":"Armando Rivera","lastName":"Moreno","suffix":""},{"id":228833108,"identity":"4c682341-4c9a-4601-8ae6-1c45e2957160","order_by":4,"name":"Jesus Efren Ospina Noreña","email":"","orcid":"","institution":"Universidad Nacional de Colombia - Sede Bogotá: Universidad Nacional de Colombia","correspondingAuthor":false,"prefix":"","firstName":"Jesus","middleName":"Efren Ospina","lastName":"Noreña","suffix":""}],"badges":[],"createdAt":"2023-08-07 22:27:25","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-3243407/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-3243407/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1007/s12518-024-00579-2","type":"published","date":"2024-07-22T16:16:15+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":42378114,"identity":"dd113dea-0fa8-4316-a7c2-edc8e9c3edf4","added_by":"auto","created_at":"2023-08-30 16:12:18","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":171259,"visible":true,"origin":"","legend":"\u003cp\u003eSpatial distribution of the response generated by simulation using a Gaussian variogram (first simulation-sim1).\u003c/p\u003e","description":"","filename":"Figure1.png","url":"https://assets-eu.researchsquare.com/files/rs-3243407/v1/3caa11f463e28c557bda7a65.png"},{"id":42379086,"identity":"61d25b04-6fb7-496f-bb92-1b0207040e13","added_by":"auto","created_at":"2023-08-30 16:20:18","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":85629,"visible":true,"origin":"","legend":"\u003cp\u003eSpatial distribution of the response generated in the first simulation for different spatial lags.\u003c/p\u003e","description":"","filename":"Figure2.png","url":"https://assets-eu.researchsquare.com/files/rs-3243407/v1/a952f3d1189e93e078764199.png"},{"id":42378108,"identity":"312add5c-3979-4b9c-ba08-eafe1bf8e08f","added_by":"auto","created_at":"2023-08-30 16:12:17","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":71059,"visible":true,"origin":"","legend":"\u003cp\u003eDistribution of the theoretical F statistic in the presence of spatial lag in 1000 simulated variance analysis.\u003c/p\u003e","description":"","filename":"Figure3.png","url":"https://assets-eu.researchsquare.com/files/rs-3243407/v1/5793ca604f27f2dbad6a056e.png"},{"id":42378110,"identity":"994637a4-1e18-4f91-b3a5-0dde460cb9c9","added_by":"auto","created_at":"2023-08-30 16:12:18","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":98268,"visible":true,"origin":"","legend":"\u003cp\u003eDistribution of the relationship between the probability of the theoretical F statistic and the probability associated with the quantile in the presence of spatial lag and its relation to the Moran index.\u003c/p\u003e","description":"","filename":"Figure4.png","url":"https://assets-eu.researchsquare.com/files/rs-3243407/v1/dbbb60cdcb5b1c5a05750d4c.png"},{"id":42379085,"identity":"5f9dbb32-160d-48b2-9fa8-d1b74490804b","added_by":"auto","created_at":"2023-08-30 16:20:17","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":104947,"visible":true,"origin":"","legend":"\u003cp\u003eDistribution of the relationship between the probability of the theoretical F statistic and the probability associated with the quantiles in the presence of spatial lag and their relationship with the spatial lag coefficient.\u003c/p\u003e","description":"","filename":"Figure5.png","url":"https://assets-eu.researchsquare.com/files/rs-3243407/v1/644351d83bbd35351026a1f4.png"},{"id":42378112,"identity":"80bdb0f6-2ded-4340-9b56-24cf1b584233","added_by":"auto","created_at":"2023-08-30 16:12:18","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":155234,"visible":true,"origin":"","legend":"\u003cp\u003eDistribution of the O\u003csub\u003eT\u003c/sub\u003e statistic in relation to the quantile of the theoretical F distribution in each ordering according to the Moran index\u003c/p\u003e","description":"","filename":"Figure6.png","url":"https://assets-eu.researchsquare.com/files/rs-3243407/v1/e25fdfd68a3aaf21fbca599d.png"},{"id":42378113,"identity":"aeb650f3-168a-4c0d-88a0-4ed984cf3565","added_by":"auto","created_at":"2023-08-30 16:12:18","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":181545,"visible":true,"origin":"","legend":"\u003cp\u003eSpatial lag coefficient in relation to the F value of ANOVA treatments and the Moran index.\u003c/p\u003e","description":"","filename":"Figure7.png","url":"https://assets-eu.researchsquare.com/files/rs-3243407/v1/d2f98ef61a0c3b4917771a61.png"},{"id":61596251,"identity":"3494bec6-8f50-4cf7-a2b2-af29a0cf4ccf","added_by":"auto","created_at":"2024-08-01 17:26:01","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":1176640,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-3243407/v1/228de4f6-d50f-4320-a8e7-9c8401e524db.pdf"}],"financialInterests":"","formattedTitle":"The effect of spatial lag on modeling geomatic covariates using analysis of variance","fulltext":[{"header":"Introduction","content":"\u003cp\u003eAnalysis of variance (ANOVA) is a very useful tool in the biological and industrial context, and it has undoubtedly gained ground in modern fields such as geomatics, where data captured with remote or proximal sensors predominate. ANOVA has been one of the standard tools used to analyze experimental data for decades; however, currently there are so many application modalities that it is not enough to appear in the methodology that applied an ANOVA (AOV, ANAVAR, ANDEVA) to have complete clarity of what is behind the method and its proper use, since aspects such as balancing, blocking, the presence of covariates, temporal measurements, geo-referenced information, discrete or categorical response(s) of univariate or multivariate nature generate their own analytic strategies for estimating parameters of the associated model (Christensen, \u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e2019\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eCurrently, remote sensing data capture technologies, and drone images or special cameras allow the creation of spectral indexes (Adak et al., \u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e2021\u003c/span\u003e) that are associated with data that usually have very different structure to the experimental design, obligating users to do some kind of aggregation to experimental plots or other units in order to be able to treat this information as a source of additional variation in the ANOVA or simply as the same response. These data that involve some form of georeferencing, may also incorporate some spatial dependency structure (and many other temporal and space-times). Regarding different aspects such as analysis of vegetation cover, climatic and meteorological variables, among others, it is important to highlight that these vary due to many factors, including the time scale, therefore, it is necessary to develop evaluations with updated information from two or more periods, which allow comparisons to identify temporary variations of the object of study, before the attribution of the cause of said change. This process is known as multitemporal evaluation, which consists of comparing the object in the same geographical area in different periods, in order to know its variation, for this there are several evaluation methods (Gao et al., \u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e2012\u003c/span\u003e). If these aspects are overlooked in the inferential process, it could invalidate its application, especially if they are analyzed with the methods usually studied in basic courses of experimental design (Kenny, \u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e1995\u003c/span\u003e; Kenny \u0026amp; Judd, \u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e1986\u003c/span\u003e; Zimmerman et al., \u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e2007\u003c/span\u003e; Zimmerman \u0026amp; Zumbo, \u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e2010\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eIt is surprising that many users of the technique incompletely report aspects for consideration prior to the ANOVA (experimental design) in the ANOVA\u0026rsquo;s own results; and they report insufficient details of analysis after ANOVA (interactions or multiple comparisons) to be certain of the scope and interpretation of the results (Acutis et al., \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2012\u003c/span\u003e). This allows its statistical software package to set the tone for substantive data analysis decision-making usurping in some way the preponderant role of the users, perhaps assuming that by the fame of the software or its recognition it can be taken for granted that everything developed there is directly applicable or interpretable.\u003c/p\u003e \u003cp\u003eSeveral documents highlight the misuse of ANOVA in two-way designs due to heterogeneous variance (Rowell \u0026amp; Walters, \u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e2006\u003c/span\u003e) highlighting defective data (Pearce, \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2006\u003c/span\u003e). Defective data includes the non-specification of the type of ANOVA (Abubakar et al., \u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e2022\u003c/span\u003e), confusion in the choice of the type of effects (fixed or random) (Bennington \u0026amp; Thayne, \u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e1994\u003c/span\u003e), assumption too quickly of the weak argument that not rejecting the null hypothesis implies a strong conclusion (Rong, \u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e2000\u003c/span\u003e), inappropriate use of the terms associated with the error in the construction of the statistical F or the type of sums of squares (Li \u0026amp; Lomax, \u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2011\u003c/span\u003e). Such defective data even includes the belief, supported in recognized literature that only a couple of assumptions are necessary to make use of the technique (Acutis et al., \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2012\u003c/span\u003e). This ignores the supposed key to independence (Lindman, \u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e1992\u003c/span\u003e) in which it is closely linked to the reality of the data captured with sophisticated sensors or equipment related to agriculture, such as performance monitors, spectrometric or other techniques for capturing data from radar or satellite images, or images from spectral or hyperspectral cameras, etc. All this information is usually converted into indices that usually relate to the physical-chemical and microbiological properties of soils or water as well as the vegetation, specifically in the development of the crop. In turn, this information is quite often incorporated into the structure of the experimental design in order to associate the data with each observation or experimental unit and thus, using ANOVA, evaluate the effect of the treatments (Atik \u0026amp; Akdemir, \u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e2022\u003c/span\u003e; Barbosa et al., \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e2020\u003c/span\u003e; Ding et al., \u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e2020\u003c/span\u003e; Firozjaei et al., \u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e2020\u003c/span\u003e; Stoy et al., \u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e2022\u003c/span\u003e; Volcani et al., \u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e2005\u003c/span\u003e; Ya\u0026rsquo;acob et al., \u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e2014\u003c/span\u003e). In the documents cited above, how SIMPLE (Simple, Informative, Meaningful, Powerful, Logical, Effective,) (Mcintosh, \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e2015\u003c/span\u003e) were these analyses in general? In the case of data which in many cases have been shown to have spatial or temporal dependency, shouldn\u0026acute;t an adequate description be made of the essential aspects of the ANOVA to involve the dependency structure present in the data? Knowing the weaknesses of ANOVA in certain contexts and of certain types of data can very surely enhance its benefits.\u003c/p\u003e \u003cp\u003eA form of spatial dependence can arise with remote sensor data or similar technologies, such as spatial lag or overlap (Shukl \u0026amp; Subrahmanyam, \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e1999\u003c/span\u003e), in order to find data associated with experiments where some design model was adjusted for the comparison of treatments, using as response(s) or some performance indicator evaluated in field or some estimate generated by some type of sensor (Wang et al., \u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). Some of these indices are added to each polygon that make up the experimental plot to finally incorporate them into the model, sometimes as a response and sometimes as covariates (if treatments do not affect this covariate, or a multiple-slope model has been specified). The literature points out how these indices show some kind of spatial dependence (dos Santos et al., \u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e2021\u003c/span\u003e; H. Zhang et al., \u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e2011\u003c/span\u003e); however, the rituality of traditional analysis mainly by knowledge of the methods that usually do not treat spatial dependence allows many excellent data captures, invalidating the work of field or laboratory tests because of the analytic technique (Gotway \u0026amp; Cressie, \u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e1990\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eWith all the reasons and documentation presented, it is important to understand how ANOVA can be invalidated, at least in the context of spatial lag; but the same results extend to other situations that may have the same effect, such as edge effects, competition, poor specification of a model, etc. (Christensen \u0026amp; Bedrick,1997). Using conditional simulations, a Gaussian variogram was generated with fully defined parameters, and 1000 simulations were obtained to generate spatial dependence. The final goal was to show the invalidity of the ANOVA when spatial lag was presented using the Moran index and the spatial lag coefficient, although the other assumptions usually found as the most important in the context of the ANOVA are fulfilled.\u003c/p\u003e"},{"header":"Methodology","content":"\u003cp\u003eIn the linear models treated in experiment design, Gauss-Markov dominates, written as \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\mathbf{y} = \\mathbf{X}\\varvec{\\beta } + \\varvec{\\epsilon }\\)\u003c/span\u003e\u003c/span\u003e ; \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(E\\left(\\varvec{\\epsilon }\\right)=0\\)\u003c/span\u003e\u003c/span\u003e and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(Cov\\left(\\varvec{\\epsilon }\\right)={}^{2}\\varvec{V}\\)\u003c/span\u003e\u003c/span\u003e, where \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\varvec{V}\\)\u003c/span\u003e\u003c/span\u003e is a known matrix representing some correlation structure between residuals. The different models depend on the assumptions made about \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\varvec{V}\\)\u003c/span\u003e\u003c/span\u003e, and it is precisely the simplest one dealt with in the design where it is assumed that \u003cstrong\u003eV\u0026thinsp;=\u0026thinsp;I\u003c/strong\u003e (an identity matrix), that is, it assumes independence and identical distribution, something that cannot be assumed so quickly when data is associated with information from remote sensing technologies. To deal with this limitation of the model, an extended model by spatial lag is used written as:\u003c/p\u003e\n\u003cdiv id=\"Equ1\" class=\"Equation\"\u003e\n\u003cdiv id=\"FileID_Equ1\" class=\"mathdisplay\"\u003e$$\\mathbf{y} = \\mathbf{X}\\varvec{\\beta } + {{\\alpha }}_{\\text{s}}\\mathbf{W} \\mathbf{X}\\varvec{\\beta } + \\varvec{\\epsilon }$$\u003c/div\u003e\n\u003cdiv class=\"EquationNumber\"\u003e1\u003c/div\u003e\n\u003c/div\u003e\n\u003cp\u003ewhere \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\mathbf{y}\\)\u003c/span\u003e\u003c/span\u003e is a random vector of length \u003cem\u003en\u003c/em\u003e denoting the response of \u003cem\u003en\u003c/em\u003e experimental units where the response associated with the spectral index, or any other index could usually be involved. \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\mathbf{X}\\)\u003c/span\u003e\u003c/span\u003e is the known design matrix of dimension \u003cem\u003en\u003c/em\u003e\u0026times;\u003cem\u003ep\u003c/em\u003e and incomplete range (as a design model), \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\varvec{\\beta }\\)\u003c/span\u003e\u003c/span\u003e is a vector of unknown parameters of length \u003cem\u003ep\u003c/em\u003e, which, in the case that develops consists of the effect of blocks and treatments, \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({{\\alpha }}_{\\text{s}}\\)\u003c/span\u003e\u003c/span\u003eis the coefficient associated with the effect of spatial lag (overlap) and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\mathbf{W}\\)\u003c/span\u003e\u003c/span\u003e = (\u003cem\u003ew\u003c/em\u003e\u003csub\u003eij\u003c/sub\u003e ) is the matrix of spatial weights of dimension \u003cem\u003en\u003c/em\u003e\u0026times;\u003cem\u003en\u003c/em\u003e, where \u003cem\u003ew\u003c/em\u003e\u003csub\u003eij\u003c/sub\u003e represents the effect of unit \u003cem\u003ej\u003c/em\u003e on unit \u003cem\u003ei\u003c/em\u003e. A conditional simulation was generated from a Gaussian variogram (Gr\u0026auml;ler et al, 2019) with parameters as described in the following R script:\u003c/p\u003e\n\u003cp\u003e\u003cimg src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAA5cAAACpCAYAAAC/DcmlAAAAAXNSR0IArs4c6QAAAARnQU1BAACxjwv8YQUAAAAJcEhZcwAAFiUAABYlAUlSJPAAADNzSURBVHhe7d2/ztzG+fZx6z0MQzCE2MegwvjZhQo7B6AiVqXKgFMLTuNSjQPXCaBKlZxCB5CoUOEELnQMcSAIgk5D7170XvKt2/eQ3B1y/34/wGD32YfLJYczw7lnyN0bbzc+AAAAAACgw//bPgIAAAAAsDeCSwAAAABAN4JLAAAAAEA3gksAAAAAQDeCSwAAAABAN4JLAAAAAEA3gksAAAAAQDeCSwAAAABAN4JLAAAAAEA3gksAAAAAQDeCSwAAAABAN4JLAAAAAEA3gksAAAAAQDeCSwAAAABAN4JLAAAAAEA3gksAAAAAQDeCSwAAAABAN4JLAAAAAEC3G283ts/3duPGje0zAAAAAMC5WCAcfIeZSwAAAABAt0VnLpeMegEAAAAA61gjhmPmEgAAAADQjeASAAAAANCN4BIAAAAA0I3gEgAAAADQjeASAAAAANCN4BIAAAAA0I3gEgAAAADQ7aSCy3/84x/D761U6Y9//ON2qTYtM7b8J5988u7/1+5///sfebHhMvGf//xn+woAAACAfTBzeYUUUH388cfbvyCfffYZgw4AAABAh5MKLv/0pz998Pbt23fp3//+9/Y/8/zzn/989149x+8psPzll1+G5z/++OOQV8fgWeo///nP21eWNXf9//3vf9/lh+g9mtUFAAAAsBtmLq+ILhVWIPXll18OQaWCeXzwwR/+8IchP7755pvh7y+++GJ4BAAAADDfWQeXY/do/vWvf90uVdM9dvk9Wl/me/I0AxbvU4wp07LVctUsmu8T9T2i+T1KnknTPvm1altjfuT/a3//9a9/DZfDzpnVjfevKmkb/Fo1sxfvZ3XyPkX+31dffTX8/fe///299yjl+x+rY6VU8f92Wb/97W9/GwJvBeBT5QcAAADA+6525lL32GUKSKqgTRSkzLlPUUGWlq3o9VbQosBPgc+Yb7/9dvvsgw++++677bPf+DVtZ56VfPjw4XuPY7Qd2p5I68yviQPueGmpaflWIDeX8qs6VqK8XpoCTHn06NHwCAAAAGCesw4u8z2a33///fY/82j5/L4qaIt06eScz1Kw5e3y5Zby/Pnz7bOaL1nV/ZAVr0vrj4GbnjvAqwJIB4ZTl8J6dtXb4dQKrJ8+fTo8Ol9iqgJO/8/7V73v//7v/4b/mY+Tk+/F1frzLKqX2WX9kS6R9exlNUMLAAAAoHa1M5cKUDwTqEcHT2NBhYIcz2zpPQ5WIn1BjF5TkGL37t3bPhunQMiXrMbAubWuJ0+ebJ+9H1DmANL7o6Bpimdd86Wz2q+xmdsqz7TdY4HcHMrnOGMrWufYtvRyfr948WJ4BAAAADDtaoPLDz/8cPvsV/FLXN68ebN99r67d+9un61jThAaAysHggrsPDPpmc2otT+ZZ0LnBKHmwE+f73saz92tW7e2zwAAAADMdbXB5ZriF+8ote4Z3FecpdQ9or40VebOko6JM6VzaHY1BqTe76W+FEfBc8xPpeqSWwAAAADHQ3C51boUdlcKfP7yl78Mz31Z666/1zklXvb6+PHjd5+nAK+6DDXP0q4h/saoA01tV2+AqXtANVPr3+R0WvOy2JcvX26fAQAAAJiL4HIrfhPqvvcJxp8aWTqgzHz5a9zu+/fvb5+9zzORcdmKg9Bnz54Nj5GC77mzhQo0vf9TX2A0Rpfp6tJf7evUFxEtyQMNt2/fHh4BAAAATLva4DLOfMXfY6zuWZwrzn7+/PPP22fvB336YpwlVJe/jgVg3q/WT62IglDli4LIGCgryGvNFCrvqtlJf9lQ6xLbmzdvDo+tn22R169fD4852NXnOdBt3U86Z/0VHUMF4drfXS8PBgAAAK7ZVc9c+v69OKPnb4Pdx507d7bPfr0ktLo/UM/12liQN0f+xtSpoNjB6NRPregyW1FQ5u3XPaO61NWfl4OuuK9ODuoePHgwPGZxdji/18GqZw6dZ06+DFi0bXFwwOasv+Kg+uuvvx4eAQAAAMxzlcFl9RuSCpx0L18PfXNqtW4Ffvn1V69ebZ/tLwZArSDOFGxpOxSoVcGYabl8Sa/el3+axPR6Fdg6P8dm/2LQXdF7q2Pi+y6tNRs8tf5MgaVnLfPPnwAAAAAYd2PTSe+LqDY0GyQLrAo7+OSTT4YASrOKreAv83tEAe8u9zLqOCvwWurS3lOhS2E9KyvKHy6JBQAAwCVbI4bjC33OlGbZHCS2vsinosDQgdRXX331rlCZ7q+sLtn1bOelXS6qYDsGlqpcBJYAAADA7pi5PCN5hk32nUmM64rHTYGlgs7Kvp91yjyTq0uB432aAAAAwCVbI4YjuDwjObhcK9jTLGX+2RLdV9nzZUcAAAAATgfBJQAAAACg2xoxHPdcbuleQ2Xw2DepAktSWVOZG/tZGt1bq2Vy0iz2JXC9y0mXK+MwlNfVMdCxuQbe37GfKML1cvnAaaCvBpw+gkvgwNSJ1ckxX3qc6eTp3wvNXrx4sX2Gc+aOUu/v3gIALlseCFxjEDauX+lSBrJxWASXuEjqrLtx1Ozf0uKM2y4jqAos//KXvwzP9RMyugyh+jkYNegOPvWFQ1ouptu3bw//O3f6EqW4X/4G5EvmwQWlzz77bPvq8ei+7XgMVC7P0b4zGt5vftsWeJ/OQ26r5qRLm/33/p/ClTTajnx+1N96fYkA0Pua6bs9ruUqFiyH4BI4EAW8Dix1Upjz26Rq2KufRuHnUs6TBxc8sKAvygIAYEoeaPZA4A8//DA89vjiiy+GR52TvH79FrqcwiAozstJBpdxVkhJM09+bYlZqDwa11qnL0HIo0KtEXIt7xEuv9eXu/neuTgCtuvyfq01OugZkXMYZVLeaVudtG/xtZzneXmnmD/i1+PPqejS0vgepZhHuTzElLfDxyk2tpphzO+rLnP87rvvhkf97EkrOPT7/a3AHpmMKe+zxNkwp2obxP+XXNfiup0vOjYxj5wnzgv93/xazrdWndlXVR6WGDWv8jGnns/R7JhO2nMGFk5ddQzycY9aeVu9J5dLp0quk1P10e1oTlPt5pw65nrieqT/x+XHys4u+bMm56f3IfJxqf5n+l9sEzKXm1b75ONz6P0W7/vYPnr7qu33eyPvb3zuchCPucX993PnZyxPzh+/pnVXfMzGyl6LzlMONHLAEYMQp2PM/juPnMbk+qhUlTMfp7HzcJWfsfw4tcr5rpS/ud/gb/DvrSsqI9pH7W/8VQBdVeUB0KX2A1diU2C7aTULrert999//259VdoU9O2S+9l07Mv1On355ZfbJd++3VS04bVNpdu+8iuvIy4rWl5Jr8d1bhrj9/7W+/dZXtvh1yp6PW/TKfI+jKWY59X/Y1J+WfX/KjlPc15XKXKZmEpxm8SfM3V88nqqpG2IcvmJqaov/l9V1+K6Xd60jriMPi+/Zs6fuXUm82fmfYz8mVXK+b6rqfZHScssxfnYu91LcnlyHanE/Mip2pdqOadcRseOQVUuXObGUtymXHadxvZ3bh2L5bf1OVX+7FqH1+Y8zWXdr4/xvre4LdByWWxzDk2fO5bMx3WsnEc+trk86P35b/FyU8vHY+PjUpVhf/5SvB3HOEYSzyWtelPJ+RlTPl+N1UenmP8ut6001rb0iHnRw3lTnduW+gycLpfTJZ3UzKVGX3zZ4KZAa0+HtGnMhteW4NHtTWV6t3591lI2jcwwcq716jNEs2j6DO/H69evh0fZZXmNWm0q+PA8jyL5b8+O7cOjdVPJo6j7iO/VvvgYbE6O21ffF/czLu+U+fVYZuKxdtK9flFcRsckijMavj8tlplqu/J9lD/99NPweP/+/eGxxe/3Nihf4nqV4m+bKn+qezO9/5q1bc3IqK5tTibv3uPyl2kd3kdtjz5Pr8X39I6c7srb7OR64Xy2aiS5Sh6F9sxiTrE8Xfu9ea7DuWzG9iuKMyqxjDrdunVr+9/fxHKp5Pqm9+eyluvkVH3UyHz8n8tOyz51TMvpdS3j5bVP8vjx4+HRtD9ev5eN6fPPPx/+J1q2Kr9VatX7OZ49ezY8qo3wenTctV+x7duH2l6VHa0rb+PTp0+HxwcPHgyPh+Iymsv0klQetM8uB6onKh/Oz1evXg2PpuX1P5czLa/3ah3y8uXL4VG+/vrr4fHJkyfDo7ls+TMlzpaOpdbM7anQfikpP3y8dPwkzyqqDis/c9vgfMmXlerKEv3feZ3LhVI+D1TLeHt+/vnn4VGWrMM+3lP9iiluUz/66KPhMfrwww+HxzV+Ux2X66SCS59YVOFj518dg9i525cDFVX4OPWvz+o9YUZelztN6nTpM27evDn8nU8iuyzvRiR3UBRUar9ivp0iNfAWLw1Uw+WGuKVq3NSAV1+IM5feq3XE8hCD+KW48V76i3hal9pqv3zifPjw4fCYqU7FE6TyoHUC8bFyh8P10WX2zZs3w+Mh6HhlvSfXMWo3HDBVn31NVI5dh3NZUflxvXHnzp1bUUctX9YluZOmv/Nratem2oe17FvHtL+xbbp79+72WZvbiainfduX9tN1XAOy6uTquPvc1KsVDPl+5KqcrCWWUQfVtmR9V35qv9yB137q2LrzHoNFcT/I/QCVf9UL500sKy5b8fwqDprmlL1zo/zQ8Yllxf2i58+fD4/mOpxvR1B+aj35uO9K21CdO13O1+CAWXVyqTbCZS06ZF3E5Tip4NINwqeffjo8Ls1B2poVXg1VPvneu3dv++z3dl3ejYhOhj656MSvjkzvfnm0birFQGwXcRRubvCm/VUeifbRI3pVJ+yUtYK2Xh5ZrTp8rketz557QvKMVLTUyWxt2veqDOeUgxlTmXVg6bzO8j0/eWZfr3lg69x5EKFVf+/cuTM8uqPsn8xRGTrXTso+dUxtVmt/87JaznVM7xsrL1q2Kr9Vitsby6dSbIv1WXotUx33cfYVP7nt93tjUn6p05tfj6pgyPu85kBRJZbpNctobjOn9jMHhGPnd223j1UsO8rfvF9q66ryktNa56ylVDOrLvN5212Hc5lU0v9abfsa9q3DkQbvdF5Su7FvfwxY00l+oY9H8paWRwbPlUfLPdPr0d9zGp3c5SSuE0UOcNwJyx35feUTjkeyT5kDbAffmevRIU+ch5IDOgeAS1IH3J1q5WFVZtWRy7MFuWMtl/LTMfES/YpnZXoHf/T+mIdKxyjHh6pj6iDGq3NUnrXPS1yaWH3Ri8q181Wf1RosiDM9cft6qfy4TXcw5FmnQw9cxdtUzpln55yPxwrWT83ctqi3zZJ8a5Fv81qSAkutV/Vx6cCyqgtL5Auuz0kGlxjnIPLRo0dDxVdndomZgUPcc7kvNaIezYuBpva9Z3t8/4koaPdntDpb+1rj/hUf71bH1iPyrY7xOVJ51/HSZUw+Vkqtjq/yPZbdVsr36MTAMl8OGekeT+Vv3JZcdnw53CXwZVOtDoevDunZX9Vn5anyLebrMcrxIeuYgirvq8uzPje2HS7/c5JnJzVQFtu2al/y5YIWA1MFofm4x212Ul6ojc6vZ75CR8GQttXbif1olkt576uafKtMDtZ9zptKa5yz1tZql0R5kctkTD1tlsqv8kziOqvyvE8dthhYLjkI433Pt2yJ27hzLA84npMKLl3AfSlVlL+oYx++P6x3BjPfJ3Joyid1YNVJ8G8TjV1KeyrijHS+x0ENbqsDlynQjMvue7+EPtMji8rP1qWRSxgr2z3cqc0nIfGXCJzCSWGpOuOBhN57ZMYoL2Ng2bo0SVQW8yVY+fLyQ8/ErMl1uDWz71sb3NZ6xlaDQGMdP1Pee7DsVPLtGHXMQZs+W23dnLxr0Xpi26a2KJZPpYo6sg743En2+WYJMRjy7Noxrr7xgEnvZaDxctRj8aWzOk4qN2veAnSKfH6N5VTl3fWoqsNLcPltDdIsQeVL/ZU12kb3HzVhkS31pUG4LicVXPpb8fIIqUZP1eEY45GesZO874/RumIjo5OoO5OR1+XLT2XOthyCK7oaTDWcYx3guda+59KNvGi7HSjoWPj1TMvouObZyXhMqmMeb0xvHS+PyEnsLKsRj39Xl4rEQFnLTnX+4ij9kvxFIiq/cRt8IhJfLnUIa9cZdwBzkO7LYvNx0PJVGc7JnW+VRbcFGh1eol4tRfs2p51bk+qwrxzI26A66nrj/NTynslVHa86d3HW2HUtDx5oGbUZEuttFAPfqfq4i7XrmNZTHU/llfdZ+ejHqvxWqafsaj+1bx50U/KAZm6Lezj40Xq1fu9ni/LEdSDOqvZwPmkbYln0Z2UeOIltud7nNuiYHJy73FTBuo5lVV5y6g22TcfJx6yq//tSPY9lUWXWxyAPtrucqQ632qCqzXB5VH6ObbvzOy7jOiRxQmOfOhz3bdf+l9oW53+LPscBeMxTtU0+d7cCWtcTpaXqJC7ApgB302oWWlXz94U2HZrh8fvid3jEy20qyPaVmteTk1/X55s+Ky+n5NfjsqLPjp/v5f69/Y0jPfr9suvymf439v9T5H2aStY6XjE5v7JqWSfnWfW/KlXlSq9VyyppuzMv39reaNPID8tWn5u16oxStR3+3xRvQ1yHP8tymfXfOfn1XGeyqf2eUx6U5uRxZSwvnao8nauVPzlV2++8UZpTLvblPBjLQ29HlTZB+Xap31TLOcV9ifs4llrlyHUsp7FjNmd/x8pFXPdY+W39T/kV15eT26pDcj5G8diM5ZXeO7eOxHVW5SbT53r5VhnYR+sYxPbG4jbEFOt25LJj/izvr/PAeebP1Oviz4vlQH+39t+ft2T+RN7+ucfY26M0Vm7mauW/U2u74nZUyfmdtd4X83dq3U6tz5jSKp9OcVuy2CaOiXUxp7HjFo/H2HbgdPn4Lenk7rnU7NmmgG7/+tWm0L/3W189NOqzaXy2f/1K668uK9UIX7XsWt9muytvm2cJzoFGyDaN0fav32xOnOXspY6X/teyKcPNEXodqyl6f/W5em98vVrXriO7HunOMyC9VGeqPNqckPaeZd7X2nWmKg/6W8cxvn4pX9LRcqyZS1N+53ZadKyrEW4tr/JYibOUHtXP9Fp8vVX3lpptydasY8qvqk0UvX7o9l1lS8cxHy8dG7/WmgHaldapcqS2dtdL/TyrtAR9dt7f1rHV+aZa9lTOw76q6ZBXrMwVr/jppTzPx0F1tFUfVYdb9UzlvVWe9L6qjxBpmXze07apzYqvt666OAWttld50+pjZUvWSZy3G5vC9PvStCNPty+wqiZdtqBLDFRhdz0JXSrluxquQwcQa3E5UkO+Vifx2FyORZ0qnZSAXajzr3uKLqXe43ppkE3tvYKCXYIznSt2fc+10KWJOn+e0jnU5701+4g4Lurk+Vojhju5mcvqfg6NkLpDTmD5K1/bfg5f5BP5+nxdyx/Fa/Uv+UsI1PB6BkT3jFR5AbSo/mgk+dzqPVDxF6/s8kU+bi+P8eU/p05BnM4rp3YO1RfF5Jk9XA7qJH5HM5e9tJqFVvXe9eE5/bi9R+FaxWvblc7x+nYdw7gPOen4XwsdP+3ztZdrzLMJKofy8v0R7sEDlpLPAZugY/ufaT4H0ma+L+an0inRee6azuvXhjp5/tZoN07usliNzPvbGqOxa+KvRcybc76kUrOU/lbJ6JIu8QUA/J5mOfzNl7T5y3AfTJbohwG4HmtcFns291wCAAAAAJaxRgx3cvdcAjh9/u2snm+9jb+/FZNm6K+N7jXXvuueqWPRFQXahnj/c6QZp3icnFrL+/7qnI79TbfANXB97fld0kPU4diumNvDVttySnQO9Pb73FXtk8RlcZlO4Vx+CgguAeBI3AnxD1VjP+6wOS35Uz/uLKy1fiyrFRApXXuH7xTdvHlzeNS3BtutW7eGx3O4FSpuo39qpdqnTOXxHILnS5XbhqljkduVqQEWfQmplvOAw7UhuARwFPqqfF2G4aT7iK+JAkvfe6ZOiPLgmF/jrnu4tQ2te7n1Td3xeLV+M87022hxed03vzTPBGTKzyVO6upwVIG/1q/jB1yyQ9Rhi531jz76aPusrTWIcMx66UCz9Xue+r/y8cft73Pquye4kuOwWucMf3t/RWUqfxeM6kK1vO4h1zH2wILed43nCoJLADgwdYwcWCpIu9TfdV2bf8pCXwzjDrA7btUXw+1D6/O6lfxTQqf4I/X4TT5uSvwG3+lxIBZnAD3z5xnMTAFZq36rXT30lQUKJBxMiPelFTh6oE7vUZDSc+ky9pPbBh+/fHWDypLP1bFN8U/rtGY8dU73uegYZfLYCC6BM6QGzaNmfu5GUY/6uzUK5/85TY2cev1OU/LyrcZ3Tb6Uceyk7WWymH9KGnX0ax6B1KP/9nPno0dG4/LZw4cPh0edfDQ7kHkdXme+NHNsVi5vv1M+uXm7c3I5OnXKA3XM1CmI3ziqjptP/K38n0uzuPm3lRWguFN4aC4H3i+Vj3jsqnJRHeeqzrvc+PjnejzWOarKp1+rtikvP9UGnSPtt/ZN+ennSua/q/rmYxHTWFmO61eq1hnl5ZWq43QIMajMqhlMlUvXPQ3MxQBBKQZ5lsubkl7LvFyV135fpVV+x/ZNPKioqyOOkf9uG5wXOZ/yNrlc6nU/9/k9lqn8vnhOjCm3KXEdVRn29k2V7zGePc4eP348PL58+XJ4tKdPnw6POqfEc4HOOSprmvFstY1a3oORVXm7aJtM7qbVLLQqADP4NzI3Dd67+qe0CVZ+97f596haaXPC3i75m03jWS7rlN9TLaOk7Z3ifdJ29tJ2+bMr/n/+zUhvQys5P53POf+1vphn1X77OGi5Fm+flmltU5VPY9uvbY1yWXGa+zua3o9qHytxn5bgvK+2d9dt25XWrXRocZ+9DTHlMpHLZ0w5b7zO1rqVKlNtRN6m1vJLlQsf+9j2HYO3Ix8D/R3rad7vXeqwtOqxU37P2PHNxypbug7vKu7rXGPlM7cdPlZV2fF7luZ9qo7t2uJnt/JJx9xcdpxPTrkMxn3J/8spc71Rip8d69MaWsfB+VLVjdhmjvE64v6cEm2b0pKYuQTOmEY8Nw3Wu9ExXX6xaSTf3Q/36tWr4VHu378/PG4aT7Ui79Km4RtezyNr+lvrlmr5zCOY+vy8vEb38mjmLjxKOpU8gqzRyU0nbXhejUJ7NDJeJqfl/PurcfuVXy3Kf+W980k38WsbvN/V5a4///zz8Pj1118Pj2O0Xm2Tjqe3x/vl2U/T6Gm1/U6ff/758D/zpVlOLkPHUh3PKvl4erS4mt3wpXZrXG7szx8rF2tTOVP58rFzfa/k+qi/xWUl07olvsd1Ptcl1UuV0bgtSi6jmdsUlbW8vF6P69+1zp8qtRFuc/23ymXcb4ttkF53/viY5Rku1QFfshfz1Mtnem8uO3H53Kbsqjo+VcrlaC7PLrX2ryXmpZLr7vPnz4fHY/Js2LNnz4bHY3A/IuaT67DPlZGW1zFwu6MyGMt4nsnz/5ximc/9Al3J43ORb3sQXQatchuvUlE5qspXlaa4bN27d294NG9rdR+tzz15tjPzuf7FixfD4zUguATOmBp4BVJu5HRC0MnKDaEbPTXCaiT1//yj5e6Aq1MTG3qdQCR3XLW8O5um9+n9OonES0fEjbaDqkNxMO3Pj9TBcgfDvFzeX+VXXtaUnwpQ4+VPzl91euNJ1HxMPv300+Fxik7G8dJZr78VHEh1mU4+LpfC92dFU5ej7cud+dzJOTSVuxg4+4tX8iXW2sZ83KfKgdatdUWtgZBHjx4Nj7ljrMuJtZ5Ieac2xXUmcl7+9NNPw+MSdJzcscwd2EOKZcXtptsal1PXV9/HqzYolmEdM3e4YwDojn/OUy1fBWB+bz5eWl7rGGtTToHL/FQZjvSe3B7kAOLYVC50rqja7UPQ56vOx3zy+bPi83wMuGIZj22TlsttpT4ntw+RyrL+rzzRgJQHvqtz+RI0kKWyrzpW3aYiuQxJde6puH+2ZPt26ggugTOWT7KtE4JnMO/cuTM8Zm7oX79+PTz6JKcTRauxjfw+dR7jaKGSv3hhanRvjE42OvlNpXxSE5004knbo+a5g+GO1Zz9tZzfysfqJBTt0oFwZzTS+v16XJdedxCs/yvv950hOIbqeFZpl47lkhSgKF+VYjk7hrGOX69WGyHxSgiVPXX+lB9TZV7evHkzPKqe5TaiKs/71PkWtUHHuPdbcmCufW21MR6Mqv7vwai4v25T55YHv1fbkI+B279d2qesOj5V2rcOV4N1kvflnNo9Ofbse5whzKrzdj53tgZfezggVZ9CSZ+R64XKUVW+qtSiwFWDzQos86DXUuYGoZeE4BK4Aj5BVJcQijuH7jy6IzhX7HS29HRa9uWTXry0R7MDrQ6eO7nnSifkOGPh2ZtTv3Swhwc2oqXLmgJLBSgqT8cOLHelbc+d7yXNLVvVccqWyFvV69ip9JUICp56vghkbS6zrTbIs0QxwNq1nLeCs2jXtv8c5EusW980i/XE/FfyYEaL+iTxXJZnP5egwFKBq9qIqcCyqmtz2rRrRXAJXAF/pXsrCHTD2Qo+54r3/eTU+v3EOfa9/+rBgwfDoy/fU0dbHazee4t6zJnlmTLWSYyjuT45a/lTDzCr41klz0rkAZHIHeQl9jkGlmt0cNak2TptuzpPsS4ek/IxbktMMbhc6p5LBZsOME/hHrsWl+dW3XaZXmIArMp7p2rQba7q+FRp35lFX2GT3+9t92BipM/TzFTcR5eHU3FuA1a7cD3OfYOxy2JFfRLfTyzVlQcqB7FcjaUsBpZjZd71rRp08bmn9ZM5do1BKMElcAUcNLY6Vx5F9OUb1Sj5GF+y5SDuVKjD5ns3FCQ4qKwuy9JJRMtVI5RLftmCT0T73oOq/RBt71Sg6kBzbN/OlS/NqsrckydPhsfey0eVXwrOlH/nFliq46V6PXYfUQ+XvdYMRH799u3bw6M6dKi5I+s6Hrm9iMH02ABLpRWcnQvXZwUdc9oy36vXG0yunV9qm+e05+dGx0iBvcrdrpec+lJd543ak6WuPNB61A7NaRt9WbvPKZHPPXfv3h0eW1w/85fqXTKCS+AKKMhwA+0TrrmzohOAG1qd5NzRicvrZKFRQDX4kd6n5fV6ayYhf+4ueu6/codEgaX2vxrdFp/M8v0n1f728IloTiCuz80jtr6kK9/LpQ5QlfeerZVT7rxUx7NKHhiIZS6WLeWDA5hqEEGUpx7RHuukuizs+kUSridKUzNra3GHJg8o5frcw8FKLKPe90xlz8tX/5e4bfvW+SoQcJ0Zu580lokquDsED35pe+Ox0T6pky7+0h9xZ1X/i8srH+Osj8XgrMonvdaz79XxqVKrXk7R+1yGVPenzinOkzyQ53zIZccDf7G+KxCp8nIpPg5j9z2ay2er/pwaz/blgSbts1+rZvWU52rX/YWFHtxVOY/lU+WhKl9VMr1f61H5mRPw+nytc0qsMyp72kaVx6nzqs/1HmC7BgSXwJXwCVONZDxJOfDIMzPV8mqQ1ZjqMfMJQOuL63fq7cjuyx0Zn8x8qWzm/c/bLw5Il7gx351sfU7Vwcu03Xl7WifGKu/dsdYo7aVxmYtl1B3B3tkKlVfXDeVhzFOnU+arCXL5ceAtczroY1xn4mdonW4nJH6jZLw0Pm6T0xJXCCiozuuVVp05JTl48va7TKsdijMtreV1jKtBNC3v17VOL++0ZhC1FJUh7avkc1ks2+LBBAUTcTnXa7eXboerOqP3rtl2Os9P7RtslxDLqvNTKZYzPY8DcA7+VE597tY508eg92oUDzRU50ql3B7qs317SawzLmtTt/s4UFY9nQpCLwnBJXAl1NDHETzTiVqv54ZPy+cOuhr4VmOq92s9VadGDWvPPZe9vE1TDXzOH+fN0twp18lqbKZAn++OkGkfqhlanYhbAZVeP/WO9T5c5jLlWezYHNOxZi61/7nsuDzH13sGfZT/c8un6fOrzrq2bex9c1U/8aPP22XdMSA+NLWTVf6og1tdmq3lleeR8rh1CZ7WkY+ZqZ04lXozRsey1dYpLxyUqM2LXwojOhcof+LrnuXXvufl9fdabafbBm3TLvmuunIulNd5e1X+4usujzoXejA0l3UdAx1bLXvoNrU6t7otHaNBCwXKUtXdS3ZjkzndPSdF8bLAqgDg5OiSOY1m68Q2FpzuQicej+DmjrU6/Hotv47zok6QLne7to6FaN/diTwXmmVQZ5C+DNY01vZP8SXoxxysxTS3f6IBCg94nKI1YjhmLgFgq7pEUB1OBZYaNV0qsBSdbDxKrpOQGnh9Fi6DRuF1XC/xcjfTIEhVZn0/UnUVwynTvVHnts04H6ovaucdWE7N8md6v85FvZeGYj1q+3SMHVhqxvOUA8u1EFwCwIZO3PkeHqU1L2vRSUejhXRoL4vKki7v0uWN53CZ4b70hR2qH7nO+H6kc5qx9YzQNc4y4/DU7u86+6hZTp0rrjFYOTdq+3WML7n9H8NlsQCw5cviomNcmqrghMticQ40Up+/SEUzMly2BwCnb40YjuASAAAAAK7MGjEcl8ViNbqhWYVWSTMxme5J8v+Vxr418xj8u2dzfi7iEDSrpu3hvrzL5bpwLVzHTq3uY3++5+jQ3+h4LDo/uN5W57ldaR1eHwCcI4JLIHEQpxvnAQA4Bp2HfC8oAJyLiwkuGek7P7rRWdPwTkve+BxHk6tvAG2J99zpviFtFzfP4xxd20y3r4SgM35c+kIctZvXcq+wv5RLaYlvk/Zvt/qbpDXIeS2zwAAuAzOXwJYCUgeW+hppvpACAHAMDlr1pV46H+0ySAoAx3RywaVnmyTfk5dH73xvh5eXuLxS1SDn9baWOzW+P0mPc/fBeeS8i/dBOlX3O1XLjeWRt82pNVtSrddpjGdhqhTvc/Fr/h0pqX5eotrn7777bnjU7xJNjUDn7VFg6teWuEczzry21ul7c3K9EG9LPA5eXscx3tfj/POx8XHedXlxOWjlgZaN61ibt8fHW8+dqnyL4rJKMS/Fx2jqeLuu5jIX68LUtuTy0HpPPE4eKKl+KiLn/9z19/CxdxorA3E/pt7jfNRPf4hmevL78jGau/64XCzn5noW/9dal2l7l87bMd4ecVl0qvbJ8rJOOS9by7XWHfffx87rdBnpyZ8l67xStR/Vcq11e59UHvzc64z1rlVePAOs81huQwDgJL1dgFaz0Krerev7779/99zp448/3i71q2+++eZ3y+SkZaJqvU55/afmyy+/LLc7pry/c/JoE0xtl3779pdffimXccrrF+VbtWxMWu/Usi3VsjHFbar+X6W4z/Ljjz8OryuPp0wdB62rx1gZVdL/xceqKrdeh5cVL5/LhPYnv7bP8qJ81d/VNsX1HYqPVStPq+PtfahSzE8vF1/z58R9dNmKvL5WinJex5S333k8lbSc7bL+ffk4tFKsj86vVsrmtD+xTu66/pincTv9es6jqgyYPztuz9q87a06EMuvtZZVynW7VV9a9VzvV8plIh+XmNe78Hpb+1CV6dY+OMX6ItUyVZsnrl+5nuX9rY6DedlWngLAvtwGLWmRtS25YV6XUmxs1ai2Gm+Zsw2xkxBP7vH1sQZ+irbP6xlL+35GPhn7RJNPUvFEmE9ozsN4Mo0n8bgP8XW/pjS2fst5kU/OFpep5H3O69HnVCfcmCfV/zPvRywXlbjeKObD1DrGxOMS1xM7Si4/Lrc+ppGXj2XNyyu5gxWPk3g/tOyuy5uXieVHvE1x2bXF8uN9kLhveXuq1+Pykf6O6/XnxdeUR/EYeZl83Lx+pUjvz2WqVQ6jqgxUdlm/1zmV4r7F98Q8dT4oxbKiz9Y2RTH/c7ky152Y95V91l/lh8t51iorom3Lx31t3halWBa8T9X2eN9aeTHG6815bF63kmg5/63P8/v9OXH5seR9i+UqloV4XGI5FG9DVVe0TXn5yOttHdfW/vnz/P5WfpmWaX0GAOzL7dGSTvaey00D/MG33367/evXLwno/YKAH374YXjcNOjvfUmLLoHcNPrD8/wD6qdK+aM8Ee3L5iQ6PJcXL15sn71PyzgP45fp+It0dInO5kQ3PFcexS/Y0edZXH/88WznoehzNifC7V/70WVC8RtbtW35clV9jvOhhy9Jun379vDY8vjx4+Ex7qtoGzadg+1f+3vy5MnwqHXFMqq6oGOyFN9P6ku5fHxv3bo1PL5582Z4tF2W//rrr4dH74upbqkMLvGlF7tS3sV7aLUNsc6YL33V/sXt1HPnf7w0TWU8tksqrzp2sdyqbH3xxRfvnvt/z549Gx5NdbGishXLguS/e6y9/kePHg2POU91PKpjoM/Odbp1vPaxz/r1HtdvXXapcqL2KLaLFtcVLyFVudGxd/04NLVZ8Zw65xi/fv16++w3S5UNt6FuQ5S/OufcvHlz+PvVq1fD477m1vno5cuX22e/0TbFcrsvn1O9f2o7dDy87tZlsablVeamlgOAYzvZ4HKpE1jkRvnTTz8dHqMYSO3beKuT6YBtLMUT/L58grI7d+5sn7VPyvfv398+q8X35fu04j2MXi7mk058MQ+XEINYdTyWOMG3xABhjAODpffVnKf37t0bHtfgTnI0Vt92Xf7u3bvDYxx4cCd7qgyuparzFgNjdy5V3mP5V/LAU+xwK9j2gIyDzgcPHgyP/ltlJgfh6uSuWZ7XoraratNyivXJ+bNGm35IDkh1PD1Q0tonl3MPRokHW1w/Dq3VZvn4RA8fPhweXQ+q+w57VOeLVpu373l1bp0X11nfnz91b+Y+8nHfdZBhjW0CgDVc1bfFusPz4YcfDo+ZTniSTzzXohq1rXi5a80nl5M1zA1yT5mCJgekDirdyT71AGPOwFIchPGgjt73888/vxc06m+vzx3daiZoDgWqMdBVWtJa65+Tny15ezyws5R91h+vWMizn5HKudoJrdN5oMDlXAYVtP0x6HTQpdRzTE+VjokCVLft2nfvb5x9BgBMu6rg0iN/raDIJ9NW8DlF6/cJaSzlb51cwtzAcC5d7pVHhp3GOlXnilHhZXkWQkGlLwdc8rLetSmIqMq+Upwh+eijj4ZHzbI/f/78XbCpIEJ/u63pmenWZZj6JtS8TUvZZf3+dtSp1FOf/BkqL3F7lKdL6Fm/v5FWfKlzi2emnj59+q7N9zdSnwMHXEox0FQAdsgA85Dn1ThL6kBTs7enEGBewsAjgOtwVcFlnFHIfAmbnMPIchbv3xq7HGjM559/vn32/uVcLTEIry6t6hXvf9TI+ZodGh/z1v2qpg6H9rXalnwP3T7cKd93hst8n9uxKJjyzI0vETzW5YC7cHCY7xdtcRnVbKb21XVP61FnUG2NO6niy9nndhTVqXVgvsal2GuvP7alc+qvlvElp0vcPpD1rF9BuCgI1/vVDoxdLuryrs9TfVQ5WCOPD8GBpq9ImGonL4HqqO+p/emnn4bHY1J5Uxk6x/4JgOtyMcGlO3BjI4y+r0In+7icOhweke6ZXYmjnmNpiU6TttedNXVwHNz1dGDiJYvqcLozFWmE2B0qneRix9mva7u0XG/AqfXH2QR9VhwEEH1mNUsS70mN9/61xJm2MZ6tyLMWU/vrPFGq8tUc3MR7XEXvUbmN3MnQ58Z8WSLvl+CZG22LjuNUp0h10nk01mlfU7xftNoGHcc4Q+J9cjDvuqcgU/ut2ctYVvx//S+ux5elZr4EV+uJ4ra5Hcg8qzo20LDr+tV2VW1aTjF4dh3O+an9VTsTeaY3v+4gWFoDLx7s0nKtPNl3/f6/g3B/UYzKSTyXRLH90vGec4+dZ1WVlpiJ24fbqioP3Zbme/7XdIjzqs4h1XH0+cD3TB+Lt21qthwATsHFBJdudPMXccTAQyd7B49xOQdIelxjtHwt2l5tfwye5sw4jon3FKkz5TxyyuLnaTu0TAw4e8Vv+xMF1XF74r5HOcCO71HKHTctr+3WPucANvIlweosxvWJR/V7Ol4qf86/uH5tV/XFOu68xnyRnkGSpcSZymN9kc+u1EZ4tsLlOaaqbOs1B9Dm8qfjljumXr8GC7xeHb/q+HomNNfFWO71+VUg7FnVXFaV3FntWf9c8Ytw4mdovTHPJNbbuGwcbNHzakBJx87HR4/x/d7+fdavIEuvaVvj+cHHMb43i+X+HGbuo5yHSqJ8yO3rJdBxzPurMivH7he4jK35RW8AsJSLCS7VgckdlYpOEjGAMnXs5l6qdgrciYrUiew96ev9GgVu5aXyKd5zqeWr/FRwU23jPrQ97shVWpejKj924UA5zgpXtD2R9jO/1kPlMOad11+Nniv4rpY9BZ650Tbt+kU+x5wp0LYqD2O+msp17mh6YMuzzub358vUtf5cnvV3dS+z6lcuxz7G8fWqvCr/q7oZ9ax/Lm1H/oyx9rbKe70/vp7XZ3Pa8F3Wr/32a/n46Dh6QMCBV+Zyr+WUD7vwzPOhaTuVFxWV/zzgdwlUbqpznl5r5cWheKBDZaj3/A4Ah3Bj03B2t5w+sR67Eb50ujTSI6nq/OzaWcE0zWj68lN1LHbpSPn4cGx+5Y55FZC1+D0KtnYNSIFToxlizTopyJ8bGOjKCQ1w7fIeXB6XHVGbeE6D3wDOxxox3FV9oQ8wRUGQLyn1ZXz5XpzqEkEFpVpeASmB5a88o7fL5YD6Zk0hsMS500CJggMFBrsEifoyqV3fg8uhcqPzjgNLnVMILAGcE4JLIFGAqRGc6jIpnfire/E821ld2nhNFIg7TzSDu8vlgHqv8nHqUk7glGnwSeVfAaI8fPhweJxDg1RqX5b45mmcP52HLvEyZACXjeASaNBJXSf3OIumQMkzm5E6klqWWcvfKLDcJdhWPisPmbHBpdj18m4PbNGOXC8de5UBJQA4R9xzCQAAAABXZo0YjplLAAAAAEA3gksAAAAAQDeCSwAAAABAN4JLAAAAAEA3gksAAAAAQDeCSwAAAABAN4JLAAAAAEA3gksAAAAAQDeCSwAAAABAN4JLAAAAAEA3gksAAAAAQDeCSwAAAABAN4JLAAAAAEA3gksAAAAAQDeCSwAAAABAN4JLAAAAAEA3gksAAAAAQDeCSwAAAABAN4JLAAAAAEA3gksAAAAAQDeCSwAAAABAN4JLAAAAAEA3gksAAAAAQDeCSwAAAABAN4JLAAAAAEA3gksAAAAAQDeCSwAAAABAN4JLAAAAAEA3gksAAAAAQDeCSwAAAABAN4JLAAAAAEA3gksAAAAAQDeCSwAAAABAN4JLAAAAAEA3gksAAAAAQDeCSwAAAABAtxtvN7bP93bjxo3tMwAAAADAuVggHHyHmUsAAAAAQLdFZi4BAAAAANeNmUsAAAAAQDeCSwAAAABAN4JLAAAAAEA3gksAAAAAQDeCSwAAAABAN4JLAAAAAEA3gksAAAAAQDeCSwAAAABAN4JLAAAAAEA3gksAAAAAQDeCSwAAAABAN4JLAAAAAEA3gksAAAAAQDeCSwAAAABAN4JLAAAAAEA3gksAAAAAQDeCSwAAAABAN4JLAAAAAEA3gksAAAAAQDeCSwAAAABAN4JLAAAAAECnDz74/w3a5Dm71emwAAAAAElFTkSuQmCC\" alt=\"\" /\u003e\u003c/p\u003e\n\u003cp\u003eThe 1000 simulations were assigned to a 10m\u0026times;6m grid for 60 polygon-associated cells representing the experimental units. With this matrix and their respective latitude and longitude coordinates artificially created for the grid (per row[y] and column[x]), three treatments were assigned in two blocks to associate the ANOVA to the simple factorial design in generalized and random complete blocks without interaction. The random generation of this data set ensured that in the 1000 runs of the ANOVA there was no effect of the treatments nor that the blocking was efficient in not distorting the distribution of the test statistic (F) as regards the no centrality parameter. The image generated for the grid with treatment-block-repeat assignment (T.B.R) is shown in Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003e.\u003c/p\u003e\n\u003cp\u003eEach cell associated with each polygon was labeled with the respective treatment (T:{1,2,3}), blocking ratio (B:{1,2}), and repetition (R:{1,2,...,10}). To generate the spatial lag, the cells were spatially varied (exchanging) by two in two in the combinations of T.B.R within the grid to modify the distance between these units and change the pattern of the neighborhood and, therefore, of the spatial dependence.\u003c/p\u003e\n\u003cp\u003eIn Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003e, each image with different ordering (ord0[data from the original variogram with high spatial dependence] up to ord60 [data with all grid elements exchanged and thus with low spatial dependence]) represents the same T.B.R response.\u003c/p\u003e\n\u003cp\u003eFor all systems, the Moran index was calculated (J. Zhang et al., \u003cspan class=\"CitationRef\"\u003e2014\u003c/span\u003e) and applied to the residual design model once the ANOVA ran in all simulations, thus creating different scenarios of spatial dependence. This index was calculated with the \u003cem\u003eape\u003c/em\u003e library of R (Paradis \u0026amp; Schliep, \u003cspan class=\"CitationRef\"\u003e2019\u003c/span\u003e).\u003c/p\u003e\n\u003cp\u003eIn each simulation, the ANOVA F statistic was extracted, normality and homoscedasticity of residues were tested for treatments. The same O\u003csub\u003eT\u003c/sub\u003e lag statistic was calculated based on the Chi-square distribution (with a degree of freedom) that appears in Darghan et al. (\u003cspan class=\"CitationRef\"\u003e2012\u003c/span\u003e) and that is annexed in this document in a script of R for use, as well as the coefficient associated with the spatial lag (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({{\\alpha }}_{\\text{s}}\\)\u003c/span\u003e\u003c/span\u003e) which was obtained from:\u003c/p\u003e\n\u003cdiv id=\"Equa\" class=\"Equation\"\u003e\n\u003cdiv id=\"FileID_Equa\" class=\"mathdisplay\"\u003e$${ \\widehat{{\\alpha }}}_{\\text{s}}=\\frac{\\mathbf{y}\\mathbf{{\\prime }}{\\mathbf{M}}_{0}\\varvec{W}(\\mathbf{I}-{\\mathbf{M}}_{0})\\mathbf{y}}{\\varvec{y}\\varvec{{\\prime }}{\\mathbf{M}}_{0}\\varvec{W}\\varvec{W}{\\mathbf{M}}_{0}\\varvec{y}} \\left(2\\right)$$\u003c/div\u003e\n\u003c/div\u003e\n\u003cp\u003ewhere \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({\\mathbf{M}}_{0}\\)\u003c/span\u003e\u003c/span\u003e is associated with the design matrix having already repaired the model in the Christensen (\u003cspan class=\"CitationRef\"\u003e2011\u003c/span\u003e) style, with \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({\\mathbf{M}}_{0}={\\varvec{X}}_{0}{\\left({\\varvec{X}}_{0}\\varvec{{\\prime }}{\\varvec{X}}_{0}\\right)}^{-1}{\\varvec{X}}_{0}^{\\varvec{{\\prime }}}\\)\u003c/span\u003e\u003c/span\u003e where \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({\\varvec{X}}_{0}=\\varvec{X}\\varvec{Z}\\)\u003c/span\u003e\u003c/span\u003e being \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\varvec{Z}\\)\u003c/span\u003e\u003c/span\u003e generated with the \"\u003cem\u003emz\u003c/em\u003e\" function of R available in (GitHub1: supplementary material). This matrix depends on the selected design, and without loss of generality, it will be enough to modify the design matrix (in balanced cases) and thus obtain this special matrix for the new design. The equation in (2) also depends on the weight matrix and in this case the inverses of the Euclidean distances were used with type C standardization according to R \u003cem\u003espdep\u003c/em\u003e library (Bivand, \u003cspan class=\"CitationRef\"\u003e2022\u003c/span\u003e). For all spatial variations all these measurements were calculated in addition to the ANOVA (AOV function of R of the library \u003cem\u003estats\u003c/em\u003e (R Core Team, \u003cspan class=\"CitationRef\"\u003e2022\u003c/span\u003e)), where the residuals of the model were extracted to test normality (function Shapiro.test of the library \u003cem\u003estats\u003c/em\u003e of R) and homoscedasticity test in the treatments (Bartlett.test function of R \u003cem\u003estats\u003c/em\u003e). All the information formed a dimension matrix (31000\u0026times;7) that appears available in github2 (supplementary material- 31 ordering scenarios for 1000 simulations and seven responses). With all these data diagrams that were made that related the p value of the Moran index and p value of the ANOVA statistic F with the statistic O\u003csub\u003eT\u003c/sub\u003e, the spatial lag coefficient in each treatment was arranged to demonstrate how this order was related to the Moran index and the lag coefficient, thus highlighting the little or lack of usefulness of the ANOVA in the presence of a dependence or spatial lag.\u003c/p\u003e"},{"header":"Simulation Results","content":"\u003cp\u003eFigure \u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e is one of the most important in this document because it shows the effect of the spatial lag in the distributional behavior of the F statistic. The monotonous continuous line represents the theoretical distribution of the F statistic (central) with 2 and 56 degrees of freedom (for the established design). The dotted vertical line represents the quantile value for a 95% confidence level and the x-axis represents the values of the F obtained in the 1000 simulations. The histogram describes the observed distributional behavior of the F statistic in the presence of spatial lag. It is evident how the theoretical curve F cannot be used to obtain the p value of the ANOVA even though practically all the generated F values do not reject the null hypothesis of the null effect of treatments. Even in the presence of spatial lag, both scenarios lead to no rejection, but the p value, usually generated from the ANOVA is not the measure that represents what is occurring. This would also apply if the rejection of the null hypothesis was even due to a difference in treatment means.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eFigure \u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e shows the probabilities associated with theoretical F and simulated F in the ord0, ord30 and ord60 systems to indicate the spatial dependence according to the Moran index. Already in the last order the spatial lag has been lost and therefore the p value of Moran is greater than 5% in almost all simulations. The dotted line represents the points where both probabilities coincide and happens when the statistic F in both cases tends to zero. It is also evident that it is more likely not to reject the null hypothesis of equality of means in spatial lag scenarios, making it easier to obtain false negatives in the presence of this effect.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eFigure \u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e is like Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e except that the coefficient associated with the spatial lag for the same orders (ord0, ord30 and ord60) is highlighted here. The estimated coefficient of lag is near zero in spatial variations of the greatest occurrence of re-orders.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eFigure \u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e shows the O\u003csub\u003eT\u003c/sub\u003e statistic for spatial lag for the null hypothesis Hoα\u003csub\u003es\u003c/sub\u003e\u0026thinsp;=\u0026thinsp;0. This hypothesis is tested with one degree of freedom. Also shown is a vertical line (for the 95% confidence level) in different systems to delimit in the x-axis the quantum of the theoretical F. The crossing of these lines yields four regions: A) [spatial lag and no treatment effect], B) [spatial lag and no treatment effect], C) [spatial lag and no treatment effect] and D) [spatial lag and no treatment effect]. From orders higher than or equal to 36 (ord36) practically all the data fall in region B of absence of spatial lag and null effect of the treatments. In the ord0, where the data with Moran index indicate spatial dependence, for a fraction of these data there is a spatial lag (B) and in the other there is none (A). To better understand these results, Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e describes the proportion of data in each region.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eDistribution of the proportion of data by region in Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"5\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/th\u003e \u003cth align=\"left\" colspan=\"4\" nameend=\"c5\" namest=\"c2\"\u003e \u003cp\u003eRegion\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eOrder\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u003cb\u003eA\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cb\u003eB\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e\u003cb\u003eC\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u003cb\u003eD\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eOrd0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e49.7\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e50.2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e0.1\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eOrd12\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e37.9\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e62\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e0.1\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eOrd24\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e30.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e69.8\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e0.0\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eOrd36\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e2.7\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e97.2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e0.1\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eOrd48\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e2.8\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e97.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e0.1\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eOrd60\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e1.3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e98.6\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e0.1\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003eFrom Table \u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e it is evident that from the Ord36 most of the data falls into B and at most 10% between regions C and D. This clearly shows as the spatial lag is lost; the spatial dependence is lost even though the ANOVA shows no evidence in the data to reject the null hypothesis of null effect.\u003c/p\u003e \u003cp\u003eFinally, Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003e separates by ordering the behavior of the estimated spatial lag coefficient \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({\\widehat{\\alpha }}_{s}\\)\u003c/span\u003e\u003c/span\u003e where an average at zero and greater variability is clearly perceived when there is greater spatial dependence and thus spatial lag. The clearest points indicate the cases with no spatial dependence when the estimated values of \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({\\widehat{\\alpha }}_{s}\\)\u003c/span\u003e\u003c/span\u003eoscillate approximately between \u0026minus;\u0026thinsp;0.01 and 0.01, maintaining a similar variability since the ord36, when the Moran index already suggests spatial independence.\u003c/p\u003e \u003cp\u003eAlthough ANOVA yields in practice results associated with not rejecting the hypothesis of the null effect of treatments, the p value of the test is not appropriate as it does not reflect the probability of obtaining test results for the F statistic, at least as extreme as the result observed, on the assumption that the null hypothesis is correct.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eSpatial dependence has a direct impact on ANOVA validity. This violation of the assumption of independence by spatial lag (although it can be violated in many ways (Christensen, \u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e1997\u003c/span\u003e) may have several important implications in the ANOVA, highlighting here the underestimation of the p-value of Exhibit F that affects the validity of the associated hypothesis tests and thus the interpretation of results, since the spatial dependence introduces a structure in the data that is not taken into account in the conventional ANOVA that is usually taught in basic courses of experimental design.\u003c/p\u003e \u003cp\u003eAlthough there are different methods for detecting spatial dependency, surely the Moran index or some of its modifications are the most used and are already programmed in several statistical packages. In the case of the lag coefficient, it could also be a convenient measure adapted to each experimental design using the design matrix. As this represents the ratio of two independent quadratic forms, it is possible not only to obtain their moments but even their distributional behavior, making this a measure of alternative interest to judge spatial lag in the context of experimental design, provided that a weight matrix is defined ensuring symmetry in the matrices associated with the quadratic forms.\u003c/p\u003e"},{"header":"Conclusions","content":"\u003cp\u003eSpatial dependence represents a significant violation of the assumption of independence of the residues in the ANOVA. This assumption assumes that observations are independent of each other, that is, that the value of an observation is not related to the value of another nearby observation in space. When spatial dependence exists, close observations tend to be more like each other than more distant observations, meaning that the residues of neighboring observations are correlated and do not meet the assumption of independence. Since the observations taken by remote or proximal sensors of physic-chemical or microbiological properties of soils and vegetation can be grouped in space by evolutionary processes, degradation, crop or soil management or other agents can generate similar predictor variables. This similarity generates a structure in the correlation of observations that manifests as spatial dependence. If this dependence is not considered and in this case is caused by spatial lag in treatments, there is a risk of obtaining incorrect estimates of the p value associated with the ANOVA F test statistic; therefore, this leads to inappropriate inferences and, therefore, to interpretations of little or no value. Spatial autocorrelation can lead to overestimation or underestimation of the significance of treatment effects that can lead to erroneous acceptance or incorrect rejection of associated hypotheses.\u003c/p\u003e \u003cp\u003eAddressing spatial dependence in the analysis of variance by Moran's index or by the lag coefficient for a specified experimental design can help to identify the problem and to take the necessary measures using the appropriate procedures or models for an estimate of some response or an appropriate comparison of treatments.\u003c/p\u003e \u003cp\u003eSpatial dependence represents a critical violation of the independence assumption in the ANOVA and should not be omitted in analyses where studies have proven some form of dependence, as is the case in many of the phenomena that are modeled in the context of geomatics.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eConflict of Interest\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors declare that there is no conflict of interest.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eData Availability\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe data as well as the R scripts to generate the simulations can be found at:\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eGitHub 1\u003c/strong\u003e:\u0026nbsp;https://github.com/CarlosRivera1212/rezago_aov\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eGitHub 2\u003c/strong\u003e:\u0026nbsp;https://github.com/darklley/Rezago\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eAbubakar, B., Uthman, Y. A., Jatau, A. I., Danbatta, A., Nuhu, H. N., \u0026amp; Mustapha, M. (2022). Misuse of analysis of variance in African biomedical journals: a call for more vigilance. \u003cem\u003eBulletin of the National Research Centre 2022 46:1\u003c/em\u003e, \u003cem\u003e46\u003c/em\u003e(1), 1\u0026ndash;8. https://doi.org/10.1186/S42269-022-00924-8\u003c/li\u003e\n\u003cli\u003eAcutis, M., Scaglia, B., \u0026amp; Confalonieri, R. (2012). Perfunctory analysis of variance in agronomy, and its consequences in experimental results interpretation. \u003cem\u003eEuropean Journal of Agronomy\u003c/em\u003e, \u003cem\u003e43\u003c/em\u003e, 129\u0026ndash;135. https://doi.org/10.1016/J.EJA.2012.06.006\u003c/li\u003e\n\u003cli\u003eAdak, S., Bandyopadhyay, K., Sahoo, R., Mridha, N., Shrivastava, M., \u0026amp; Purakayastha, T. (2021). Prediction of Wheat Yield Using Spectral Reflectance Indices Under Different Tillage, Residue and Nitrogen Management Practices. \u003cem\u003eCurrent Science\u003c/em\u003e, \u003cem\u003e121\u003c/em\u003e, 402. https://doi.org/10.18520/cs/v121/i3/402-413\u003c/li\u003e\n\u003cli\u003eAtik, M., \u0026amp; Akdemir, B. (2022). Spatial Variability of NDVI and SPAD for Variable- and Fixed-Rate Fertilization in Olive Orchards. \u003cem\u003eErwerbs-Obstbau\u003c/em\u003e, \u003cem\u003e64\u003c/em\u003e(4), 647\u0026ndash;655. https://doi.org/10.1007/S10341-022-00748-W/METRICS\u003c/li\u003e\n\u003cli\u003eBarbosa, J. de A., de Faria, R. T., Coelho, A. P., Dalri, A. B., \u0026amp; Palaretti, L. F. (2020). Nitrogen fertilization management in white oat using spectral indices. \u003cem\u003ePesquisa Agropecu\u0026aacute;ria Tropical\u003c/em\u003e, \u003cem\u003e50\u003c/em\u003e, 1\u0026ndash;9. https://doi.org/10.1590/1983-40632020V5064924\u003c/li\u003e\n\u003cli\u003eBennington, C. C., \u0026amp; Thayne, W. V. (1994). Use and misuse of mixed model analysis of variance in ecological studies. \u003cem\u003eEcology\u003c/em\u003e, \u003cem\u003e75\u003c/em\u003e(3), 717\u0026ndash;722. https://doi.org/10.2307/1941729\u003c/li\u003e\n\u003cli\u003eBivand, R. (2022). R Packages for Analyzing Spatial Data: A Comparative Case Study with Areal Data. \u003cem\u003eGeographical Analysis\u003c/em\u003e, \u003cem\u003e54\u003c/em\u003e(3), 488\u0026ndash;518. https://doi.org/10.1111/GEAN.12319\u003c/li\u003e\n\u003cli\u003eChristensen, R.,Bedrick,E.J. (1997). Testing the independence assumption in linear models. Journal of the American Statistical Association, 92, 1006\u0026ndash;1016.\u003c/li\u003e\n\u003cli\u003eChristensen, R. (2011). \u003cem\u003ePlane Answers to Complex Questions\u003c/em\u003e. https://doi.org/10.1007/978-1-4419-9816-3\u003c/li\u003e\n\u003cli\u003eChristensen, R. (2019). \u003cem\u003eAdvanced Linear Modeling\u003c/em\u003e. https://doi.org/10.1007/978-3-030-29164-8\u003c/li\u003e\n\u003cli\u003eDarghan, A., Surendra, P., \u0026amp; Goit\u0026iacute;a, A. (2012). An application of the overlapping test in a two-way classification reparameterized model. \u003cem\u003eRevista de La Facultad de Agronomia\u003c/em\u003e, \u003cem\u003e29\u003c/em\u003e, 124\u0026ndash;137.\u003c/li\u003e\n\u003cli\u003eDing, Y., Jiang, Y., Yu, H., Yang, C., Wu, X., Sun, G., Fu, X., \u0026amp; Dou, X. (2020). Measurement Method for Height-Independent Vegetation Indices Based on an Active Light Source. \u003cem\u003eSensors 2020, Vol. 20, Page 1830\u003c/em\u003e, \u003cem\u003e20\u003c/em\u003e(7), 1830. https://doi.org/10.3390/S20071830\u003c/li\u003e\n\u003cli\u003edos Santos, R. A., Filgueiras, R., Mantovani, E. C., Fernandes-Filho, E. I., Almeida, T. S., Venancio, L. P., \u0026amp; da Silva, A. C. B. (2021). Surface reflectance calculation and predictive models of biophysical parameters of maize crop from RG-NIR sensor on board a UAV. \u003cem\u003ePrecision Agriculture\u003c/em\u003e, \u003cem\u003e22\u003c/em\u003e(5), 1535\u0026ndash;1558. https://doi.org/10.1007/S11119-021-09795-X/METRICS\u003c/li\u003e\n\u003cli\u003eFirozjaei, M. K., Fathololoumi, S., Mijani, N., Kiavarz, M., Qureshi, S., Homaee, M., \u0026amp; Alavipanah, S. K. (2020). Evaluating the Spectral Indices Efficiency to Quantify Daytime Surface Anthropogenic Heat Island Intensity: An Intercontinental Methodology. \u003cem\u003eRemote Sensing 2020, Vol. 12, Page 2854\u003c/em\u003e, \u003cem\u003e12\u003c/em\u003e(17), 2854. https://doi.org/10.3390/RS12172854\u003c/li\u003e\n\u003cli\u003eGao, J., Li, X., and Brierley, G. (2012). Topographic influence on wetland distribution and change in Maduo County, Qinghai-Tibet Plateau, China. Journal of Mountain Science, 9(3):362{371.\u003c/li\u003e\n\u003cli\u003eGotway, C. A., \u0026amp; Cressie, N. A. C. (1990). A Spatial Analysis of Variance Applied to Soil-Water Infiltration. \u003cem\u003eWater Resources Research\u003c/em\u003e, \u003cem\u003e26\u003c/em\u003e(11). https://doi.org/10.1029/WR026i011p02695\u003c/li\u003e\n\u003cli\u003eGr\u0026auml;ler,B.,Pebesma,E.,Heuvelink,G.(2016). Spatio-Temporal interpolation using gstat. The R Journal, 8(1), pp. 204-218.\u003c/li\u003e\n\u003cli\u003eKenny, D. A. (1995). The effect of nonindependence on significance testing in dyadic research. \u003cem\u003ePersonal Relationships\u003c/em\u003e, \u003cem\u003e2\u003c/em\u003e(1), 67\u0026ndash;75. https://doi.org/10.1111/J.1475-6811.1995.TB00078.X\u003c/li\u003e\n\u003cli\u003eKenny, D. A., \u0026amp; Judd, C. M. (1986). Consequences of Violating the Independence Assumption in Analysis of Variance. \u003cem\u003ePsychological Bulletin\u003c/em\u003e, \u003cem\u003e99\u003c/em\u003e(3), 422\u0026ndash;431. https://doi.org/10.1037/0033-2909.99.3.422\u003c/li\u003e\n\u003cli\u003eLi, J., \u0026amp; Lomax, R. G. (2011). Analysis of Variance: What Is Your Statistical Software Actually Doing? \u003cem\u003eHttp://Dx.Doi.Org/10.1080/00220973.2010.481684\u003c/em\u003e, \u003cem\u003e79\u003c/em\u003e(3), 279\u0026ndash;294. https://doi.org/10.1080/00220973.2010.481684\u003c/li\u003e\n\u003cli\u003eLindman, H. R. (1992). \u003cem\u003eAnalysis of Variance in Experimental Design\u003c/em\u003e. https://doi.org/10.1007/978-1-4613-9722-9\u003c/li\u003e\n\u003cli\u003eMcintosh, M. S. (2015). Can Analysis of Variance Be More Significant? \u003cem\u003eAgronomy Journal\u003c/em\u003e, \u003cem\u003e107\u003c/em\u003e(2), 706\u0026ndash;717. https://doi.org/10.2134/AGRONJ14.0177\u003c/li\u003e\n\u003cli\u003eParadis, E., \u0026amp; Schliep, K. (2019). ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R. \u003cem\u003eBioinformatics\u003c/em\u003e, \u003cem\u003e35\u003c/em\u003e(3), 526\u0026ndash;528. https://doi.org/10.1093/BIOINFORMATICS/BTY633\u003c/li\u003e\n\u003cli\u003ePearce, S. C. (2006). Defective data in the analysis of variance. \u003cem\u003eHttp://Dx.Doi.Org/10.1080/02664768600000022\u003c/em\u003e, \u003cem\u003e13\u003c/em\u003e(2), 139\u0026ndash;147. https://doi.org/10.1080/02664768600000022\u003c/li\u003e\n\u003cli\u003eR Core Team. (2022). \u003cem\u003eR: A Language and Environment for Statistical Computing\u003c/em\u003e. https://www.R-project.org/\u003c/li\u003e\n\u003cli\u003eRong, Y. (2000). Statistical Methods and Pitfalls in Environmental Data Analysis. \u003cem\u003eEnvironmental Forensics\u003c/em\u003e, \u003cem\u003e1\u003c/em\u003e(4), 213\u0026ndash;220. https://doi.org/10.1006/ENFO.2000.0022\u003c/li\u003e\n\u003cli\u003eRowell, J. G., \u0026amp; Walters, D. E. (2006). A critical appraisal of the analysis of variance for two-way classifications. \u003cem\u003eHttp://Dx.Doi.Org/10.1080/02664768500000021\u003c/em\u003e, \u003cem\u003e12\u003c/em\u003e(2), 157\u0026ndash;169. https://doi.org/10.1080/02664768500000021\u003c/li\u003e\n\u003cli\u003eShukl, G. K., \u0026amp; Subrahmanyam, G. S. V. (1999). A note on an exact test and confidence interval for competition and overlap effects. \u003cem\u003eBiometrics\u003c/em\u003e, \u003cem\u003e55\u003c/em\u003e(1), 273\u0026ndash;276. https://doi.org/10.1111/J.0006-341X.1999.00273.X\u003c/li\u003e\n\u003cli\u003eStoy, P. C., Khan, A. M., Wipf, A., Silverman, N., \u0026amp; Powell, S. L. (2022). The spatial variability of NDVI within a wheat field: Information content and implications for yield and grain protein monitoring. \u003cem\u003ePLOS ONE\u003c/em\u003e, \u003cem\u003e17\u003c/em\u003e(3), e0265243. https://doi.org/10.1371/JOURNAL.PONE.0265243\u003c/li\u003e\n\u003cli\u003eVolcani, A., Karnieli, A., \u0026amp; Svoray, T. (2005). The use of remote sensing and GIS for spatio-temporal analysis of the physiological state of a semi-arid forest with respect to drought years. \u003cem\u003eForest Ecology and Management\u003c/em\u003e, \u003cem\u003e215\u003c/em\u003e(1\u0026ndash;3), 239\u0026ndash;250. https://doi.org/10.1016/J.FORECO.2005.05.063\u003c/li\u003e\n\u003cli\u003eWang, C., Nie, S., Xi, X., Luo, S., Sun, X., Zhang, J., Lin, X., Baghdadi, N., Gloaguen, R., \u0026amp; Thenkabail, P. S. (2016). Estimating the Biomass of Maize with Hyperspectral and LiDAR Data. \u003cem\u003eRemote Sensing 2017, Vol. 9, Page 11\u003c/em\u003e, \u003cem\u003e9\u003c/em\u003e(1), 11. https://doi.org/10.3390/RS9010011\u003c/li\u003e\n\u003cli\u003eYa\u0026rsquo;acob, N., Ismail, N. S., Mustafa, N., \u0026amp; Yusof, A. L. (2014). Investigation of flood pattern using ANOVA statistic and remote sensing in Malaysia. \u003cem\u003eIOP Conference Series: Earth and Environmental Science\u003c/em\u003e, \u003cem\u003e20\u003c/em\u003e(1), 012030. https://doi.org/10.1088/1755-1315/20/1/012030\u003c/li\u003e\n\u003cli\u003eZhang, H., Lan, Y., Lacey, R., Hoffmann, W. C., \u0026amp; Westbrook, J. K. (2011). Spatial Analysis of NDVI Readings with Different Sampling Densities. \u003cem\u003eTransactions of the ASABE\u003c/em\u003e, \u003cem\u003e54\u003c/em\u003e(1), 349\u0026ndash;354. https://doi.org/10.13031/2013.36244\u003c/li\u003e\n\u003cli\u003eZhang, J., Atkinson, P., \u0026amp; Goodchild, M. F. (2014). Scale in spatial information and analysis. \u003cem\u003eScale in Spatial Information and Analysis\u003c/em\u003e, 1\u0026ndash;347.\u003c/li\u003e\n\u003cli\u003eZimmerman, D. W., Williams, R. H., \u0026amp; Zurabo, B. D. (2007). Effect of nonindependence of sample observations on some parametric and nonparametric statistical tests. \u003cem\u003eHttp://Dx.Doi.Org/10.1080/03610919308813123\u003c/em\u003e, \u003cem\u003e22\u003c/em\u003e(3), 779\u0026ndash;789. https://doi.org/10.1080/03610919308813123\u003c/li\u003e\n\u003cli\u003eZimmerman, D. W., \u0026amp; Zumbo, B. D. (2010). Correction for Nonindependence of Sample Observations in ANOVA F Tests. \u003cem\u003eHttp://Dx.Doi.Org/10.1080/00220973.1992.9943874\u003c/em\u003e, \u003cem\u003e60\u003c/em\u003e(4), 367\u0026ndash;381. https://doi.org/10.1080/00220973.1992.9943874\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":true,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"applied-geomatics","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"agmj","sideBox":"Learn more about [Applied Geomatics](http://link.springer.com/journal/12518)","snPcode":"12518","submissionUrl":"https://submission.nature.com/new-submission/12518/3","title":"Applied Geomatics","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false},"keywords":"spatial lag, spatial overlapping, spectral index, sensing","lastPublishedDoi":"10.21203/rs.3.rs-3243407/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-3243407/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eIn recent years, statistical methods involving spatial considerations have been developed, for example, those incorporating data with some type of georeferencing. The descriptive part of geographic information systems currently provides many visualization and analytic tools; however, the latter is still quite limited. In this sense, research of a spatial nature is seen as combining non-spatial statistical methods for inferential treatment that can certainly invalidate the excellent capture work with advanced tools such as those observed every day in the geomatic context. This prompted the current document, drawing attention to how geomatic information analyzed with statistical methods that imply independence in modeled observations can be invalid. The Moran index is compared with a proposal for a spatial lag coefficient in the context of experimental design so that users of variance analysis do not apply this well-known procedure in a ritualistic way, perhaps revising some assumptions and perhaps ignoring more important ones. The distortion of the p value generated from the analysis of variance is clear in the presence of spatial dependence. In this case it is associated with the lag or spatial overlap. The methodology is simple to adopt in other experimental designs with the simple consideration of the design matrix and its reparameterization and the choice of the appropriate weight matrix. This will allow users to reconsider the traditional method of analysis and incorporate some methodology to support spatial dependency structures.\u003c/p\u003e","manuscriptTitle":"The effect of spatial lag on modeling geomatic covariates using analysis of variance","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2023-08-30 16:12:13","doi":"10.21203/rs.3.rs-3243407/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"reviewerAgreed","content":"","date":"2023-09-01T09:55:34+00:00","index":0,"fulltext":""},{"type":"reviewersInvited","content":"","date":"2023-08-25T09:20:46+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2023-08-11T05:44:38+00:00","index":"","fulltext":""},{"type":"submitted","content":"Applied Geomatics","date":"2023-08-09T17:21:17+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"applied-geomatics","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"agmj","sideBox":"Learn more about [Applied Geomatics](http://link.springer.com/journal/12518)","snPcode":"12518","submissionUrl":"https://submission.nature.com/new-submission/12518/3","title":"Applied Geomatics","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false}}],"origin":"","ownerIdentity":"6a1c119d-f93f-4feb-81d5-e2e27ca23765","owner":[],"postedDate":"August 30th, 2023","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[],"tags":[],"updatedAt":"2024-08-01T17:08:27+00:00","versionOfRecord":{"articleIdentity":"rs-3243407","link":"https://doi.org/10.1007/s12518-024-00579-2","journal":{"identity":"applied-geomatics","isVorOnly":false,"title":"Applied Geomatics"},"publishedOn":"2024-07-22 16:16:15","publishedOnDateReadable":"July 22nd, 2024"},"versionCreatedAt":"2023-08-30 16:12:13","video":"","vorDoi":"10.1007/s12518-024-00579-2","vorDoiUrl":"https://doi.org/10.1007/s12518-024-00579-2","workflowStages":[]},"version":"v1","identity":"rs-3243407","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-3243407","identity":"rs-3243407","version":["v1"]},"buildId":"J0_U0BvcaRcwD8yVFaRlm","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}