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(GeneName=row.names(HG),HG)
write.table(GeneResultsNetworkScreening,file="HubGene.txt",row.names=F,sep="\t")