MarronWand              package:nor1mix              R Documentation

_M_a_r_r_o_n-_W_a_n_d _D_e_n_s_i_t_i_e_s _a_s '_n_o_r_M_i_x' _O_b_j_e_c_t_s

_D_e_s_c_r_i_p_t_i_o_n:

     The fifteen density examples used in Marron and Wand (1992)'s
     simulation study have been used in quite a few subsequent studies,
     can all be written as normal mixtures and are provided here for
     convenience and didactical examples of normal mixtures. Number 16
     has been added by Jansen et al.

_U_s_a_g_e:

      MW.nm1  # Gaussian
      MW.nm2  # Skewed
      MW.nm3  # Str Skew
      MW.nm4  # Kurtotic
      MW.nm5  # Outlier
      MW.nm6  # Biomdal
      MW.nm7  # Separated (bimodal)
      MW.nm8  # Aymmettric Bimodal
      MW.nm9  # Trimodal
      MW.nm10 # Claw
      MW.nm11 # Double Claw
      MW.nm12 # Asymmetric Claw
      MW.nm13 # Asymm. Double Claw
      MW.nm14 # Smooth   Comb
      MW.nm15 # Discrete Comb
      MW.nm16 # Distant Bimodal

_A_u_t_h_o_r(_s):

     Martin Maechler

_S_o_u_r_c_e:

     The first part  are translated from Steve Marron's Matlab code at
     <URL:
     http://www.stat.unc.edu/postscript/papers/marron/parameters/nmpar.
     m>

_R_e_f_e_r_e_n_c_e_s:

     Marron, S. and Wand, M. (1992) Exact Mean Integrated Squared
     Error; _Annals of Statistcs_ *20*, 712-736.

     For number 16,
      Janssen, Marron, Verb..., Sarle (1995) ....

_E_x_a_m_p_l_e_s:

     MW.nm10
     plot(MW.nm14)

     ## These are defined as norMix() calls in  ../R/zMarrWand-dens.R
     require("nor1mix")
     ppos <- which("package:nor1mix" == search())
     nms <- ls(pat="^MW.nm", pos = ppos)
     nms <- nms[order(as.numeric(substring(nms,6)))]
     for(n in nms) {
        cat("\n",n,":\n"); print(get(n, pos = ppos))
     }

     ## Plot all of them:
     op <- par(mfrow=c(4,4), mgp = c(1.2, 0.5, 0), tcl = -0.2,
               mar = .1 + c(2,2,2,1), oma = c(0,0,3,0))
     for(n in nms) { plot(get(n, pos = ppos))}
     mtext("The Marron-Wand Densities", outer = TRUE,
           font = 2, cex = 1.6)
     par(op)

