f_plot_linERR_wf: plot likelihood function from wf

In 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_wf(object, formula, data, id_name, doses, times, covars, 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
doses	sub data set of grouped doses
times	sub data set of times relatives to doses
covars	sub data set of the covars that will be involved in the model (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(AgeAtEntry,age_at_event,outcome) ~ lin(dose_cum) + strata(sex)

# fit the model
fit1 <- f_fit_linERR_wf(formula1,data=cohort_wf,id_name="id",doses=cohort_wf[,45:79],
                        times=cohort_wf[,10:44],covars=cohort_wf[,c("sex","country")],
                        lag=2,exclusion_done = FALSE)
                        
# plot the partial loglikelihood function
f_plot_linERR_wf(fit1,formula1,data=cohort_wf,id_name="id",doses=cohort_wf[,45:79],
                 times=cohort_wf[,10:44],covars=cohort_wf[,c("sex","country")],
                 lag=2,exclusion_done = FALSE)

rERR documentation built on May 2, 2019, 11:10 a.m.