plot_nnm_reconstruction: MM Algorithm - Plotting the Reconstruction

In gettingtothebottom: Learning Optimization and Machine Learning for Statistics

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)

gettingtothebottom documentation built on May 29, 2017, 8:28 p.m.