gdata: Dataset for RSA

In RSAgeo: Resampling-Based Analysis of Geostatistical Data

Description

Simulated date matrix with response variable distributed with 'powered exponential' covariance matrix, with φ=25, τ=1, κ=1, σ=1.

Usage

data(gdata)

Format

A sample dataset with 1000 rows and 4 columns. Each row represents one observation. The first two column of the matrix consists of x1 and x2, where x1 and x2 provide the 2D location for each obsevation. The third column of the matrix gives the y value, which is the response value. Covariate is included in Column 4. This is the default form for the variable Data in function RSA.

Details

####code used to generate gdata####
require("geoR")
require("RandomFields")
DataNum=1000;
gData=grf(DataNum,grid="irreg",DataNum,DataNum,xlims=c(0,100),ylims=c(0,100),nsim=1,mean=0, cov.mode="powered.exponential",cov.par=c(1,25),nugget=1,kappa=1)
l=gData$coords
y=gData$data
x=rnorm(DataNum)
y2=y+.5+x
gdata=cbind(l,y2,x)

Examples

data(gdata)
## This is an illustration of the function. 
## In practice, Total_Iteration is recommended to be set to 2500 or higher.
## N_subset is recommemded to be set to 300 or higher.
RSA(gdata,N_subset=50,Stepscale=40,Total_Iteration=100,Warm=20)

RSAgeo documentation built on May 30, 2017, 8:10 a.m.