Слитый лассо с пакетом Genlasso в R - PullRequest
Новая
       9

Слитый лассо с пакетом Genlasso в R

0 голосов
/ 23 мая 2018 
n = 150

g <- matrix(0,n,Nrofreps)
X <- array(0, dim = c(Nrofreps,n,1))
  e <- rnorm(100+n,1)
  a <- c(rep(0,100+0.3*n), rep(5,0.4*n), rep(5,0.3*n))
  beta <- c(rep(0.0,100+0.7*n), rep(0.0,0.3*n))

  y <- rep(0,101+n)

  for (z in 1:(100+n)){
    y[z+1]<- a[z]+beta[z]*y[z]+e[z]
  }
  y2 <- y[101:(100+n)]
  g[,1]=y2
  X[1,,1]=rep(5,n)
 a1 = fusedlasso1d(g[,1],X = as.matrix(X[1,,]), minlam = 3, gamma = 0.1)

Пока это мой код.В будущем я хочу расширить его и быть в состоянии сделать больше перерисовок.По этой причине некоторые переменные имеют более высокие размеры, чем необходимо.Этот код ошибки: 

Error: (mak <- max(k)) <= m - 1 is not TRUE.

Я не знаю, почему это или как это исправить, пожалуйста, помогите

...