Группа LARS в R

Question

Я хотел бы проверить групповой алгоритм LARS в R.

Для этого я сначала создал следующий набор данных

set.seed(21021995)
n = 100                                             
mu = c(0,0,0,0)                         
Sigma = matrix(c(1,0.3,0.4,0.7,0.3,1,0.1,0.2,0.4,0.1,1,0.5,0.7,0.2,0.5,1), 
nrow=4, ncol=4) 
observations = mvrnorm(n, mu, Sigma)            
x1 = observations[,1]                     
x2 = observations[,2]                    
x3 = observations[,3]                     
x4 = observations[,4]                     
epsilon = rnorm(100,mean=0,sd=0.1)        
yvector = -0.5*x1+4*x2+2*x2^2-2*x3-x3^2-3*x3^3+3*x4+0.5*x4^2+epsilon   
dataframe = as.data.frame(cbind(yvector,x1,x2,x2^2,x3,x3^2,x3^3,x4,x4^2))
colnames(dataframe) = c("yvector","x1","x2","x2^2","x3","x3^2", 
"x3^3","x4","x4^2")
assign = c(1, 2, 2, 3, 3, 3, 4, 4)

Теперь я хочу применить group LARS следующим образом (используя библиотеку grpreg)

LARS=rgrplars(yvector ~ x1+x2+x2^2+x3+x3^2+x3^3+x4+x4^2,
         data=dataframe,assign=assign,sMax=15,robust=FALSE)

Но здесь я всегда получаю ошибку

Error in lmrob.S(x, y, control = control, mf = mf) : 
NA/NaN/Inf in foreign function call (arg 1)

Может кто-то помочь мне с этим?