# NOT RUN {
# identify predictor variables used in mixture
mixvars=names(lwqs_data)[5:9]
# run model
model=lwqs(data=lwqs_data,
timevar="time",
wqs_parms=list(formula=out ~ wqs,
data = lwqs_data,
mix_name=mixvars,
b1_constr = TRUE,
b1_pos=TRUE,
b = 5,
q = 5,
validation = 0,
family = "gaussian",
seed = 1),
outcome="out",
ID="ID")
# use extract_weights to access time-varying predictor weights
timeweights=extract_weights(model)
# }
Run the code above in your browser using DataLab