nunosempere.github.io/rat/eamentalhealth

EA Mental Health Exploratory Plots

The Code to produce them was:

# EA Mental Health - Clean

# Dependencies:

  # install.packages("tidyr")
  # install.packages("scales")
  # install.packages("tidyverse")
  # install.packages("viridisLite")
  # install.packages("viridis")
  # install.packages("magrittr")

library(stringi)
library(magrittr)
library(tidyr)
library(scales)
library(tidyverse)
library(viridisLite)
library(viridis)
library(forcats)

setwd("/home/nuno/Documents/Effective Altruism/EA Mental Health/")
setwd("/home/nuno/Documents/Effective Altruism/EA Mental Health/2")

A <- read.csv("Responses.csv", header=TRUE, sep=",", stringsAsFactors = FALSE)
  ## Note: if you don't select stringsAsFactors = FALSE, you can't manipulate elements as strings (for example, with grepl) unless you first process your data. For example: as.character(A[,column of interest]) -> A[,column of interest]

# Q1:How involved are you with the effective altruism community
colnames(A)[3]
(ggplot(data = A, aes(x = How.involved.are.you.in.the.Effective.Altruism.Community.)) +
  geom_bar(
      aes(y = (..count..)/sum(..count..)), 
      fill="royalblue4"
      )+
  labs(
    title="Q1: How involved are you in the Effective Altruism community?",
    subtitle="n =303", 
    x=NULL, 
    y="Frequency (%)", 
    caption="@EA Mental Health Survey"
    )+
  scale_y_continuous(
    labels = percent_format(accuracy = 1)
    )+
  coord_flip()+
  theme(
    legend.title = element_blank(), 
    plot.subtitle = element_text(hjust = 0.5), 
    plot.title = element_text(hjust = 0.5),
    legend.position="bottom")
) %>%
  ggsave("Q1.png", plot =., units="in", width=10, height=5, dpi=800)

## Q2: Do you attend EA Meetings
colnames(A)[4]

### This adds a line:
for( i in c(1:length(A[,4]))){
  if(A[,4][i]=="No, but I regularly participate in an EA online group"){
    A[,4][i] = "No, but I regularly participate\n in an EA online group"
  }
}

ifelse(A[,4]=="", "Did not answer", A[,4]) -> A[,4]

(ggplot(data = A, aes(x = Do.you.attend.EA.meetings.)) +
  geom_bar(
    aes(
      y = (..count..)/sum(..count..)), 
      fill="royalblue4"
      )+
  coord_flip()+
  labs(
    title="Q2:Do you attend EA Meetings?", 
    subtitle="n =303", 
    x=NULL, 
    y="Frequency (%)",
    caption="@EA Mental Health Survey"
    )+
  scale_y_continuous(
    labels = percent_format(accuracy = 1)
    )+
  theme(
    legend.title = element_blank(), 
    plot.subtitle = element_text(hjust = 0.5), 
    plot.title = element_text(hjust = 0.5),
    legend.position="bottom")
) %>%
  ggsave("Q2.png", plot =., units="in", width=10, height=5, dpi=800)

## Q3: How much impact do EA ideas have on your life
ini = 2
i=ini+3
colnames(A)[i] %>%
  gsub("\\."," ",.) %>%
  c("Q",i-ini,": ",.,"?") %>%
  paste(., collapse="") %>%
  gsub(" \\?","?",.)  -> titulo

ifelse(A[,i]=="", "Did not answer", A[,i]) -> A[,i]

(ggplot(data = A, aes(x = get(colnames(A)[i]))) +
  geom_bar(
    aes(y = (..count..)/sum(..count..)), 
    fill="royalblue4")+
  coord_flip()+
  labs(
    title=titulo, 
    subtitle="n =303", 
    x=NULL, 
    y="Frequency (%)", 
    caption="@EA Mental Health Survey"
    )+
  scale_y_continuous(
    labels = percent_format(accuracy = 1)
    )+
  theme(legend.title = element_blank(), plot.subtitle = element_text(hjust = 0.5), plot.title = element_text(hjust = 0.5),legend.position="bottom")
) %>%
ggsave("Q3.png", plot =., units="in", width=10, height=5, dpi=800)


## Q4: Do you donate part of your income to GiveWell recommended charities?
ini = 2
i=ini+4

colnames(A)[i] %>%
  gsub("\\."," ",.) %>%
  c("Q",i-ini,": ",.,"?") %>%
  paste(., collapse="") %>%
  gsub(" \\?","?",.)  -> titulo

ifelse(A[,i]=="", "Did not answer", A[,i]) -> A[,i]

(ggplot(data = A, aes(x = get(colnames(A)[i]))) +
  geom_bar(
    aes(y = (..count..)/sum(..count..)),
    fill="royalblue4"
    )+
  labs(
    title=titulo, 
    subtitle="n =303", 
    x=NULL, 
    y="Frequency (%)", 
    caption="@EA Mental Health Survey"
    )+
  scale_y_continuous(
    labels = percent_format(accuracy = 1)
    )+
  theme(
    legend.title = element_blank(), 
    plot.subtitle = element_text(hjust = 0.5), 
    plot.title = element_text(hjust = 0.5),
    legend.position="bottom"
    ) 
  )%>%
  ggsave("Q4.png", plot =., units="in", width=10, height=5, dpi=800)

  ## I just discovered piping! Elegant!

