rhosp (version 1.10)

estimParam: compute the parametric 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 mu is the estimator of maximum likelihood of the random variable \(Zi / Xi=x\). the risk constant R is given by the equation (star) \(R=mu\).

Usage

estimParam(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

the best estimation of mu : the value of mu which maximize the loglikelihood

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)

Examples

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


# }