-
Notifications
You must be signed in to change notification settings - Fork 1
/
Prescribing for anxiety depressions psychosis.R
37 lines (29 loc) · 1.72 KB
/
Prescribing for anxiety depressions psychosis.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
# ScotPHO indicator: population prescribed drugs for anxiety/depression/psychosis
# Part 1 - Prepare basefile
# Part 2 - Run analysis functions
###############################################.
## Packages/Filepaths/Functions ----
###############################################.
source("1.indicator_analysis.R") #Normal indicator functions
source("2.deprivation_analysis.R") # deprivation function
###############################################.
## Part 1 - Prepare basefile ----
###############################################.
#Data comes from prescribing team
presc_anx <- read_csv(file=paste0(data_folder, 'Received Data/prescrib_anx_depression_ScotPHO Indicator by Datazone 2010 to 2021.csv')) %>%
group_by(datazone2011, year) %>% summarise(numerator = sum(patients)) %>%
ungroup() %>% rename(datazone = datazone2011)
saveRDS(presc_anx, file=paste0(data_folder, 'Prepared Data/prescriptions_anxiety_raw.rds'))
saveRDS(presc_anx, file=paste0(data_folder, 'Prepared Data/prescriptions_anxiety_depr_raw.rds'))
###############################################.
## Part 2 - Run analysis functions ----
###############################################.
analyze_first(filename = "prescriptions_anxiety", geography = "datazone11", pop= "DZ11_pop_allages",
measure = "percent", yearstart = 2010, yearend = 2021, time_agg = 1)
analyze_second(filename = "prescriptions_anxiety", measure = "percent", time_agg = 1,
ind_id = 20401, year_type = "financial")
#Deprivation analysis function
analyze_deprivation(filename="prescriptions_anxiety_depr", measure="percent", time_agg=1,
yearstart= 2014, yearend=2021, year_type = "financial",
pop= "depr_pop_allages", ind_id = 20401)
##END