## Q5: Which of these conditions have you been formally diagnosed with?
i=ini+5

colnames(A)[i] %>%
  gsub("\\."," ",.) %>%
  c("Q",i-ini,": ",.,"?") %>%
  paste(., collapse="") %>%
  gsub(" \\?","?",.)  -> titulo

### This question will require some processing:
List_of_Mental_Illnesses = c("Depression","Anxiety", "Obsessive Compulsive Disorder", "Disordered Eating", "Alcoholism", "Drug Addiction", "Borderline Personality Disorder", "Bipolar Disorder", "Autism", "ADHD", "Schizophrenia")

MI=list()
for( j in List_of_Mental_Illnesses){
  MI$count = c(MI$count, 
                sum(
                  grepl(j, A[,i])
                )
              )
  MI$type = c(MI$type, j)
}

MI$frequency = 100*MI$count / l
as.data.frame(MI) -> MI

(ggplot(data = MI, aes(x =type, y = frequency)) +
  geom_histogram(
    stat="identity", 
    fill="royalblue4")+
  scale_y_continuous(
    labels=function(x) paste0(x,"%"),
    limits = c(0, 33)
    )+
  scale_x_discrete(
    limits=
      List_of_Mental_Illnesses[
          order(MI$count)# c(1,2,10,8,3,4,9,7,11,6,5) ## This orders the ilnesses by frequency
        ]
    )+
  coord_flip()+    
  labs(
    title=titulo, 
    subtitle="n =303", 
    x=NULL, 
    y="Frequency (%)", 
    caption="@EA Mental Health Survey"
    )+
  theme(
    legend.title = element_blank(), 
    plot.subtitle = element_text(hjust = 0.5), 
    plot.title = element_text(hjust = 0.5),
    legend.position="bottom"
    ) 
  ) %>%
    ggsave("Q5c.png", plot =., units="in", width=10, height=5, dpi=800)

## Q6: Which of these conditions do you think you may have  but have never been formally diagnosed with?
  ## I guess we can just reuse the code from Q5

i=ini+6
  
colnames(A)[i] %>%
  gsub("\\."," ",.) %>%
  c("Q",i-ini,": ",.,"?") %>%
  paste(., collapse="") %>%
  gsub(" \\?","?",.)  -> titulo

titulo= "Q6: Which of these conditions do you think you may have\nbut have never been formally diagnosed with?"

List_of_Mental_Illnesses = c("Depression","Anxiety", "Obsessive Compulsive Disorder", "Disordered Eating", "Alcoholism", "Drug Addiction", "Borderline Personality Disorder", "Bipolar Disorder", "Autism", "ADHD", "Schizophrenia")

MI2=list()
for( j in List_of_Mental_Illnesses){
  MI2$count = c(MI2$count, 
               sum(
                 grepl(j, A[,i])
               )
  )
  MI2$type = c(MI2$type, j)
}

MI2$frequency = 100*MI2$count / l
as.data.frame(MI2) -> MI2

(ggplot(data = MI2, aes(x =type, y = frequency)) +
    geom_histogram(
      stat="identity", 
      fill="royalblue4")+
    scale_y_continuous(
      labels=function(x) paste0(x,"%"),
      limits = c(0, 33)
    )+
    scale_x_discrete(
      limits=
        List_of_Mental_Illnesses[
            order(MI2$count)
            #c(2,1,10,9,4,3,8,5,7,6,11) ## This orders the ilnesses by frequency
          ]
    )+
    coord_flip()+    
    labs(
      title=titulo, 
      subtitle="n =303", 
      x=NULL, 
      y="Frequency (%)", 
      caption="@EA Mental Health Survey"
    )+
    theme(
      legend.title = element_blank(), 
      plot.subtitle = element_text(hjust = 0.5), 
      plot.title = element_text(hjust = 0.5),
      legend.position="bottom"
    ) 
) %>%
  ggsave("Q6c.png", plot =., units="in", width=10, height=5, dpi=800)


## Q5 and Q6 together. Both (formally diagnosed or not) count.
i=ini+5
colnames(A)[i] %>%
  gsub("\\."," ",.) %>%
  c("Q",i-ini,": ",.,"?") %>%
  paste(., collapse="") %>%
  gsub(" \\?","?",.)  -> titulo1

colnames(A)[i+1] %>%
  gsub("\\."," ",.) %>%
  c("Q",i-ini+1,": ",.,"?") %>%
  paste(., collapse="") %>%
  gsub(" \\?","?",.)  -> titulo2

paste(c(titulo1, titulo2), collapse="\n+\n")-> titulo

MI$diagnosed = "Diagnosed"
MI2$diagnosed = "Intuited"

as.data.frame(rbind(MI, MI2)) -> MI3

(ggplot(data = MI3, aes(x =type, y = frequency, fill=diagnosed)) +
    geom_histogram(
      stat="identity",
      position= position_stack(reverse = TRUE)
      )+
    scale_fill_manual(
      values=c("royalblue3","steelblue1"))+
    scale_x_discrete(
      limits=
        List_of_Mental_Illnesses[
          order(MI$count +MI2$count ) ## This orders the ilnesses by frequency
          ]
    )+
    scale_y_continuous(
      labels=function(x) paste0(x,"%"))+
    coord_flip()+
    labs(
      title=titulo, 
      subtitle="n =303", 
      x=NULL, 
      y="Frequency (%)", 
      caption="@EA Mental Health Survey"
      )+
    theme(
      legend.title = element_blank(), 
      plot.subtitle = element_text(hjust = 0.5), 
      plot.title = element_text(hjust = 0.5),
      legend.position="bottom"
      ) 
  ) %>%
    ggsave("Q5and6.png", plot =., units="in", width=10, height=5, dpi=800)

