Srovnání dvou metod měření (Passing-Bablok, Deming)
Stručně
Dvě metody měření změřily stejné vzorky. Srovnání metod v otázce „Shodují se hodnocení?“ (u přesně dvou číselných sloupců, vedle Bland-Altmanova diagramu) ukazuje, zda se liší pevnou hodnotou (konstanta není 0), nebo činitelem (sklon není 1). Oba postupy, Passing-Bablok i Deming, považují obě měření za zatížená chybou; obyčejná regrese to nedělá a dala by příliš plochý sklon.
Passing-Bablok
- Počítá jako balíček R mcr (verze 1.3.3.1): medián sklonů všech dvojic bodů, posunutý podle pravidla Passinga a Babloka, úhlová míra, analytický interval. Je odolná vůči odlehlým hodnotám a nevyžaduje normální rozdělení.
- Je zablokováno u záporných hodnot (jako v mcr; posunutí dat by změnilo konstantu a neprovádí se) a u nekladné souvislosti. Potřebuje spojité hodnoty: s méně než 20 různými hodnotami se objeví poznámka, interval je hrubý. Porovnává každou dvojici s každou jinou, takže je povoleno nejvýše 4000 vzorků.
- Cusum test linearity ukazuje jen testovou statistiku h a graf reziduí, bez p-hodnoty: vzorec p-hodnoty není v originále zdokumentován. Nenápadná hodnota není důkazem linearity.
Deming a poměr λ
- λ je rozptyl chyby měření na vodorovné ose dělený rozptylem chyby na svislé ose. Výchozí je λ = 1: obě metody měří se stejně velkou chybou. Směr je v každém programu jiný; Statable následuje balíček R mcr.
- λ lze určit z opakovaných měření (dvě měření na vzorek a metodu), nikoli z dat samotného srovnání. Chybné λ zkresluje přímku; tabulka ukazuje sklon pro λ = 0,5, 1 a 2, aby bylo vidět, jak moc λ záleží. Interval pochází z jackknife podle Linneta (1990).
Čtení výsledku
- Neobsahuje-li interval konstanty 0, je tu náznak konstantní odchylky; neobsahuje-li interval sklonu 1, náznak proporcionální. Je to jen náznak, nikoli důkaz, že jsou metody rovnocenné nebo se liší.
- Obrázek ukazuje vzorky, přímku identity y = x (čárkovaně) a obě přímky. Přímka platí jen v rozsahu naměřených hodnot. Bland-Altman ukazuje rozdíly vůči průměru, Passing-Bablok a Deming přímku; obojí patří k sobě.
- S váhou případu je srovnání metod zablokováno: každý vzorek se počítá jednou.
Pro odborníky
Passing-Bablok krok za krokem sleduje calc_PaBa z mcr (tolerance rozdílu 1e-12, posun od úhlů do −π/4, pravidlo pořadí, C zaokrouhleno na celé číslo) a shoduje se s mcr i s referenčními sadami programu „Calanas“ z testu balíčku na 1e-9. Deming: Linnet (1990), jackknife s t(n − 2), jako mcr. Zdroje: Passing a Bablok (1983), J Clin Chem Clin Biochem 21, 709-720; Linnet (1990), Statistics in Medicine 9, 1463-1473.