nGenes=ncol(datExpr) nSample=nrow(datExpr) #选基因 datME=moduleEigengenes(datExpr,moduleColors,trapErrors=FALSE)$eigengenes HG=networkScreeningGS(datExpr,datME,t(geneTraitSignificance),oddPower = 3,blockSize = 1000,minimumSampleSize = 4, addGS = TRUE) GeneResultsNetworkScreening=data.frame(Gene...
回答于 2018-04-03 19:28