## Q7: During the past 14 days  how many hours did you miss from work because of mental health problems?

i=ini+7

read.csv("horus.csv", header=TRUE, sep=",", stringsAsFactors = FALSE) %>%
  .[complete.cases(sapply(.,as.integer)),] %>%
  .[order(as.integer(.$hours_worked)),] -> horus

horus$person2 = c(1:dim(horus)[1])

### We take the opportunity to calculate the average hours worked and missed.
colnames(horus)

    mean(as.integer(horus$hours_worked))
      # the mean is 59.58182 ~ 59.5 hours / 2 weeks
    sd(as.integer(horus$hours_worked))
      # 31.19616
    mean(as.integer(horus$hours_missed_processed))
      # 5.749091
    sd(as.integer(horus$hours_missed_processed))
      # 14.44113
    mean(as.integer(horus$hours_missed_processed_other_reasons))
      # 6.578182
    sd(as.integer(horus$hours_missed_processed_other_reasons))
      # 13.7492

### Problem: not everyone has answered each of the three questions, and the above only takes into account those who have. If we take into account 
    read.csv("horus.csv", header=TRUE, sep=",", stringsAsFactors = FALSE)-> horus2
    mean(as.integer(horus2$hours_worked), na.rm=TRUE)
    # the mean is 58.31707 ~ 58.3 hours / 2 weeks
    sd(as.integer(horus2$hours_worked), na.rm=TRUE)
    # 31.74407
    mean(as.integer(horus2$hours_missed_processed), na.rm=TRUE)
    # 6.190141
    sd(as.integer(horus2$hours_missed_processed), na.rm=TRUE)
    # 15.05628
    mean(as.integer(horus2$hours_missed_processed_other_reasons), na.rm=TRUE)
    # 6.758007
    sd(as.integer(horus2$hours_missed_processed_other_reasons), na.rm=TRUE)
    # 14.06154

### However, what we're interested in is in the ratio. The following plot ignores partial nonrespondents:

### We continue modifying the list further, to arrive at the plot we want

gather(horus,"type of hour", "hours", c(2:4)) -> horus
# Gather transforms:
    #  Person, hour_worked, hours_missed, hours_missed_non_mental ilness
    #person hours_worked hours_missed_processed hours_missed_processed_other_reasons
    #1      1           40                      0                                    0
    #> horus[c(1:10),]
    #person hours_worked hours_missed_processed hours_missed_processed_other_reasons
    #1       1           40                      0                                    0
    #2       2           80                      0                                    0
    #3       3          100                      2                                   10
    #4       4           80                      0                                    0

## Into:
    #person type of hour hours
    #1      1 hours_worked    40
    #2      2 hours_worked    80
    #3      3 hours_worked   100
    #4      4 hours_worked    80
    #5      5 hours_worked     0

horus$hours = as.integer(horus$hours)

(ggplot(data = horus, aes(x =person2, y = hours, fill=get("type of hour"))) +
  geom_histogram(stat="identity") +
  scale_fill_viridis_d(
    begin=1, 
    end=.2, 
    labels=
      c(
        " Work hours missed for mental health reasons",
        " Work hours missed for other reasons", 
        " Total hours of work" 
      )
    )+
  labs(
    title="Q7-9: Mental health and work hours", 
    subtitle="n =303", 
    x="Persons, ordered by hours worked", 
    y="Hours", 
    caption="@EA Mental Health Survey")+
  theme(
    legend.title = element_blank(), 
    plot.subtitle = element_text(hjust = 0.5), 
    plot.title = element_text(hjust = 0.5),
    legend.position="bottom", 
    legend.direction = "vertical")
)%>%
    ggsave("Q7-9c.png", plot =., units="in", width=12, height=8, dpi=800)

## Q10: During the past 14 days  how much did mental health problems affect your productivity while you were working?

ini=2
i=ini+10

colnames(A)[i] %>%
  gsub("\\."," ",.) %>%
  c("Q",i-ini,": ",.,"?") %>%
  paste(., collapse="") %>%
  gsub(" \\?","?",.)  -> titulo

titulo ="Q10: During the past 14 days, how much did mental health problems \naffect your productivity while you were working?"

(ggplot(data = A, aes(x = get(colnames(A)[i]))) +
  geom_bar(
    aes(y = (..count..)/sum(..count..)), 
    fill="royalblue4"
    )+
  labs(
    title=titulo, 
    subtitle="n =303", 
    x=NULL, y="Frequency (%)", 
    caption="@EA Mental Health Survey"
    )+
  scale_y_continuous(
    labels = percent_format(accuracy = 1)
    )+
  theme(
    legend.title = element_blank(), 
    plot.subtitle = element_text(hjust = 0.5), 
    plot.title = element_text(hjust = 0.5),
    legend.position="bottom") 
)%>%
  ggsave("Q10.png", plot =., units="in", width=10, height=5, dpi=800)

## [1] "Q11: During the past 14 days  how much did mental health problems affect your ability to do your regular daily activities  other than work at a job?"

ini = 2
i=ini+11

