home *** CD-ROM | disk | FTP | other *** search
/ Big Green CD 8 / BGCD_8_Dev.iso / NEXTSTEP / UNIX / Educational / R-0.49-MI / R-0.49-I / help / bootstrap / abcpar < prev    next >
Encoding:
Text File  |  1997-09-13  |  2.1 KB  |  70 lines
  1.     
  2.     _P_a_r_a_m_e_t_r_i_c  _A_B_C  _c_o_n_f_i_d_e_n_c_e _l_i_m_i_t_s
  3.     
  4.          abcpar(x, tt, S, etahat, mu, n=rep(1,length(x)),lambda=0.001,
  5.           alpha=c(0.025, 0.05, 0.1, 0.16))
  6.     
  7.     _A_r_g_u_m_e_n_t_s:
  8.     
  9.                 x: vector of data
  10.     
  11.                tt: function of expectation parameter mu defining
  12.                    the parameter of interest
  13.     
  14.                 S: maximum likelihood estimate of the  covari-
  15.                    ance matrix of x
  16.     
  17.            etahat: maximum likelihood estimate of the natural
  18.                    parameter eta
  19.     
  20.                mu: function giving expectation of x in terms of
  21.                    eta
  22.     
  23.                 n: optional argument containing denominators for
  24.                    binomial (vector of length length(x))
  25.     
  26.            lambda: optional argument specifying step size for
  27.                    finite difference calculation
  28.     
  29.             alpha: optional argument specifying confidence lev-
  30.                    els desired
  31.     
  32.     _V_a_l_u_e_s:
  33.     
  34.          list with the following components
  35.     
  36.         call: the call to abcpar
  37.     
  38.       limits: The nominal confidence level, ABC point, quadratic
  39.               ABC point, and standard normal point.
  40.     
  41.        stats: list consisting of  observed value of tt,
  42.               estimated standard error and estimated bias
  43.     
  44.     constants: list consisting of a=acceleration constant,
  45.               z0=bias adjustment, cq=curvature component
  46.     
  47.     _R_e_f_e_r_e_n_c_e_s:
  48.     
  49.          Efron, B, and DiCiccio, T. (1992) More accurate confi-
  50.          dence intervals in exponential families. Bimometrika
  51.          79, pages 231-245.
  52.     
  53.          Efron, B. and Tibshirani, R. (1993) An Introduction to
  54.          the Bootstrap.  Chapman and Hall, New York, London.
  55.     
  56.     _E_x_a_m_p_l_e_s:
  57.     
  58.          # binomial
  59.          # x is a p-vector of successes, n is a p-vector of
  60.          #  number of trials
  61.          S <- matrix(0,nrow=p,ncol=p)
  62.          S[row(S)==col(S)] <- x*(1-x/n)
  63.          mu <- function(eta,n)n/(1+exp(eta))
  64.          etahat <- log(x/(n-x))
  65.          #suppose p=2 and we are interested in mu2-mu1
  66.          tt <- function(mu)mu[2]-mu[1]
  67.          x <- c(2,4); n <- c(12,12)
  68.          a <- abcpar(x, tt, S, etahat,n)
  69.     
  70.