Nothing
##
matlogic <- function(x)
{
n=nrow(x)
m=ncol(x)
nblev <- apply(x,2,function(u){nlevels(as.factor(u))})
## Keep names....
rownames <- rownames(x)
colnames <- colnames(x)
i <- 0
colnamesnew <- c(apply(x,2,function(u){ i<<- i+1;paste(colnames[i],levels(as.factor(u)),sep=".")}),recursive=TRUE)
k <- sum(nblev)
res <- as.integer(matrix(0,ncol=k,nrow=n))
x <- c(x,recursive=TRUE)
result <- .C(C_matind,
as.integer(nblev),
as.integer(x),
res=res,
as.integer(n),
as.integer(m),
as.integer(k),
PACKAGE="amap")
result <- matrix(result$res,ncol=k)
rownames(result) <- rownames
colnames(result) <- colnamesnew
result
}
burt <- function(x)
{
ind <- matlogic(x)
t(ind) %*% ind
}
## x: table de burt, ou table
afc <- function (x)
{
f <- as.matrix(x/sum(x))
fi <- apply(f,1,sum)
fj <- apply(f,2,sum)
## Dr = diag(fi)
## Dc = diag(fj)
f <- (1/fi) * t(t(f)/fj)
acp(f,wI=fi,wV=fj,center=TRUE,reduce=FALSE)
}
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.