# \donttest{
data("apparelTrans")
# Fit pnbd standard model, no covariates
clv.data.apparel <- clvdata(apparelTrans, time.unit="w",
estimation.split=52, date.format="ymd")
pnbd.apparel <- pnbd(clv.data.apparel)
# summary about model fit
summary(pnbd.apparel)
# Add static covariate data
data("apparelStaticCov")
data.apparel.cov <-
SetStaticCovariates(clv.data.apparel,
data.cov.life = apparelStaticCov,
names.cov.life = "Gender",
data.cov.trans = apparelStaticCov,
names.cov.trans = "Gender",
name.id = "Id")
# fit model with covariates and regualization
pnbd.apparel.cov <- pnbd(data.apparel.cov,
reg.lambdas = c(life=2, trans=4))
# additional summary about covariate parameters
# and used regularization
summary(pnbd.apparel.cov)
# }
Run the code above in your browser using DataLab