estimParam {rhosp}R Documentation

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 estimatoçr 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)

Author(s)

Christophe Dutang and Julie Barthes

Examples


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


[Package rhosp version 1.04 Index]