https://github.com/cran/IQCC
Tip revision: ef4202fc61e884dce81758e0f31c6f1ae429e523 authored by Emanuel P. Barbosa on 10 March 2010, 00:00:00 UTC
version 0.6
version 0.6
Tip revision: ef4202f
T2.1.R
T2.1 <- function(estat, m, n)
{
t2 <- vector()
if(n == 1)
{
for (i in 1:m)
{
T2 <- (t(estat[[3]][i, ]) %*% solve(estat[[2]]) %*% (estat[[3]][i, ]))
t2 <- c(t2, T2)
}
}
if(n > 1)
{
for (i in 1:m)
{
T2 <- n * (t(estat[[3]][i, ] - estat[[1]]) %*% solve(estat[[2]]) %*% (estat[[3]][i, ] - estat[[1]]))
t2 <- c(t2, T2)
}
}
return(t2)
}