estimNonParam: compute the parametric estimation

In rhosp: Side Effect Risks in Hospital : Simulation and Estimation

Description

estimParam compute the parametric estimation of the first model on the input file containing the stay duration Xi and the side effect reporting Zi. this estimator of R (and not mu) is the smallest root of the equation (*) : P(Zi=0)E[exp(-RXi)/Zi=0]=1

Usage

estimNonParam(fileName, toplot = TRUE, header = TRUE, ks = FALSE, DV = FALSE)

Arguments

fileName	the file in which the simulated data will be stored
toplot	a logical variable to plot the result of this estimation
header	a logical for : has the input file an header
ks	a logical for : do you want the Kolmogorv Smirnov test
DV	a logical for : do you want to calculate the De Vielder estimation

Value

a list of the following components

CR	the CR risk constant calculated with the value of R
R	the risk constant estimated by this estimation
T	the vector of observations of the random variable T
lambdaHat	the best estimation of lambda : the ESBVM of lambda
muHat	not available with this estimation
lambdaEmp	the estimation of the lambda of the De Vielder method (only available if DV=TRUE and if there is no problem with De Vielder method)
muEmp	the estimation of the mu of the De Vielder method (only available if DV=TRUE)

Author(s)

Christophe Dutang and Julie Barthes

Examples

#res<-estimDV("data.rda",TRUE,TRUE)

rhosp documentation built on May 1, 2019, 11:27 p.m.