https://github.com/cran/IQCC
Raw File
Tip revision: ef4202fc61e884dce81758e0f31c6f1ae429e523 authored by Emanuel P. Barbosa on 10 March 2010, 00:00:00 UTC
version 0.6
Tip revision: ef4202f
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)
}
back to top