doldlogspline           package:polspline           R Documentation

_L_o_g_s_p_l_i_n_e _D_e_n_s_i_t_y _E_s_t_i_m_a_t_i_o_n - _1_9_9_2 _v_e_r_s_i_o_n

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

     Probability density function ('doldlogspline'), distribution
     function ('poldlogspline'), quantiles  ('qoldlogspline'), and 
     random samples ('roldlogspline') from a logspline density that was
     fitted using the 1992 knot deletion algorithm ('oldlogspline'). 
     The 1997 algorithm using knot deletion and addition is available
     using the 'logspline' function.

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

     doldlogspline(q, fit) 
     poldlogspline(q, fit) 
     qoldlogspline(p, fit) 
     roldlogspline(n, fit) 

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

       q: vector of quantiles. Missing values (NAs) are allowed.  

       p: vector of probabilities. Missing values (NAs) are allowed.  

       n: sample size. If 'length(n)' is larger than 1, then 
          'length(n)' random values are returned. 

     fit: 'oldlogspline' object, typically the result of
          'oldlogspline'. 

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

     Elements of  'q' or 'p' that are missing will cause the 
     corresponding elements of the result to be missing.

_V_a_l_u_e:

     Densities ('doldlogspline'), probabilities ('poldlogspline'),
     quantiles ('qoldlogspline'),  or a random sample ('roldlogspline')
     from an 'oldlogspline' density that was fitted using knot
     deletion.

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

     Charles Kooperberg clk@fhcrc.org.

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

     Charles Kooperberg and Charles J. Stone.  Logspline density
     estimation for censored data (1992). _Journal of Computational and
     Graphical Statistics_, *1*, 301-328.

     Charles J. Stone, Mark Hansen, Charles Kooperberg, and Young K.
     Truong. The use of polynomial splines and their tensor products in
     extended linear modeling (with discussion) (1997).  _Annals of
     Statistics_, *25*, 1371-1470.

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

     'logspline', 'oldlogspline', 'plot.oldlogspline',
     'summary.oldlogspline'

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

     x <- rnorm(100)
     fit <- oldlogspline(x)
     qq <- qoldlogspline((1:99)/100, fit)
     plot(qnorm((1:99)/100), qq)                  # qq plot of the fitted density
     pp <- poldlogspline((-250:250)/100, fit)
     plot((-250:250)/100, pp, type = "l")
     lines((-250:250)/100, pnorm((-250:250)/100)) # asses the fit of the distribution
     dd <- doldlogspline((-250:250)/100, fit)
     plot((-250:250)/100, dd, type = "l")
     lines((-250:250)/100, dnorm((-250:250)/100)) # asses the fit of the density
     rr <- roldlogspline(100, fit)                # random sample from fit

