Nothing
## ----load,message=FALSE-------------------------------------------------------
library(dejaVu)
## -----------------------------------------------------------------------------
summary(simData)
## -----------------------------------------------------------------------------
#analyse observed data
library(MASS)
obsMod <- glm.nb(y~z+offset(log(fupTime)), data=simData)
summary(obsMod)
## -----------------------------------------------------------------------------
covar.df <- data.frame(id=1:dim(simData)[1], z=simData$z)
dejaData <- MakeDejaData(covar.df, arm="z", Id="id")
## -----------------------------------------------------------------------------
eventTimes <- expandEventCount(count = simData$y, time = simData$fupTime)
## -----------------------------------------------------------------------------
obsData <-
ImportSim(
dejaData,
event.times = eventTimes,
status = "dropout",
study.time = 1,
censored.time = simData$fupTime,
allow.beyond.study = T
)
## -----------------------------------------------------------------------------
obsFit <- Simfit(obsData,equal.dispersion=TRUE)
## -----------------------------------------------------------------------------
imputed.data.sets <- Impute(fit = obsFit,impute.mechanism = weighted_j2r(trt.weight=0),
N=10)
## -----------------------------------------------------------------------------
fitted <- Simfit(imputed.data.sets,family="negbin")
summary(fitted)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.