sommer (version 3.2)

plot.variogram.MMERM: Default variogram.MMERM plotting

Description

Takes a fitted variogram.MMERM object produced by variogram.MMERM() and plots the associated sample variogram using an wireframe plot (package lattice).

Usage

# S3 method for variogram.MMERM
plot(x, stnd=TRUE, ...)

Arguments

x

an object of class variogram.MMERM as produced by variogram.MMERM().

stnd

argument for ploting the residuals to know if they should be standarized.

…

further arguments passed to or from other methods. Not yet implemented.

Details

This function as well as function variogram.MMERM() can only be used for regular two dimensional data.

References

Gilmour, A.R., Cullis, B.R., and Verbyla, A.P. (1997). Accounting for Natural and Extraneous Variation in the Analysis of Field Experiments. Journal of Agricultural, Biological, and Environmental Statistics, 2, 269 - 293.

Stefanova, K.T., Smith, A.B., and Cullis, B.R. (2009). Enhanced Diagnostics for the Spatial Analysis of Field Trials. Journal of Agricultural, Biological, and Environmental Statistics, 14, 392 - 410.

See Also

mmer2,mmer, variogram.MMERM

Examples

# NOT RUN {
data(CPdata)
#### look at the data and fit the model
head(CPpheno)
mix1 <- mmer2(Yield~1,
              random=~ id
                      + Rowf + Colf,
              rcov=~units,
              data=CPpheno)
summary(mix1)

# compute the variogram
vm0 <- variogram(mix1, xcoor = "Row", ycoor = "Col")
# plot the variogram
plot(vm0$F1)
# }