gistfile1.txt

pan_bpi_YN <-pancreas%>%dplyr::select(Patientnr,BPI_YN_WKB,BPI_YN_WK_1,BPI_YN_WK_2,BPI_YN_WK_4,BPI_YN_WK_5,BPI_YN_WK_7,BPI_YN_WK_11,BPI_YN_WK_19,BPI_YN_WK_31);dim(pan_bpi_YN);head(pan_bpi_YN)
pan_bpi_2  <-pancreas%>%dplyr::select(Patientnr,BPI_2_WKB,BPI_2_WK_1,BPI_2_WK_2,BPI_2_WK_4,BPI_2_WK_5,BPI_2_WK_7,BPI_2_WK_11,BPI_2_WK_19,BPI_2_WK_31);dim(pan_bpi_2)
pan_bpi_3  <-pancreas%>%dplyr::select(Patientnr,BPI_3_WKB,BPI_3_WK_1,BPI_3_WK_2,BPI_3_WK_4,BPI_3_WK_5,BPI_3_WK_7,BPI_3_WK_11,BPI_3_WK_19,BPI_3_WK_31);dim(pan_bpi_3)
pan_bpi_4  <-pancreas%>%dplyr::select(Patientnr,BPI_4_WKB,BPI_4_WK_1,BPI_4_WK_2,BPI_4_WK_4,BPI_4_WK_5,BPI_4_WK_7,BPI_4_WK_11,BPI_4_WK_19,BPI_4_WK_31);dim(pan_bpi_4)
pan_bpi_6  <-pancreas%>%dplyr::select(Patientnr,BPI_6_WKB,BPI_6_WK_1,BPI_6_WK_2,BPI_6_WK_4,BPI_6_WK_5,BPI_6_WK_7,BPI_6_WK_11,BPI_6_WK_19,BPI_6_WK_31);dim(pan_bpi_6)
pan_bpi_7A <-pancreas%>%dplyr::select(Patientnr,BPI_7A_WKB,BPI_7A_WK_1,BPI_7A_WK_2,BPI_7A_WK_4,BPI_7A_WK_5,BPI_7A_WK_7,BPI_7A_WK_11,BPI_7A_WK_19,BPI_7A_WK_31);dim(pan_bpi_7A)
pan_bpi_7B <-pancreas%>%dplyr::select(Patientnr,BPI_7B_WKB,BPI_7B_WK_1,BPI_7B_WK_2,BPI_7B_WK_4,BPI_7B_WK_5,BPI_7B_WK_7,BPI_7B_WK_11,BPI_7B_WK_19,BPI_7B_WK_31);dim(pan_bpi_7B)
pan_bpi_7C <-pancreas%>%dplyr::select(Patientnr,BPI_7C_WKB,BPI_7C_WK_1,BPI_7C_WK_2,BPI_7C_WK_4,BPI_7C_WK_5,BPI_7C_WK_7,BPI_7C_WK_11,BPI_7C_WK_19,BPI_7C_WK_31);dim(pan_bpi_7C)
pan_bpi_7D <-pancreas%>%dplyr::select(Patientnr,BPI_7D_WKB,BPI_7D_WK_1,BPI_7D_WK_2,BPI_7D_WK_4,BPI_7D_WK_5,BPI_7D_WK_7,BPI_7D_WK_11,BPI_7D_WK_19,BPI_7D_WK_31);dim(pan_bpi_7D)
pan_bpi_7E <-pancreas%>%dplyr::select(Patientnr,BPI_7E_WKB,BPI_7E_WK_1,BPI_7E_WK_2,BPI_7E_WK_4,BPI_7E_WK_5,BPI_7E_WK_7,BPI_7E_WK_11,BPI_7E_WK_19,BPI_7E_WK_31);dim(pan_bpi_7E)
pan_bpi_7F <-pancreas%>%dplyr::select(Patientnr,BPI_7F_WKB,BPI_7F_WK_1,BPI_7F_WK_2,BPI_7F_WK_4,BPI_7F_WK_5,BPI_7F_WK_7,BPI_7F_WK_11,BPI_7F_WK_19,BPI_7F_WK_31);dim(pan_bpi_7F)
pan_bpi_7G <-pancreas%>%dplyr::select(Patientnr,BPI_7G_WKB,BPI_7G_WK_1,BPI_7G_WK_2,BPI_7G_WK_4,BPI_7G_WK_5,BPI_7G_WK_7,BPI_7G_WK_11,BPI_7G_WK_19,BPI_7G_WK_31);dim(pan_bpi_7G)
pan_bpi_YN_melt<-reshape2::melt(pan_bpi_YN, id=c("Patientnr"));dim(pan_bpi_YN_melt);head(pan_bpi_YN_melt)

