spmodel <- list(model = 'blocks',
coord = globulus[, c('x','y')],
id = 'bl')
fm1 <- remlf90(fixed = phe_X ~ gg,
spatial = spmodel,
data = globulus)
#> No specification of initial variances.
#> Using default value of 1 for all variance components.
#> See ?breedR.getOption.
plot(fm1, "spatial")
#> Error in x$call$spatial$model: object of type 'symbol' is not subsettable