## Function to calculate angles

angle <- function(x,y) { 
  calculo <- sum(x*y)/(sqrt(sum(x^2))*sqrt(sum(y^2)))
  angleprevio<-acos(calculo)
  angle<-angleprevio*180/pi
  return(angle)  
}