Nothing
computeCausOrder <- function(G)
# Copyright (c) 2013 Jonas Peters [peters@stat.math.ethz.ch]
# All rights reserved. See the file COPYING for license terms.
{
p <- dim(G)[2]
remaining <- 1:p
causOrder <- rep(NA,p)
for(i in 1:(p-1))
{
root <- min(which(colSums(G) == 0))
causOrder[i] <- remaining[root]
remaining <- remaining[-root]
G <- G[-root,-root]
}
causOrder[p] <- remaining[1]
return(causOrder)
}
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.