Результаты проведенных клинико-биохимических исследований явились основанием для разработки способа прогнозирования течения диффузного остеопороза у больных множественной миеломой. Нами предложен способ прогнозирования течения данного заболевания, основанный на учете биохимических показателей костного ремоделирования и минеральной плотности костной ткани. Прогнозирование осуществлялось с помощью построения аналитической модели при помощи программы Mathcad, версия 11.0b.
Развитие и характер течения диффузного остеопороза у больных множественной миеломой, степень снижения минеральной плотности костной ткани обусловливают биохимические показатели костного ремоделирования. Каждый из показателей, определяющий течение диффузного остеопороза у больных множественной миеломой, имеет определенное значение для степени снижения минеральной плотности костной ткани, то есть определенную информативность.
Основываясь на физиологической сущности рассматриваемого явления (минеральная плотность костной ткани пациентов), предложена переменная Y(в нашем случае ультразвуковая денситометрия (УЗД), являющаяся непрерывной функцией ряда независимо действующих факторов (показателей, обусловливающих активность процесса костного ремоделирования), вследствие чего использовалась полиномиальная модель.
Аппроксимация лабораторных данных проводилась при помощи полиномов 2-й степени методом наименьших квадратов. Степень полинома лимитировалась числом пациентов, включенных в исследование, и числом включенных факторов. Общий вид модели:
где Yu = прогнозируемый параметр;
β – коэффициент;
m – степень полинома;
Х – независимо влияющий фактор;
εu – погрешность модели.
Выбор компонентов модели осуществлялся по степени тесноты связи вариационных рядов УЗД после 12 месяцев лечения и всех исходных параметров. Включенные факторы были центрированы.
Для построения аналитической модели использовались значения из вариационных рядов общей лактатдегидрогеназы, костной щелочной фосфатазы, тартратрезистентной кислой фосфатазы и уровня кальция сыворотки крови до начала терапии.
Зависимость УЗД после 12 месяцев терапии от выбранных факторов в общем виде представлена в виде следующего уравнения:
где Y – УЗД;
A – уровень общей лактатдегидрогеназы;
B – уровень костной щелочной фосфатазы;
C – уровень тартрат–резистентной кислой фосфатазы;
D – уровень кальция сыворотки крови;
k – коэффициенты;
ε – погрешность модели
Точность аппроксимации данных регрессионной моделью оценивалась с помощью метода исследования остатков. Коэффициенты аналитической модели представлены в таблице 1.
Все больные множественной миеломой с диффузным остеопорозом были разделены нами на 4 группы (А1, В1,B2,C).
Группа А - 20 человек, получали на фоне стандартной схемы лечения (ПХТ по программе М2) антирезорбтивный препарат миакальцик (назальный спрей) в дозе 200 МЕ в сутки в течение двух месяцев с двухмесячным перерывом в интермиттирующем режиме в течение года.
Группа В1 - 18 человек, получали на фоне основной схемы лечения антирезорбтивный препарат - алендронат (фосамакс) в дозе 10 мг в сутки ежедневно утром натощак в течение года.
Группа B2 - 18 человек, получали на фоне стандартной схемы лечения - алендронат (фосамакс) в дозе 70 мг в неделю однократно утром натощак в течение года.
Группа C- контрольная группа, 20 человек, получали стандартную схему лечения множественной миеломы.
Выбранная нами модель была проверена на больных множественной миеломой с диффузным остеопорозом, вошедших в наше исследование. Для всех больных с известными исходами ультразвуковой денситометрии (УЗД) после 12 месяцев лечения был осуществлен прогноз, т.е. рассчитаны ожидаемый уровень УЗД после 12 месяцев лечения и отклонение (остаток) рассчитанного. Остаток был равен разнице ожидаемого и реального уровней УЗД.
Результаты компьютерной обработки позволяют сделать заключение по каждой группе больных множественной миеломой с диффузным остеопорозом.
В группе А ожидаемый уровень УЗД после 12 месяцев лечения мало отличался от реального значения, в результате чего остаток был слишком мал (остаток Rmax=0,25; остаток Rmin=0,066; SD|R| = 0,072).
В группе В1 получены аналогичные результаты. Rmaxсоставил 0,167, а Rminсоставил 0,045; SD|R| = 0,051.
В группе B2ожидаемый уровень УЗД после 12 месяцев лечения практически совпал с результатами в группе В1 и также мало отличался от реального значения (Rmax=0,16; Rmin=0,032; SD|R| = 0,039).
В группе С получены следующие значения: ожидаемый и реальный уровни УЗД после 12 месяцев лечения практически совпадали полностью, а остаток был слишком мал (Rmax=1,8*10-11; Rmin=5,3*10-12; SD|R| = 5,2*10-12).
Таким образом, на основании полученных данных мы пришли к выводу, что построенная аналитическая модель обладает достаточной точностью и позволяет прогнозировать изменение минеральной плотности костной ткани при динамическом наблюдении больных множественной миеломой с диффузным остеопорозом в зависимости от исходных показателей метаболизма костной ткани.