Calculate SSM parameters with bootstrapped confidence intervals for a variety of different analysis types. Depending on what arguments are supplied, either mean-based or correlation-based analyses will be performed, one or more groups will be used to stratify the data, and contrasts between groups or measures will be calculated.
ssm_analyze(.data, scales, angles, measures = NULL, grouping = NULL, contrast = "none", boots = 2000, interval = 0.95, listwise = TRUE)
| .data | Required. A data frame containing at least circumplex scales. |
|---|---|
| scales | Required. The variable names or column numbers for the
variables in |
| angles | Required. A numeric vector containing the angular displacement
of each circumplex scale included in |
| measures | Optional. The variable names or column numbers for one or
more variables in |
| grouping | Optional. The variable name or column number for the variable
in |
| contrast | Optional. A string indicating what type of contrast to run. Current options are "none" for no contrast, "model" to find SSM parameters for the difference scores, or "test" to find the difference between the SSM parameters. Note that only two groups or measures can be contrasted at a time (default = "none"). |
| boots | Optional. A single positive integer indicating how many bootstrap resamples to use when estimating the confidence intervals (default = 2000). |
| interval | Optional. A single positive number between 0 and 1 (exclusive) that indicates what confidence level to use when estimating the confidence intervals (default = 0.95). |
| listwise | Optional. A logical indicating whether missing values should be handled by listwise deletion (TRUE) or pairwise deletion (FALSE). Note that pairwise deletion may result in different missing data patterns in each bootstrap resample and is slower to compute (default = TRUE). |
A list containing the results and description of the analysis.
A tibble with the SSM parameter estimates
A list with the number of bootstrap resamples (boots), the confidence interval percentage level (interval), and the angular displacement of scales (angles)
A language object containing the function call that created this object
A tibble containing the mean scale scores
A string indicating what type of SSM analysis was done
Other ssm functions: ssm_append,
ssm_plot, ssm_table
# Load example data data("jz2017") # Single-group mean-based SSM ssm_analyze(jz2017, scales = PA:NO, angles = octants())#> Call: #> ssm_analyze(.data = jz2017, scales = PA:NO, angles = octants()) #> #> Profile [All]: #> Estimate Lower CI Upper CI #> Elevation 0.917 0.889 0.946 #> X-Value 0.351 0.326 0.379 #> Y-Value -0.252 -0.280 -0.222 #> Amplitude 0.432 0.405 0.462 #> Displacement 324.292 321.063 327.900 #> Model Fit 0.878 #># Single-group correlation-based SSM ssm_analyze(jz2017, scales = PA:NO, angles = octants(), measures = c(NARPD, ASPD))#> Call: #> ssm_analyze(.data = jz2017, scales = PA:NO, angles = octants(), #> measures = c(NARPD, ASPD)) #> #> Profile [NARPD]: #> Estimate Lower CI Upper CI #> Elevation 0.202 0.171 0.235 #> X-Value -0.062 -0.095 -0.028 #> Y-Value 0.179 0.144 0.213 #> Amplitude 0.189 0.153 0.225 #> Displacement 108.967 99.207 118.622 #> Model Fit 0.957 #> #> Profile [ASPD]: #> Estimate Lower CI Upper CI #> Elevation 0.124 0.089 0.159 #> X-Value -0.099 -0.132 -0.063 #> Y-Value 0.203 0.168 0.240 #> Amplitude 0.226 0.190 0.264 #> Displacement 115.927 107.038 124.248 #> Model Fit 0.964 #># Multiple-group mean-based SSM ssm_analyze(jz2017, scales = PA:NO, angles = octants(), grouping = Gender)#> Call: #> ssm_analyze(.data = jz2017, scales = PA:NO, angles = octants(), #> grouping = Gender) #> #> Profile [Female]: #> Estimate Lower CI Upper CI #> Elevation 0.946 0.908 0.983 #> X-Value 0.459 0.421 0.500 #> Y-Value -0.310 -0.354 -0.268 #> Amplitude 0.554 0.510 0.600 #> Displacement 325.963 322.176 329.821 #> Model Fit 0.889 #> #> Profile [Male]: #> Estimate Lower CI Upper CI #> Elevation 0.884 0.840 0.925 #> X-Value 0.227 0.192 0.262 #> Y-Value -0.186 -0.224 -0.150 #> Amplitude 0.294 0.260 0.331 #> Displacement 320.685 313.718 327.851 #> Model Fit 0.824 #># Multiple-group mean-based SSM with contrast ssm_analyze(jz2017, scales = PA:NO, angles = octants(), grouping = Gender, contrast = "model")#> Call: #> ssm_analyze(.data = jz2017, scales = PA:NO, angles = octants(), #> grouping = Gender, contrast = "model") #> #> Contrast [Male - Female]: #> Estimate Lower CI Upper CI #> Elevation -0.062 -0.120 -0.004 #> X-Value -0.232 -0.282 -0.178 #> Y-Value 0.124 0.065 0.184 #> Amplitude 0.263 0.209 0.322 #> Displacement 151.858 140.282 164.263 #> Model Fit 0.855 #># Single-group correlation-based SSM with contrast ssm_analyze(jz2017, scales = PA:NO, angles = octants(), measures = c(NARPD, ASPD), contrast = "test")#> Call: #> ssm_analyze(.data = jz2017, scales = PA:NO, angles = octants(), #> measures = c(NARPD, ASPD), contrast = "test") #> #> Contrast [ASPD - NARPD]: #> Estimate Lower CI Upper CI #> Elevation -0.079 -0.118 -0.043 #> X-Value -0.037 -0.076 -0.001 #> Y-Value 0.024 -0.013 0.062 #> Amplitude 0.037 0.000 0.077 #> Displacement 6.960 -2.835 16.963 #> Model Fit 0.007 #># Multiple-group correlation-based SSM ssm_analyze(jz2017, scales = PA:NO, angles = octants(), measures = NARPD, grouping = Gender)#> Call: #> ssm_analyze(.data = jz2017, scales = PA:NO, angles = octants(), #> measures = NARPD, grouping = Gender) #> #> Profile [Female_NARPD]: #> Estimate Lower CI Upper CI #> Elevation 0.172 0.126 0.219 #> X-Value -0.080 -0.129 -0.036 #> Y-Value 0.202 0.151 0.250 #> Amplitude 0.217 0.164 0.272 #> Displacement 111.669 100.392 123.891 #> Model Fit 0.972 #> #> Profile [Male_NARPD]: #> Estimate Lower CI Upper CI #> Elevation 0.244 0.191 0.295 #> X-Value -0.029 -0.075 0.015 #> Y-Value 0.146 0.101 0.192 #> Amplitude 0.149 0.105 0.195 #> Displacement 101.248 84.307 118.790 #> Model Fit 0.902 #># Multiple-group correlation-based SSM with contrast ssm_analyze(jz2017, scales = PA:NO, angles = octants(), measures = NARPD, grouping = Gender, contrast = "test")#> Call: #> ssm_analyze(.data = jz2017, scales = PA:NO, angles = octants(), #> measures = NARPD, grouping = Gender, contrast = "test") #> #> Contrast [NARPD: Male - Female]: #> Estimate Lower CI Upper CI #> Elevation 0.072 0.006 0.142 #> X-Value 0.051 -0.012 0.113 #> Y-Value -0.056 -0.125 0.008 #> Amplitude -0.068 -0.139 -0.002 #> Displacement -10.421 -30.970 10.708 #> Model Fit -0.071 #>