X=read.table(commandArgs(3))
m=lm(V4~.^2,data=X)
summary(m)