nmf.gmm.inference: Statistical inference for an NMF-GMM fit (given basis)

In nmfkc: Non-Negative Matrix Factorization with Kernel Covariates

Statistical inference for an NMF-GMM fit (given basis)

Description

This function is experimental and still under development. The interface may change in future versions: argument names, defaults and the contents of the returned object are not yet stable. Code written against it today may need adjusting after an update. The rest of the package does not carry this caveat.

nmf.gmm.inference adds Wald inference on the covariate-coefficient matrix C (=\Theta) of a fitted nmf.gmm object, conditional on the estimated basis \hat X and mixture. It reports the analytic outer-product (mixture-information) standard error — the coverage-calibrated choice at K>1 — together with a wild-bootstrap standard error and percentile confidence interval, mirroring nmfre.inference / nmfkc.inference. Supported for the tied-covariance variant.

Usage

nmf.gmm.inference(object, Y, A = object$A, ...)

Arguments

object	A fitted "nmf.gmm" object (with cov = "tied").
Y	The data matrix used in the fit.
A	The covariate matrix used in the fit. Default object$A.
...	Additional arguments, named as in the other wild-bootstrap inference functions of the package: wild.B (bootstrap replicates, default 500), wild.level (confidence level, default 0.95) and seed (default 123).

Value

object with class c("nmf.gmm.inference", "nmf.gmm") and added components: coefficients (data frame with Basis, Covariate, Estimate, SE, BSE, z_value, p_value, CI_low, CI_high), C.se.outer, C.se.sandwich, C.se.boot (Q x R matrices) and Xi (= X C).

See Also

nmf.gmm, nmf.gmm.select

Examples

set.seed(1)
Y <- matrix(abs(rnorm(20 * 60)) + 1, 20, 60); A <- rbind(1, rnorm(60))
fit <- nmf.gmm(Y, A, rank = 2, K = 2)
fit <- nmf.gmm.inference(fit, Y, A, wild.B = 300)
head(fit$coefficients)

nmfkc documentation built on Oct. 30, 2026, 9:06 a.m.