log_lik: Evaluate the negative log-likelihood of an intercept-only...

In plmmr: Penalized Linear Mixed Models for Correlated Data

Evaluate the negative log-likelihood of an intercept-only Gaussian plmm model

Description

This function allows you to evaluate the negative log-likelihood of a linear mixed model under the assumption of a null model in order to estimate the variance parameter, eta.

Usage

log_lik(eta, n, s, U, y, rot_y = NULL)

Arguments

eta	The proportion of variance in the outcome that is attributable to causal SNP effects. In other words, signal-to-noise ratio.
n	The number of observations
s	The singular values of K, the realized relationship matrix
U	The left-singular vectors of the standardized design matrix
y	Continuous outcome vector.
rot_y	Optional: if y has already been rotated, then this can be supplied.

Value

the value of the log-likelihood of the PLMM, evaluated with the supplied parameters

plmmr documentation built on April 4, 2025, 12:19 a.m.