Я изо всех сил пытаюсь найти решение для представления данных для рисования графика Бланда-Альтмана (с помощью blandr / epade или других пакетов), где я могу по-разному раскрасить разные методы.
Я хотел бы сравнить измерения, сделанные на смартфоне-акселерометре, с видеоанализом. Поэтому я прочитал все данные и представил их в длинном формате.
Date Participant Rater Machine_ASVZ Machine Set Contraction_Mode Time_Video Smartphone Time_Smartphone
3 2019.01.26 CV CV A1 LEG PRESS 1 Con1 1340 1 1650
4 2019.01.26 CV CV A1 LEG PRESS 1 Ecc1 2580 1 2450
5 2019.01.26 CV CV A1 LEG PRESS 1 Con2 1160 1 1350
6 2019.01.26 CV CV A1 LEG PRESS 1 Ecc2 2300 1 2100
7 2019.01.26 CV CV A1 LEG PRESS 1 Con3 1220 1 1400
8 2019.01.26 CV CV A1 LEG PRESS 1 Ecc3 2120 1 1950
Оценка - это коэффициент двух Рейтингов, которые независимо анализируются. Time_Video и Смартфон - это коэффициент двух измерений смартфона (Time_Smartphone).
Поскольку графики Бланда-Альтмана отдают предпочтение двум векторам (по одному для каждого метода измерений), я приведу и мутирую data.frame, чтобы вычислить среднее значение для анализа видео с двумя оценщиками (Time_Video) и двух измерений на смартфоне (Time_Smartphone) :
Date Participant Machine_ASVZ Machine Set Contraction_Mode Time_Video_CV Time_Video_GRFD Time_Smartphone_1 Time_Smartphone_3 Rater_Mean Smartphone_Mean
1 2019.01.26 CV A4 LEG EXTENSION 1 Con1 1160 1140 1300 1300 1150 1300
2 2019.01.26 CV A4 LEG EXTENSION 1 Ecc1 2340 2320 2250 2300 2330 2275
3 2019.01.26 CV A4 LEG EXTENSION 1 Con2 1160 1160 1200 1150 1160 1175
4 2019.01.26 CV A4 LEG EXTENSION 1 Ecc2 2340 2360 2300 2300 2350 2300
5 2019.01.26 CV A4 LEG EXTENSION 1 Con3 1220 1200 1250 1250 1210 1250
6 2019.01.26 CV A4 LEG EXTENSION 1 Ecc3 2020 2000 2000 2000 2010 2000
>
Теперь я не знаю, как вставить уровни (Видео, Смартфон) для раскрашивания сюжета Бланда-Альтмана так, чтобы измерения видео и смартфона были различимы.
bland_altman_plot<- blandr.draw(as.numeric(unlist(dplyr::select(filter(melt_subset, Source=="Rater_Mean"),"Time"))), as.numeric(unlist(dplyr::select(filter(melt_subset, Source=="Smartphone_Mean"),"Time"))), method1name = "Method 1",
method2name = "Method 2",
plotTitle = paste("Bland-Altman-Plot ",i, sep=""),
sig.level = 0.95, LoA.mode = 1, annotate = FALSE, ciDisplay = TRUE,
ciShading = FALSE, normalLow = FALSE, normalHigh = FALSE,
lowest_y_axis = FALSE, highest_y_axis = FALSE, point_size = 2,
overlapping = FALSE, plotter = "ggplot2", x.plot.mode = "means",
y.plot.mode = "difference", plotProportionalBias = TRUE,
plotProportionalBias.se = TRUE, assume.differences.are.normal = TRUE)
Как можно переносить уровни в широком формате, чтобы их можно было раскрасить по факторам на графике, таком как график Бланда-Альтмана, с использованием пакета blandr?
