Se a resposta for gaussiana, use a rmvtnorm() par gerar resíduos dados correlacionados. Para isso crie um arquivo com as coordenadas, como é um quadrado látino, no espaço elas formarão um retângulo regular, que pode ser obtido com a expand.grid(). Obtenha a matriz de distâncias com a dist(). Escolha um modelo de dependência espacial (esponencial, gaussiano, esférico). Gere a matriz de covariância, pode usar o comando outer() para gera a matriz. Passe para rmvtnorm() como argumento Sigma=. Crie o vetor de parâmetros dos efeitos dos tratamentos, crie a matriz de delinamento com model.matriz(). Obtenha os valores esperados. Some aos resíduos. Os dados foram gerados.