Since roughly last November, I have had IBS-like symptoms but never met Rome 3 criteria until January in terms of percentage of abnormal stools (just hit 25%). I suffer mostly from stomach pain which I attribute to my colon. I am under the supervision of a GI doc and PCP and have had several tests done – everything normal. Right now, we are in the stage of trying different medications and seeing what works.

I’ve been keeping a spreadsheet of each med I take, how much cardio I do, and an AM/PM health quality index (HQI) which i define as follows…

HQI scores from 1 to 4 based on levels of my functionality vs symptoms;

Level 4, no symptoms present and I do everything I wanted to do during the time period (AM or PM)

Level 3, symptoms but able to cope in such a way that I do everything I wanted to do during the time period

Level 2, symptoms but unable to cope in such a way that I don’t get to do everything I set out to do

Level 1, require medical attention due to symptom severity so that I can do anything. (i never hit this level)

I group every BM by bristol stool scale using an android app. For each day, I create another data product from this called adverse event (which is really just a stool outside of the BSS 3-5 range).

With the metrics described above, I try to achieve the highest HQI possible while minimizing adverse events. I average over 3 and 7 day means because the medicines take time to reach therapeutic levels due to their half life.

I am currently taking a few things:

- 40mg daily of nexium in 2, 20mg doses
- 2-3 librax/day
- 1 clariten-d prn before bed for seasonal allergies
- 5k IU vit d when i can remember to do so
- 1 glass Metamucil mostly daily to help normalize stools since the meds arent completely working (and fiber has been shown to help the biome greatly)

Here is the output of the regression I ran from yesterday. Can you see the culprit giving me issues?

```
===================================================================================================================
==================== bss_adverse_events_3day ~ cardio + nexium + librax + clrtn + vitd + mtmcl ====================
===================================================================================================================
OLS Regression Results
===================================================================================
Dep. Variable: bss_adverse_events_3day R-squared: 0.306
Model: OLS Adj. R-squared: 0.236
Method: Least Squares F-statistic: 4.344
Date: Sun, 12 Mar 2017 Prob (F-statistic): 0.00107
Time: 23:31:35 Log-Likelihood: -109.56
No. Observations: 66 AIC: 233.1
Df Residuals: 59 BIC: 248.5
Df Model: 6
Covariance Type: nonrobust
==============================================================================
coef std err t P>|t| [0.025 0.975]
------------------------------------------------------------------------------
Intercept 2.8728 0.393 7.308 0.000 2.086 3.659
cardio -0.5571 0.288 -1.935 0.058 -1.133 0.019
nexium -0.8889 0.258 -3.445 0.001 -1.405 -0.373
librax 0.0720 0.174 0.414 0.681 -0.276 0.420
clrtn 1.1026 0.393 2.806 0.007 0.316 1.889
vitd -0.0406 0.072 -0.563 0.575 -0.185 0.104
mtmcl -0.5411 0.444 -1.218 0.228 -1.430 0.348
==============================================================================
Omnibus: 5.193 Durbin-Watson: 0.955
Prob(Omnibus): 0.075 Jarque-Bera (JB): 4.918
Skew: 0.668 Prob(JB): 0.0855
Kurtosis: 2.969 Cond. No. 11.7
==============================================================================
```

regression for hqi 3 day mean in future

```
========================================================================================================
==================== hqi_3daymean ~ cardio + nexium + librax + clrtn + vitd + mtmcl ====================
========================================================================================================
OLS Regression Results
==============================================================================
Dep. Variable: hqi_3daymean R-squared: 0.300
Model: OLS Adj. R-squared: 0.229
Method: Least Squares F-statistic: 4.212
Date: Sun, 12 Mar 2017 Prob (F-statistic): 0.00136
Time: 23:31:35 Log-Likelihood: -19.365
No. Observations: 66 AIC: 52.73
Df Residuals: 59 BIC: 68.06
Df Model: 6
Covariance Type: nonrobust
==============================================================================
coef std err t P>|t| [0.025 0.975]
------------------------------------------------------------------------------
Intercept 3.1519 0.100 31.446 0.000 2.951 3.352
cardio -0.0630 0.073 -0.858 0.394 -0.210 0.084
nexium 0.1930 0.066 2.933 0.005 0.061 0.325
librax 0.0441 0.044 0.993 0.325 -0.045 0.133
clrtn -0.0951 0.100 -0.949 0.347 -0.296 0.105
vitd -0.0146 0.018 -0.795 0.430 -0.051 0.022
mtmcl 0.0346 0.113 0.306 0.761 -0.192 0.261
==============================================================================
Omnibus: 1.538 Durbin-Watson: 0.431
Prob(Omnibus): 0.463 Jarque-Bera (JB): 0.875
Skew: 0.002 Prob(JB): 0.646
Kurtosis: 3.564 Cond. No. 11.7
==============================================================================
```

I have normalized the regression by dosage so you can interpret the coef as the contribution of a dose. E.g. in the first regression table a coef of 1 means that a dose of the item contributes to 1 adverse event (linearly).

So the problem is the clariten-D I believe. It creates adverse events and does nothing for my HQI. So this week I am quitting it and will evaluate for the next 2 weeks to see if helps improve my symptoms.