loglik_kappa_est: Kappa log-likelihood

View source: R/utils.R

loglik_kappa_estR Documentation

Kappa log-likelihood

Description

Applicable to a single latent field, or multiple latent fields if common smoothness is assumed

Usage

loglik_kappa_est(par, delta, D_diag, mesh, C1 = 1/(4 * pi), Q = NULL)

Arguments

par

Vector of length two containing values of log kappa and log residual variance at which to compute log likelihood

delta

Estimate of delta (subject effect or deviation)

D_diag

Diagonal values of D matrix (template standard deviations)

mesh

Object of class "templateICA_mesh" containing the triangular mesh (see make_mesh)

C1

For the unit variance case, τ^2 = C1/κ^2, where C1 = 1/(4π) when α=2, ν=1, d=2

Q

Equal to the number of ICs for the common smoothness model, or NULL for the IC-specific smoothness model

Details

Compute log-likelihood of kappa given an initial estimate of delta

Value

Value of negative log likelihood


templateICAr documentation built on Feb. 16, 2023, 8:14 p.m.