library(codeutils)
library(hplot)
library(igraph)
rundir <- pathtopath(getDBdir(),"/A_code/sizemod/R/")
dir(rundir)
checkvar <- c("readLBMdata","initiateModel",
"preliminaryfit","getproduction","readLBMctrl")
checkvar <- c("makeoutputsm")
checkvar <- c("dynamics","negLLP","do_modfit","getvcov")
checkvar <- c("do_ESS","printESS")
outmap <- getfunmap(indir=rundir,checkfun=checkvar,sortby="functions",
files="",outfile="")
# arrows imply a function is called by whichever is pointing to it
plotprep(width=11, height=10)
parset(cex=1.25)
plot(graph_from_edgelist(as.matrix(outmap$alledge)))
str1(outmap)
outmap$isolates[,"function"]
outmap$edg
str1(outmap)
funs <- outmap$functions
nC <- nrow(funs)
counts <- numeric(nC)
for (i in 1:nC) {# i = 5
counts[i] <- length(removeEmpty(unlist(strsplit(funs[i,"calls"],", "))))
}
pickMany <- which(counts > 4)
funs[pickMany,"function"]
pickESS <- which(funs[,"function"] == "do_ESS")
funs[pickESS,3:4]
outmap$isolates[,"function"]
g1 <- graph_from_literal(writectrlfile -+ writedatafile -+ initiatemodel -+
equilibrium)
library(hplot)
plotprep(width=6, height=5)
parset(cex=1.25)
plot(g1)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.