# APPLY PENALIZED LEAST SQUARES WITH CROSS VALIDATION TO DATASET #2. source("a1funcs.r") Xtrn <- as.matrix (read.table("Xtrn2.txt",header=TRUE)) ytrn <- scan("ytrn2.txt") Xtst <- as.matrix (read.table("Xtst2.txt",header=TRUE)) ytst <- scan("ytst2.txt") cat("\n") try <- c(0, 0.25, 0.5, 1, 2, 4, 8, 16, 32) cat("CV ARRAY:\n\n") cva <- cross.val.array(Xtrn,ytrn,20,try,10) print(round(sqrt(cva),2)) cat("\nBest value:",round(sqrt(min(cva)),2),"\n") cat("\n\n") cat("TEST ARRAY:\n\n") ta <- penalized.least.squares.test.array (Xtrn,ytrn,20,try,Xtst,ytst) print(round(sqrt(ta),2)) cat("\nBest value:",round(sqrt(min(ta)),2),"\n")