ddegras/sparsegfl

Hybrid Approach to Sparse Group Fused Lasso

sparsegfl's Issues

Documentation

I've tested out your spatiotemporal lasso, and it is fantastically fast. Two questions:

1. Using the below code, how does one actually plot out the results? I see the beta matrix there, in the results, after this code, but it's not conformable with the X matrix, so unlike the genlasso package... I can't exactly plot the results.
   `
   m <- 5 # number of (time-varying) observation variables
   n <- 100 # time series length
   p <- 10 # number of (time-varying) predictor variables

x <- array(rnorm(mnp),c(m,p,n)) # random coefficients
beta <- matrix(0,p,n) # random coefficients with 50% sparsity
beta[1:5,1:50] <- runif(5) # and one change point at time 51
beta[6:p,51:n] <- runif(5)
y <- matrix(,m,n)
for (i in 1:n){ y[,i] <- x[,,i] %*% beta[,i] + rnorm(m,sd=.01)}

lambda1 <- .05
lambda2 <- 1.5
alpha <- .9

result <- SGFL(x, y, lambda1, lambda2, alpha,
intercept=FALSE, verbose=TRUE)

`
2. Is there an example that uses parallelization? Does one do this with SNOW?

type definitions not known

What are the definitions of the types
mat (a 2 dimensional array?)
cube (a 3 dimensional array?)
uword (an integer?)
in the allmethods.cpp?