gistfile1.txt

##  QOL original scales
str(qol_combined)
qol_combined2<-qol_combined
qol_combined2<-ddply(qol_combined2,.(Patientnr),.drop=FALSE,transform,baselinechange=value-value[1])
qol_combined3<-qol_combined
qol_combined3<-ddply(qol_combined3,.(Patientnr),.drop=FALSE,transform,change=value-lag(value))
qol_combined$qol_baselinechange<-qol_combined2$baselinechange
qol_combined$qol_change<-qol_combined3$change
qol_combined$qol_baselinechange[qol_combined$time==0]<-NA
qol_combined$qol_change[qol_combined$time==0]<-NA
str(qol_combined)
# Check visually if worked
g1<-ggplot(qol_combined) + geom_bar(aes(as.factor(time), qol_baselinechange), position = "dodge", stat = "summary", fun.y = "mean") +facet_wrap(~variable,scales = "free") + theme_bw()
g2<-ggplot(qol_combined) + geom_bar(aes(as.factor(time), qol_change), position = "dodge", stat = "summary", fun.y = "mean") +facet_wrap(~variable, scales = "free") + theme_bw()
grid.arrange(g1,g2,ncol=2, top="QOL Changes ")