-
Notifications
You must be signed in to change notification settings - Fork 1
/
New_Mexico_SGP.R
32 lines (27 loc) · 1.02 KB
/
New_Mexico_SGP.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
###################################################################
###
### R script for 2023 SGP analyses for NM MSSA
###
###################################################################:w
### Load packages
require(SGP)
require(data.table)
### Load data
load("Data/New_Mexico_Data_LONG.Rdata")
### Parameters
parallel.config <- list(BACKEND="PARALLEL", WORKERS=list(PERCENTILES=4, BASELINE_PERCENTILES=4, PROJECTIONS=4, LAGGED_PROJECTIONS=4, SGP_SCALE_SCORE_TARGETS=4))
### Run abcSGP
New_Mexico_SGP <- abcSGP(
sgp_object = New_Mexico_Data_LONG,
steps = c("prepareSGP", "analyzeSGP", "combineSGP", "summarizeSGP", "visualizeSGP", "outputSGP"),
sgp.percentiles = TRUE,
sgp.projections = TRUE,
sgp.projections.lagged = TRUE,
sgp.percentiles.baseline = FALSE,
sgp.projections.baseline = FALSE,
sgp.projections.lagged.baseline = FALSE,
sgPlot.demo.report=TRUE,
parallel.config = parallel.config
)
### Save results
save(New_Mexico_SGP, file="Data/New_Mexico_SGP.Rdata")