genscore (version 1.0.2.2)

get_elts_exp: The R implementation to get the elements necessary for calculations for the exponential square-root setting (a=0.5, b=0.5).

Description

The R implementation to get the elements necessary for calculations for the exponential square-root setting (a=0.5, b=0.5).

Usage

get_elts_exp(
  hdx,
  hpdx,
  x,
  centered = TRUE,
  profiled_if_noncenter = TRUE,
  scale = "",
  diagonal_multiplier = 1
)

Value

A list that contains the elements necessary for estimation.

n

The sample size.

p

The dimension.

centered

The centered setting or not. Same as input.

scale

The scaling method. Same as input.

diagonal_multiplier

The diagonal multiplier. Same as input.

diagonals_with_multiplier

A vector that contains the diagonal entries of $\mathbf{\Gamma }$ after applying the multiplier.

setting

The setting "exp".

g_K

The $\mathit{g}$ vector. In the non-profiled non-centered setting, this is the $\mathit{g}$ sub-vector corresponding to $\mathbf{K}$.

Gamma_K

The $\mathbf{\Gamma }$ matrix with no diagonal multiplier. In the non-profiled non-centered setting, this is the $\mathbf{\Gamma }$ sub-matrix corresponding to $\mathbf{K}$.

g_eta

Returned in the non-profiled non-centered setting. The $\mathit{g}$ sub-vector corresponding to $\mathit{\eta }$.

Gamma_K_eta

Returned in the non-profiled non-centered setting. The $\mathbf{\Gamma }$ sub-matrix corresponding to interaction between $\mathbf{K}$ and $\mathit{\eta }$.

Gamma_eta

Returned in the non-profiled non-centered setting. The $\mathbf{\Gamma }$ sub-matrix corresponding to $\mathit{\eta }$.

t1,t2

Returned in the profiled non-centered setting, where the $\mathit{\eta }$ estimate can be retrieved from ${\mathit{t}}_{\mathbf{1}}-{\mathit{t}}_{\mathbf{2}}\hat{\mathbf{K}}$ after appropriate resizing.

Arguments

hdx

A matrix, $h(\mathbf{x})$ applied to the distance of x from the boundary of the domain, should be of the same dimension as x.

hpdx

A matrix, $h^{\prime }(\mathbf{x})$ applied to the distance of x from the boundary of the domain, should be of the same dimension as x.

x

An n by p matrix, the data matrix, where n is the sample size and p the dimension.

centered

A boolean, whether in the centered setting (assume $\mathit{\mu }=\mathit{\eta }=0$) or not. Default to TRUE.

profiled_if_noncenter

A boolean, whether in the profiled setting (${\lambda }_{\mathit{\eta }}=0$) if non-centered. Parameter ignored if centered=TRUE. Default to TRUE.

scale

A string indicating the scaling method. Returned without being checked or used in the function body. Default to "norm".

diagonal_multiplier

A number >= 1, the diagonal multiplier.

Details

For details on the returned values, please refer to get_elts_ab or get_elts.

Examples

n <- 50
p <- 30
eta <- rep(0, p)
K <- diag(p)
domain <- make_domain("R+", p=p)
x <- gen(n, setting="exp", abs=FALSE, eta=eta, K=K, domain=domain, finite_infinity=100,
       xinit=NULL, seed=2, burn_in=1000, thinning=100, verbose=FALSE)
h_hp <- get_h_hp("min_pow", 1, 3)
h_hp_dx <- h_of_dist(h_hp, x, domain) # h and h' applied to distance from x to boundary
elts <- get_elts_exp(h_hp_dx$hdx, h_hp_dx$hpdx, x, centered=TRUE, scale="norm", diag=1.5)
elts <- get_elts_exp(h_hp_dx$hdx, h_hp_dx$hpdx, x, centered=FALSE, profiled_if_noncenter=TRUE,
      scale="norm", diag=1.7)
elts <- get_elts_exp(h_hp_dx$hdx, h_hp_dx$hpdx, x, centered=FALSE, profiled_if_noncenter=FALSE,
      scale="norm", diag=1.7)