R/QR.lasso.admm.R

Defines functions QR.lasso.admm

Documented in QR.lasso.admm

QR.lasso.admm = function(X,y,tau,lambda,rho,beta,maxit)
{

if(missing(maxit)){
maxit = 200
}

n=dim(X)[1]
x=cbind(rep(1,n),X)

if(missing(beta)){
	beta=solve(t(x)%*%x,t(x)%*%y)
}

if(missing(rho)){
	rho=0.4
}

if(missing(lambda)){
	lambda=1
}


betah=QRPADMMCPP(X,y,beta,maxit,tau,rho,lambda)
beta=matrix(betah[-1])
b=betah[1]
return(list(beta=beta,b=b))
}

Try the cqrReg package in your browser

Any scripts or data that you put into this service are public.

cqrReg documentation built on June 7, 2022, 9:06 a.m.