plot.ppr                package:stats                R Documentation

_P_l_o_t _R_i_d_g_e _F_u_n_c_t_i_o_n_s _f_o_r _P_r_o_j_e_c_t_i_o_n _P_u_r_s_u_i_t _R_e_g_r_e_s_s_i_o_n _F_i_t

_D_e_s_c_r_i_p_t_i_o_n:

     Plot ridge functions for projection pursuit regression fit.

_U_s_a_g_e:

     ## S3 method for class 'ppr':
     plot(x, ask, type = "o", ...)

_A_r_g_u_m_e_n_t_s:

       x: A fit of class '"ppr"' as produced by a call to 'ppr'. 

     ask: the graphics parameter 'ask': see 'par' for details. If set
          to 'TRUE' will ask between the plot of each cross-section. 

    type: the type of line to draw 

     ...: further graphical parameters 

_V_a_l_u_e:

     None

_S_i_d_e _E_f_f_e_c_t_s:

     A series of plots are drawn on the current graphical device, one
     for each term in the fit.

_S_e_e _A_l_s_o:

     'ppr', 'par'

_E_x_a_m_p_l_e_s:

     require(graphics)

     with(rock, {
     area1 <- area/10000; peri1 <- peri/10000
     par(mfrow=c(3,2))# maybe: , pty="s")
     rock.ppr <- ppr(log(perm) ~ area1 + peri1 + shape,
                     data = rock, nterms = 2, max.terms = 5)
     plot(rock.ppr, main="ppr(log(perm)~ ., nterms=2, max.terms=5)")
     plot(update(rock.ppr, bass=5), main = "update(..., bass = 5)")
     plot(update(rock.ppr, sm.method="gcv", gcvpen=2),
          main = "update(..., sm.method=\"gcv\", gcvpen=2)")
     })

