semisccs | R Documentation |

The function fits the semiparametric self-controlled case series method where the age effect is left unspecified, as published in Farrington and Whitaker (2006).

semisccs(formula, indiv, astart, aend, aevent, adrug, aedrug, expogrp = list(), washout = list(), sameexpopar = list(), dataformat="stack", data)

`formula` |
model formula. The dependent variable should always be "event" e.g. event ~ itp. There is no need to specify age in the model formula. |

`indiv` |
a vector of individual identifiers of cases. |

`astart` |
a vector of ages at which the observation periods start. |

`aend` |
a vector of ages at end of observation periods. |

`aevent` |
a vector of ages at event, an individual can experience multiple events. |

`adrug` |
a list of vectors of ages at start of exposures or a list of matrices if the exposures have multiple episodes ( |

`aedrug` |
a list of vectors of ages at which exposure-related risk ends or a list of matrices if there are multiple episodes (repeat exposures in different columns) of the same exposure type. The dimension of each item of |

`expogrp` |
list of vectors of days to the start of exposure-related risk, counted from |

`washout` |
list of vectors with days to start of washout periods counted from |

`sameexpopar` |
a vector of logical values. If TRUE (the default) no dose effect is assumed: the same exposure parameters are used for multiple doses of the same exposure type, presented in |

`dataformat` |
the way the input data are assembled. It accepts "multi" or "stack" (the default), where "multi" refers to a data assembled with one row representing one event and "stack" refers to a data frame where repeated exposures of the same type are stack in one column. In the "multi" dataformat different episodes of the same type are recorded as separate columns in the dataframe. |

`data` |
a data frame containing the input data. The data should be in 'stack' or 'multi' (see |

In the standard SCCS method both age and exposure effects are modelled using step functions. However, mis-specification of age groups in the standard SCCS may lead to bias in the exposure related relative incidence estimates. In the semiparametric SCCS no age groups are pre-specified. A parameter for each day an event occurred is fitted, which means that this method is only suitable for small to medium sized data sets. An alternative for large data sets is provided by `smoothagesccs`

.

The function returns age and exposure related relative incidence estimates along with 95% confidence limits.

Yonas Ghebremichael-Weldeselassie, Heather Whitaker, Paddy Farrington.

Farrington, C. P., Whitaker, H. J. (2006). Semiparametric analysis of case series data. Applied Statistics, 55(5): 553–594.

Farrington, P., Whitaker, H., and Ghebremichael-Weldeselassie, Y. (2018). Self-controlled Case Series Studies: A modelling Guide with R. Boca Raton: Chapman & Hall/CRC Press.

`standardsccs`

, `smoothagesccs`

, `smoothexposccs`

# Example 1 # Semiparametric model for the ITP and MMR vaccine data itp.mod1 <- semisccs(event~mmr, indiv=case, astart=sta, aend=end, aevent=itp, adrug=mmr, aedrug=mmr+42, expogrp=c(0,15,29), data=itpdat) itp.mod1 # Example 2 # Data on itp and mmr vaccine # Sex and mmr interaction included itp.mod2 <- semisccs(event~factor(sex)*mmr, indiv=case, astart=sta, aend=end, aevent=itp, adrug=mmr, aedrug=mmr+42, expogrp=c(0,15,29), data=itpdat) itp.mod2

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.