# code for example 01.03.01 [chapter 1, part 2, the first example]

mydata = read.table("data010301.txt")  
X = mydata$V1   # define X
Y = mydata$V2   # define Y

myreg = lm(Y~X) # do the linear regression

anova(myreg)


# Check F(1-alpha, 1, n-2 ) = t(1-alpha/2, n-2)

alpha <- c(0.01, 0.05, 0.1)
n <- 25

qf(1-alpha, 1, n-2)
qt(1-alpha/2, n-2)^2