plotData {geoCount}R Documentation

Plot Geostatistical 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 ljing918@gmail.com

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)

[Package geoCount version 1.150120 Index]