Dear Stata Forum Users, do really need your help. I am reporting results of my randomized study on medication adjustment and face a challenge to adjust for the baseline values.
1) dataset:
I have 400 patients: 200 control/200 intervention, variables recorded: id (patient study number), r_treat (intervention or control group); nmedicine1 and nmedicine3 (number of medications received at day 3 and day 90), polypharmacy1 and polypharmacy3 - at day 3 and 90 (generated dichotome variables: if 10 and more drugs - polypharmacy==1, if 0-9 drugs polypharmacy==0)
2) Question : I compared number of persons with polypharmacy1 in the two groups and polypharmacy3 in with command:
*calculation risk ratios -table 90 days
csi 38 60 132 106, or
BUT i got a note that I should adjust for the baseline hyperpolypharmacy1 when I compare hyperpolypharmacy3.
-got an advice to calculate 2x2 table and then compare with z-test either by hand og using mcci to calculate the difference from baseline to 90 days and then again do z-test
I am really lost how to do it practically - I am sure there is a smart way to do it in STATA, Could you please help me with the commands?
3) example of my dataset:
[CODE]
* Example generated by -dataex-. To install: ssc install dataex
clear
input int id float(nmedicine1 nmedicine3) byte r_treat float(hyperpolyfarmacy1 hyperpolyfarmacy3)
9 11 12 1 1 1
11 9 8 2 0 0
12 4 3 1 0 0
17 10 7 1 1 0
20 10 8 2 1 0
39 7 6 2 0 0
53 10 8 2 1 0
54 3 3 2 0 0
55 9 7 1 0 0
56 9 6 2 0 0
57 6 6 1 0 0
59 6 6 1 0 0
65 7 9 2 0 0
68 11 14 2 1 1
69 6 6 2 0 0
74 6 9 1 0 0
79 8 7 1 0 0
82 5 5 1 0 0
93 7 7 2 0 0
97 10 10 1 1 1
101 11 10 1 1 1
104 6 5 2 0 0
105 4 7 2 0 0
Thanks alot on advance for your help! It is really appreciated! Sincerely, Natallia
1) dataset:
I have 400 patients: 200 control/200 intervention, variables recorded: id (patient study number), r_treat (intervention or control group); nmedicine1 and nmedicine3 (number of medications received at day 3 and day 90), polypharmacy1 and polypharmacy3 - at day 3 and 90 (generated dichotome variables: if 10 and more drugs - polypharmacy==1, if 0-9 drugs polypharmacy==0)
2) Question : I compared number of persons with polypharmacy1 in the two groups and polypharmacy3 in with command:
*calculation risk ratios -table 90 days
csi 38 60 132 106, or
BUT i got a note that I should adjust for the baseline hyperpolypharmacy1 when I compare hyperpolypharmacy3.
-got an advice to calculate 2x2 table and then compare with z-test either by hand og using mcci to calculate the difference from baseline to 90 days and then again do z-test
I am really lost how to do it practically - I am sure there is a smart way to do it in STATA, Could you please help me with the commands?
3) example of my dataset:
[CODE]
* Example generated by -dataex-. To install: ssc install dataex
clear
input int id float(nmedicine1 nmedicine3) byte r_treat float(hyperpolyfarmacy1 hyperpolyfarmacy3)
9 11 12 1 1 1
11 9 8 2 0 0
12 4 3 1 0 0
17 10 7 1 1 0
20 10 8 2 1 0
39 7 6 2 0 0
53 10 8 2 1 0
54 3 3 2 0 0
55 9 7 1 0 0
56 9 6 2 0 0
57 6 6 1 0 0
59 6 6 1 0 0
65 7 9 2 0 0
68 11 14 2 1 1
69 6 6 2 0 0
74 6 9 1 0 0
79 8 7 1 0 0
82 5 5 1 0 0
93 7 7 2 0 0
97 10 10 1 1 1
101 11 10 1 1 1
104 6 5 2 0 0
105 4 7 2 0 0
Thanks alot on advance for your help! It is really appreciated! Sincerely, Natallia
Comment