Displaying 1 result from an estimated 1 matches for "mard1".
Did you mean:
mardi
2007 Jun 16
1
mardia's test
...of Mardia's test, what does the line " x1 = x[x[, p] == i, -p]" mean? Thanks a lot!
function (x)
{
p = dim(x)[2]
f = p - 1
clases = length(table(x[, p]))
for (i in 1:clases) {
x1 = x[x[, p] == i, -p]
ndat = dim(x1)[1]
mo3 = mo3(x1)
mard1 = ndat * mo3/6
cat("Mardia's test for class", i, "\n")
cat("mard1=", mard1, "\n")
p1 = 1 - pchisq(mard1, df = f * (f + 1) * (f + 2)/6)
cat("pvalue for m3=", p1, "\n")
mo4 = mo4(x1)
mard2...