gistfile1.txt

## Interference
Interference_BPI2<-Interference_BPI
Interference_BPI2<-ddply(Interference_BPI2,.(Patientnr),.drop=FALSE,transform,baselinechange=Interference-Interference[1])
Interference_BPI3<-Interference_BPI
Interference_BPI3<-ddply(Interference_BPI3,.(Patientnr),.drop=FALSE,transform,change=Interference-lag(Interference))
Interference_BPI$Interference_baselinechange<-Interference_BPI2$baselinechange
Interference_BPI$Interference_change<-Interference_BPI3$change
Interference_BPI$Interference_baselinechange[Interference_BPI$time==0]<-NA
Interference_BPI$Interference_change[Interference_BPI$time==0]<-NA
Interference_BPI
# Check visually if worked
g2<-ggplot(Interference_BPI,aes(x=as.factor(time),  y=Interference, group=as.factor(Patientnr))) + geom_line() + geom_point() + facet_wrap(~Patientnr) + theme_bw()
g1<-ggplot(Interference_BPI,aes(x=as.factor(time),  y=as.factor(Patientnr),fill=Interference)) + geom_tile(colour="white",size=0.25, na.rm=TRUE)+
  labs(x="",y="")+scale_y_discrete(expand=c(0,0))+scale_fill_viridis(discrete = FALSE, option="magma", direction = -1)+  theme_grey(base_size=10)+marc_discrete
g3<-ggplot(Interference_BPI) + 
  geom_bar(aes(as.factor(time), Interference_baselinechange), 
           position = "dodge", stat = "summary", fun.y = "mean") + facet_wrap(~Patientnr)+ theme_bw()
g4<-ggplot(Interference_BPI) + 
  geom_bar(aes(as.factor(time), Interference_change), 
           position = "dodge", stat = "summary", fun.y = "mean") + facet_wrap(~Patientnr)+ theme_bw()
grid.arrange(g1,g2,g3,g4,ncol=4, top="BPI Interference Changes ")

## Relief 
BPI_6 <- bpi_combined%>%
  dplyr::select(Patientnr,variable,value,time)%>%
  dplyr::filter(variable == "BPI_6_")%>%
  dplyr::arrange(Patientnr,time)
BPI_62<-BPI_6
BPI_62<-ddply(BPI_62,.(Patientnr),.drop=FALSE,transform,baselinechange=value-value[1])
BPI_63<-BPI_6
BPI_63<-ddply(BPI_63,.(Patientnr),.drop=FALSE,transform,change=value-lag(value))
BPI_6$relief_baselinechange<-BPI_62$baselinechange
BPI_6$relief_change<-BPI_63$change
BPI_6
# Check visually if worked
g1<-ggplot(pan_bpi_6_melt, aes(as.factor(variable), as.factor(Patientnr), fill=value))+ 
  geom_tile(colour="white",size=0.25, na.rm=TRUE)+
  labs(x="",y="")+scale_y_discrete(expand=c(0,0))+
  scale_fill_viridis(discrete = FALSE, option="magma")+
  theme_grey(base_size=10)+marc_discrete
g2<-ggplot(BPI_6, aes(x=as.factor(time),  y=value, group=as.factor(Patientnr))) + geom_line() + geom_point() + facet_wrap(~Patientnr) + theme_bw()
g3<-ggplot(BPI_6) + geom_bar(aes(as.factor(time), relief_baselinechange), 
                             position = "dodge", stat = "summary", fun.y = "mean") + facet_wrap(~Patientnr)+ theme_bw()
g4<-ggplot(BPI_6) + geom_bar(aes(as.factor(time), relief_change), 
                             position = "dodge", stat = "summary", fun.y = "mean") + facet_wrap(~Patientnr)+ theme_bw()
grid.arrange(g1,g2,g3,g4,ncol=4, top="BPI Relief Changes")