gettingtothebottom (version 3.2)

plot_spect: MM Algorithm - Plotting the Spectroscopic Signal

Description

plot_spect Function for plotting the spectroscopic signal

Usage

plot_spect(n, y, X, b, nnm)

Arguments

n

Number of samples

y

Nonnegative response

X

Nonnegative design matrix

b

Nonnegative initial regression vector

nnm

NNM object from generate_nnm 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)

plot_spect(n,y,X,b,nnm)