SpiecEasi
SpiecEasi copied to clipboard
How to filter output the results
Hi, I followed the sample tutorial to analyze the intestinal data and got too many edges.What methods are available to filter the results? Thanks! This is the code I am running:
#glasso
se.gl.amgut <- spiec.easi(amgut1.filt, method='glasso', lambda.min.ratio=1e-2,
nlambda=20, pulsar.params=list(rep.num=50))
#mb
se.mb.amgut <- spiec.easi(amgut1.filt, method='mb', lambda.min.ratio=1e-2,
nlambda=20, pulsar.params=list(rep.num=50))
secor <- as.matrix(cov2cor(getOptCov(se.gl.amgut)))
sebeta <- as.matrix(symBeta(getOptBeta(se.mb.amgut), mode='maxabs'))