samEL: Training function of Sparse Additive Poisson Regression

In SAM: Sparse Additive Modelling

Training function of Sparse Additive Poisson Regression

Description

Fit a sparse additive Poisson regression model on training data.

Usage

samEL(
  X,
  y,
  p = 3,
  lambda = NULL,
  nlambda = NULL,
  lambda.min.ratio = 0.25,
  thol = 1e-05,
  max.ite = 1e+05,
  regfunc = "L1",
  dfmax = NULL,
  verbose = FALSE,
  dev.ratio.thr = NULL,
  dev.change.thr = NULL
)

Arguments

X	Numeric training matrix with n rows (samples) and d columns (features).
y	Response vector of length n; values must be non-negative integers.
p	The number of basis spline functions. The default value is 3.
lambda	Optional user-supplied regularization sequence. If provided, use a decreasing sequence; warm starts are used along the path and are usually much faster than fitting a single value.
nlambda	The number of lambda values. The default value is 20.
lambda.min.ratio	Smallest lambda as a fraction of lambda.max (the smallest value that keeps all component functions at zero). The default is 0.25.
thol	Stopping tolerance. The default value is 1e-5.
max.ite	Maximum number of iterations. The default value is 1e5.
regfunc	A string indicating the regularizer. The default value is "L1". You can also assign "MCP" or "SCAD" to it.
dfmax	Maximum number of non-zero groups allowed. When the number of non-zero groups reaches dfmax, the regularization path is terminated early. NULL (default) means no limit.
verbose	Logical; if TRUE, print iteration info for each lambda.
dev.ratio.thr	Deviance ratio threshold for early stopping. When the deviance ratio 1 - D(\lambda)/D_0 exceeds this value the regularization path is terminated early. NULL (default) disables this criterion.
dev.change.thr	Relative deviance change threshold for early stopping. When the relative change in deviance over the last few lambda steps falls below this value, the path is terminated. NULL (default) disables this criterion.

Details

The solver combines block coordinate descent, fast iterative soft-thresholding, and Newton updates. Computation is accelerated by warm starts and active-set screening.

Value

p	The number of basis spline functions used in training.
X.min	Per-feature minimums from training data (used to rescale test data).
X.ran	Per-feature ranges from training data (used to rescale test data).
lambda	Sequence of regularization parameters used in training.
w	Solution path matrix with size d*p+1 by length(lambda); each column corresponds to one regularization parameter.
df	Degrees of freedom along the solution path (number of non-zero component functions).
knots	The p-1 by d matrix. Each column contains the knots applied to the corresponding variable.
Boundary.knots	The 2 by d matrix. Each column contains the boundary points applied to the corresponding variable.
func_norm	Functional norm matrix (d by length(lambda)); each column corresponds to one regularization parameter.

See Also

SAM,plot.samEL,print.samEL,predict.samEL

Examples

## generating training data
n = 200
d = 100
X = 0.5*matrix(runif(n*d),n,d) + matrix(rep(0.5*runif(n),d),n,d)
u = exp(-2*sin(X[,1]) + X[,2]^2-1/3 + X[,3]-1/2 + exp(-X[,4])+exp(-1)-1+1)
y = rep(0,n)
for(i in 1:n) y[i] = rpois(1,u[i])

## Training
out.trn = samEL(X,y)
out.trn

## plotting solution path
plot(out.trn)

## generating testing data
nt = 1000
Xt = 0.5*matrix(runif(nt*d),nt,d) + matrix(rep(0.5*runif(nt),d),nt,d)
ut = exp(-2*sin(Xt[,1]) + Xt[,2]^2-1/3 + Xt[,3]-1/2 + exp(-Xt[,4])+exp(-1)-1+1)
yt = rep(0,nt)
for(i in 1:nt) yt[i] = rpois(1,ut[i])

## predicting response
out.tst = predict(out.trn,Xt)

SAM documentation built on March 13, 2026, 9:08 a.m.