HELP with matched case-control study

Posts: 0

HELP with matched case-control study

Hi guys... so glad I found this board.

I need help with a matched case-control study. I had an outcome of interest 'TMC' and matched 1 case of 'TMC' to 2 controls by age. Now, I want to analyse with SAS to see whether variables of interest are different (2-sided) in cases vs. controls. The data is non-parametric, non-symmetrical for the most part.

'TMC' is the condition (case=1, control=0), 'daysadmit' is the continuous variable

For the continuous variables, I used the following code:

PROC GLM DATA = descriptive ;
CLASS match tmc ;
MODEL daysadmit = tmc match(tmc);
MEANS tmc;

Is this correct? Do I need to look at the p-value for match(tmc), or do I need that nested variable in there at all? Or, is sufficient to have 'match' in the class statement and only have daysadmit=tmc as the model?
Also, is there any value to putting in this statment:
TEST H = tmc E = match(tmc);

For the categorical variables I used:

data binomial;
set descriptive;
if TMC=0 then status=2;
if TMC=1 then status=1;
if TMC=1 then censor=1;
if TMC=0 then censor=0;

proc phreg data=binomial;
model status*censor(0)= fever /ties=discrete risklimits;
strata match;

Where I made status=1 controls, status=2 cases (of toxic megacolon), censor=0 controls, censor=1 cases.

The problem here is that the numbers are strange. I want to see whether the frequency of a bunch of variables is different in the case group from the control group. The the predictor variable listed above ('fever'), 90% of the case group ('TMC') had fever while 25% of the control group had fever. Using the above program, I get a p of 0.026. With leukocytosis 'leuko' as the predictor, 80% of the case group had leukocytosis while 5 % of the control group had leukocytosis. It looks like it should be significant, but I get a p value of 0.996.

HELP! Can someone help me with what to do?
Post a Question
Discussion Stats
  • 0 replies
  • 1 in conversation