colnames(A)[i] %>%
  gsub("\\."," ",.) %>%
  c("Q",i-ini,": ",.,"?") %>%
  paste(., collapse="") %>%
  gsub("  "," ",.) %>% 
  gsub("   \\?","?",.) %>% 
  gsub(" \\?","?",.)-> titulo

titulo = "Q11: During the past 14 days, how much did mental health problems\n affect your ability to do your regular daily activities\n other than work at a job?"

(ggplot(data = A, aes(x = get(colnames(A)[i]))) +
    geom_bar(aes(
      y = (..count..)/sum(..count..)), 
      fill="royalblue4")+
    labs(
      title=titulo, 
      subtitle="n =303", 
      x=NULL, y="Frequency (%)", 
      caption="@EA Mental Health Survey"
      )+
    scale_y_continuous(
      labels = percent_format(accuracy = 1)
      )+
    theme(
      legend.title = element_blank(), 
      plot.subtitle = element_text(hjust = 0.5), 
      plot.title = element_text(hjust = 0.5),
      legend.position="bottom") 
)%>%
  ggsave(paste(c("Q", i-ini,".png"), collapse=""), plot =., units="in", width=10, height=5, dpi=800)

## Q12: Are the past 14 days representative of your average mental health?

ini = 2
i=ini+12

colnames(A)[i] %>%
  gsub("\\."," ",.) %>%
  c("Q",i-ini,": ",.,"?") %>%
  paste(., collapse="") %>%
  gsub("   \\?","?",.) %>% 
  gsub(" \\?","?",.)-> titulo

titulo

unique(A[,i])
ifelse(grepl("No", A[,i]),"No", A[,i]) -> A[,i]
ifelse(grepl("Yes", A[,i]),"Yes", A[,i]) -> A[,i]
ifelse(A[,i]=="","Did not answer", A[,i]) -> A[,i]
unique(A[,i])

(ggplot(data = A, aes(x = get(colnames(A)[i]))) +
    geom_bar(
      aes(y = (..count..)/sum(..count..)), 
      fill="royalblue4")+
    labs(
      title=titulo, 
      subtitle="n =303", 
      x=NULL, y="Frequency (%)", 
      caption="@EA Mental Health Survey"
      )+
    scale_y_continuous(
      labels = percent_format(accuracy = 1)
      )+
    theme(
      legend.title = element_blank(), 
      plot.subtitle = element_text(hjust = 0.5), 
      plot.title = element_text(hjust = 0.5),
      legend.position="bottom"
      ) 
)%>%
  ggsave(paste(c("Q", i-ini,".png"), collapse=""), plot =., units="in", width=10, height=5, dpi=800)

## Perhaps do some linear regression model? Predict hours out of of the abstract questions + how representative it is? / whether it has been better. Exploratory plots.

## [1] "Q13: If you answered  No  to the previous question  how has your mental health differed over the past 14 days?"

ini = 2
i=ini+13

colnames(A)[i] %>%
  gsub("\\."," ",.) %>%
  c("Q",i-ini,": ",.,"?") %>%
  paste(., collapse="") %>%
  gsub("   \\?","?",.) %>% 
  gsub(" \\?","?",.)-> titulo

titulo = "Q13: If you answered 'No' to the previous question,\nhow has your mental health differed over the past 14 days?"

  ## The unprocessed data is just a comment. Must be processed.
  ## Levels: Much worse, worse, better, Much better

#representative = c()
#for(j in c(1:length(A[,i]))){
#  if(A[,i][j]!=""){
#    representative=c(representative, readline(paste(c("j=",j," Answer=", A[,i][j],". Value = "),collapse="")))
#  }
#  else{
#    representative=c(representative, NA)
#  }
#}

representative2 <- read.csv(file = "Q13_how_has_your_mental_health_differed_over_the_past_14_days.csv",col.names=TRUE)

as.data.frame(representative2[!is.na(representative2)]) -> rep

(ggplot(data = rep, aes(x = get(colnames(rep)[1]))) +
    geom_bar(
      aes(y = (..count..)/sum(..count..)), 
      fill="royalblue4"
      )+
    labs(
      title=titulo, 
      subtitle="n =303",
      x=NULL, y="Frequency (%)", 
      caption="@EA Mental Health Survey")+
    scale_y_continuous(
      labels = percent_format(accuracy = 1)
      )+
    theme(
      legend.title = element_blank(), 
      plot.subtitle = element_text(hjust = 0.5), 
      plot.title = element_text(hjust = 0.5),
      legend.position="bottom") )%>%
  ggsave(paste(c("Q", i-ini,".png"), collapse=""), plot =., units="in", width=10, height=5, dpi=800)

## [1] "Q14: How challenging was it to receive the mental healthcare services you needed within the past 12 months?"

ini = 2
i=ini+14

colnames(A)[i] %>%
  gsub("\\."," ",.) %>%
  c("Q",i-ini,": ",.,"?") %>%
  paste(., collapse="") %>%
  gsub("   \\?","?",.) %>% 
  gsub(" \\?","?",.)-> titulo

titulo="Q14: How challenging was it to receive the mental healthcare services\n you needed within the past 12 months?"

ifelse(grepl("I have not", A[,i]), "I have not needed mental healthcare services\n in the past 12 months.", A[,i]) -> A[,i]
ifelse(grepl("I did not seek", A[,i]), "I did not seek out mental healthcare services \nin the past 12 months but believe \nI could have benefitted from them.", A[,i]) -> A[,i]
ifelse(A[,i]=="", "No answer", A[,i]) -> A[,i]


(ggplot(data = A, aes(x = get(colnames(A)[i]))) +
    geom_bar(aes(y = (..count..)/sum(..count..)), fill="royalblue4")+
    #coord_flip()+
    coord_flip()+
    scale_x_discrete(limits=unique(A[,i])[rev(c(8,1,5,4,3,2,6))])+
    labs(
      title=titulo, 
      subtitle="n =303",
      x=NULL, y="Frequency (%)", 
      caption="@EA Mental Health Survey")+
    scale_y_continuous(
      labels = percent_format(accuracy = 1)
    )+
    theme(
      legend.title = element_blank(), 
      plot.subtitle = element_text(hjust = 0.5), 
      plot.title = element_text(hjust = 0.5),
      legend.position="bottom") 
)%>%
  ggsave(paste(c("Q", i-ini,".png"), collapse=""), plot =., units="in", width=10, height=5, dpi=800)

## [1] "Q15: How challenging is it to find useful information on mental healthcare services?"

ini = 2
i=ini+15

colnames(A)[i] %>%
  gsub("\\."," ",.) %>%
  c("Q",i-ini,": ",.,"?") %>%
  paste(., collapse="") %>%
  gsub("   \\?","?",.) %>% 
  gsub(" \\?","?",.)-> titulo

ifelse(A[,i]=="", "No answer", A[,i]) -> A[,i]
ifelse(grepl("I have not looked for", A[,i]),"I have not looked for information\n on mental healthcare services.", A[,i]) -> A[,i]

(ggplot(data = A, aes(x = get(colnames(A)[i]))) +
    geom_bar(aes(y = (..count..)/sum(..count..)), fill="royalblue4")+
    #coord_flip()+
    coord_flip()+
    scale_x_discrete(limits=unique(A[,i])[rev(c(7,4,1,5,3,2,6))])+
    labs(
      title=titulo, 
      subtitle="n =303",
      x=NULL, y="Frequency (%)", 
      caption="@EA Mental Health Survey")+
    scale_y_continuous(
      labels = percent_format(accuracy = 1)
    )+
    theme(
      legend.title = element_blank(), 
      plot.subtitle = element_text(hjust = 0.5), 
      plot.title = element_text(hjust = 0.5),
      legend.position="bottom") 
)%>%
  ggsave(paste(c("Q", i-ini,".png"), collapse=""), plot =., units="in", width=12, height=5, dpi=800)

## [1] "Q16: Do you experience financial difficulties as a result of mental healthcare?"

ini = 2
i=ini+16

colnames(A)[i] %>%
  gsub("\\."," ",.) %>%
  c("Q",i-ini,": ",.,"?") %>%
  paste(., collapse="") %>%
  gsub("   \\?","?",.) %>% 
  gsub(" \\?","?",.)-> titulo

titulo = "Q16: Do you experience financial difficulties\nas a result of mental healthcare?"

ifelse(A[,i]=="", "No answer", A[,i]) -> A[,i]

(ggplot(data = A, aes(x = get(colnames(A)[i]))) +
    geom_bar(aes(y = (..count..)/sum(..count..)), fill="royalblue4")+
    #coord_flip()+
    coord_flip()+
    scale_x_discrete(limits=unique(A[,i])[rev(c(1,4,3,2,5))])+
    labs(
      title=titulo, 
      subtitle="n =303",
      x=NULL, y="Frequency (%)", 
      caption="@EA Mental Health Survey")+
    scale_y_continuous(
      labels = percent_format(accuracy = 1)
    )+
    theme(
      legend.title = element_blank(), 
      plot.subtitle = element_text(hjust = 0.5), 
      plot.title = element_text(hjust = 0.5),
      legend.position="bottom") 
)%>%
  ggsave(paste(c("Q", i-ini,".png"), collapse=""), plot =., units="in", width=12, height=5, dpi=800)

## [1] "I am currently receiving the mental healthcare I need"

ini = 2
i=ini+17

colnames(A)[i] %>%
  gsub("\\."," ",.) %>%
  c("Q",i-ini,": ",.,"?") %>%
  paste(., collapse="") %>%
  gsub("   \\?","?",.) %>% 
  gsub(" \\?","?",.)-> titulo
titulo
titulo="Q17: I am currently receiving the mental healthcare I need"

ifelse(A[,i]=="", "No answer", A[,i]) -> A[,i]

(ggplot(data = A, aes(x = fct_infreq(get(colnames(A)[i])))) +
    geom_bar(aes(y = (..count..)/sum(..count..)), fill="royalblue4")+
    coord_flip()+
    scale_x_discrete(limits=unique(A[,i])[rev(c(5,4, 2,1,3,6))])+
    labs(
      title=titulo, 
      subtitle="n =303",
      x=NULL, y="Frequency (%)", 
      caption="@EA Mental Health Survey")+
    scale_y_continuous(
      labels = percent_format(accuracy = 1)
    )+
    theme(
      legend.title = element_blank(), 
      plot.subtitle = element_text(hjust = 0.5), 
      plot.title = element_text(hjust = 0.5),
      legend.position="bottom") 
)%>%
  ggsave(paste(c("Q", i-ini,".png"), collapse=""), plot =., units="in", width=12, height=5, dpi=800)

