A regresszión alapuló költségmodellek egyik gyakori hátránya, hogy a betegút egyes szakaszai alatt felmerülő költségeket csupán aggregált formában modellezik, ha az egészségügyi mutatók becslésére szolgáló betegutakat csak korlátozottan, egy-egy időablakban figyeljük meg. Mindez a gazdasági modellként történő használat korlátait jelentheti, ugyanis az átalányköltség már nem veszi figyelembe a betegenkénti eltérések okozta varianciát, s jelentős eltéréseket okozhat a való- ságban tapasztalt és szimulált költségek között. A költségmodelleket a betegutak komplex leírására alkalmas Cox-féle arányos hazárd statisztikai modellel együtt fejleszthetjük tovább, ami az egyes szimulált betegút szeletek pontosabb költségtételekkel való társításával segítheti a kockázatok finanszírozási nézőpont szerinti értékelését. A Cox-féle arányos hazárd modell azon mó- dosításával, hogy az események egymás utáni bekövetkeztekor nem csak a közöttük lévő időtávot, hanem a költségkülönbségeket is vizsgáljuk, módunk van az egyes eseményeket követő költségek nagyságát és valószínűségét megbecsülni. A kiadások eloszlásának, várható értékének és szóródásának pontosabb előrejelzése fokozza az egyes befolyásoló tényezők eltéréséből származó különbségek kimutathatóságát.
Angol absztrakt:
A common problem with cost regression models applied to aggregated health expenditures is that they cannot reflect on real and particular patient-path’s events as we are restricted to use a limited time frame approach for observation. Fitting only average-like costs to patient-path pieces instead of modelling detailed cost elements imposes restrictions on the applicability of these models and may cause large differences between the observed and simulated costs. As an alternative approach we propose a procedure analogous to the Cox’s proportional hazard model that is able to describe complex patient-paths. If simulated patientpaths pieces are matched with cost components then such cost models can improve the precision of the corresponding financing model. The magnitude and probability of the cost elements occurring after certain events can be estimated by a modified Cox’s proportional hazard model that is scaled according to cost differences instead of times elapsed. More accurate prediction of the distribution, expected value and standard deviation of the costs makes it possible to describe the differences – caused by modifying factors – more accurately
[1]53:419-434.[1] Cox dR (1972): Regression models and life-tables (with discussion). J. Royal stat. soc. B, 34:187-220
[2] Kalbfleisch Jd, Prentice RL (2002): the statistical Analysis of Failure time data. 2nd Edition. Wiley.
[3] therneau tM,Grambsch PM (2000): Modeling survival data. Extending the Cox Model. springer.
[4] Ross sM (1995) stochastic Processes. Wiley.
[5] Zhao H, Cheng Y, Bang H (2011): some insight on censored cost estimators. statistics in Medicine 30:2381-2388.
[6] Lin dY, Feuer EJ, Etzioni R, Wax Y (1997): Estimating Medical Costs from Incomplete Follow-Up data. Biometrics 53:419-434.
A cikket sikeresen a könyvepolcára helyeztük!
Tisztelt Felhasználónk!
A cikket a könyvespolcára helyeztük. A későbbiekben
bármikor elérheti a cikket a könyvespolcán található listáról.
A cikk megtekintéséhez onine regisztráció szükséges!
Tisztelt Látogató!
Az Ön által megtekinteni kívánt cikk az IMEONLINE cikkadatbázisához tartozik, melynek olvasása online regisztrációhoz kötött.
A regisztrálást követően fogja tudni megtekinteni a cikk tartalmát!
A megadott cikk nem elérhető!
Tisztelt Felhasználónk!
Az Ön által megtekinteni kívánt cikk nem elérhető a rendszerben!
A megadott cikk nem elérhető!
Tisztelt Felhasználónk!
Az Ön által megtekinteni kívánt cikk nem elérhető a rendszerben!
Sikeresen szavazott a cikkre!
Tisztelt Felhasználónk!
Köszönjük a szavazatát!
A szavazás nem sikerült!
Tisztelt Felhasználónk!
Ön már szavazott az adott cikkre!
Cikk megtekintése
Tisztelt Felhasználónk!
A cikk több nyelven is elérhető! Kérjük, adja meg, hogy melyik nyelven kívánja megtekinteni az adott cikket!
Cikk megtekintésének megerősítése!
Tisztelt Felhasználónk!
Az Ön által megtekintetni kívánt cikk tartalma fizetős szolgáltatás.
A megtekinteni kívánt cikket automatikusan hozzáadjuk a könyvespolcához!
A cikket bármikor elérheti a könyvespolcok menüpontról is!
egészség-gAzdAságtAn MódsZERtAn Eseményt követő kiadások eloszlásának és várható értékének modellezése Cox-féle arányos hazárd modellel Lang Zsolt, Rakonczai Pál, Lazányi István, Lohn Eszter, Bacskai Miklós, HealthWare Kft. A regresszión alapuló költségmodellek egyik gyakori hátránya, hogy a betegút egyes szakaszai alatt felmerülő költségeket csupán aggregált formában modellezik, ha az egészségügyi mutatók becslésére szolgáló betegutakat csak korlátozottan, egy-egy időablakban figyeljük meg. Mindez a gazdasági modellként történő használat korlátait jelentheti, ugyanis az átalányköltség már nem veszi figyelembe a betegenkénti eltérések okozta varianciát, s jelentős eltéréseket okozhat a valóságban tapasztalt és szimulált költségek között. A költségmodelleket a betegutak komplex leírására alkalmas Cox-féle arányos hazárd statisztikai modellel együtt fejleszthetjük tovább, ami az egyes szimulált betegút szeletek pontosabb költségtételekkel való társításával segítheti a kockázatok finanszírozási nézőpont szerinti értékelését. A Cox-féle arányos hazárd modell azon módosításával, hogy az események egymás utáni bekövetkeztekor nem csak a közöttük lévő időtávot, hanem a költségkülönbségeket is vizsgáljuk, módunk van az egyes eseményeket követő költségek nagyságát és valószínűségét megbecsülni. A kiadások eloszlásának, várható értékének és szóródásának pontosabb előrejelzése fokozza az egyes befolyásoló tényezők eltéréséből származó különbségek kimutathatóságát. A common problem with cost regression models applied to aggregated health expenditures is that they cannot reflect on real and particular patient-path’s events as we are restricted to use a limited time frame approach for observation. Fitting only average-like costs to patient-path pieces instead of modelling detailed cost elements imposes restrictions on the applicability of these models and may cause large differences between the observed and simulated costs. As an alternative approach we propose a procedure analogous to the Cox’s proportional hazard model that is able to describe complex patient-paths. If simulated patientpaths pieces are matched with cost components then such cost models can improve the precision of the corresponding financing model. The magnitude and probability of the cost elements occurring after certain events can be estimated by a modified Cox’s proportional hazard model that is scaled according to cost differences instead of times elapsed. More accurate prediction of the distribution, expected value and standard deviation of the costs makes it possible to describe the differences – caused by modifying factors – more accurately. 34 IMe XI. évfolyAM 1. száM 2012. feBruár Bevezetés Az egyre pontosabb költségbecsléseket használó egészség-gazdaságtani modellek kiemelt szerepet tölthetnek be a finanszírozói döntések során. A szolgáltatásvásárlói és biztosítói funkciók fejlesztésében az elvárt eredményesség mérése mellett a kiadások várható értékének megfelelő közelítése is az alapvető feladatok közé tartozik. Különösen a finanszírozási kockázatok vizsgálatakor értékelődnek fel a költség-elemzések és modellek azon formái, melyek képesek az egyéni betegúthoz tartozó kiadási variancia részletesebb kifejezésére és bemutatására is. A megbetegedési eseményeket időben egymásra épülő terápiás szakaszok, különböző egészségügyi ellátások követik. Ezek az események, történések betegről betegre nagy változatosságot mutatnak, részleteikben már nehezen tipizálhatók, illetve nem sorolhatók be kevés számú homogén betegút-csoportba. A betegek mindenkori állapotát, későbbi megbetegedéseit és ezek kezelését egyéni adottságaikon, érzékenységükön, korábbi kórtörténetükön kívül véletlen körülmények is jelentős mértékben alakítják. Az egyes betegút-elemekhez költségek társíthatók, amelyek nagysága és időbeni megoszlása szintén függ a beteg pillanatnyi egészségi állapotától, kórtörténetétől, egyéni adottságaitól és nem kis mértékben a beteg körülményeit, betegségét befolyásoló véletlen hatásoktól. A betegút egyik alkotó elemét, (pl. a terápiaváltást, hospitalizációt) kiragadva vizsgálható, hogy mennyi idő múlva következik be újabb változás, mekkora egészségügyi kiadás halmozódik fel addig, illetve a változás milyen típusú esemény, a beteg állapotának javulása-e vagy romlása, esetleg újabb terápiaváltás, újabb hospitalizáció vagy a betegút egyéb irányváltása lesz. A kiinduló állapottól egy újabb potenciális állapotig tartó idő betegről betegre változik, a betegre jellemző szisztematikus és ad hoc, véletlenszerű hatások függvényében, aminek valószínűség-eloszlása Coxféle arányos hazárd modell keretében becsülhető meg. A modell megadja a túlélési függvényt, ami kifejezi annak valószínűségét, hogy az állapotváltásig tartó idő egy előre adott időküszöböt meghalad. A váltást befolyásoló körülmények, jellemzők a modell magyarázó változói, kovariánsai, ezek a váltás relatív kockázatát számszerűsítik, és ezen keresztül módosítják a túlélési függvényt [1, 2, 3]. A lehetséges kiinduló- és célállapotok, események mindegyikéhez egy-egy Cox-modellt illesztve, az eredményekből megállapítható, hogy a kiinduló helyzetből adott idő elteltével adott irányú váltás mekkora valószínűséggel várható. A Cox-modellek túlélési függvényeiből és átmenet-valószínű- egészség-gAzdAságtAn MódsZERtAn ségeiből a betegutak számítógépes szimulálását lehetővé tevő szemi-Markov folyamat állítható össze [4]. A betegútszimuláció során az aktuális állapothoz az átmenet-valószínűség- és túlélési valószínűségek szerint kisorsolható a következő állapot típusa és a váltásig tartó idő. A betegúton eltelt időben előrehaladva a beteg ellátásához, kezeléséhez társuló költségek kumulálódnak. A betegút kiinduló és célállapotokra való szakaszolása után az egyes szakaszokon illesztett Cox-modellekhez kapcsolódóan a fellépő költségek megbecsülhetők. Elterjedt módszerek a költségek várható értékének, átlagának becslésére léteznek, ezek figyelembe veszik, hogy az adatok cenzoráltak, csak egy időablakon belül állnak rendelkezésre. A költségek valószínűség-eloszlását, szórását azonban általában nem modellezik [5, 6]. A kumulálódó költségek skálájára alkalmas módon illesztett Cox-modellekkel az átalányköltségek mellett a költségek szóródása, betegenkénti eloszlása is vizsgálható. A módszerek és eszközök szemléltetéséhez, a COPd betegek gyógyszeres terápiás kezeléseire és egyéb egészségügyi ellátási eseményeire vonatkozóan bemutatunk egy időés költségalakulást egyaránt leíró, számítógépes szimulációs betegút-modellt. A Cox-modelleket az Országos Egészségbiztosítási Pénztár adatbázisaiban fellelhető COPd-vel kezelt betegek adataira illesztettük retrospektív módon a 2008-2010-es évek közötti időszakban. A szimulált betegutak adataiból kiszámítottuk az események előtti és utáni időtartamok, költségek eloszlását, átlagát, szórását. A modell lehetőséget teremt a betegút állapotai közötti kapcsolatok feltérképezésére, támogatást nyújthat a betegellátási és finanszírozási szerkezet tervezett módosításaira vonatkozó hatástanulmányokhoz. AnyAg és Módszer A COPd-adatok forrása az OEP finanszírozási adatbázisa, s az adott diagnózis mellett 2008.01.01. – 2010.09.30. közötti időszakban követett, 626.667 páciens longitudinális adatainak leválogatását jelenti. A betegek által szedett gyógyszereket, gyógyszer-kombinációkat hét fő terápiás csoportba osztottuk (1. táblázat). A betegeket gyógyszerforgalmi adataik alapján időben dinamikusan terápiás csoportokhoz, lépcsőkhöz rendeltük, illetve ha a beteg által szedett gyógyszerek összetétele változott, akkor ennek megfelelően átsoroltuk. A betegutakat a terápiaváltásokon kívül a beteg állapotával összefüggő más események is tagolják, ezek felsorolása és értelmezése a 2. táblázatban található. 2. táblázat A betegutakat tagoló főbb ellátási események tapasztalatunk szerint a COPd-betegekre hét fő terápiás csoporttal és hat kiemelt eseménnyel megfelelően illeszkedő, a költségek alakulását is jól bemutató betegút-modell állítható össze. Az új gyógyszerek megjelenésére és a finanszírozási szabályok változására tekintettel a modelleket két egymást átfedő kétéves időablakra, 2008-2009-re és 2009-2010-re illesztettük egy 33.980 fős random betegmintára. Egy terápiaváltás vagy esemény bekövetkezése után az újabb eseményig, terápiaváltásig tartó idő szisztematikus és véletlen hatásoktól függ, ennek valószínűség-eloszlását Cox-féle arányos hazárd modellel becsültük meg a ténylegesen megfigyelt COPd betegút-szakaszok adataiból. Elképzelhető, hogy az új esemény vagy terápiaváltás a betegek egy részénél már nem megfigyelhető, mert az időablakon kívül következik be, azaz az adatok cenzoráltak. A Cox-féle túléléselemzés statisztikai módszerei képesek a cenzorálás okozta torzító hatást kiküszöbölni. A túléléselemzés eredménye az ún. túlélési függvény, ami adott t időponthoz megadja annak valószínűségét, hogy az eseményig, terápiaváltásig tartó idő legalább . Ebből kiszámítható az ún. hazárdfüggvény, ami adott t időpontban megadja az időegység alatt bekövetkező új események, váltások számának a kockázatnak kitett betegek számára vetített arányát, az ún. hazárdot (1. ábra). 1. ábra A túlélési és hazárdfüggvények összehasonlítása 1. táblázat Az elemzésben szereplő terápiás csoportok A túlélés elemzés során alkalmazott legegyszerűbb, ún. Kaplan-Meier módszer nem vesz figyelembe kovariánsokat, magyarázó változókat, egyetlen túlélési függvényt (és hazárdfüggvényt) illeszt a teljes populációra. A medián túlélési időnél a túlélési valószínűség½. A medián értéke és stan- IMe XI. évfolyAM 1. száM 2012. feBruár 35 egészség-gAzdAságtAn MódsZERtAn dard hibája kiszámítható a túlélési függvényre kapott becslésekből. Amennyiben a modellben kovariánsok is szerepelnek, akkor a Cox-féle arányos hazárd modellhez jutunk. Minden kovariánshoz tartozik relatív hazárd; ennyiszeresre változik a hazárd, ha a kovariáns értéke 1-gyel nő (pl. 1 évvel idősebb lesz a beteg). A kovariánsok lehetnek az eltelt időtől függetlenek (pl. a beteg neme), vagy függhetnek is az időtől (pl. a terápia típusa, dózisa, terápiaváltások száma adott időpontig). A cenzorált adatokból történő túlélés elemzési becslések torzítatlanságához szükséges, hogy a vizsgált eseményig, váltásig tartó idő ne függjön össze, ne korreláljon a cenzorálásig tartó idővel. néha ez a feltétel nem teljesül, pl. ha az esemény a kórházi felvétel, a cenzorálás pedig a beteg halálának időpontja. Ebben az esetben pozitív korreláció várható, ui. egészségesebb betegeknél rendszerint mind a hospitalizációig, mind pedig a halálozásig tartó idő megnő. Az általunk elemzett COPd betegek körében a mortalitás alacsony, ezért ez a torzító körülmény elhanyagolható. A vizsgált eseményig, terápiaváltásig a költségek minden beteg esetén növekednek, ezért az eltelt idő helyett a felhalmozódott költség függvényében is modellezhetjük az esemény bekövetkezésének valószínűségét. Így jutunk el a költségskálás túlélési függvény és a költségskálás hazárdfüggvény fogalmához. A túlélési függvény most adott költségküszöbhöz rendeli annak valószínűségét, hogy az eseményig kumulálódó összes költség legalább z. A hazárd adott költségszinthez megadja az egységnyi költségnövekedés során bekövetkező események számának a kockázatnak kitett betegek számához viszonyított arányát. A költségskálás túlélés elemzési modellekbe szintén illeszthetőek kovariánsok, többek között az eltelt idő is. A hozzájuk tartozó relatív költséghazárdok módosítják a költségskálás túlélési függvényt és a hazárdfüggvényt. Betegeket külön-külön követve valamely eseményig, a költség-idő összefüggés monoton növekvő, ugyanakkor betegről betegre más. A legegyszerűbb – ritkán előforduló, inkább elméleti – esetben, amikor minden betegre ugyanaz a költség-idő kapcsolat, ez a közös trend könnyen megkapható: számítsuk ki a költségskálás és az időskálás KaplanMeier túlélési függvényeket és minden időpontnak feleltessük meg azt a költséget, amihez ugyanaz a túlélési valószínűség tartozik. A betegeket pontosabban jellemző költség-idő trendeket kapunk, ha a költség- és időskálás Cox-modellekben a betegek költség- és időfüggő vagy állandó jellemzőit kovariánsként figyelembe vesszük. Ezek az ún. fix hatásos kovariánsos modellek azonban a betegek egy-egy rétegére továbbra is állandó trendeket feltételeznek. Egyedi, betegspecifikus összefüggéseket úgy nyerhetünk, ha betegspecifikus, random hatásos kovariánsokat, ún. frailty tagokat is beillesztünk a modellbe. Ezek egyéni érzékenységeket, exacerbációkat jelenítenek meg a költség-idő trendben. [3] A szimulált betegút-szakaszokhoz tartozó egyedi költségeket a most ismertetett idő- és költségskálás frailty modellek idő-költség összefüggéseiből kaphatjuk meg. Minden ál- 36 IMe XI. évfolyAM 1. száM 2012. feBruár lapotváltás után először kisorsoljuk a beteg következő eseményéig, váltásáig szükséges időt, majd a beteg egyéni érzékenységét. A következő eseményig kumulálódó költséget a betegspecifikus költségtrendből számítjuk ki, az eltelt idő behelyettesítésével. Eseményt, terápiaváltást követően adott időpontig halmozódó költségek elemzése nem csak közvetve, költségekkel társított számítógépes betegút-szimulációval lehetséges, hanem megfelelően beállított költségskálás Coxmodellel is. A költségek általában bizonyos események bekövetkezéséhez köthetők. Jellemzően hospitalizációhoz, műtéthez, gyógyszeres terápia megkezdéséhez kapcsolódnak gyógyítási költségek. Ezek az esemény, a váltás megtörténte után, bizonyos késéssel jelentkeznek. Annak érdekében, hogy a költségeknek az eltelt időtől való függését modellezhessük, bevezetünk egy mesterséges eseményt a vizsgált eseményt követő t-edik időpontban és költségskálás Cox-modellt illesztünk a mesterséges eseményig halmozódó költségre. A költségskálás túlélési függvény a t-edik időpontig kumulálódott költség valószínűségeloszlását adja meg. Ha a t időpontot véletlenszerűen választjuk – mondjuk egyenletes eloszlással egyéves időintervallumban –, akkor t értéke betegről betegre más lesz, t folytonos kovariánsként illeszthető be a költségskálás Cox-modellbe. Előfordulhat, hogy a kisorsolt t idő elteltével a beteget már nem figyeljük meg (terápiát váltott, meghalt, vagy már előbb cenzorálódott). Ebben az esetben is a t időponthoz tartozó mesterséges eseményt rendeljük a beteghez azzal, hogy ennek az eseménynek a szempontjából a beteg cenzorált. A t idő függvényében megjeleníthetjük a költségtúlélési függvényt, a kumulált költségek várható értékét, mediánját, alsó vagy felső kvartilisét (4. és 5. ábra). A modellben az eltelt időn kívül természetesen további kovariánsok is szerepelhetnek (pl. nem, életkor, korábbi és jelen anamnézis stb.). Ez a megközelítés egyszerűbben kivitelezhető, mint az idő- és költségskálás frailty modellek, azonban azáltal, hogy nem különíti el, nem modellezi szeparáltan a betegspecifikus, egyénenként különböző költség-idő összefüggéseket, az eljárás csak viszonylag alacsony cenzorálási arány mellett ad megfelelő eredményt. A betegenként eltérő sebességgel kumulálódó költségek ui. általában pozitív korrelációt mutatnak a cenzorálásig halmozódó költséggel [6]. A COPd költségadatokon ezért az események, terápiaváltások után csak aránylag rövid idő elteltéig felgyülemlő költségekre alkalmaztuk ezt a módszert. Retrospektív követéses adatok cenzorálási arányait egyébként könnyen csökkenthetjük, ha a feltétlenül szükséges követési időnél valamivel hosszabbat választunk. Összességében kijelenthető, hogy korlátainak tudatosítása mellett ezt az eljárást érdemes használni egyszerűsége és az eredmények szemléletessége miatt. A Cox-modellek illesztését, a betegutak szimulációját, a költségek statisztikai elemzését támogató számítógépes programokat az R 14.1-es statisztikai szoftvercsomaggal készítettük el. egészség-gAzdAságtAn MódsZERtAn eredMények Lévén ez egy metodikai közlemény, a COPd adatok elemzését, az eredményeket, a levont következtetéseket csak részlegesen, a módszerek illusztrálásához és megértéséhez szükséges mértékben mutatjuk be. Úgy gondoljuk, hogy a legfontosabb eredmény maga a betegutak és költségek modellezését lehetővé tevő matematikai statisztikai modell és az erre épülő számítógépes betegút-szimulációt végző szoftver. A Cox-modellek illesztését, a túlélési függvények és a magyarázó változókhoz társuló relatív hazárdok kiszámítását követően 10000 egymástól független betegutat generáltunk számítógépes szimulációval a 2009 és 2010-es évekre. A betegutak induló terápiáit, a betegek demográfiai és morbiditási jellemzőit az adott időszak kezdetének OEP COPd-vel kezelt betegek adataiból kinyerhető arányokra állítottuk be. Az eredményeket összevetettük a valódi COPdadatok negyedévre képzett átlagaival és gyakoriságaival. A szimuláció során a Cox-modellek paramétereit a pontosabb egyezés érdekében egy erre a célra kifejlesztett optimalizáló program segítségével tovább finomítottuk, amire azért van szükség, mert a Cox-eljárás nem teljesen a betegút-szimuláció szempontjából optimalizálja a becsléseket. A Coxmodellek paraméterei ugyanis csak a becslések standard hibáinak megfelelő mértékben pontosak. A szimuláció a betegutak generálása során a véletlen hibával terhelt becsléseket többször is felhasználja, a hibák ez által szisztematikussá válnak, hatásuk felerősödik, elvesztik véletlen természetüket, mintegy „beleégnek” a betegutakba. Ha az illeszkedés a gyakorlati követelményeknek még mindig nem felel meg, akkor szükség lehet a Cox-modellek finomítására, bővítésére újabb magyarázó változókkal, interakciós hatásokkal. Ilyen esetben részletesen össze kell hasonlítani a modell szerinti és a valódi betegáramlási, terápiaváltási gyakoriságokat és ezek időbeli változását. Az adjusztálást követően a finomított paraméterekkel újrafuttattuk a betegút-szimulációt és ellenőriztük a valódi COPd adatokkal való egyezés mértékét. A szimulált összes beteg számát arányosítással egyenlővé tettük a COPd adatbázis vizsgált időszakban előforduló betegeinek számával. Első lépésben a terápiák adott időszakbeli gyakoriságait és a terápiás napokat hasonlítottuk össze. Az adott negyedév kezdetekor a terápiás lépcsőn tartózkodó betegek arányai mellett vizsgálatra került az aktuális terápia hossza is (2. ábra). Az egyezés összességében kielégítő, nagyobb eltérés a terápiamentes esetek és a szteroid alapú terápiás lépcső (3) gyakoriságaiban mutatkozik. Az adott időszakra jellemző, megszakítás nélküli terápiahosszak valamivel gyengébb illeszkedést mutatnak, mint a betegszám-gyakoriságok. A kombinált szteroid terápiák (3) esetében az egy beteg esetében megfigyelhető, osztatlan terápiás időszakok hoszszabbnak, míg a második lépcső (2, 21, 22) tekintetében valamivel rövidebbnek bizonyultak a valóságban tapasztalhatónál. 2. ábra Betegek és terápiás napok terápiatípusok közötti megoszlása a megfigyelt negyedéves állapotokban Az ellenőrzés második lépésében a szimulált betegutakhoz tartozó költségek terápiák közötti megoszlását hasonlítottuk össze a valódi COPd költségekkel. A költségek pontossága függ a betegszám- és terápianap-megoszlások pontosságától is, ezért a költségskálás modellek elsődleges Cox-becsléseinek standard hibái általában nagyobbak, mint az időskálás modelleké. Itt is szükség van a becsült paraméterek szimuláció alapján történő finomítására. A finomítás után a betegutakhoz szimulált költségek terápiák közötti megoszlása jó egyezést ad a valós költségmegoszlással, a statisztikai feldolgozás alapját adó, közel 34 ezer fős OEP-minta kiadási szerkezetének vizsgálatakor (3. ábra). 3. ábra Éves terápiás költségek adott terápiás lépcső mellett kezelt betegeknél (finomítás után) A betegség lefolyásától függetlenül a betegek egy részének terápia tartása nem az elvártnak megfelelő. Erre mutat rá, hogy az adott időszak kezdetén terápiához nem rendelhető betegek jelentős költségeket halmoznak fel, ami nem magyarázható csak az újabb betegek belépésével. A betegek egyéb jellemzőitől, így nemétől és különösen életkorától is jelentős mértékben függenek a terápiás költségek, miután a költségek relatív kockázatai jelentős különbségeket mutatnak. A terápiára kerülés utáni költségek időtől való függésére kapott eredmények közül, a xantin származékok eseteiben megfigyelhető költségtúlélési függvényeit mutatja a 4. ábra. IMe XI. évfolyAM 1. száM 2012. feBruár 37 egészség-gAzdAságtAn MódsZERtAn 4. ábra Költség túlélési görbe alakulása az elhúzódó hatású xantin származékokra váltást követően 5. ábra A medián költség alakulása a terápiára váltásokat követően eltelt napok függvényében Látható, hogy 30 nappal a terápiára kerülés után a betegek 10%-ának költsége nagyobb 4 ezer Ft-nál, 60 napnál pedig a betegek 10%-ának 12 ezer Ft-nál nagyobb a költsége. Ezek az eredmények 20 évesnél fiatalabb férfiakra vonatkoznak, a nőkre és más korcsoportokra az eredmények a Cox-modell relatív hazárdjai szerint módosulnak. jellemzők melletti költségvarianciájának kimutatása hasznos segítséget nyújt a valóságot egyre akkurátusabban leképző költség-hatékonysági és finanszírozói modellek fejlesztéséhez. A terápiára kerülés utáni mediánköltségek az eltelt idő függvényében terápiánként különböző, jellegzetes mintázatot mutatnak, amit a 41-60 év közötti nőbetegek eseteiben figyelhetünk meg az 5. ábrán. A harmadik terápiás lépcső betegeire eső gyógyszerkiadások időben lineárisan nőnek, míg a második lépcsős terápiáké eleinte alacsonyak, majd kb. 40 nap elteltével ugrásszerű növekedésnek indulnak. A költségnövekedések időbeli alakulása terápiák között öszszehasonlítható, s részben visszaadja a kezelés és a betegegyüttműködés mintáit is. A betegség irreverzibilis lefolyása ellenére a beteg-compliance nem nevezhető optimálisnak a vizsgált betegkörben, illetve a harmadik terápiás lépcsőhöz rendelt betegek sokáig ingadoznak a tartós szteroid alapú terápiákra való állás előtt. A megfigyelés idejének elteltével azonban a betegek folyamatosan veszik fel kiadási tételeiket, ami a lineáris jellegű becslést adja. A második lépcsőre állított betegek új, magasabb költségű gyógyszerei a ráállítás pillanatától kezdve jelentős finanszírozási terhet jelentenek, amihez az első hónap végi terápiás kontrollt követő dózismódosítások, illetve a megjelenő rövid hatású szerek járulékos terheinek tendenciái társulnak változó mértékben. A terápiákat követő gyógyszerköltségekről a Cox-modellekből közvetlenül megállapítható, hogy az egymást követő terápiás lépcsők költségei hogyan viszonyulnak egymáshoz az időben. A kiegészítő és rohamterápiák jelentette kiadási többletek mellett a különbségek idővel el is mosódhatnak, illetve a LAMA (22) és különösen a LABA/ICs terápiák (3) költségeinek betegek közötti szóródása is igen magas. Érdemes vizsgálni, hogy mely betegcsoportokon hogyan alakulnak a gyógyszerköltségek adott időszaki mediánjai, amit az exacerbációk és szövődményekhez kapcsolódó kiadások is jelentősen befolyásolhatnak. Mindezek adott terápia melletti megjelenésének, elkerülésének, illetve adott beteg- Kialakítottunk egy terápiaváltásokat és gyógyítási eseményeket időben és költségskálán leíró betegút modellt. Ennek keretében vizsgálható, hogy a terápiaváltások és események bekövetkezésének valószínűségét, a váltásig tartó idő eloszlását, várható értékét, szórását a betegek jellemzői, korábbi eseményei, terápiái milyen mértékben módosítják. Az illeszkedés szimulációs kísérletekkel ellenőrizhető. Az illesztett, megfelelő módon kalibrált és ellenőrzött számítógépes modellel korlátlan számú betegút generálható, melyek segítségével képet kaphatunk a vizsgált terápiás időtartamok és költségek várható értékeiről, szórásáról, a kiugró értékek gyakoriságáról. Az eszközzel hatástanulmányok készíthetők új terápiák várható forgalmáról, költségeiről, a meglévő terápiás készítmények forgalmával való kölcsönhatásáról is. Jelenleg a költségekre vonatkozó modellek egy adott esemény után a rákövetkező eseményig adódó költségekre irányulnak. Ez a betegút-szimulációt jól támogatja, ugyanakkor nem rendeli a költségeket egyértelműen eseményekhez. Modellünk leírja, hogy az eseményekhez, terápiákhoz milyen kockázatok, időtartamok, költségek társulnak, az oksági kapcsolatokat azonban nem tárja fel, annak ellenére, hogy a társulás statisztikai szignifikanciája kiértékelhető. A jövőben a kauzalitás irányába is tovább kívánunk lépni. A diszkrét eseményvezérelt modellek fejlesztése jelentős nehézségekkel járhat, de szerepük felértékelődik, amikor a finanszírozási döntések során egyre nagyobb figyelem szentelődik a beteg és ellátási jellemzők kiadási szerkezetre tett hatásaira, illetve azok ez irányba történő lehetséges befolyásolására. Az idővel és a megelőző betegút eseményeivel korrigált pillanatnyi hazárdok eseményekre és az ezzel kapcsolatban felmerülő költségekre tett hatásainak vizsgálata nem csak a teljesítmény-alapú finanszírozási és 38 IMe XI. évfolyAM 1. száM 2012. feBruár MegBeszélés egészség-gAzdAságtAn MódsZERtAn betegség-menedzsment programok értékelésénél, de a kockázatok hatékony kezelését igénylő alap- és járulékos biztosítási csomagok kialakításánál is alapvető fontosságúvá válhat a közeljövőben. IROdALOMJEGYZÉK [1] Cox dR (1972): Regression models and life-tables (with discussion). J. Royal stat. soc. B, 34:187-220 [2] Kalbfleisch Jd, Prentice RL (2002): the statistical Analysis of Failure time data. 2nd Edition. Wiley. [3] therneau tM,Grambsch PM (2000): Modeling survival data. Extending the Cox Model. springer. [4] Ross sM (1995) stochastic Processes. Wiley. [5] Zhao H, Cheng Y, Bang H (2011): some insight on censored cost estimators. statistics in Medicine 30:23812388. [6] Lin dY, Feuer EJ, Etzioni R, Wax Y (1997): Estimating Medical Costs from Incomplete Follow-Up data. Biometrics 53:419-434. A sZERZŐK BEMUtAtÁsA lang zsolt jelenleg a Healthware tanácsadó Kft. biostatisztikusaként dolgozik. 1982-ben szerzett okleveles matematikus végzettséget az Eötvös Loránd tudományegyetemen. 1983-1989 között az MtA számítástechnikai és Automatizálási Kutató Intézet Geofizikai Osztályának, 1990-2008-ig az Abacom Kft., 1994-2008 között az OEP Közgazdasági és Biztosításpolitikai Főosztály statisztikai Osztályának munkatársa. 2001-2008 között a semmelweis Egyetemen és az Országos Környezetegészségügyi Intézet Környezetepidemiológiai Osztályán dolgozik. 20072008-ban a szent István Egyetem Állatorvos-tudományi Kar Biomatematikai és számítástechnikai tanszékének munkatársa. 1992-től a nemzetközi Klinikai Biostatisztikai társaság Magyar tagozatának tagja. rakonczai Pál alkalmazott matematikus, 2004-ben szerzett oklevelet az Eötvös Loránd tudományegyetem természettudományi Karán. 2004-től a Matematikai Intézet Valószinűségelmélet és statisztika tanszék doktorandusz hallgatója, 2007-től egyetemi tanársegéd. Főbb éredeklődési területei: többdimenziós összefüggési struktúrák, extrém érték elmélet, lineáris modellek. lazányi István zsolt Jelenleg a HealthWare tanácsadó Kft. fejlesztő informatikusa. 2003-ban végzett a Budapesti Műszaki és Gazdaságtudományi Egyetem műszaki informatika szakán. 2003-2006 között BMGE Irányítástechnika és Informatika tanszékének Phd képzésében vett részt. Phd fokozatát 2009-ben szerezte meg. dr. Bacskai Miklós a HealthWare tanácsadó Kft. ügyvezető igazgatója. diplomáit a debreceni Orvostudományi Egyetem Általános Orvostudományi Karán (1999), és a szegedi tudományegyetem Közgazdaságtudományi Kará- nak orvos-szakközgazdász szakán (2000) szerezte. 2000től az Országos Egészségbiztosítási Pénztár Gyógyszerügyi Főosztályán dolgozott, 2001-2002 között elemzési osztályvezető, 2002-2004 között gyógyszertámogatási osztályvezető munkakörben. szakmai gyakorlatait farmakoökonómia és információs rendszerszervezés területén szerezte. XIII. Outsourcing Konferencia Időpont: 2012. november 7. (szerda) Helyszín: Best Western Hotel Hungaria (1074 Budapest, Rákóczi út 90.) IMe XI. évfolyAM 1. száM 2012. feBruár 39