# Try out knn and linear regression on the prostate cancer data. trn = read.table("prostate-train",head=F) tst = read.table("prostate-test",head=F) plot (trn[,9], trn[,1], xlab="Log of Prostate-Specific Antigen", ylab="Log of Prostate Cancer Volume") m = lm(trn[,1] ~ trn[,9]) print(summary(m)) abline(m) b0 = coef(m)[1] b1 = coef(m)[2] points (tst[,9], b0+tst[,9]*b1, col="red", pch=5) source("knn.r") pred_5nn = knn (5, matrix(trn[,9],ncol=1), trn[,1], matrix(tst[,9],ncol=1)) points (tst[,9], pred_5nn, col="green", pch=5)