R : Copyright 2002, The R Development Core Team Version 1.5.1 (2002-06-17) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type `license()' or `licence()' for distribution details. R is a collaborative project with many contributors. Type `contributors()' for more information. Type `demo()' for some demos, `help()' for on-line help, or `help.start()' for a HTML browser interface to help. Type `q()' to quit R. > invisible(options(echo = TRUE)) > library(rpart) > > target.lrn <- read.table("../lrn/num/472.dat",header=T,colClasses="numeric") > target.val <- read.table("../val/num/472.dat",header=T,colClasses="numeric") > target.tst <- read.table("../tst/num/472.dat",header=T,colClasses="numeric") > > > y.lrn <- target.lrn[,1] > y.val <- target.val[,1] > y.tst <- target.tst[,1] > y <- c(y.lrn,y.val,y.tst) > > n.lrn <- length(y.lrn) > n.val <- length(y.val) > n.tst <- length(y.tst) > n <- length(y) > > rm(target.lrn,target.val,target.tst) > > wts <- mat.or.vec(n,1) ; for (i in 1:n.lrn) wts[i]=1 > idx.lrn <- 1:n.lrn > idx.val <- (n.lrn+1):(n.lrn+n.val) > idx.tst <- (n.lrn+n.val+1):n > > mod <- read.table("../cty_mod.txt", + header=F,colClasses="character",col.names=c("file","feature","type")) > > n.mod <- length(mod$file) > > first.time <- TRUE > > for (i in 1:n.mod) { + + fn.lrn <- paste("../lrn/",mod$type[i],"/",mod$file[i],".dat",sep="") + fn.val <- paste("../val/",mod$type[i],"/",mod$file[i],".dat",sep="") + fn.tst <- paste("../tst/",mod$type[i],"/",mod$file[i],".dat",sep="") + print(mod$feature[i]) + + if (mod$type[i]=="chr") { + + f.lrn <- read.table(fn.lrn, + header=T,colClasses="character",blank.lines.skip=F) + f.val <- read.table(fn.val, + header=T,colClasses="character",blank.lines.skip=F) + f.tst <- read.table(fn.tst, + header=T,colClasses="character",blank.lines.skip=F) + + f <- c(f.lrn[,1],f.val[,1],f.tst[,1]) + + if (mod$feature[i]=="STATE") { + f[f=="AS"|f=="DC"|f=="DE"|f=="MA"|f=="ME"|f=="NH"] <- "S1" + f[f=="OH"|f=="RI"|f=="VI"|f=="WV"] <- "S1" + f[f=="AA"|f=="AE"|f=="AP"|f=="CT"|f=="GU"|f=="MD"] <- "S2" + f[f=="NJ"|f=="NY"|f=="PA"|f=="PA"|f=="VA"|f=="VT"] <- "S2" + f[f=="WY"] <- "S2" + f[f=="AK"|f=="UT"|f=="MS"] <- "S3" + f[f=="NE"|f=="ND"] <- "S4" + f[f=="SD"|f=="SC"] <- "S5" + } + + f <- as.factor(f) + + n.lev <- nlevels(f) + f.name <- levels(f) + if(n.lev==2) f.name <- c(mod$feature[i],mod$feature[i]) + print(paste(" nlevels = ",n.lev)) + + f <- model.matrix(y ~ f - 1) # Note: Intercept removed. + f <- f[,2:n.lev] # Note: First dummy deleted. + f.name <- f.name[2:n.lev] # Note: First name deleted. + + } else { + + f.lrn<-read.table(fn.lrn, + header=T,colClasses="numeric",blank.lines.skip=F) + f.val<-read.table(fn.val, + header=T,colClasses="numeric",blank.lines.skip=F) + f.tst<-read.table(fn.tst, + header=T,colClasses="numeric",blank.lines.skip=F) + + f <- c(f.lrn[,1],f.val[,1],f.tst[,1]) + + f[is.na(f)] <- 0 + + f.name <- mod$feature[i] + + } + + if (first.time) { + X <- f + X.names <- f.name + first.time <- FALSE + } else { + X <- cbind(prev.X,f) + X.names <- c(prev.X.names,f.name) + } + + prev.X <- X + prev.X.names <- X.names + + } [1] "LASTGIFT" [1] "PEPSTRFL" [1] " nlevels = 2" [1] "STATE" [1] " nlevels = 33" [1] "RECP3" [1] " nlevels = 2" [1] "DOB" [1] "MAILCODE" [1] " nlevels = 2" [1] "MHUC2" [1] "LASTDATE" [1] "MINRAMNT" > > rm(prev.X,prev.X.names) > rm(f.lrn,f.val,f) > > dimnames(X) <- list(NULL,X.names) > > ctrl <- rpart.control(minsplit=50,maxsurrogate=0,maxdepth=15,cp=0.0001) > tree <- rpart(y~X,weights=wts,method="anova",control=ctrl) > > source("psopts.r") > postscript(file="cty_tree_0001.eps") > par.default <- par(no.readonly=TRUE) > par(mar=c(0,0,0,0)) > plot(tree) > #text(tree) > par(par.default) > dev.off() null device 1 > > print(ctrl) $minsplit [1] 50 $minbucket [1] 17 $cp [1] 1e-04 $maxcompete [1] 4 $maxsurrogate [1] 0 $usesurrogate [1] 2 $surrogatestyle [1] 0 $maxdepth [1] 15 $xval [1] 10 > print(tree$frame) var n wt dev yval complexity ncompete 1 X.LASTGIFT 95412 66900 1.344638e+06 0.7887979 3.446296e-03 4 2 X.LASTGIFT 94735 66431 1.163124e+06 0.7666839 1.014759e-03 4 4 X.LASTDATE 87823 61590 8.612719e+05 0.7265038 6.310128e-04 4 8 X.CA 84505 59297 8.117689e+05 0.7034229 3.088016e-04 4 16 X.DOB 69620 48887 6.425555e+05 0.6671457 3.088016e-04 4 32 X.DOB 21836 15295 1.821454e+05 0.5705459 2.599639e-04 4 64 5053 3524 2.205468e+04 0.4016742 4.837636e-05 0 65 X.LASTDATE 16783 11771 1.599601e+05 0.6211027 2.599639e-04 4 130 10159 7101 6.929827e+04 0.5355584 3.769154e-05 0 131 X.GA 6624 4670 9.053087e+04 0.7511777 2.599639e-04 4 262 6225 4386 4.735445e+04 0.7081623 9.626605e-05 0 263 X.LASTGIFT 399 284 4.304297e+04 1.4154930 2.599639e-04 4 526 350 249 2.861181e+03 0.7108434 4.039607e-05 0 527 49 35 3.917857e+04 6.4285714 1.000000e-04 0 33 X.MINRAMNT 47784 33592 4.602024e+05 0.7111291 3.088016e-04 4 66 X.DOB 47713 33539 4.559886e+05 0.7077805 3.088016e-04 4 132 X.PEPSTRFL 47668 33507 4.461139e+05 0.7054720 3.088016e-04 4 264 X.S2 24941 17529 2.371947e+05 0.6392835 3.088016e-04 4 528 X.MINRAMNT 24796 17429 2.328205e+05 0.6331975 3.088016e-04 4 1056 X.DOB 23498 16532 1.903260e+05 0.6152311 1.200120e-04 4 2112 23355 16435 1.872127e+05 0.6099179 4.377387e-05 0 2113 X.LASTGIFT 143 97 3.034227e+03 1.5154639 1.200120e-04 4 4226 91 60 5.341833e+02 0.6166667 1.131242e-05 0 4227 X.MINRAMNT 52 37 2.372973e+03 2.9729730 1.200120e-04 3 8454 27 17 0.000000e+00 0.0000000 1.000000e-04 0 8455 25 20 2.095000e+03 5.5000000 1.000000e-04 0 1057 X.LASTDATE 1298 897 4.239086e+04 0.9643255 3.088016e-04 4 2114 X.MO 1151 804 1.780634e+04 0.7835821 1.008615e-04 4 4228 1121 780 1.603010e+04 0.7115385 7.359709e-05 0 4229 30 24 1.640625e+03 3.1250000 1.000000e-04 0 2115 X.FL 147 93 2.433118e+04 2.5268817 3.088016e-04 4 4230 126 82 2.336890e+03 1.0365854 4.362976e-05 0 4231 21 11 2.045455e+04 13.6363636 1.000000e-04 0 529 X.MHUC2 145 100 4.261000e+03 1.7000000 1.089537e-04 4 1058 117 79 1.533544e+03 1.0759494 9.922053e-05 0 1059 28 21 2.580952e+03 4.0476190 1.000000e-04 0 265 X.MINRAMNT 22727 15978 2.087582e+05 0.7780855 3.088016e-04 4 530 X.LASTGIFT 20594 14518 1.714173e+05 0.7337960 3.088016e-04 4 1060 X.HI 19096 13456 1.326969e+05 0.7033479 1.566053e-04 4 2120 X.LASTGIFT 19061 13425 1.293410e+05 0.6984171 1.566053e-04 4 4240 X.MO 10081 7099 4.952138e+04 0.6121989 1.566053e-04 4 8480 9668 6797 3.745602e+04 0.5933500 4.143485e-05 0 8481 X.LASTDATE 413 302 1.200860e+04 1.0364238 1.566053e-04 4 16962 338 255 1.891349e+03 0.7019608 2.472518e-05 0 16963 X.DOB 75 47 9.933957e+03 2.8510638 1.566053e-04 3 33926 57 36 2.838889e+02 0.9444444 1.000000e-04 0 33927 18 11 9.090909e+03 9.0909091 1.000000e-04 0 4241 8980 6326 7.970765e+04 0.7951707 4.843901e-05 0 2121 35 31 3.214194e+03 2.8387097 1.000000e-04 0 1061 X.LA 1498 1062 3.854981e+04 1.1195857 3.088016e-04 4 2122 X.AZ 1476 1047 2.781499e+04 0.9971347 1.609494e-04 4 4244 X.WA 1436 1015 2.384201e+04 0.9231527 1.609494e-04 4 8488 X.LASTGIFT 1363 963 1.965675e+04 0.8348910 1.079469e-04 4 16976 X.FL 1204 848 1.408390e+04 0.6945755 1.079469e-04 4 33952 1066 756 9.843286e+03 0.5476190 1.000000e-04 0 33953 138 92 4.090120e+03 1.9021739 1.000000e-04 0 16977 159 115 5.433043e+03 1.8695652 8.303197e-05 0 8489 X.DOB 73 52 4.038827e+03 2.5576923 1.609494e-04 3 16978 38 29 9.655172e+01 0.3448276 1.000000e-04 0 16979 35 23 3.621217e+03 5.3478261 1.000000e-04 0 4245 40 32 3.791219e+03 3.3437500 1.000000e-04 0 2123 22 15 9.623333e+03 9.6666667 1.000000e-04 0 531 X.KS 2133 1460 3.702930e+04 1.2184932 2.104615e-04 4 1062 X.S4 2099 1441 3.527336e+04 1.1679389 1.988769e-04 4 2124 X.MHUC2 2071 1422 3.203215e+04 1.1181435 1.967935e-04 4 4248 X.DOB 1870 1283 2.482399e+04 0.9968823 1.529680e-04 4 8496 674 443 5.205530e+03 0.6591422 9.097673e-05 0 8497 X.DOB 1196 840 1.954128e+04 1.1750000 1.529680e-04 4 16994 935 646 1.159241e+04 0.8916409 7.610955e-05 0 16995 X.NC 261 194 7.724273e+03 2.1185567 1.529680e-04 4 33990 244 182 6.086066e+03 1.7912088 1.000000e-04 0 33991 17 12 1.322917e+03 7.0833333 1.000000e-04 0 4249 X.DOB 201 139 7.015165e+03 2.2374101 1.967935e-04 4 8498 155 105 2.477962e+03 1.3523810 6.525073e-05 0 8499 46 34 4.200971e+03 4.9705882 1.000000e-04 0 2125 28 19 2.973789e+03 4.8947368 1.000000e-04 0 1063 34 19 1.472947e+03 5.0526316 1.000000e-04 0 133 45 32 9.687500e+03 3.1250000 1.000000e-04 0 67 X.LASTGIFT 71 53 3.975472e+03 2.8301887 3.088016e-04 4 134 54 38 1.585526e+02 0.6578947 6.360547e-06 0 135 17 15 3.183333e+03 8.3333333 1.000000e-04 0 17 X.LASTGIFT 14885 10410 1.688469e+05 0.8737867 1.684229e-04 4 34 4264 2972 1.766347e+04 0.6404509 2.715973e-05 0 35 X.DOB 10621 7438 1.509570e+05 0.9670207 1.487544e-04 4 70 3662 2580 4.259179e+04 0.7980620 4.588091e-05 0 71 X.DOB 6959 4858 1.082524e+05 1.0567517 1.487544e-04 4 142 1210 842 1.156165e+04 0.6235154 9.586147e-05 0 143 X.RECP3 5749 4016 9.649962e+04 1.1475847 1.487544e-04 4 286 X.MHUC2 5669 3966 9.323081e+04 1.1292738 1.487544e-04 4 572 X.MINRAMNT 5073 3550 8.137875e+04 1.0756901 1.487544e-04 4 1144 X.LASTDATE 4848 3391 6.704332e+04 1.0371572 1.127426e-04 4 2288 3342 2349 4.389636e+04 0.9429544 3.344251e-05 0 2289 X.DOB 1506 1042 2.307912e+04 1.2495202 1.127426e-04 4 4578 1481 1023 2.122008e+04 1.1847507 3.291781e-05 0 4579 25 19 1.623684e+03 4.7368421 1.000000e-04 0 1145 X.LASTDATE 225 159 1.422302e+04 1.8974843 1.487544e-04 2 2290 123 89 1.323091e+03 0.6820225 1.505037e-05 0 2291 X.LASTDATE 102 70 1.260127e+04 3.4428571 1.487544e-04 2 4582 48 35 9.671429e+02 1.2857143 1.000000e-04 0 4583 X.DOB 54 35 1.130840e+04 5.6000000 1.487544e-04 1 9166 30 17 7.821176e+02 2.4117647 1.000000e-04 0 9167 24 18 1.019028e+04 8.6111111 1.000000e-04 0 573 X.LASTGIFT 596 416 1.175488e+04 1.5865385 1.487544e-04 4 1146 356 245 5.627984e+03 0.9918367 8.762257e-05 0 1147 X.LASTGIFT 240 171 5.916105e+03 2.4385965 1.464345e-04 4 2294 X.LASTDATE 211 152 4.398467e+03 2.0592105 1.023842e-04 4 4588 X.DOB 175 121 2.598959e+03 1.6033058 1.023842e-04 4 9176 94 67 5.261194e+02 0.5970149 1.457171e-05 0 9177 81 54 1.920815e+03 2.8518519 6.613374e-05 0 4589 36 31 1.676194e+03 3.8387097 1.000000e-04 0 2295 29 19 1.320737e+03 5.4736842 1.000000e-04 0 287 X.DOB 80 50 3.162000e+03 2.6000000 1.487544e-04 4 574 X.LASTDATE 59 37 1.191892e+03 0.9459459 1.055242e-04 4 1148 41 30 0.000000e+00 0.0000000 1.000000e-04 0 1149 18 7 1.050000e+03 5.0000000 1.000000e-04 0 575 21 13 1.580769e+03 7.3076923 1.000000e-04 0 9 X.LASTGIFT 3318 2293 4.865447e+04 1.3233755 3.211386e-04 4 18 1858 1265 1.233696e+04 0.9509881 7.324178e-05 0 19 X.RECP3 1460 1028 3.592622e+04 1.7816148 3.211386e-04 4 38 X.IL 1141 805 2.233470e+04 1.4248447 1.917314e-04 4 76 X.CA 1066 755 1.789554e+04 1.2966887 1.649161e-04 4 152 X.DOB 823 575 1.064208e+04 1.0400000 1.117416e-04 4 304 234 154 3.441558e+02 0.1948052 6.180293e-06 0 305 589 421 1.014767e+04 1.3491686 8.767326e-05 0 153 X.MINRAMNT 243 180 7.094550e+03 2.1166667 1.649161e-04 4 306 218 162 4.566179e+03 1.6975309 8.474370e-05 0 307 25 18 2.243778e+03 5.8888889 1.000000e-04 0 77 X.DOB 75 50 4.239520e+03 3.3600000 1.917314e-04 4 154 X.LASTDATE 51 35 1.147143e+03 1.7142857 1.147415e-04 4 308 29 21 0.000000e+00 0.0000000 1.000000e-04 0 309 22 14 9.928571e+02 4.2857143 1.000000e-04 0 155 24 15 2.776400e+03 7.2000000 1.000000e-04 0 39 X.LASTGIFT 319 223 1.311917e+04 3.0695067 2.874122e-04 4 78 X.LASTGIFT 182 127 4.534256e+03 2.1535433 1.530425e-04 4 156 164 114 3.085546e+03 1.7236842 9.829945e-05 0 157 18 13 1.242923e+03 5.9230769 1.000000e-04 0 79 X.CA 137 96 8.337406e+03 4.2812500 2.874122e-04 4 158 X.LASTDATE 106 72 4.778653e+03 2.9305556 1.710704e-04 4 316 68 48 1.594667e+03 1.6666667 7.082807e-05 0 317 38 24 2.953958e+03 5.4583333 1.000000e-04 0 159 31 24 3.033333e+03 8.3333333 1.000000e-04 0 5 X.ID 6912 4841 3.004873e+05 1.2778785 7.653015e-04 4 10 X.LASTDATE 6881 4821 2.942413e+05 1.2497843 7.653015e-04 4 20 X.CA 6648 4660 2.712720e+05 1.1721481 7.238545e-04 4 40 X.LASTDATE 5117 3601 1.872529e+05 1.0041655 7.201795e-04 4 80 X.CO 4193 2943 9.760063e+04 0.8579681 4.303385e-04 4 160 X.LASTDATE 4064 2848 8.699592e+04 0.7998596 2.711977e-04 4 320 X.KY 986 689 6.382824e+03 0.3236575 2.523923e-04 4 640 964 676 3.520371e+03 0.2263314 2.654250e-05 0 641 22 13 2.523077e+03 5.3846154 1.000000e-04 0 321 X.OR 3078 2159 8.040699e+04 0.9518296 2.711977e-04 4 642 X.HI 2982 2085 7.337324e+04 0.8978417 2.711977e-04 4 1284 X.MI 2958 2066 7.036598e+04 0.8697967 2.711977e-04 4 2568 X.GA 2720 1900 5.314348e+04 0.7805263 1.997897e-04 4 5136 X.LASTGIFT 2579 1793 4.393865e+04 0.7133296 1.012011e-04 4 10272 X.DOB 695 480 6.203125e+03 0.3125000 1.012011e-04 4 20544 678 467 3.728587e+03 0.2141328 6.885537e-05 0 20545 17 13 2.307692e+03 3.8461538 1.000000e-04 0 10273 1884 1313 3.763021e+04 0.8598629 6.982776e-05 0 5137 X.MHUC2 141 107 9.061065e+03 1.9065421 1.997897e-04 4 10274 96 72 1.657653e+03 0.5694444 9.628687e-05 0 10275 45 35 7.009886e+03 4.6571429 1.000000e-04 0 2569 X.LASTGIFT 238 166 1.703405e+04 1.8915663 2.711977e-04 4 5138 X.MINRAMNT 218 154 6.388909e+03 1.1818182 1.500927e-04 4 10276 X.MINRAMNT 193 139 4.063755e+03 0.8057554 1.101654e-04 4 20552 141 102 7.439608e+02 0.3137255 1.791867e-05 0 20553 X.DOB 52 37 3.227027e+03 2.1621622 1.101654e-04 4 41106 28 20 0.000000e+00 0.0000000 1.000000e-04 0 41107 24 17 3.023529e+03 4.7058824 1.000000e-04 0 10277 25 15 2.123333e+03 4.6666667 1.000000e-04 0 5139 20 12 9.572000e+03 11.0000000 1.000000e-04 0 1285 24 19 2.828947e+03 3.9473684 1.000000e-04 0 643 X.PEPSTRFL 96 74 6.856446e+03 2.4729730 2.711977e-04 4 1286 49 34 0.000000e+00 0.0000000 1.000000e-04 0 1287 47 40 6.471775e+03 4.5750000 1.000000e-04 0 161 X.LASTDATE 129 95 1.030680e+04 2.6000000 4.303385e-04 4 322 X.LASTGIFT 110 82 4.763720e+03 1.4024390 1.769179e-04 4 644 80 59 3.932203e+02 0.3389831 2.470579e-05 0 645 30 23 4.132609e+03 4.1304348 1.000000e-04 0 323 19 13 4.683692e+03 10.1538462 1.000000e-04 0 81 X.LASTDATE 924 658 8.930806e+04 1.6580547 7.201795e-04 4 162 X.FL 644 455 2.142831e+04 0.9230769 3.671891e-04 4 324 X.LASTGIFT 563 398 7.452244e+03 0.5904523 1.541469e-04 4 648 X.DOB 528 375 5.448149e+03 0.4346667 1.541469e-04 4 1296 317 221 6.221719e+02 0.1131222 3.365134e-06 0 1297 X.DOB 211 154 4.770338e+03 0.8961039 1.541469e-04 4 2594 190 142 1.824648e+03 0.4225352 2.688776e-05 0 2595 21 12 2.537000e+03 6.5000000 1.000000e-04 0 649 35 23 1.846609e+03 3.1304348 1.000000e-04 0 325 X.DOB 81 57 1.362456e+04 3.2456140 3.671891e-04 4 650 49 36 6.076389e+02 0.6944444 1.000000e-04 0 651 32 21 1.238095e+04 7.6190476 1.000000e-04 0 163 X.MINRAMNT 280 203 6.708306e+04 3.3054187 7.201795e-04 4 326 X.DOB 198 142 9.507218e+03 1.3732394 2.054519e-04 4 652 117 78 6.169872e+02 0.3205128 3.277501e-05 0 653 X.DOB 81 64 8.698438e+03 2.6562500 2.054519e-04 4 1306 62 50 2.802000e+03 1.4000000 9.275901e-05 0 1307 19 14 5.535714e+03 7.1428571 1.000000e-04 0 327 X.DOB 82 61 5.581164e+04 7.8032787 7.201795e-04 3 654 34 24 8.625000e+02 1.2500000 1.000000e-04 0 655 48 37 5.324989e+04 12.0540541 1.000000e-04 0 41 X.MHUC2 1531 1059 8.357193e+04 1.7433522 7.238545e-04 4 82 X.MHUC2 1502 1041 7.750283e+04 1.6130740 7.238545e-04 4 164 X.LASTGIFT 1327 916 4.960097e+04 1.3681332 4.411132e-04 4 328 X.LASTDATE 571 381 9.291276e+03 0.7401575 2.332239e-04 4 656 432 287 3.367115e+03 0.4041812 2.297669e-05 0 657 X.LASTDATE 139 94 5.792851e+03 1.7659574 2.332239e-04 4 1314 69 44 0.000000e+00 0.0000000 1.000000e-04 0 1315 X.PEPSTRFL 70 50 5.534880e+03 3.3200000 2.332239e-04 4 2630 47 37 2.432432e+03 1.3513514 1.000000e-04 0 2631 23 13 2.550923e+03 8.9230769 1.000000e-04 0 329 X.LASTDATE 756 535 4.005244e+04 1.8153458 4.411132e-04 4 658 X.LASTGIFT 281 194 4.068521e+03 0.7268041 1.008299e-04 4 1316 262 181 2.875249e+03 0.5027624 8.278859e-05 0 1317 19 13 1.057692e+03 3.8461538 1.000000e-04 0 659 X.MHUC2 475 341 3.562327e+04 2.4346334 4.411132e-04 4 1318 X.DOB 320 235 1.956444e+04 1.6638298 4.411132e-04 4 2636 124 90 0.000000e+00 0.0000000 1.000000e-04 0 2637 X.LASTDATE 196 145 1.916065e+04 2.6965517 4.411132e-04 4 5274 X.PEPSTRFL 137 102 5.696755e+03 1.2843137 1.939406e-04 4 10548 84 62 0.000000e+00 0.0000000 1.000000e-04 0 10549 53 40 5.435975e+03 3.2750000 3.652386e-05 0 5275 X.DOB 59 43 1.277791e+04 6.0465116 4.411132e-04 4 10550 38 27 2.407407e+03 1.8518519 1.000000e-04 0 10551 21 16 9.093750e+03 13.1250000 1.000000e-04 0 1319 X.LASTDATE 155 106 1.560966e+04 4.1434906 4.411132e-04 4 2638 44 30 0.000000e+00 0.0000000 1.000000e-04 0 2639 X.MINRAMNT 111 76 1.489129e+04 5.7790789 3.480662e-04 4 5278 X.DOB 68 47 6.110638e+03 3.8297872 1.603548e-04 3 10556 29 18 3.777778e+02 1.1111111 1.000000e-04 0 10557 39 29 5.517241e+03 5.5172414 1.000000e-04 0 5279 43 29 8.312634e+03 8.9382759 1.000000e-04 0 165 X.DOB 175 125 2.744419e+04 3.4080000 7.238545e-04 4 330 X.DOB 114 80 3.703550e+03 1.0750000 1.604273e-04 4 660 76 56 0.000000e+00 0.0000000 1.000000e-04 0 661 38 24 3.487833e+03 3.5833333 1.000000e-04 0 331 X.MINRAMNT 61 45 2.253111e+04 7.5555556 7.238545e-04 4 662 35 23 1.171739e+03 1.5217391 1.000000e-04 0 663 26 22 1.964659e+04 13.8636364 1.000000e-04 0 83 29 18 5.029611e+03 9.2777778 1.000000e-04 0 21 X.MHUC2 233 161 2.212825e+04 3.4968944 7.653015e-04 4 42 X.DOB 216 151 1.863270e+04 2.9072848 7.653015e-04 4 84 X.LASTGIFT 143 99 5.627293e+03 1.3131313 1.404597e-04 4 168 74 52 0.000000e+00 0.0000000 1.000000e-04 0 169 X.DOB 69 47 5.438426e+03 2.7659574 1.380199e-04 4 338 24 16 0.000000e+00 0.0000000 1.000000e-04 0 339 45 31 5.252839e+03 4.1935484 1.000000e-04 0 85 X.LASTDATE 73 52 1.227483e+04 5.9423077 7.653015e-04 4 170 X.DOB 53 39 3.209692e+03 2.5384615 1.967328e-04 4 340 28 20 0.000000e+00 0.0000000 1.000000e-04 0 341 25 19 2.945158e+03 5.2105263 1.000000e-04 0 171 20 13 7.257692e+03 16.1538462 1.000000e-04 0 43 17 10 2.650400e+03 12.4000000 1.000000e-04 0 11 31 20 5.324950e+03 8.0500000 1.000000e-04 0 3 X.PEPSTRFL 677 469 1.768801e+05 3.9211087 3.155550e-03 4 6 X.MHUC2 428 296 3.762162e+04 1.6216216 9.874031e-04 4 12 X.MINRAMNT 385 275 1.907964e+04 1.0363636 5.528384e-04 4 24 X.LASTGIFT 366 264 9.245360e+03 0.7007576 4.430923e-04 4 48 311 223 8.428251e+02 0.1793722 1.994797e-05 0 49 X.LASTDATE 55 41 8.012195e+03 3.5365854 4.430923e-04 3 98 38 25 0.000000e+00 0.0000000 1.000000e-04 0 99 17 16 7.210938e+03 9.0625000 1.000000e-04 0 25 19 11 9.090909e+03 9.0909091 1.000000e-04 0 13 43 21 1.721429e+04 9.2857143 1.000000e-04 0 7 X.LASTDATE 249 173 1.350154e+05 7.8554913 2.738128e-03 4 14 X.DOB 210 148 7.440576e+04 5.9594595 1.202532e-03 4 28 X.DOB 179 124 5.045209e+04 4.6532258 1.202532e-03 4 56 74 46 0.000000e+00 0.0000000 1.000000e-04 0 57 X.LASTDATE 105 78 4.886868e+04 7.3974359 1.202532e-03 4 114 X.LASTGIFT 50 37 1.001081e+04 3.2432432 6.841260e-04 3 228 33 26 0.000000e+00 0.0000000 1.000000e-04 0 229 17 11 9.090909e+03 10.9090909 1.000000e-04 0 115 X.LASTGIFT 55 41 3.764312e+04 11.1463415 1.202532e-03 4 230 17 13 0.000000e+00 0.0000000 1.000000e-04 0 231 38 28 3.527811e+04 16.3214286 1.000000e-04 0 29 31 24 2.264896e+04 12.7083333 1.000000e-04 0 15 39 25 5.692784e+04 19.0800000 1.000000e-04 0 nsurrogate 1 0 2 0 4 0 8 0 16 0 32 0 64 0 65 0 130 0 131 0 262 0 263 0 526 0 527 0 33 0 66 0 132 0 264 0 528 0 1056 0 2112 0 2113 0 4226 0 4227 0 8454 0 8455 0 1057 0 2114 0 4228 0 4229 0 2115 0 4230 0 4231 0 529 0 1058 0 1059 0 265 0 530 0 1060 0 2120 0 4240 0 8480 0 8481 0 16962 0 16963 0 33926 0 33927 0 4241 0 2121 0 1061 0 2122 0 4244 0 8488 0 16976 0 33952 0 33953 0 16977 0 8489 0 16978 0 16979 0 4245 0 2123 0 531 0 1062 0 2124 0 4248 0 8496 0 8497 0 16994 0 16995 0 33990 0 33991 0 4249 0 8498 0 8499 0 2125 0 1063 0 133 0 67 0 134 0 135 0 17 0 34 0 35 0 70 0 71 0 142 0 143 0 286 0 572 0 1144 0 2288 0 2289 0 4578 0 4579 0 1145 0 2290 0 2291 0 4582 0 4583 0 9166 0 9167 0 573 0 1146 0 1147 0 2294 0 4588 0 9176 0 9177 0 4589 0 2295 0 287 0 574 0 1148 0 1149 0 575 0 9 0 18 0 19 0 38 0 76 0 152 0 304 0 305 0 153 0 306 0 307 0 77 0 154 0 308 0 309 0 155 0 39 0 78 0 156 0 157 0 79 0 158 0 316 0 317 0 159 0 5 0 10 0 20 0 40 0 80 0 160 0 320 0 640 0 641 0 321 0 642 0 1284 0 2568 0 5136 0 10272 0 20544 0 20545 0 10273 0 5137 0 10274 0 10275 0 2569 0 5138 0 10276 0 20552 0 20553 0 41106 0 41107 0 10277 0 5139 0 1285 0 643 0 1286 0 1287 0 161 0 322 0 644 0 645 0 323 0 81 0 162 0 324 0 648 0 1296 0 1297 0 2594 0 2595 0 649 0 325 0 650 0 651 0 163 0 326 0 652 0 653 0 1306 0 1307 0 327 0 654 0 655 0 41 0 82 0 164 0 328 0 656 0 657 0 1314 0 1315 0 2630 0 2631 0 329 0 658 0 1316 0 1317 0 659 0 1318 0 2636 0 2637 0 5274 0 10548 0 10549 0 5275 0 10550 0 10551 0 1319 0 2638 0 2639 0 5278 0 10556 0 10557 0 5279 0 165 0 330 0 660 0 661 0 331 0 662 0 663 0 83 0 21 0 42 0 84 0 168 0 169 0 338 0 339 0 85 0 170 0 340 0 341 0 171 0 43 0 11 0 3 0 6 0 12 0 24 0 48 0 49 0 98 0 99 0 25 0 13 0 7 0 14 0 28 0 56 0 57 0 114 0 228 0 229 0 115 0 230 0 231 0 29 0 15 0 > summary.rpart(tree,file="cty_tree_0001.out") > > yhat <- tree$frame[tree$where,5] > ehat <- y - yhat > mse.lrn <- sum(ehat[idx.lrn]^2)/n.lrn > mse.val <- sum(ehat[idx.val]^2)/n.val > mse.tst <- sum(ehat[idx.tst]^2)/n.tst > > print(paste(" mse.lrn = ",mse.lrn)) [1] " mse.lrn = 19.0171539150594" > print(paste(" mse.val = ",mse.val)) [1] " mse.val = 19.6427237028794" > print(paste(" mse.tst = ",mse.tst)) [1] " mse.tst = 18.9090330593579" > > x0 <- c(0,n.lrn) > y0 <- c(0,sum(y.lrn-0.68)) > > x1 <- (1:n.lrn) > y1 <- yhat[1:n.lrn] > y1 <- y1-0.68 > y1 <- sort(y1) > y1 <- y1[n.lrn:1] > y1 <- cumsum(y1) > > idx <- 1:n.lrn > > print(paste("maximum profit in learning sample is ", max(y1))) [1] "maximum profit in learning sample is 12806.7" > print(paste("maximum occurs at ", idx[y1==max(y1)])) [1] "maximum occurs at 24572" > > idx <- seq(1,n.lrn,length=200) > > x1 <- x1[idx] > y1 <- y1[idx] > > source("psopts.r"); > postscript(file="cty_lif_0001.eps"); > > plot(x=c(x0,x1),y=c(y0,y1),ylab="dollars",xlab="size of mailing",type="n") > lines(x=x0,y=y0,col="green") > lines(x=x1,y=y1,col="red") > > dev.off() null device 1 > proc.time() [1] 2486.81 8.80 2498.87 0.00 0.00 >