getCSC: Function to return Jacobson (et al.) Clinically Significant...
In cpsyctc/CECPfuns: Package of Utility Functions for Psychological Therapies, Mental Health and Well-being Work (Created by Chris Evans and Clara Paz)

Description Usage Arguments Value Background History/development log Author(s) References See Also Examples

Fairly trivial function to compute CSC, mainly useful if using group_by() piping to find CSCs for various groups/subsets in your data.

1	getCSC(formula1, data)

`formula1`	defines variables to use as score ~ grp where score has scores and grp has grouping
`data`	is the data, typically cur_data() when called in tidyverse pipe

A single numeric value for the CSC

Like the CSC, the RCI comes out of the classic paper Jacobson, Follette & Revenstorf (1984). The thrust of the paper was about trying to bridge the gap between (quantitative) researchers who then, and still, tend to think in terms of aggregated data about change in therapy, and clinicians who tend to think about individual client's change.

The authors came up with two indices: the RCI, and the one here: the CSC or Clinically Significant Change (cut-off).

The authors defined three methods to set a criterion between "clinical" and "non-clinical" scores. Generally we prefer "help-seeking" and "non-help-seeking" to "clinical" and "non-clinical" to avoid automatically jumping into a disease/medical model.

For a measure cued so that higher scores indicate more problems the methods were as follows.

Method a: the score 2 SD below the mean in the help-seeking group.
Method b: the score 2 SD above the mean in the non-help-seeking group.
Method c: the score between the means of the two groups defined so that, for Gaussian distributions of scores, the cutting point would misclassify the same proportions of each group. This is often spoken of as the point halfway between the two means but the definition is actually more subtle than that and will only be halfway between the means where the two groups have the same score SD. Where the groups have different SD the formula for method c is this.

\loadmathjax \mjdeqn\frac

SD_HS*M_notHS + SD_notHS*M_HSSD_HS+SD_notHS

(with SD for Standard Deviation (doh!) and M for Mean)

Started before 5.iv.21

Chris Evans

Jacobson, N. S., Follette, W. C., & Revenstorf, D. (1984). Psychotherapy outcome research: Methods for reporting variability and evaluating clinical significance. Behavior Therapy, 15, 336–352.
Evans, C., Margison, F., & Barkham, M. (1998). The contribution of reliable and clinically significant change methods to evidence-based mental health. Evidence Based Mental Health, 1, 70–72. https://doi.org/0.1136/ebmh.1.3.70

getBootCICSC for CSC and bootstrap CI around it

Other RCSC functions: classifyScoresVectorByRCI(), getBootCICSC(), getRCIfromSDandAlpha()

## Not run: 
# example uses very crude simulated data but illustrates utility of being able to pipe grouped
# data into getCSC()
library(tidyverse)
set.seed(12345)
n <- 250 # total sample size
scores <- rnorm(n) # get some random numbers
# now random allocation to help-seeking or non-help-seeking samples
grp <- sample(c("HS", "notHS"), n, replace = TRUE)
# and random allocation by gender
gender <- sample(c("F", "M"), n, replace = TRUE)
list(scores = scores,
    grp = grp,
    gender = gender) %>%
 as_tibble() %>%
 mutate(scores = if_else(gender == "F", scores + .13, scores), # make women .13 higher scoring
        ### then make help-seeking group score 1.1 higher
        scores = if_else(grp == "HS", scores + 1.1, scores)) -> tibDat
### now get CSC overall
tibDat %>%
 summarise(CSC = getCSC(scores ~ grp, cur_data()))

### get CSC by gender
tibDat %>%
 group_by(gender) %>%
 summarise(CSC = getCSC(scores ~ grp, cur_data()))

## End(Not run)