I'm investigating the association between a (life time mean) vitamin level and the occurrence of a disease (yes/no).

I am using a 4600 patients data set (after cleaned) among 160 developed the disease and the rest haven't.

I have done a logistic regression where my dependent variable was the disease status (yes/no) and the independent variables were those I have adjusted to using stepwise approach and clinical literature. My variable of interest I have coded to dichotomic variable for each patient 1 insufficient vs 0 sufficient using accepted level. Then I have found the same for deficient vs sufficient and then for severe Def vs suf. Every time using the same model and always vs sufficient levels.

The odds ratios were always significant and the as the deficiency was stronger the odds ratio were getting stronger (the expB was bigger each time).

=====

Since it's quite convincing there's some sort of trend I wanted to show a dose response. How would you do so as for analysis and as for graphically?

=====

The only idea I came up with was this: having 75 and above are normal making new dichotomic variables instead for three levels (deficient vs sufficient, insufficient vs sufficient and svr deficient vs insufficient) to do so for a lot of levels saying every 5 nmol interval, for each one writing the odds ratio as appears in the same adjusted log reg and then to show a scatter plot y for odds ratio x for vit level and since the odds ratio are increasing gradually it would like nicely. I can also calculate regression line r square and so on. How does it sound? Any other ideas?

====

Thank you

David