[R] help with mars output

Russell Brown russ.brown at rogers.com
Sun May 4 17:24:37 CEST 2008


I need help in interpreting some output of polymars.

The returned model by my command is:

m$model
  pred1    knot1 pred2 knot2       coefs           SE
1     0       NA     0    NA  1.32222163 0.0007806758
2     1       NA     0    NA -0.10904285 0.0006735827
3     1 1.193575     0    NA  1.25396217 0.0205092438
4     1 1.205400     0    NA -0.18020096 0.0261931173
5     1 1.230725     0    NA  0.04709508 0.0118559207
6     1 1.285400     0    NA -0.02593593 0.0097733958
7     1 1.317275     0    NA  0.03516571 0.0177043207


I'm trying to confirm the generated output of mars by manually entering the basis functions into my dataset, What I don't understand is that in the knot1 column there are to rows containing NA.  In the polymars literature I was reading this means that these basis functions are linear.

BF1 = MAX(0, A2 - 1.193575)
BF2 =MAX(0, 1.193575 -  A2)
BF3 = MAX(0, A2 - 1.2054)
BF5 = MAX(0, A2 - 1.230725)
BF7 = MAX(0, A2 - 1.2854)
BF11 = MAX(0, A2 - 1.317275)

RESPONSE = 1.32222163 - 0.10904285 * C2 + 1.25396217 * D2 - 0.18020096 *         E2 + 0.04709508 * F2 - 0.02593593 * G2 + 0.03516571 * H2

On the attached, the CLOSE column is my source data, the ESTIMATE is the column generated by polymars, the PREDICTION column is my formula.

Can I have some help with this please.







More information about the R-help mailing list