R中的genlasso程序包-对偶坐标超出了边界吗？

Question

我试图在R中使用软件包genlasso解决我的问题。但是，当我绘制对偶图时，某些对偶坐标似乎在边界之外，这没有任何意义。我不知道为什么。

我不知道如何在此处上传图片，但是以下代码可以生成该图

library(genlasso)
set.seed(95)
x = matrix(c(rnorm(3,2,1),rnorm(3,8,1)),ncol=1)
no.exp = nrow(x)
dim=ncol(x)/2
D = matrix(0,nrow=dim*no.exp*(no.exp-1),ncol=2*dim*no.exp)
Dindex=1
Drow=1-2*dim
repeat{
   j=Dindex+2*dim
   repeat{
       Drow=Drow+2*dim
       D[Drow:(Drow+2*dim-1),Dindex:(Dindex+2*dim-1)]=diag(rep(1,2*dim))
       D[Drow:(Drow+2*dim-1),j:(j+2*dim-1)]=-diag(rep(1,2*dim))

       if(j+2*dim-1==2*dim*no.exp){break}
       j=j+2*dim
   }
 Dindex=Dindex+2*dim
 if(Dindex==2*dim*no.exp+1-2*dim){break}
}
model = genlasso(x,D=D)
plot(model,type='dual')