trun.q {gamlss.tr} | R Documentation |
Creates a function to produce the inverse of a truncated cumulative density function generated from a current GAMLSS family distribution
trun.q(par, family = "NO", type = c("left", "right", "both"), ...)
par |
a vector with one (for left or right truncation) or two elements for both |
family |
a gamlss.family object, which is used to define the distribution and the link functions of the various parameters.
The distribution families supported by gamlss() can be found in gamlss.family .
Functions such as BI() (binomial) produce a family object. |
type |
whether left , right or in both sides truncation is required, (left is the default) |
... |
for extra arguments |
Returns a q family function
Mikis Stasinopoulos d.stasinopoulos@londonmet.ac.uk and Bob Rigby r.rigby@londonmet.ac.uk
Rigby, R. A. and Stasinopoulos D. M. (2005). Generalized additive models for location, scale and shape,(with discussion), Appl. Statist., 54, part 3, pp 507-554.
Stasinopoulos D. M., Rigby R.A. and Akantziliotou C. (2003) Instructions on how to use the GAMLSS package in R. Accompanying documentation in the current GAMLSS help files, (see also http://www.gamlss.com/).
trun.d
, trun.q
, trun.r
, gen.trun
# continuous #---------------------------------------------------------------------------------------- # left test1<-trun.q(par=c(0), family="TF", type="left") test1(.6) qTF(pTF(0)+0.6*(1-pTF(0))) #---------------------------------------------------------------------------------------- # right test2 <- trun.q(par=c(10), family="BCT", type="right") test2(.6) qBCT(0.6*pBCT(10)) #---------------------------------------------------------------------------------------- # both test3<-trun.q(par=c(-3,3), family="TF", type="both") test3(.6) qTF(0.6*(pTF(3)-pTF(-3))+pTF(-3)) #---------------------------------------------------------------------------------------- #---------------------------------------------------------------------------------------- # FOR DISCRETE DISTRIBUTIONS # trucated q function # left test4<-trun.q(par=c(0), family="PO", type="left") test4(.6) qPO(pPO(0)+0.6*(1-pPO(0))) #---------------------------------------------------------------------------------------- # right test5 <- trun.q(par=c(10), family="NBI", type="right") test5(.6) qNBI(0.6*pNBI(10)) test5(.6, mu=10, sigma=2) qNBI(0.6*pNBI(10, mu=10, sigma=2), mu=10, sigma=2) #---------------------------------------------------------------------------------------- # both test6<-trun.q(par=c(0,10), family="NBI", type="both") test6(.6) qNBI(0.6*(pNBI(10)-pNBI(0))+pNBI(0)) #----------------------------------------------------------------------------------------