具有着色级别的Bland-Altman的Data.Frame格式

时间:2019-06-25 10:04:11

标签: r

我正在努力寻找一种数据表示的解决方案,以绘制Bland-Altman图(带有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和Smartphone是两个Smartphone(Time_Smartphone)测量的因素。

由于Bland-Altman图偏爱两个向量(每个方法测量一个),因此我对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图进行着色的级别(视频,智能手机),以便使视频和智能手机的测量结果可区分。

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软件包在诸如Bland-Altman图的图中按因子进行着色?

enter image description here

1 个答案:

答案 0 :(得分:0)

几乎在已经发布的两种情况下,我手动解决了绘图问题。您可以使用首选的绘图功能(blandaltmanleh或blandr)绘制第一组数据,然后添加点。如果您对自制的BA线不满意,则最后一种选择是计算所有数据的BA统计信息,手动绘制边距和线,然后以两步/颜色绘制点。