> Sys.time(); [1] "2014-10-21 10:57:43 EEST" > summary(m2, digits = 2, style = "pmax") Call: earth(formula=fmla, data=Dm, degree=2) ELTYY = 7.1 - 0.25 * pmax(0, 27 - DURURA) + 0.17 * pmax(0, ACT4 - -0.61) + 0.23 * pmax(0, ACT5 - 0.13) - 0.37 * pmax(0, 0.13 - ACT5) + 0.13 * pmax(0, 8 - SOSVERK) + 0.028 * pmax(0, 25 - SAIPV) + 0.00019 * pmax(0, ELAKETUb - 1563) - 7.9 * pmax(0, 1.7 - JOB2) * pmax(0, JOB4 - 2.3) - 1.2 * pmax(0, 1.7 - JOB2) * pmax(0, ACT2 - 1.5) + 0.11 * pmax(0, 1.7 - JOB2) * pmax(0, 1.5 - ACT2) + 3 * pmax(0, JOB2 - 1.7) * pmax(0, ACT3 - -0.031) - 0.0056 * pmax(0, JOB2 - 1.7) * pmax(0, 1850 - ELAKETUb) + 0.012 * pmax(0, JOB4 - 2.1) * pmax(0, 1563 - ELAKETUb) + 0.29 * pmax(0, ACT2 - -0.72) * pmax(0, 0.13 - ACT5) + 2 * pmax(0, ACT2 - 1.9) * pmax(0, 8 - SOSVERK) - 0.062 * pmax(0, 1.9 - ACT2) * pmax(0, 8 - SOSVERK) + 0.37 * pmax(0, 0.13 - ACT5) * pmax(0, 0.055 - ACT6) - 0.62 * pmax(0, 0.13 - ACT5) * pmax(0, ASUM01 - 0) - 0.011 * pmax(0, ACT6 - 1.7) * pmax(0, 1563 - ELAKETUb) - 0.00047 * pmax(0, 1.7 - ACT6) * pmax(0, 1563 - ELAKETUb) - 0.086 * pmax(0, 8 - SOSVERK) * pmax(0, 1 - LAHTO1) + 0.31 * pmax(0, 8 - SOSVERK) * pmax(0, -0.16 - LAHTO2) Selected 23 of 42 terms, and 14 of 21 predictors Importance: DURURA, SOSVERK, LAHTO1, SAIPV, ACT5, ASUM01, JOB2, ACT2, ... Number of terms at each degree of interaction: 1 7 15 GCV 1.2 RSS 1406 GRSq 0.22 RSq 0.28 > evimp(m2, trim=FALSE); nsubsets gcv rss DURURA 22 100.0 100.0 SOSVERK 21 80.5 84.7 LAHTO1 21 80.5 84.7 SAIPV 20 65.6 73.5 ACT5 19 50.1 62.6 ASUM01 19 50.1 62.6 JOB2 18 44.5 58.4 ACT2 18 44.5 58.4 ELAKETUb 15 38.9 52.1 ACT4 14 34.3 48.6 ACT3 13 30.4 45.4 LAHTO2 11 24.4 40.0 ACT6 10 20.7 37.0 JOB4 3 12.3 20.5 MIES-unused 0 0.0 0.0 IKATUTK-unused 0 0.0 0.0 EDUCI-unused 0 0.0 0.0 TUTKP-unused 0 0.0 0.0 JOB1-unused 0 0.0 0.0 JOB3-unused 0 0.0 0.0 ACT1-unused 0 0.0 0.0 > sink();