bpi_gYN<-ggplot(pan_bpi_YN_melt, aes(as.factor(variable), as.factor(Patientnr), fill=as.factor(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 = TRUE, option="magma", direction = -1 )+
  theme_grey(base_size=10)+marc_discrete
pan_bpi_2_melt<-reshape2::melt(pan_bpi_2, id=c("Patientnr"));dim(pan_bpi_2_melt)
bpi_g2<-ggplot(pan_bpi_2_melt, aes(as.factor(variable), as.factor(Patientnr), fill=as.factor(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 = TRUE, option="magma", direction = -1 )+
  theme_grey(base_size=10)+marc_discrete
pan_bpi_3_melt<-reshape2::melt(pan_bpi_3, id=c("Patientnr"));dim(pan_bpi_3_melt)
bpi_g3<-ggplot(pan_bpi_3_melt, aes(as.factor(variable), as.factor(Patientnr), fill=as.factor(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 = TRUE, option="magma", direction = -1 )+
  theme_grey(base_size=10)+marc_discrete
pan_bpi_4_melt<-reshape2::melt(pan_bpi_4, id=c("Patientnr"));dim(pan_bpi_4_melt)
bpi_g4<-ggplot(pan_bpi_4_melt, aes(as.factor(variable), as.factor(Patientnr), fill=as.factor(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 = TRUE, option="magma", direction = -1 )+
  theme_grey(base_size=10)+marc_discrete
pan_bpi_6_melt<-reshape2::melt(pan_bpi_6, id=c("Patientnr"));dim(pan_bpi_6_melt)
bpi_g6<-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", direction = -1 )+
  theme_grey(base_size=10)+marc_discrete
pan_bpi_7A_melt<-reshape2::melt(pan_bpi_7A, id=c("Patientnr"));dim(pan_bpi_7A_melt)
bpi_g7A<-ggplot(pan_bpi_7A_melt, aes(as.factor(variable), as.factor(Patientnr), fill=as.factor(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 = TRUE, option="magma", direction = -1 )+
  theme_grey(base_size=10)+marc_discrete
pan_bpi_7B_melt<-reshape2::melt(pan_bpi_7B, id=c("Patientnr"));dim(pan_bpi_7B_melt)
bpi_g7B<-ggplot(pan_bpi_7B_melt, aes(as.factor(variable), as.factor(Patientnr), fill=as.factor(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 = TRUE, option="magma", direction = -1 )+
  theme_grey(base_size=10)+marc_discrete
pan_bpi_7C_melt<-reshape2::melt(pan_bpi_7C, id=c("Patientnr"));dim(pan_bpi_7C_melt)
bpi_g7C<-ggplot(pan_bpi_7C_melt, aes(as.factor(variable), as.factor(Patientnr), fill=as.factor(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 = TRUE, option="magma", direction = -1 )+
  theme_grey(base_size=10)+marc_discrete
pan_bpi_7D_melt<-reshape2::melt(pan_bpi_7D, id=c("Patientnr"));dim(pan_bpi_7D_melt)
bpi_g7D<-ggplot(pan_bpi_7D_melt, aes(as.factor(variable), as.factor(Patientnr), fill=as.factor(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 = TRUE, option="magma", direction = -1 )+
  theme_grey(base_size=10)+marc_discrete
pan_bpi_7E_melt<-reshape2::melt(pan_bpi_7E, id=c("Patientnr"));dim(pan_bpi_7E_melt)
bpi_g7E<-ggplot(pan_bpi_7E_melt, aes(as.factor(variable), as.factor(Patientnr), fill=as.factor(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 = TRUE, option="magma", direction = -1 )+
  theme_grey(base_size=10)+marc_discrete
pan_bpi_7F_melt<-reshape2::melt(pan_bpi_7F, id=c("Patientnr"));dim(pan_bpi_7F_melt)
bpi_g7F<-ggplot(pan_bpi_7F_melt, aes(as.factor(variable), as.factor(Patientnr), fill=as.factor(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 = TRUE, option="magma", direction = -1 )+
  theme_grey(base_size=10)+marc_discrete
pan_bpi_7G_melt<-reshape2::melt(pan_bpi_7G, id=c("Patientnr"));dim(pan_bpi_7G_melt)
bpi_g7G<-ggplot(pan_bpi_7G_melt, aes(as.factor(variable), as.factor(Patientnr), fill=as.factor(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 = TRUE, option="magma", direction = -1 )+
  theme_grey(base_size=10)+marc_discrete
grid.arrange(bpi_gYN,bpi_g2,bpi_g3,bpi_g4,bpi_g6,bpi_g7A,
             bpi_g7B,bpi_g7C,bpi_g7D,bpi_g7E,bpi_g7F,bpi_g7G,ncol=3)