R/plot.nonparasccs.R

In SCCS: The Self-Controlled Case Series Method

plot.nonparasccs <- function(x, type="l", conf.int=0.95, ...) {
  
   fit <- x
   z <- qnorm((1 + conf.int)/2, 0, 1)
   
   rho <- fit$exposure
   #lci <- fit$exposure - z*fit$se
   #uci <- fit$exposure + z*fit$se
   
   lci <- rho*exp(-z*fit$se*(1/rho))
   uci <- rho*exp(z*fit$se*(1/rho))
   
   
   par(mfrow=c(1,2))
 
   #rho <- fit$exposure
   #lci <- rho*exp(-z*fit$se*(1/rho))
   #uci <- rho*exp(z*fit$se*(1/rho))
   
  # Age related relative incidence function
  plot(fit$ageaxis, fit$age, type=type, ylab = "relative incidence", xlab = "age (days)")
  
  # Exposure related relative incidence function
  plot(fit$timesinceexpo, fit$exposure, type=type, ylab= "relative incidence", xlab= "days since start of risk period", ylim = c(0, (max(uci)+2)))
  #plot(fit$timesinceexpo, fit$exposure, type=type, ylab= "Relative incidence", xlab= "Days since start of risk period")
  #plot(fit$timesbiceexpo, fit$exposure, type=type, ylab= "Relative incidence", xlab= "Days since start of risk period", ylim = c((min(fit$lci)-0.5), (max(fit$uci)+2)))
  lines(fit$timesinceexpo, lci, lty=2)
  lines(fit$timesinceexpo, uci, lty=2)
}