# NOT RUN {
# identify predictor variables used in mixture
mixvars=names(lwqs_data)[5:9]
# run model. Note for example run-time only 1 bootstrap (b=1) is used. Set b to be >50
model=lwqs(data=lwqs_data,
timevar="time",
wqs_parms=list(formula=out ~ wqs,
data = lwqs_data,
mix_name=mixvars,
b1_constr = TRUE,
b1_pos=FALSE,
b = 5,
q = 5,
validation = 0,
family = "gaussian",
seed = 1),
outcome="out",
ID="ID")
# use extract_mixture to access time-varying wqs index
mixtime=extract_mixture(model)
# }
Run the code above in your browser using DataLab