Receiving = list()
for( j in unique(A[,i])){
  Receiving$answer = c(Receiving$answer,sum(A[,i]==j)/303)
  Receiving$category = c(Receiving$category,j)
}

Receiving = as.data.frame(Receiving)
factor(Receiving$category, levels = Receiving$category) -> Receiving$category
Receiving = Receiving[rev(c(1,2,5,6,4,3)),]
factor(Receiving$category, levels = Receiving$category) -> Receiving$category

(ggplot(data = Receiving, aes(x="", y = answer, fill = category)) +
    geom_bar(stat="identity")+
    labs(
        title=titulo, 
        subtitle="n =303", 
        x=NULL, 
        y=NULL, 
        caption="@EA Mental Health Survey"
        )+
    scale_fill_viridis_d(
        #limits=unique(A[,i])[rev(c(6,1,2,5,4,3))]
        )+
    coord_polar(
        "y", 
        start=0) +
    theme(
        axis.text = element_blank(),
        axis.ticks = element_blank(),
        panel.grid  = element_blank())+
    theme(
        legend.title = element_blank(), 
        plot.subtitle = element_text(hjust = 0.5), 
        plot.title = element_text(hjust = 0.5)
        ) 
) %>%
  ggsave(paste(c("Q", i-ini,"-7",".png"), collapse=""), plot =., units="in", width=8, height=6, dpi=800)


## "Q18: How satisfied are you\nwith the mental healthcare you ve received?"

ini = 2
i=ini+18

colnames(A)[i] %>%
  gsub("\\."," ",.) %>%
  c("Q",i-ini,": ",.,"?") %>%
  paste(., collapse="") %>%
  gsub("   \\?","?",.) %>% 
  gsub(" \\?","?",.)-> titulo

titulo="Q18: How satisfied are you\nwith the mental healthcare you ve received?"

ifelse(A[,i]=="", "No answer", A[,i]) -> A[,i]
ifelse(A[,i]=="I have not received mental healthcare services.", "I have not received mental\nhealthcare services.", A[,i]) -> A[,i]

(ggplot(data = A, aes(x = get(colnames(A)[i]))) +
    geom_bar(aes(y = (..count..)/sum(..count..)), fill="royalblue4")+
    coord_flip()+
    scale_x_discrete(
      limits=unique(A[,i])[rev(c(5,1,4,3,2,6))]
      )+
    labs(
      title=titulo, 
      subtitle="n =303",
      x=NULL, y="Frequency (%)", 
      caption="@EA Mental Health Survey")+
    scale_y_continuous(
      labels = percent_format(accuracy = 1)
    )+
    theme(
      legend.title = element_blank(), 
      plot.subtitle = element_text(hjust = 0.5), 
      plot.title = element_text(hjust = 0.5),
      legend.position="bottom") 
)%>%
  ggsave(paste(c("Q", i-ini,".png"), collapse=""), plot =., units="in", width=12, height=5, dpi=800)

## [1] "Q19: Which of the following treatments have you tried?"

ini = 2
i=ini+19

colnames(A)[i] %>%
  gsub("\\."," ",.) %>%
  c("Q",i-ini,": ",.,"?") %>%
  paste(., collapse="") %>%
  gsub("   \\?","?",.) %>% 
  gsub(" \\?","?",.)-> titulo

titulo = "Which treatments have you tried?"

## The following is somewhat ugly, but I didn't want to take the time to write short code.
Treatments = list()

Treatments$num = c(Treatments$num, sum(grepl("CBT", A[,i])))
Treatments$name = c(Treatments$name, "CBT")

Treatments$num = c(Treatments$num, sum(grepl("Group therapy", A[,i])))
Treatments$name = c(Treatments$name, "Group therapy")

Treatments$num = c(Treatments$num, sum(grepl("upport", A[,i])) + sum(grepl("connection",A[,i])))
Treatments$name = c(Treatments$name, "Group / friend support")

Treatments$num = c(Treatments$num, sum(grepl("Inpatient", A[,i])))
Treatments$name = c(Treatments$name, "Inpatient mental health treatment")

Treatments$num = c(Treatments$num, sum(grepl("None", A[,i])) + sum(A[,i]==""))
Treatments$name = c(Treatments$name, "None")

