Displaying 1 result from an estimated 1 matches for "pcwise".
Did you mean:
cwise
2011 Mar 18
1
PREDICTIONS from a PIECEWISE LINEAR (mixed) MODEL: THEY AIN'T LINEAR BETWEEN BREAK POINTS!!
...rt,sex,age,
nadir,pharmac2,hbs,hcv,resp2,
hivbl,switch)
# I run a model specifing the function
(t-t*)+ = 0 if t<=t*
AND
(t-t*)+ = t if t>t*
with the following syntax, emulating MJ Crawley (The R book 2007, Wiley),
# rt+rt*(rt>15)
I run the model, is a mixed effect one.......
PCWISE<-lme(chol~rt+rt*(rt>15)
+sex+age+cd4+cd4:rt+
cd4:rt*(rt>15)+nadir+hivbl+
pharmac2:rt+resp2:rt+
resp2:rt*(rt>15)+hbs+hcv+
switch,data=frame3,
random= ~rt|id,
na.action=na.omit)
# I prepare a fictious data.frame for the predictions , in fact lme predict
values inside the clusters
, id...