doubleLoop

MyData <- data.frame(Year=seq(2000,2009))

# year variables
dim <- length(MyData$Year)
yMat <- matrix(NA, nrow=dim, ncol=dim)
for(i in 1:dim) {
      yMat[i,] <- rbind(t(MyData)) # transpose
}
yMat

# dummy variables
dMat <- matrix(0, nrow=dim, ncol=dim)
for(j in 1:dim) {
      for(i in 1:dim) {
            if(i > j)
                  dMat[i,j]=1
      }
}
dMat

# linear variables
lMat <- matrix(0, nrow=dim, ncol=dim)
for(j in 1:dim) {
      for(i in 1:dim) {
            if(i > j)
                  lMat[i,j]=lMat[i-1,j]+1
      }
}
lMat

# quadratic/cube variables
qMat <- lMat^2
qMat
cMat <- lMat^3
cMat