Preparar la tabla (data.frame) de datos

Y <- c(9,12,11,14,13,
       10,6,9,9,10,
       12,14,11,13,11,
       9,8,11,7,8)

prog <- c(rep(1,5),
          rep(2,5),
          rep(3,5),
          rep(4,5))

prog <- as.factor(prog)

datos <- data.frame(prog,Y)

datos
##    prog  Y
## 1     1  9
## 2     1 12
## 3     1 11
## 4     1 14
## 5     1 13
## 6     2 10
## 7     2  6
## 8     2  9
## 9     2  9
## 10    2 10
## 11    3 12
## 12    3 14
## 13    3 11
## 14    3 13
## 15    3 11
## 16    4  9
## 17    4  8
## 18    4 11
## 19    4  7
## 20    4  8