f_plot_linERR_ef: plot likelihood function from ef

In fbr600/rERR: Excess Relative Risk Models

Description

plot the partial log likelihood function in the case of one dimension in the linear part

Usage

f_plot_linERR_ef(object, formula, data, id_name, dose_name, time_name,
  covars_names, lag, exclusion_done = F)

Arguments

object	an rERR class object
formula	Surv(entry_time,exit_time,outcome)~loglin(loglin_var1,..,loglin_varn)+ lin(lin_var1,..,lin_varm)+strata(strat_var1,...strat_varp)
data	data set returned from f_to_model_data
id_name	name of variable containing the names of subjects
dose_name	name of variable containing the doses at each time
time_name	name of the time variable
covars_names	a character vector with the names of the variables used as covariates in the fomrula (adjustments and stratification)
lag	latency period
exclusion_done	a logical indicating wheather the exclusion is already done or not

Value

rERR object with the estimation

Examples

# set the formulas for the models
formula1  <- Surv(entry_age,exit_age,outcome) ~ lin(dose_cum) + strata(sex)

# fit the model
fit1 <- f_fit_linERR_ef(formula1,data=cohort_ef,id_name="id",dose_name="dose",
                        time_name="age",covars_names=c("sex"),lag=2,exclusion_done=TRUE)
                        
# plot the partial loglikelihood function
f_plot_linERR_ef(fit1,formula1,data=cohort_ef,id_name="id",dose_name="dose",
                        time_name="age",covars_names=c("sex"),lag=2,exclusion_done=TRUE)

fbr600/rERR documentation built on May 7, 2019, 9:36 a.m.