xy = read.table("data020603.txt")

y = xy[,1];
X = xy[,2];
D = xy[,3];

DX = D*X;

r = lm(y~X+D+DX)
anova(r)