Перед тем, как оценивать эффективность терапии, хорошо бы оценить, что контрольная и экспериментальная группа не различаются до этапа терапии. В противном случае нельзя будет делать осмысленную интерпретацию изменений в терапевтической группе от терапии. Судя по результатам, группы значимо различаются только по шкале “нереагирование”.
surveys <- groups_homogeneity[, names(.SD), .SDcols = c(5:25)]
groups_diff <- cbind(
survey_scale = surveys,
groups_homogeneity[`прошли` == 1, list(median_t = unlist(lapply(.SD, function(x) median(x, na.rm = T)))), .SDcols = surveys],
groups_homogeneity[`прошли` == 0, list(median_c = unlist(lapply(.SD, function(x) median(x, na.rm = T)))), .SDcols = surveys],
groups_homogeneity[,
list(pvalue = unlist(
lapply(.SD,
function(x)
pvalue(wilcox_test(x ~ groups_homogeneity[, as.factor(`прошли`)], distribution="exact"))))),
.SDcols = surveys])graph1 <- melt(groups_diff[, list(survey_scale, median_t, median_c)], id.vars = "survey_scale", measure.vars = c("median_t", "median_c"))
graph1 <- graph1[, survey_scale := factor(survey_scale, levels = groups_diff[, survey_scale])]
p1 <- plot_ly(graph1[order(-survey_scale)],
x = ~survey_scale, y = ~value, color = ~variable, type = "bar") %>%
layout(title = "Медианные значения по шкалам<br>в экспериментальной и контрольной группах<br>до терапии",
xaxis = list(title = ""),
yaxis = list(title = ""),
margin = list(b = 200))
htmltools::tagList(p1)## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levelsСудя по результам тестов, терапия оказала значимое влияние на тревожность пациентов.
treatment_anxiety <- treatment_anxiety[, period := factor(period, levels = c("before", "after"))]
wilcoxsign_test(`Личностная тревога` ~ period, treatment_anxiety, distribution="exact")##
## Exact Wilcoxon-Pratt Signed-Rank Test
##
## data: y by x (pos, neg)
## stratified by block
## Z = 5.907, p-value = 2.842e-14
## alternative hypothesis: true mu is not equal to 0wilcoxsign_test(`Ситуационная тревога` ~ period, treatment_anxiety, distribution="exact")##
## Exact Wilcoxon-Pratt Signed-Rank Test
##
## data: y by x (pos, neg)
## stratified by block
## Z = 5.9078, p-value = 2.842e-14
## alternative hypothesis: true mu is not equal to 0graph2 <- melt(treatment_anxiety[, list(period, `Личностная тревога`, `Ситуационная тревога`)], id.vars = "period")
plot_ly(graph2, x = ~variable, y = ~value, color = ~period, type = 'box') %>%
layout(title = "Распередение уровня тревоги до и после терапии",
boxmode = "group",
xaxis = list(title = ""),
yaxis = list(title = ""))## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels## Warning: 'layout' objects don't have these attributes: 'boxmode'
## Valid attributes include:
## 'font', 'title', 'titlefont', 'autosize', 'width', 'height', 'margin', 'paper_bgcolor', 'plot_bgcolor', 'separators', 'hidesources', 'smith', 'showlegend', 'dragmode', 'hovermode', 'xaxis', 'yaxis', 'scene', 'geo', 'legend', 'annotations', 'shapes', 'images', 'updatemenus', 'ternary', 'mapbox', 'radialaxis', 'angularaxis', 'direction', 'orientation', 'barmode', 'bargap', 'mapType'Судя по результам тестов, терапия оказала значимое влияние на показатели осознанности у пациентов. К данным по параметру “Нереагирование” следует относиться с осторожностью, так как экспериментальная и контрольная группа по этой шкале значимо различались
treatment_mindfulness <- treatment_mindfulness[, period := factor(period, levels = c("before", "after"))]
lapply(treatment_mindfulness[, .SD, .SDcols = c(5:9)],
function(x)
wilcoxsign_test(x ~ treatment_mindfulness[, period], distribution="exact"))## $Наблюдение
##
## Exact Wilcoxon-Pratt Signed-Rank Test
##
## data: y by x (pos, neg)
## stratified by block
## Z = 5.9078, p-value = 2.842e-14
## alternative hypothesis: true mu is not equal to 0
##
##
## $Описание
##
## Exact Wilcoxon-Pratt Signed-Rank Test
##
## data: y by x (pos, neg)
## stratified by block
## Z = 5.9098, p-value = 2.842e-14
## alternative hypothesis: true mu is not equal to 0
##
##
## $`Осознанность действий`
##
## Exact Wilcoxon-Pratt Signed-Rank Test
##
## data: y by x (pos, neg)
## stratified by block
## Z = 5.9079, p-value = 2.842e-14
## alternative hypothesis: true mu is not equal to 0
##
##
## $Безоценочность
##
## Exact Wilcoxon-Pratt Signed-Rank Test
##
## data: y by x (pos, neg)
## stratified by block
## Z = 5.9081, p-value = 2.842e-14
## alternative hypothesis: true mu is not equal to 0
##
##
## $Нереагирование
##
## Exact Wilcoxon-Pratt Signed-Rank Test
##
## data: y by x (pos, neg)
## stratified by block
## Z = 5.9146, p-value = 2.842e-14
## alternative hypothesis: true mu is not equal to 0graph3 <- melt(treatment_mindfulness[, list(period, `Наблюдение`, `Описание`, `Осознанность действий`, `Безоценочность`, `Нереагирование`)], id.vars = "period")
plot_ly(graph3, x = ~variable, y = ~value, color = ~period, type = 'box') %>%
layout(title = "Распередение показателей осознанности<br>до и после терапии",
boxmode = "group",
xaxis = list(title = ""),
yaxis = list(title = ""),
margin = list(b = 200))## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels## Warning: 'layout' objects don't have these attributes: 'boxmode'
## Valid attributes include:
## 'font', 'title', 'titlefont', 'autosize', 'width', 'height', 'margin', 'paper_bgcolor', 'plot_bgcolor', 'separators', 'hidesources', 'smith', 'showlegend', 'dragmode', 'hovermode', 'xaxis', 'yaxis', 'scene', 'geo', 'legend', 'annotations', 'shapes', 'images', 'updatemenus', 'ternary', 'mapbox', 'radialaxis', 'angularaxis', 'direction', 'orientation', 'barmode', 'bargap', 'mapType'Судя по результам тестов, терапия оказала значимое влияния на эмоциональные схемы пациентов.
treatment_lihy <- treatment_lihy[, period := factor(period, levels = c("before", "after"))]
lapply(treatment_lihy[, .SD, .SDcols = c(5:18)],
function(x)
wilcoxsign_test(x ~ treatment_lihy[, period], distribution="exact"))## $`Инвалидация эмоций другими`
##
## Exact Wilcoxon-Pratt Signed-Rank Test
##
## data: y by x (pos, neg)
## stratified by block
## Z = 5.6074, p-value = 5.457e-12
## alternative hypothesis: true mu is not equal to 0
##
##
## $`Недостаточная осмысленность эмоций`
##
## Exact Wilcoxon-Pratt Signed-Rank Test
##
## data: y by x (pos, neg)
## stratified by block
## Z = 4.8931, p-value = 4.563e-08
## alternative hypothesis: true mu is not equal to 0
##
##
## $`Чувство вины за собственные эмоции`
##
## Exact Wilcoxon-Pratt Signed-Rank Test
##
## data: y by x (pos, neg)
## stratified by block
## Z = 4.1244, p-value = 1.173e-05
## alternative hypothesis: true mu is not equal to 0
##
##
## $`Упрощенное представление об эмоциях`
##
## Exact Wilcoxon-Pratt Signed-Rank Test
##
## data: y by x (pos, neg)
## stratified by block
## Z = 5.6106, p-value = 3.638e-12
## alternative hypothesis: true mu is not equal to 0
##
##
## $`Обесценивание эмоций`
##
## Exact Wilcoxon-Pratt Signed-Rank Test
##
## data: y by x (pos, neg)
## stratified by block
## Z = 5.3248, p-value = 5.057e-10
## alternative hypothesis: true mu is not equal to 0
##
##
## $`Страх потери контроля при переживании сильных эмоций`
##
## Exact Wilcoxon-Pratt Signed-Rank Test
##
## data: y by x (pos, neg)
## stratified by block
## Z = 4.7204, p-value = 1.894e-07
## alternative hypothesis: true mu is not equal to 0
##
##
## $`Эмоциональное оцепенение`
##
## Exact Wilcoxon-Pratt Signed-Rank Test
##
## data: y by x (pos, neg)
## stratified by block
## Z = 5.2324, p-value = 1.808e-09
## alternative hypothesis: true mu is not equal to 0
##
##
## $`Склонность к рационализации чувств`
##
## Exact Wilcoxon-Pratt Signed-Rank Test
##
## data: y by x (pos, neg)
## stratified by block
## Z = 5.6252, p-value = 3.638e-12
## alternative hypothesis: true mu is not equal to 0
##
##
## $`Прогнозируемая длительность эмоций`
##
## Exact Wilcoxon-Pratt Signed-Rank Test
##
## data: y by x (pos, neg)
## stratified by block
## Z = 5.4319, p-value = 1.601e-10
## alternative hypothesis: true mu is not equal to 0
##
##
## $`Недостаточная согласованность собственных эмоций с эмоциями других`
##
## Exact Wilcoxon-Pratt Signed-Rank Test
##
## data: y by x (pos, neg)
## stratified by block
## Z = 5.2059, p-value = 1.95e-09
## alternative hypothesis: true mu is not equal to 0
##
##
## $`Ингибирование собственных эмоций`
##
## Exact Wilcoxon-Pratt Signed-Rank Test
##
## data: y by x (pos, neg)
## stratified by block
## Z = 5.6243, p-value = 3.183e-12
## alternative hypothesis: true mu is not equal to 0
##
##
## $`Склонность к руминациям`
##
## Exact Wilcoxon-Pratt Signed-Rank Test
##
## data: y by x (pos, neg)
## stratified by block
## Z = 4.9833, p-value = 1.838e-08
## alternative hypothesis: true mu is not equal to 0
##
##
## $`Низкая эмоциональная экспрессивность`
##
## Exact Wilcoxon-Pratt Signed-Rank Test
##
## data: y by x (pos, neg)
## stratified by block
## Z = 5.5531, p-value = 1e-11
## alternative hypothesis: true mu is not equal to 0
##
##
## $`Обвинение других`
##
## Exact Wilcoxon-Pratt Signed-Rank Test
##
## data: y by x (pos, neg)
## stratified by block
## Z = 5.3156, p-value = 6.257e-10
## alternative hypothesis: true mu is not equal to 0graph4 <- melt(treatment_lihy[, .SD, .SDcols = c(4:18)], id.vars = "period")
plot_ly(graph4, x = ~variable, y = ~value, color = ~period, type = 'box') %>%
layout(title = "Распередение выраженности эмоциональных схем<br>до и после терапии",
boxmode = "group",
xaxis = list(title = ""),
yaxis = list(title = ""),
margin = list(b = 200))## Warning: Ignoring 56 observations## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels## Warning: 'layout' objects don't have these attributes: 'boxmode'
## Valid attributes include:
## 'font', 'title', 'titlefont', 'autosize', 'width', 'height', 'margin', 'paper_bgcolor', 'plot_bgcolor', 'separators', 'hidesources', 'smith', 'showlegend', 'dragmode', 'hovermode', 'xaxis', 'yaxis', 'scene', 'geo', 'legend', 'annotations', 'shapes', 'images', 'updatemenus', 'ternary', 'mapbox', 'radialaxis', 'angularaxis', 'direction', 'orientation', 'barmode', 'bargap', 'mapType'delta_activity <-
cbind(treatment_meta,
do.call(cbind, lapply(names(treatment_anxiety)[5:6],
function(x)
treatment_anxiety[period == "after", get(x)] - treatment_anxiety[period == "before", get(x)])),
do.call(cbind, lapply(names(treatment_mindfulness)[5:9],
function(x)
treatment_mindfulness[period == "after", get(x)] - treatment_mindfulness[period == "before", get(x)])),
do.call(cbind, lapply(names(treatment_lihy)[5:18],
function(x)
treatment_lihy[period == "after", get(x)] - treatment_lihy[period == "before", get(x)])))
setnames(delta_activity, c(names(treatment_meta), names(treatment_anxiety)[5:6], names(treatment_mindfulness)[5:9], names(treatment_lihy)[5:18]))
delta_mean <- delta_activity[, list(mean_delta = lapply(.SD, function(x) round(mean(x, na.rm = T), 3)),
median_delta = lapply(.SD, function(x) round(median(x, na.rm = T), 3))), .SDcols = c(8:28)]
graph5 <- melt(delta_activity[, .SD, .SDcols = c(8:28)])## Warning in melt.data.table(delta_activity[, .SD, .SDcols = c(8:28)]): To
## be consistent with reshape2's melt, id.vars and measure.vars are internally
## guessed when both are 'NULL'. All non-numeric/integer/logical type columns
## are conisdered id.vars, which in this case are columns []. Consider
## providing at least one of 'id' or 'measure' vars in future.## Warning in melt.data.table(delta_activity[, .SD, .SDcols = c(8:28)]):
## 'measure.vars' [Личностная тревога, Ситуационная тревога, Наблюдение,
## Описание, Осознанность действий, Безоценочность, Нереагирование,
## Инвалидация эмоций другими, Недостаточная осмысленность эмоций, Чувство
## вины за собственные эмоции, Упрощенное представление об эмоциях,
## Обесценивание эмоций, Страх потери контроля при переживании сильных
## эмоций, Эмоциональное оцепенение, Склонность к рационализации чувств,
## Прогнозируемая длительность эмоций, Недостаточная согласованность
## собственных эмоций с эмоциями других, Ингибирование собственных эмоций,
## Склонность к руминациям, Низкая эмоциональная экспрессивность, Обвинение
## других] are not all of the same type. By order of hierarchy, the molten
## data value column will be of type 'double'. All measure variables not of
## type 'double' will be coerced to. Check DETAILS in ?melt.data.table for
## more on coercion.plot_ly(graph5, x = ~variable, y = ~value, type = "box") %>%
layout(title = "Распередение разницы в параметрах<br>, вызванная терапией (после терапии - до терапии)",
boxmode = "group",
xaxis = list(title = ""),
yaxis = list(title = ""),
margin = list(b = 200))## Warning: Ignoring 28 observations## Warning: 'layout' objects don't have these attributes: 'boxmode'
## Valid attributes include:
## 'font', 'title', 'titlefont', 'autosize', 'width', 'height', 'margin', 'paper_bgcolor', 'plot_bgcolor', 'separators', 'hidesources', 'smith', 'showlegend', 'dragmode', 'hovermode', 'xaxis', 'yaxis', 'scene', 'geo', 'legend', 'annotations', 'shapes', 'images', 'updatemenus', 'ternary', 'mapbox', 'radialaxis', 'angularaxis', 'direction', 'orientation', 'barmode', 'bargap', 'mapType'Пробы не оказывают значимого влияния на изменения в значениях шкал.
predictors_trials <-
data.table(
cbind(survey_scale = names(delta_activity)[8:28],
do.call(rbind, lapply(names(delta_activity)[8:28], function(x) round(summary(lm(get(x) ~ `Кол-во проб`, delta_activity))$coefficients[2, ], 3)))))
kable(cbind(predictors_trials, delta_mean))| survey_scale | Estimate | Std. Error | t value | Pr(>|t|) | mean_delta | median_delta |
|---|---|---|---|---|---|---|
| Личностная тревога | -0.105 | 0.132 | -0.799 | 0.433 | -1.217 | -1 |
| Ситуационная тревога | -0.042 | 0.317 | -0.131 | 0.897 | -2.304 | -3 |
| Наблюдение | -0.006 | 0.174 | -0.033 | 0.974 | 0 | 0 |
| Описание | 0.098 | 0.091 | 1.071 | 0.296 | -0.565 | 1 |
| Осознанность действий | -0.061 | 0.121 | -0.502 | 0.621 | -1.087 | 0 |
| Безоценочность | 0.049 | 0.14 | 0.35 | 0.73 | 0.261 | 1 |
| Нереагирование | -0.045 | 0.114 | -0.389 | 0.701 | 1.565 | 2 |
| Инвалидация эмоций другими | 0.017 | 0.025 | 0.649 | 0.524 | 0.048 | 0 |
| Недостаточная осмысленность эмоций | -0.001 | 0.029 | -0.052 | 0.959 | 0.286 | 0 |
| Чувство вины за собственные эмоции | 0.019 | 0.031 | 0.623 | 0.541 | 0.048 | 0 |
| Упрощенное представление об эмоциях | 0.027 | 0.037 | 0.745 | 0.465 | 0.452 | 0.5 |
| Обесценивание эмоций | -0.025 | 0.033 | -0.762 | 0.455 | -0.19 | 0 |
| Страх потери контроля при переживании сильных эмоций | -0.068 | 0.029 | -2.383 | 0.028 | 0.238 | 0 |
| Эмоциональное оцепенение | -0.033 | 0.031 | -1.091 | 0.289 | -0.119 | 0 |
| Склонность к рационализации чувств | 0.014 | 0.033 | 0.426 | 0.675 | 0.119 | 0 |
| Прогнозируемая длительность эмоций | -0.017 | 0.022 | -0.775 | 0.448 | 0.286 | 0 |
| Недостаточная согласованность собственных эмоций с эмоциями других | 0.013 | 0.029 | 0.447 | 0.66 | -0.262 | 0 |
| Ингибирование собственных эмоций | 0.02 | 0.035 | 0.584 | 0.566 | 0.262 | 0 |
| Склонность к руминациям | -0.041 | 0.03 | -1.356 | 0.191 | 0.024 | 0 |
| Низкая эмоциональная экспрессивность | 0.032 | 0.031 | 1.033 | 0.315 | -0.333 | -0.5 |
| Обвинение других | -0.014 | 0.029 | -0.49 | 0.629 | -0.024 | 0 |
predictors_time <-
data.table(
cbind(survey_scale = names(delta_activity)[8:28],
do.call(rbind, lapply(names(delta_activity)[8:28], function(x) round(summary(lm(get(x) ~ `Общее время`, delta_activity))$coefficients[2, ], 3)))))
kable(cbind(predictors_time, delta_mean))| survey_scale | Estimate | Std. Error | t value | Pr(>|t|) | mean_delta | median_delta |
|---|---|---|---|---|---|---|
| Личностная тревога | -0.007 | 0.009 | -0.767 | 0.452 | -1.217 | -1 |
| Ситуационная тревога | -0.009 | 0.021 | -0.429 | 0.672 | -2.304 | -3 |
| Наблюдение | -0.004 | 0.011 | -0.327 | 0.747 | 0 | 0 |
| Описание | 0.007 | 0.006 | 1.233 | 0.231 | -0.565 | 1 |
| Осознанность действий | -0.001 | 0.008 | -0.136 | 0.893 | -1.087 | 0 |
| Безоценочность | 0.009 | 0.009 | 1.05 | 0.306 | 0.261 | 1 |
| Нереагирование | -0.005 | 0.007 | -0.674 | 0.507 | 1.565 | 2 |
| Инвалидация эмоций другими | 0.001 | 0.002 | 0.612 | 0.548 | 0.048 | 0 |
| Недостаточная осмысленность эмоций | 0 | 0.002 | -0.144 | 0.887 | 0.286 | 0 |
| Чувство вины за собственные эмоции | 0 | 0.002 | 0.174 | 0.864 | 0.048 | 0 |
| Упрощенное представление об эмоциях | 0.002 | 0.002 | 0.878 | 0.391 | 0.452 | 0.5 |
| Обесценивание эмоций | -0.001 | 0.002 | -0.679 | 0.506 | -0.19 | 0 |
| Страх потери контроля при переживании сильных эмоций | -0.005 | 0.002 | -2.919 | 0.009 | 0.238 | 0 |
| Эмоциональное оцепенение | -0.002 | 0.002 | -0.99 | 0.335 | -0.119 | 0 |
| Склонность к рационализации чувств | 0.001 | 0.002 | 0.354 | 0.727 | 0.119 | 0 |
| Прогнозируемая длительность эмоций | 0 | 0.001 | -0.31 | 0.76 | 0.286 | 0 |
| Недостаточная согласованность собственных эмоций с эмоциями других | 0.001 | 0.002 | 0.769 | 0.451 | -0.262 | 0 |
| Ингибирование собственных эмоций | 0.001 | 0.002 | 0.338 | 0.739 | 0.262 | 0 |
| Склонность к руминациям | -0.002 | 0.002 | -1.093 | 0.288 | 0.024 | 0 |
| Низкая эмоциональная экспрессивность | 0.001 | 0.002 | 0.707 | 0.488 | -0.333 | -0.5 |
| Обвинение других | -0.001 | 0.002 | -0.491 | 0.629 | -0.024 | 0 |
predictors_ind <-
data.table(
cbind(survey_scale = names(delta_activity)[8:28],
do.call(rbind, lapply(names(delta_activity)[8:28], function(x) round(summary(lm(get(x) ~ `Кол-во инд`, delta_activity))$coefficients[2, ], 3)))))
kable(cbind(predictors_ind, delta_mean))| survey_scale | Estimate | Std. Error | t value | Pr(>|t|) | mean_delta | median_delta |
|---|---|---|---|---|---|---|
| Личностная тревога | 2.086 | 1.006 | 2.074 | 0.051 | -1.217 | -1 |
| Ситуационная тревога | 0.167 | 2.622 | 0.064 | 0.95 | -2.304 | -3 |
| Наблюдение | 0.51 | 1.437 | 0.355 | 0.726 | 0 | 0 |
| Описание | 0.983 | 0.743 | 1.323 | 0.2 | -0.565 | 1 |
| Осознанность действий | 1.197 | 0.972 | 1.231 | 0.232 | -1.087 | 0 |
| Безоценочность | -0.532 | 1.154 | -0.461 | 0.649 | 0.261 | 1 |
| Нереагирование | -1.606 | 0.882 | -1.822 | 0.083 | 1.565 | 2 |
| Инвалидация эмоций другими | -0.262 | 0.197 | -1.329 | 0.199 | 0.048 | 0 |
| Недостаточная осмысленность эмоций | 0.05 | 0.229 | 0.217 | 0.831 | 0.286 | 0 |
| Чувство вины за собственные эмоции | 0.086 | 0.249 | 0.344 | 0.735 | 0.048 | 0 |
| Упрощенное представление об эмоциях | 0.059 | 0.297 | 0.2 | 0.844 | 0.452 | 0.5 |
| Обесценивание эмоций | 0.238 | 0.267 | 0.891 | 0.384 | -0.19 | 0 |
| Страх потери контроля при переживании сильных эмоций | 0.225 | 0.256 | 0.88 | 0.39 | 0.238 | 0 |
| Эмоциональное оцепенение | -0.272 | 0.245 | -1.11 | 0.281 | -0.119 | 0 |
| Склонность к рационализации чувств | -0.308 | 0.255 | -1.207 | 0.242 | 0.119 | 0 |
| Прогнозируемая длительность эмоций | -0.095 | 0.18 | -0.529 | 0.603 | 0.286 | 0 |
| Недостаточная согласованность собственных эмоций с эмоциями других | 0.138 | 0.234 | 0.589 | 0.562 | -0.262 | 0 |
| Ингибирование собственных эмоций | 0.268 | 0.274 | 0.978 | 0.34 | 0.262 | 0 |
| Склонность к руминациям | 0.362 | 0.241 | 1.503 | 0.149 | 0.024 | 0 |
| Низкая эмоциональная экспрессивность | 0.503 | 0.227 | 2.219 | 0.039 | -0.333 | -0.5 |
| Обвинение других | 0.305 | 0.225 | 1.355 | 0.191 | -0.024 | 0 |
predictors_group <-
data.table(
cbind(survey_scale = names(delta_activity)[8:28],
do.call(rbind, lapply(names(delta_activity)[8:28], function(x) round(summary(lm(get(x) ~ `Кол-во групп`, delta_activity))$coefficients[2, ], 3)))))
kable(cbind(predictors_group, delta_mean))| survey_scale | Estimate | Std. Error | t value | Pr(>|t|) | mean_delta | median_delta |
|---|---|---|---|---|---|---|
| Личностная тревога | 0.781 | 1.477 | 0.529 | 0.602 | -1.217 | -1 |
| Ситуационная тревога | -3.116 | 3.464 | -0.9 | 0.379 | -2.304 | -3 |
| Наблюдение | 0.924 | 1.93 | 0.479 | 0.637 | 0 | 0 |
| Описание | 1.312 | 1.001 | 1.311 | 0.204 | -0.565 | 1 |
| Осознанность действий | 1.75 | 1.301 | 1.345 | 0.193 | -1.087 | 0 |
| Безоценочность | -2.17 | 1.488 | -1.458 | 0.159 | 0.261 | 1 |
| Нереагирование | -0.388 | 1.274 | -0.305 | 0.764 | 1.565 | 2 |
| Инвалидация эмоций другими | -0.485 | 0.259 | -1.87 | 0.077 | 0.048 | 0 |
| Недостаточная осмысленность эмоций | 0.07 | 0.313 | 0.222 | 0.827 | 0.286 | 0 |
| Чувство вины за собственные эмоции | -0.106 | 0.34 | -0.31 | 0.76 | 0.048 | 0 |
| Упрощенное представление об эмоциях | 0.755 | 0.367 | 2.056 | 0.054 | 0.452 | 0.5 |
| Обесценивание эмоций | 0.26 | 0.367 | 0.71 | 0.486 | -0.19 | 0 |
| Страх потери контроля при переживании сильных эмоций | 0.175 | 0.354 | 0.495 | 0.626 | 0.238 | 0 |
| Эмоциональное оцепенение | -0.196 | 0.343 | -0.572 | 0.574 | -0.119 | 0 |
| Склонность к рационализации чувств | 0.358 | 0.352 | 1.017 | 0.322 | 0.119 | 0 |
| Прогнозируемая длительность эмоций | -0.093 | 0.247 | -0.376 | 0.711 | 0.286 | 0 |
| Недостаточная согласованность собственных эмоций с эмоциями других | -0.149 | 0.321 | -0.465 | 0.647 | -0.262 | 0 |
| Ингибирование собственных эмоций | -0.067 | 0.383 | -0.175 | 0.863 | 0.262 | 0 |
| Склонность к руминациям | 0.407 | 0.335 | 1.215 | 0.239 | 0.024 | 0 |
| Низкая эмоциональная экспрессивность | -0.126 | 0.346 | -0.365 | 0.719 | -0.333 | -0.5 |
| Обвинение других | 0.351 | 0.312 | 1.123 | 0.275 | -0.024 | 0 |
Посмотрим, различаются ли показатели после терпии у экспериментальной группы с показателями контрольной группы. Значимых различий не наблюдается, из чего можно сделать вывод, что терапия оказала значимый эффект, который, тем не менее, оказался весьма слабым и недостаточным для сильного улучшения состояния по сравнению с контрольной группой.
groups_diff2 <- cbind(
survey_scale = surveys,
groups_homogeneity_a[`прошли` == 1, list(median_t = unlist(lapply(.SD, function(x) median(x, na.rm = T)))), .SDcols = surveys],
groups_homogeneity_a[`прошли` == 0, list(median_c = unlist(lapply(.SD, function(x) median(x, na.rm = T)))), .SDcols = surveys],
groups_homogeneity_a[,
list(pvalue = unlist(
lapply(.SD,
function(x)
pvalue(wilcox_test(x ~ groups_homogeneity_a[, as.factor(`прошли`)], distribution="exact"))))),
.SDcols = surveys])
graph6 <- melt(groups_diff2[, list(survey_scale, median_t, median_c)], id.vars = "survey_scale", measure.vars = c("median_t", "median_c"))
graph6 <- graph6[, survey_scale := factor(survey_scale, levels = groups_diff2[, survey_scale])]
plot_ly(graph6[order(-survey_scale)],
x = ~survey_scale, y = ~value, color = ~variable, type = "bar") %>%
layout(title = "Медианные значения по шкалам<br>в экспериментальной и контрольной группах<br>после терапии",
xaxis = list(title = ""),
yaxis = list(title = ""),
margin = list(b = 200))## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels