Hi guys,
I'm trying to conduct anova test to assess if the interaction terms are significant for each one of these dummy variables in terms of treatment (drug vs placebo) with Bonferroni correction for multiple comparisons.
However, I have WARNINGS: LSMEAN differences for ttr are not estimable.
What am I doing wrong? Any hints or suggestions appreciated.
Thanks in advance.
data have;
input id score ttr gender marriage drug stage age cvd geography race statin surgery;
cards;
1 4 1 0 0 0 1 0 0 0 1 1 1
2 5 1 0 0 1 1 0 0 1 1 1 0
3 8 1 1 0 0 1 1 1 1 1 0 0
4 9 1 1 0 1 1 1 1 1 1 0 0
5 14 1 1 1 0 1 1 1 1 1 0 0
6 2 1 0 1 1 0 0 0 0 0 0 0
7 15 1 1 0 0 0 1 0 0 0 1 1
8 6 1 0 0 1 0 0 1 0 0 0 1
9 9 1 1 1 1 0 1 0 0 0 0 1
10 10 1 1 1 1 1 1 1 0 1 0 0
11 0 2 0 1 0 0 1 0 1 0 0 0
12 1 2 0 0 0 1 1 1 1 1 1 0
13 4 2 0 0 0 0 0 0 0 0 1 0
14 -2 2 1 0 0 1 0 1 0 1 1 1
15 6 2 0 1 0 0 0 0 1 0 1 0
16 -4 2 1 1 1 0 0 0 0 1 0 1
17 3 2 1 1 1 0 0 0 1 0 0 0
18 4 2 1 0 1 0 0 0 0 1 1 1
19 7 2 0 0 0 1 0 0 1 0 1 0
20 9 2 0 0 0 1 0 0 0 0 0 1
;
proc glm data=have;
class ttr gender marriage drug stage age cvd geography race statin surgery;
model score = ttr ttr*gender
ttr*marriage
ttr*drug
ttr*stage
ttr*age
ttr*cvd
ttr*geography
ttr*race
ttr*statin
ttr*surgery;
lsmeans ttr ttr*gender
ttr*marriage
ttr*drug
ttr*stage
ttr*age
ttr*cvd
ttr*geography
ttr*race
ttr*statin
ttr*surgery/adjust=bon;
run;
Log:
WARNING: LSMEAN differences for ttr are not estimable; no table of differences will be produced.
WARNING: LSMEAN differences for ttr*stage are not estimable; no table of differences will be
produced.
WARNING: LSMEAN differences for ttr*age are not estimable; no table of differences will be
produced.
WARNING: LSMEAN differences for ttr*race are not estimable; no table of differences will be
produced.
