gmGeostats (version 0.11.3)

anaForward: Forward gaussian anamorphosis forward transformation to multivariate gaussian scores

Description

Forward gaussian anamorphosis forward transformation to multivariate gaussian scores

Usage

anaForward(
  x,
  Y,
  sigma0,
  sigma1 = 1 + sigma0,
  steps = 30,
  plt = FALSE,
  sphere = TRUE,
  weights = NULL
)

Value

a matrix with the gaussian scores; same dimensions of x

Arguments

x

points to be transformed (a matrix)

Y

node points defining the transformation (another matrix, same nr. of columns as x)

sigma0

starting spread of the kernels

sigma1

final spread of the kernels

steps

number of steps to linearize the transform (default 30 is good)

plt

boolean, do you want to get a plot of the transformation?

sphere

boolean, should the data be pre-Y-spherified first? defaults to true

weights

vector of weights for all computations, length must be equal to number of rows of x

Author

K. Gerald van den Boogaart, Raimon Tolosana-Delgado

See Also

ana() for defining a function that carries over the transformation (by means of a closure), anaBackward() for the explicit back-transformation, sphTrans() for defining a function that carries over the spherification of the data

Examples

data("jura", package="gstat")
Y = jura.pred[,c(10,12,13)]
plot(compositions::acomp(Y))
Ylr = compositions::alr(Y)
plot(Ylr)
z = anaForward(x=Ylr, Y=Ylr, sigma0=0.1)
plot(z, asp=1)
shapiro.test(z[,1])
shapiro.test(z[,2])