data(genotype2_ExploredModels)
head(genotype2_ExploredModels)
plot(genotype2_ExploredModels[, c("dim", "logLik")],
col = "blue", xlab = "Dimension", ylab = "Log-likelihood")
# Data-driven calibration of the penalty
dimJump.R(genotype2_ExploredModels, h = 5, N=1000, header=T)
Run the code above in your browser using DataLab