looping a looped function over factors

Question

Let\'s say I have a 3 variables that were collected at 3 sites

Var1 <- c(1:20)
Var2 <- c(41:60)
Var3<-c(81:100)
Site<-rep(letters[1:3],20)