plotData: Plot Geostatistical Data

In geoCount: Analysis and Modeling for Geostatistical Count Data

Description

This function plots geostatistical data for up to three data sets.

Usage

plotData(Y=NULL, loc=NULL, Yp=NULL, locp=NULL, Yt=NULL, loct=NULL, bdry=NULL, 
         cols=1:3, pchs=1:3, size=c(0.3, 2.7), ...)

Arguments

Y, Yp, Yt	the vector of response variables
loc, locp, loct	n \times 2 matrix that indicates the coordinates of locations
bdry	a list containing the coordinates of boundaries
cols	the colors used for different sets of response variables
pchs	the shapes used for different sets of response variables
size	the minimum and maximum of the sizes
...	other parameters that control the plotting

Author(s)

Liang Jing [e-mail]

See Also

locCircle, locGrid, locSquad, simData.

Examples

## Not run: 
loc1 <- locGrid(1, 1, 10, 10)
loc2 <- locCircle(0.72, 60)
loc3 <- locSquad(0.38, 10)
loc <- rbind(as.matrix(loc1) , loc2, loc3); plot(loc)
dat <- simData(loc, cov.par = c(1, 0.2, 1))
Y <- dat$data
plotData(Y[1:nrow(loc1)], loc1, 
         Y[(nrow(loc1)+1):(nrow(loc1)+nrow(loc2))], loc2, 
         Y[(length(Y)-nrow(loc3)+1):length(Y)], loc3, 
         xlab="x", ylab="y", pchs = c(1, 16, 15)
         )
# plot boundaries
data(TexasCounty_boundary)
plotData(bdry = TexasCounty.boundary)
# plot data with the boundary
data(Rongelap)
str(Rongelap)
plotData(bdry = Rongelap$borders, Y = Rongelap$data, loc = Rongelap$coords)

## End(Not run)

geoCount documentation built on May 2, 2019, 12:38 a.m.