MJacobs1985

pan_qol_YN_melt<-reshape2::melt(pan_qol_YN, id=c("Patientnr"));dim(pan_qol_YN_melt);head(pan_qol_YN_melt)  # 30*9=270
ggplot(pan_qol_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

MJacobs1985

pan_qol_YN<-pancreas%>%dplyr::select(Patientnr, QOL_YN_WKB, QOL_YN_WK_1, QOL_YN_WK_2, QOL_YN_WK_4, QOL_YN_WK_5, QOL_YN_WK_7, QOL_YN_WK_11, QOL_YN_WK_19, QOL_YN_WK_31);dim(pan_qol_YN);head(pan_qol_YN)
pan_qol_1 <-pancreas%>%dplyr::select(Patientnr,  QOL_1_WKB,  QOL_1_WK_1,  QOL_1_WK_2,  QOL_1_WK_4,  QOL_1_WK_5,  QOL_1_WK_7,  QOL_1_WK_11,  QOL_1_WK_19,  QOL_1_WK_31);dim(pan_qol_1)
pan_qol_2 <-pancreas%>%dplyr::select(Patientnr,  QOL_2_WKB,  QOL_2_WK_1,  QOL_2_WK_2,  QOL_2_WK_4,  QOL_2_WK_5,  QOL_2_WK_7,  QOL_2_WK_11,  QOL_2_WK_19,  QOL_2_WK_31);dim(pan_qol_2)
pan_qol_3 <-pancreas%>%dplyr::select(Patientnr,  QOL_3_WKB,  QOL_3_WK_1,  QOL_3_WK_2,  QOL_3_WK_4,  QOL_3_WK_5,  QOL_3_WK_7,  QOL_3_WK_11,  QOL_3_WK_19,  QOL_3_WK_31);dim(pan_qol_3)
pan_qol_4 <-pancreas%>%dplyr::select(Patientnr,  QOL_4_WKB,  QOL_4_WK_1,  QOL_4_WK_2,  QOL_4_WK_4,  QOL_4_WK_5,  QOL_4_WK_7,  QOL_4_WK_11,  QOL_4_WK_19,  QOL_4_WK_31);dim(pan_qol_4)
pan_qol_5 <-pancreas%>%dplyr::select(Patientnr,  QOL_5_WKB,  QOL_5_WK_1,  QOL_5_WK_2,  QOL_5_WK_

MJacobs1985

pancreas$Performance_status_at_consultation<-factor(pancreas$Performance_status_at_consultation)
pancreas$Performance_status_at_consultation
pancreas$Gender<-factor(pancreas$Gender,
                        levels = c("1", "2"),
                        labels = c("Female", "Male"))
pancreas$Gender
pancreas$Disease_stage_at_consultation<-factor(pancreas$Disease_stage_at_consultation,
levels = c("1", "2"),
labels = c("Metastasized", "Non-Metastasized"))
pancreas$Disease_stage_at_consultation

MJacobs1985

pancreas <- read_delim("pancreas.csv", delim = ";", 
escape_double = FALSE, trim_ws = TRUE)
dim(pancreas)
View(pancreas)
str(pancreas)
skim(pancreas) 
DataExplorer::plot_missing(pancreas)

MJacobs1985

rm(list = ls()) 
#### LIBRARIES ####
library(readr)
library(skimr)
library(DataExplorer)
library(tidyverse)
library(reshape2)
library(psych)
library(paletteer)
library(viridis)

MJacobs1985

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_7

MJacobs1985

#### IMPORT LIBRARIES ####
rm(list = ls())
library(brms)
library(dplyr) 
library(ggplot2)
library(rethinking)
library(dmetar)
library(meta)
library(tidybayes)
library(ggridges)

MJacobs1985

library(dplyr)
library(ggplot2)
library(sjPlot)
library(sjmisc)

p_grid                   <- seq(from=0 , to=1   , length.out=100)
prior                    <- dbinom(60, size=100 , prob=p_grid)
likelihood               <- dbinom(50, size=100 , prob=p_grid)
unstd.posterior          <- likelihood * prior
posterior                <- unstd.posterior / sum(unstd.posterior)

MJacobs1985

rm(list = ls())
library(brms)
library(dplyr)
library(ggplot2)
library(grid)
library(gridExtra)
library(bayesplot)
zinb <- read.csv("http://stats.idre.ucla.edu/stat/data/fish.csv")
zinb$camper <- factor(zinb$camper, labels = c("no", "yes"))
head(zinb)

MJacobs1985

rm(list = ls())
library(ggplot)
library(dplyr)
library(readr)
library(ggthemes)
library(lubridate)
library(zoo)
library(gganimate)
library(cowplot)
library(gridExtra)