chi3: Theoretical triple-wise chi (chi_q) for a (Inverted)...

In Jbrich95/sdfExtreme: Spatial Deformation for Non-Stationary Extremal Dependence

Theoretical triple-wise χ (χ_q) for a (Inverted) Brown-Resnick process

Description

Calculates the theoretical triple-wise χ(s_i,s_j,s_k) or χ_q(s_i,s_j,s_k) for a given 3 by 3 matrix of semivariogram. values.

Usage

chi3(v_H)

chi3q(v_H, q)

Arguments

v_H	3 by 3 matrix of semivariogram. values.
q	Exceedance threshold for χ_q(s_i,s_j,s_k).

Value

Theoretical χ(s_i,s_j,s_k) or χ_q(s_i,s_j,s_k) measure for the corresponding matrix of semivariogram values.

Examples

data(Aus_Heat)
data(Aus_Heat_Output)
Z<-Aus_Heat$Temp.
Gcoords<-Aus_Heat$coords
sdf<-Aus_Heat_Output$sdf
likD.MSP<-Aus_Heat_Output$likD.MSP
likD.IMSP<-Aus_Heat_Output$likD.IMSP

Dcoords<-returnDcoord(sdf$par,Gcoords,sdf$m.ind,sphere.dis=TRUE)

 ind.triple<-c(1,2,3) ##Denotes indices of triple for which triple-wise chi calculated.

 #MSP
 v_H<-(rdist.earth(Dcoords[ind.triple,],miles=F)/likD.MSP$par[2])^likD.MSP$par[1]
 print(chi3(v_H))

  #IMSP - likD.IMSP only contains range parameter as smoothing parameter is 2.
 v_H<-(rdist.earth(Dcoords[ind.triple,],miles=F)/likD.IMSP$par[1])^2
 print(chi3q(v_H,q=0.95))

Jbrich95/sdfExtreme documentation built on March 24, 2022, 11:15 a.m.