rERR (version 0.1)

f_plot_linERR_wf: plot likelihood function from wf

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

# NOT RUN {
# 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)
# }