### R代码来自Vignette源'Yeastexample.rnw'############################################################# ###代码块数字1：设置##############################################################＃##代码块编号2：makeclustergraph ######################################################.character））ccclust = conncomp（cg1）########################################] ccl2 = cclens [cclens> 1] ####################################################### ###代码块编号4：CreateSubg ###########################################################ccl2 [[5]，litg）sg2 =子图（ccl2 [[1]]，litg）############################################################### ###代码块数字5：子图###################################################################################### ### code chunk number 6: sG2 ################################################### if (require("Rgraphviz")) plot(sG2, "neato") ################################################### ### code chunk number 7: intersect ################################################### commonG = intersection(litG, cg1) commonG ################################################### ### code chunk number 8: edgePerm ################################################### ePerm = function (g1, g2, B=500) { ans = rep(NA, B) n1 = nodes(g1) for(i in 1:B) { nodes(g1) = sample(n1) ans[i] = numEdges(intersection(g1, g2)) } return(ans) } set.seed(123) ##takes a long time #nPdist = ePerm(litG, cg1) data(nPdist) max(nPdist)