FeedbackTS (version 1.2)

krige.test: Test of spatial variation in feedback and change-in-feedback

Description

Randomization test to investigate spatial variation in a kriged index.

Usage

krige.test(krige.output, subregion, alternative, nb.rand, subregion.coverage=0.8)

Arguments

krige.output

a list of items corresponding to the output of the krige function.

subregion

a list of two vectors containing the lat/long coordinates of the vertices of a polygon. The polygon defines a subregion where one supposes variation in the predicted index krige.output$krige$predict. [object Object],[object Object]

Value

• an object from the KT.output class

item

• alternative
• nb.rand
• subregion.coverage

Details

The criterion to accept a random translation is a minimum coverage of subregion by the translated grid: the ratio between the number of nodes of the translated grid that are contained in subregion and the number of nodes of the original grid that are contained in subregion must be equal to or greater than subregion.coverage.

References

Soubeyrand, S., Morris, C. E. and Bigg, E. K. (in press). Analysis of fragmented time directionality in time series to elucidate feedbacks in climate data. Environmental Modelling and Software.

See Also

krige, KT.output-class, rain.feedback.stats

Examples

#### load data of feedback and change-in-feedback indices in 88 sites across Australia
data(rain.feedback.stats)

#### spatial coordinates of the 88 sites
coord=rain.feedback.stats[,3:4]

#### map of feedback index computed from the whole data series
stat1=rain.feedback.stats[["Feedback.whole.period"]]

#### variogram analysis and kriging of feedback index
par(mfrow=c(2,3), mar=c(5.1,4.1,4.1,2.1))
kr1=krige(coordinates=coord,statistic=stat1,
   grid=list(x=seq(110,155,0.25),y=seq(-45,-11,0.25),border="Australia",
      proj="+proj=lcc +lat_1=-18 +lat_2=-36 +lat0=-25 +lon_0=140",degrees=TRUE),
   variog.param=list(keep.distance=2/3,nb.bin=15,cov.pars=c(0.0008,10^6),
      nugget=0.0001,fix.nugget=FALSE,coordinates.scaling=10^5,statistic.scaling=0.01),
   krige.param=list(type.krige="ok",trend.d="1st"),plots=TRUE)

#### test spatial variation in feedback index
## computer intensive stage (un-comment the following command lines)
kt1=krige.test(krige.output=kr1,subregion=list(x=c(138,152,152,138),y=-c(40,40,33,33)),
   alternative="greater", nb.rand=2000)

## plot test output	
par(mfrow=c(1,2), mar=c(5.1,4.1,4.1,2.1))	
plot(kt1,digits=list(predict=3,pvalue=3),breaks=12)