r test

rt <-read.table("C:/Users/Ju/Desktop/anova/anova.txt",head=TRUE)
lm.sol <-lm(Y~X1+X2+X3,data=rt)
summary(lm.sol)

"C:/Users/Ju/Desktop/anova/anova.txt" contents:

X1 X2 X3 Y
1 12 123 531
2 13 124 539
3 14 125 547
4 15 126 555
5 16 127 563
6 17 128 571
7 18 129 579
8 19 130 587
9 20 131 595
10 21 132 603
11 22 133 611
12 23 123 222