Nothing
standtab <-
function(x, std = T, ctr = T){
ntot <- sum(x)
if(length(dim(x)) != 3){
stop("standtab: this is not a 3-way array!\n")
}
if(ntot == 0){
stop("This array has for sum zero!")
}
x <- x/ntot
pii <- margI(x)
pj <- margJ(x)
pk <- margK(x)
ui <- rep(1, length(pii))
uj <- rep(1, length(pj))
y <- rep(1, length(x))
dim(y) <- dim(x)
y <- if(ctr) prod3(y, pii, pj, pk) else 0
x <- x - y
if(std){
prod3(x, 1/sqrt(pii), 1/sqrt(pj), 1/sqrt(pk))
} else {
prod3(x, 1/pii, 1/pj, 1/pk)
}
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.