R/penalty_D1_mcp.R

Defines functions penalty_D1_mcp

## File Name: penalty_D1_mcp.R
## File Version: 0.03


penalty_D1_mcp <- function(x, lambda, eps, a=2.7)
{
    x <- abs(x)
    res <- ifelse( x < a*lambda, lambda*sqrt(x^2 + eps) - x^2 / (2*a), .5*a*lambda^2)
    return(res)
}
alexanderrobitzsch/sirt documentation built on March 18, 2024, 1:29 p.m.