T2.2.R
T2.2 <- function(datum2, estat, n)
{
if(n == 1)
{
T2 <- (t(datum2 - estat[[1]]) %*% solve(estat[[2]]) %*% (datum2 - estat[[1]]))
}
if(n > 1)
{
media <- colMeans(datum2)
T2 <- n * (t(media - estat[[1]]) %*% solve(estat[[2]]) %*% (media - estat[[1]]))
}
return(T2)
}