# NOT RUN {
set.seed(14)
beta <- c(3, 2, -1.6, -4)
noise <- 5
simData <- SimData(N=100,beta=beta, noise=noise, corr=FALSE)
stepPenal<- StepPenal(Data=simData, lamda=1.2, w=0.7)
(coeffP <- stepPenal$coeffP)
findROC(simData, coeff=coeffP)
# }
Run the code above in your browser using DataLab