Я установил конкретную формулу c в MICE.
Но результат иногда не принимается.
DBP должен быть выше, чем SBP, ---> случай 1
и MBP не является "DBP + (SBP - DBP) / 3" ---> случай 2
Как я могу сделать этот случай?
- использовать вмененные данные как таковые.
- адаптировать данные к формуле.
- et.al.
imp = mice.MICEData(mdf[['DBP', 'MBP', 'SBP']])
fml = 'MBP ~ DBP + np.divide((DBP - SBP), 3)'
mMice = mice.MICE(fml, sm.OLS, imp)
static_results = mMice.fit()
|MEASUREMENT_SOURCE_VALUE|DBP |MBP |SBP |
|1 |41.0|40.0|38.0| case 1
|7 |41.0|40.0|38.0| case 1
|13 |43.0|42.0|38.0| case 1
|24 |41.0|50.0|70.0| case 2
|53 |22.0|29.0|38.0| case 2