plot.scaleboot           package:scaleboot           R Documentation

_P_l_o_t _D_i_a_g_n_o_s_t_i_c_s _f_o_r _M_u_l_t_i_s_c_a_l_e _B_o_o_t_s_t_r_a_p

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

     'plot' method for class '"scaleboot"'.

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

     ## S3 method for class 'scaleboot':
     plot(x, models = names(x$fi),
          k=NULL, s=NULL, sp=NULL, xval = c("square", "inverse","sigma"),
          yval = c("psi", "zvalue", "pvalue"), xlab = NULL,
          ylab = NULL,log.xy = "", xlim = NULL, ylim = NULL,
          add = F, length.x = 300, main=NULL,
          col =1:6, lty = 1:5, lwd = par("lwd"), ex.pch=2:7,
          pch = 1, cex = 1, pt.col = col[1],pt.lwd = lwd[1],
          legend.x = NULL, inset = 0.1, ...)

     ## S3 method for class 'summary.scaleboot':
     plot(x, models="best",
          k=x$parex$k,s=x$parex$s,sp=x$parex$sp, ...)

     ## S3 method for class 'scalebootv':
     plot(x,models=attr(x,"models"),...)

     ## S3 method for class 'summary.scalebootv':
     plot(x, models="best",...)

     ## S3 method for class 'scaleboot':
     lines(x,z,models=names(x$fi), k=NULL,s=NULL,sp=NULL,
           length.x=z$length.x, col=z$col,lty=z$lty,lwd=z$lwd,... )

     sblegend(x="topright",y=NULL,z,inset=0.1,...)

_A_r_g_u_m_e_n_t_s:

       x: an object used to select a method.  For 'sblegend', 'x' is a
          numeric or character such as "lefttop" or "righttop", which
          is passed to 'legend'. 

  models: character vector of model names. If numeric,
          'names(object$fi)[models]' is used for each '"scaleboot"'
          object.

       k: k for extrapolation.

       s: s for extrapolation.

      sp: sp for extrapolation.

    xval: specifies x-axis. "square" for sigma^2, "inverse" for
          1/sigma, "sigma" for sigma.

    yval: specifies y-axis. "zvalue" for psi(sigma^2|beta)/sigma or
          'qnorm(1-bp[i])', "pvalue" for 1-Phi(psi(sigma^2|beta)/sigma)
          or 'bp[i]', "psi" for psi(sigma^2|beta) or
          'sqrt(sa[i])*qnorm(1-bp[i])'.

    xlab: label for x-axis.

    ylab: label for y-axis.

  log.xy: character to specify log-scale. "", "x", "y", or "xy".

    xlim: range for x-axis.

    ylim: range for y-axis.

     add: logical for adding another plot.

length.x: the number of segments to draw curves.

    main: for title.

     col: color for model curves.

     lty: lty for model curves.

     lwd: lwd for model curves.

  ex.pch: pch for extrapolation.

     pch: pch for bp points.

     cex: cex for bp points.

  pt.col: col for bp points.

  pt.lwd: lwd for bp points.

legend.x: passed to sblegend as the first argument.

     ...: further arguments passed to or from other methods.

       z: output from previous 'plot.scaleboot'.

       y: numeric passed to 'legend'.

   inset: inset distance from the margins, which is passed to 'legend'.

_D_e_t_a_i_l_s:

     The 'plot' method plots bootstrap probabilities and calls the
     'lines' method, which draws fitted curves for models.

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

     Hidetoshi Shimodaira

_S_e_e _A_l_s_o:

     'sbfit'.

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

     data(mam15)
     ## a single plot
     a <- mam15.relltest[["t4"]] # an object of class "scaleboot"
     plot(a,legend="topleft") # x=sigma^2, y=psi
     plot(a,xval="inverse",yval="zvalue",
          legend="topleft") # x=1/sigma, y=z-value
     plot(a,xval="sigma",log="x",yval="pvalue",
          legend="topleft") # x=log(sigma), y=probability
     ## multiple plots
     b <- mam15.relltest[1:15] # an object of class "scalebootv"
     plot(b) # x=sigma^2, y=psi