Treatments$num = c(Treatments$num, sum(grepl("Meditation", A[,i])))
Treatments$name = c(Treatments$name, "Meditation,mindfulness techniques,\nalternative therapies
")

Treatments$num = c(Treatments$num, sum(grepl("Prescribed", A[,i])))
Treatments$name = c(Treatments$name, "Prescribed medication")

Treatments$num = c(Treatments$num, sum(grepl("sychotherapy", A[,i])) + sum(grepl("ounselling", A[,i])) + sum(grepl("psychology service", A[,i])))
Treatments$name = c(Treatments$name, "Psychotherapy")

Treatments$num = c(Treatments$num, sum(grepl("Self-medication", A[,i])) + sum(grepl("GABA", A[,i])) + sum(grepl("St John", A[,i])))
Treatments$name = c(Treatments$name, "Self-medication")

Treatments$num = c(Treatments$num, sum(grepl("books", A[,i])) + sum(grepl("Bibliotherapy", A[,i])))
Treatments$name = c(Treatments$name, "Books")

Treatments$num = c(Treatments$num, sum(grepl("elf-organisation", A[,i])) + sum(grepl("stick to routine", A[,i])) + sum(grepl("self-discipline", A[,i])))
Treatments$name = c(Treatments$name, "Self-discipline")

Treatments$num = c(Treatments$num, sum(grepl("ports", A[,i])) + sum(grepl("running", A[,i])) + sum(grepl("physical activity", A[,i])) + sum(grepl("exercise", A[,i])))
Treatments$name = c(Treatments$name, "Exercise")

Treatments$num = c(Treatments$num, sum(grepl("CBD", A[,i])) + sum(grepl("CFAR", A[,i])) + sum(grepl("EMDR", A[,i])) + sum(grepl("gratitude", A[,i])) + sum(grepl("diet", A[,i])) + sum(grepl("free phone consultant", A[,i])) + sum(grepl("PRI", A[,i])) + sum(grepl("practice sleep hygene", A[,i])) + sum(grepl("Sun lamp", A[,i])) + sum(grepl("Talking Space", A[,i])) + sum(grepl("Neuro-linguistic programming", A[,i])) )
Treatments$name = c(Treatments$name, "Others")

as.data.frame(Treatments) ->Treatments
Treatments[order(Treatments$num),] -> Treatments
factor(Treatments$name, levels = Treatments$name) -> Treatments$name

View(Treatments)

(ggplot(data = Treatments, aes(x = name, y = num)) +
    geom_col(fill="royalblue4")+
    coord_flip()+
    labs(
      title=titulo, 
      subtitle="n =303",
      x=NULL, y="Frequency (%)", 
      caption="@EA Mental Health Survey")+
    theme(
      legend.title = element_blank(), 
      plot.subtitle = element_text(hjust = 0.5), 
      plot.title = element_text(hjust = 0.5),
      legend.position="bottom") 
)%>%
  ggsave(paste(c("Q", i-ini,"-2.png"), collapse=""), plot =., units="in", width=9, height=6, dpi=800)

## [1] "Q20: Any other thoughts you d like to share on mental healthcare access?"

i=ini +20
colnames(A)[i] %>%
  gsub("\\."," ",.) %>%
  c("Q",i-ini,": ",.,"?") %>%
  paste(., collapse="") %>%
  gsub("   \\?","?",.) %>% 
  gsub(" \\?","?",.)-> titulo

## [1] "Q21: I believe that offering mental health resources to its members is an effective use of the EA Community's resources"
i=ini +21
colnames(A)[i] %>%
  gsub("\\."," ",.) %>%
  c("Q",i-ini,": ",.,"?") %>%
  paste(., collapse="") %>%
  gsub("   \\?","?",.) %>% 
  gsub(" \\?","?",.)-> titulo
titulo = "Q21: I believe that offering mental health resources to its members is an effective use of the EA Community's resources"

ifelse(A[,i]=="", "No answer", A[,i]) -> A[,i]

(ggplot(data = A, aes(x = get(colnames(A)[i]))) +
    geom_bar(aes(y = (..count..)/sum(..count..)), fill="royalblue4")+
    coord_flip()+
    scale_x_discrete(
      limits=unique(A[,i])[rev(c(3,1,2,4,6,5))]
      )+
    labs(
      title=titulo, 
      subtitle="n =303",
      x=NULL, y="Frequency (%)", 
      caption="@EA Mental Health Survey")+
    scale_y_continuous(
      labels = percent_format(accuracy = 1)
    )+
    theme(
      legend.title = element_blank(), 
      plot.subtitle = element_text(hjust = 0.5), 
      plot.title = element_text(hjust = 0.5),
      legend.position="bottom") 
)%>%
  ggsave(paste(c("Q", i-ini,".png"), collapse=""), plot =., units="in", width=12, height=5, dpi=800)

## [1] "Q22: I believe that offering mental health resources to effective altruists is NOT likely to be one of the most effective interventions available"

i=ini +22
colnames(A)[i] %>%
  gsub("\\."," ",.) %>%
  c("Q",i-ini,": ",.,"?") %>%
  paste(., collapse="") %>%
  gsub("   \\?","?",.) %>% 
  gsub(" \\?","?",.)-> titulo

titulo = "Q22: I believe that offering mental health resources to effective altruists\nis NOT likely to be one of the most effective interventions available"

ifelse(A[,i]=="", "No answer", A[,i]) -> A[,i]

(ggplot(data = A, aes(x = get(colnames(A)[i]))) +
    geom_bar(aes(y = (..count..)/sum(..count..)), fill="royalblue4")+
    coord_flip()+
    scale_x_discrete(limits=unique(A[,i])[rev(c(6,3,1,2,5,4))])+
    labs(
      title=titulo, 
      subtitle="n =303",
      x=NULL, y="Frequency (%)", 
      caption="@EA Mental Health Survey")+
    scale_y_continuous(
      labels = percent_format(accuracy = 1)
    )+
    theme(
      legend.title = element_blank(), 
      plot.subtitle = element_text(hjust = 0.5), 
      plot.title = element_text(hjust = 0.5),
      legend.position="bottom") 
)%>%
  ggsave(paste(c("Q", i-ini,".png"), collapse=""), plot =., units="in", width=12, height=5, dpi=800)

## Q23: Please explain your answers to the questions above

i=ini +23
colnames(A)[i] %>%
  gsub("\\."," ",.) %>%
  c("Q",i-ini,": ",.,"?") %>%
  paste(., collapse="") %>%
  gsub("   \\?","?",.) %>% 
  gsub(" \\?","?",.)-> titulo

## [1] "Q24: Do you think you could personally benefit from EA community mental health resources?"

i=ini +24
colnames(A)[i] %>%
  gsub("\\."," ",.) %>%
  c("Q",i-ini,": ",.,"?") %>%
  paste(., collapse="") %>%
  gsub("   \\?","?",.) %>% 
  gsub(" \\?","?",.)-> titulo

ifelse(A[,i]=="", "No answer", A[,i]) -> A[,i]
ifelse(is.na(A[,i]), "No answer", A[,i]) -> A[,i]

(ggplot(data = A, aes(x = get(colnames(A)[i]))) +
    geom_bar(aes(y = (..count..)/sum(..count..)), fill="royalblue4")+
    scale_x_discrete(limits=unique(A[,i])[rev(c(3,2,1,4))])+
    labs(
      title=titulo, 
      subtitle="n =303",
      x=NULL, y="Frequency (%)", 
      caption="@EA Mental Health Survey")+
    scale_y_continuous(
      labels = percent_format(accuracy = 1)
    )+
    theme(
      legend.title = element_blank(), 
      plot.subtitle = element_text(hjust = 0.5), 
      plot.title = element_text(hjust = 0.5),
      legend.position="bottom") 
)%>%  
  ggsave(paste(c("Q", i-ini,".png"), collapse=""), plot =., units="in", width=12, height=5, dpi=800)


## [1] "Q25: Which  if any  of the following resources do you think you could potentially benefit from?"
ini=2
i=ini +25
colnames(A)[i] %>%
  gsub("\\."," ",.) %>%
  c("Q",i-ini,": ",.,"?") %>%
  paste(., collapse="") %>%
  gsub("   \\?","?",.) %>% 
  gsub(" \\?","?",.)-> titulo

ifelse(A[,i]=="", "No answer", A[,i]) -> A[,i]
ifelse(is.na(A[,i]), "No answer", A[,i]) -> A[,i]

Answers = c("No answer",
"Non-clinical mental healthcare supporter/advisor",
"Employee assistance programs",
"Mental healthcare navigators",
"EA-knowledgeable therapists",
"Funds for treatment",
"Other")

s=0
t=as.vector(A[,i])
for(j in c(1:303)){
  if(!(t[j]=="" | t[j]%in%Answers)){
    cat("\n",j+1," - ",t[j], ":")
    s=s+as.integer(readline())  
  }
} ## s=31

You_could_benefit=list()
for(j in Answers){
  You_could_benefit$thing = c(You_could_benefit$thing, j)
  You_could_benefit$count = c(You_could_benefit$count, sum(grepl(j, A[,i])))
}
as.data.frame(You_could_benefit) -> You_could_benefit
You_could_benefit$count[7] = 31

    ##  158  -  Someone answered they could benefit from heroin. Heroin :

(ggplot(data = You_could_benefit, aes(x = thing, y = count/303 )) +
    geom_col(fill="royalblue4")+
    #coord_flip()+
    labs(
      title=titulo, 
      subtitle="n =303", 
      x=NULL, 
      y="Frequency (%)", 
      caption="@EA Mental Health Survey"
      )+
    coord_flip()+
    scale_x_discrete(
      #limits=unique(A[,i])[rev(c(3,2,1,4))]
      )+
    scale_y_continuous(
      labels = percent_format(accuracy = 1)
      )+
    theme(
      legend.title = element_blank(),
      plot.subtitle = element_text(hjust = 0.5), 
      plot.title = element_text(hjust = 0.5),
      legend.position="bottom") 
)%>%
  ggsave(paste(c("Q", i-ini,".png"), collapse=""), plot =., units="in", width=12, height=5, dpi=800)



## [1] "Q26: Any other thoughts you d like to share on EA sponsored mental health resources?"

ini=2
i=ini +26
colnames(A)[i] %>%
  gsub("\\."," ",.) %>%
  c("Q",i-ini,": ",.,"?") %>%
  paste(., collapse="") %>%
  gsub("   \\?","?",.) %>% 
  gsub(" \\?","?",.)-> titulo


## Country

ini=2
i=ini +27
colnames(A)[i] %>%
  gsub("\\."," ",.) %>%
  c("Q",i-ini,": ",.,"?") %>%
  paste(., collapse="") %>%
  gsub("   \\?","?",.) %>% 
  gsub(" \\?","?",.)-> titulo

## Age

ini=2
i=ini +28
colnames(A)[i] %>%
  gsub("\\."," ",.) %>%
  c("Q",i-ini,": ",.,"?") %>%
  paste(., collapse="") %>%
  gsub("   \\?","?",.) %>% 
  gsub(" \\?","?",.)-> titulo

(ggplot(data=A) +
    geom_bar(
      aes(x=Age),
      fill="royalblue4"
      )+
    labs(
      title=titulo,
      subtitle="n =303", 
      x="Age", 
      y="Relative frequency", 
      caption="@EA Mental Health Survey"
    )+
    theme(
      legend.title = element_blank(),
      plot.subtitle = element_text(hjust = 0.5), 
      plot.title = element_text(hjust = 0.5),
      axis.text.y = element_blank(),
      legend.position="bottom") 
) %>%
  ggsave(paste(c("Q", i-ini,".png"), collapse=""), plot =., units="in", width=12, height=5, dpi=800)