gettingtothebottom (version 3.2)

plot_nnm_reconstruction: MM Algorithm - Plotting the Reconstruction

Description

plot_nnm_reconstruction Function for plotting the nnm_sol reconstruction

Usage

plot_nnm_reconstruction(nnm, X, nnm_sol)

Arguments

nnm

NNM object from generate_nnm function

X

Nonnegative design matrix

nnm_sol

Solution object from nnls_mm function

Examples

# Setup mixture example
n <- 1e3
p <- 10
nnm <- generate_nnm(n,p)

set.seed(12345)
X <- nnm$X
b <- double(p)
nComponents <- 3
k <- sample(1:p,nComponents,replace=FALSE)
b[k] <- matrix(runif(nComponents),ncol=1)
y <- X%*%b + 0.25*matrix(abs(rnorm(n)),n,1)

# Obtain solution to mixture problem
nnm_sol <- nnls_mm(y,X,runif(p))

# Plot the reconstruction
plot_nnm_reconstruction(nnm,X,nnm_sol)