Dear colleagues,
Question: When using graphical methods to assess agreement between two scales (concord/batplot commands) Is there any command to generate a list of the observations that lie outside the limits of agreement?
I have a cohort of patients (n=366) who had two questionnaires administered, the Barthel Index a measure of functional independence, and the EQ-5D-3L utiltiy score a health-related quality of life score.
The data has been log transformed and converted to z scores.
I am attempting to measure the agreement between the two instruments including using graphical methods to display agreement, using concord and batplot commands. The commands work well and graphs are shown below.
I would now like to investigate the outliers, the observations that lie outside the limits of agreement. Is there any command to generate a list of the observations that lie outside the limits of agreement? I am currently labelling outliers and recording study id manually, which is a little time consuming and inelegant.
Grateful for any suggestions,
concord z_w7 z_y7, loa mlabel(record_id)
batplot z_w7 z_y7 , info moptions(mlabp(9)) sc(jitter(4)) valabel (record_id)
Question: When using graphical methods to assess agreement between two scales (concord/batplot commands) Is there any command to generate a list of the observations that lie outside the limits of agreement?
I have a cohort of patients (n=366) who had two questionnaires administered, the Barthel Index a measure of functional independence, and the EQ-5D-3L utiltiy score a health-related quality of life score.
The data has been log transformed and converted to z scores.
I am attempting to measure the agreement between the two instruments including using graphical methods to display agreement, using concord and batplot commands. The commands work well and graphs are shown below.
I would now like to investigate the outliers, the observations that lie outside the limits of agreement. Is there any command to generate a list of the observations that lie outside the limits of agreement? I am currently labelling outliers and recording study id manually, which is a little time consuming and inelegant.
Grateful for any suggestions,
concord z_w7 z_y7, loa mlabel(record_id)
batplot z_w7 z_y7 , info moptions(mlabp(9)) sc(jitter(4)) valabel (record_id)
Comment