Netcompare Error after permutation about netcomi HOT 7 CLOSED

Hey Xiao,

Can't see what causes these errors at the moment and can't reproduce them with my data. Could you please send me a reproducible example with either your data (can be anonymized, or only a subset) or any public data via email? Then I will have a closer look into the issue.

Thanks a lot,
Stefanie

from netcomi.

@stefpeschel I just sent you an email with the data. Let me know if you need more files. Thank you!

from netcomi.

@stefpeschel, I have the same problem with my data.
I'm calculating my permutations for all samples (ITS data) - it is ok, but when I subset the samples by sample type (bulk soil, rhizosphere, shoots and roots), I have this error only for bulk soil and rhizosphere (but only when I set the permuation number for 1000, when I was only running sample code with 3 permutations, everyting was fine). Here is the code:
`ps.rhizo<-subset_samples(ps, Type == "rhizosphere")
permutations <- 1000

flist <- filterfun(kOverA(5, 2e-05))
ent.logi <- filter_taxa(ps.rhizo, flist)
ent.trim <- filter_taxa(ps.rhizo, flist, TRUE)
identical(ent.trim, prune_taxa(ent.logi, ps.rhizo))
identical(sum(ent.logi), ntaxa(ent.trim))
ps.rhizo <- filter_taxa(ps.rhizo, flist, TRUE)

agglomerated <- tax_glom(ps.rhizo,"Species")
agglomerated_healthy = subset_samples(agglomerated, health == "healthy")
agglomerated_unhealthy = subset_samples(agglomerated, health == "unhealthy")

net <- netConstruct(data = agglomerated_healthy,
data2 = agglomerated_unhealthy,
filtTax = "highestVar",
filtTaxPar = list(highestVar = 50),
measure = "spring",
measurePar = list(nlambda=100,rep.num=10),
normMethod = "none",
zeroMethod = "none",
sparsMethod = "none",
dissFunc = "signed",
verbose = 3,
seed = 420)

analyzed <- netAnalyze(net,
centrLCC = TRUE,
avDissIgnoreInf = TRUE,
sPathNorm = FALSE,
clustMethod = "cluster_fast_greedy",
hubPar = c("degree", "between", "closeness"),
hubQuant = 0.9,
lnormFit = TRUE,
normDeg = FALSE,
normBetw = FALSE,
normClose = FALSE,
normEigen = FALSE,
weightDeg = FALSE)

compared <- netCompare(analyzed, permTest = TRUE, nPerm = permutations,
storeAssoPerm = TRUE,
fileStoreAssoPerm ="Asso_perm_comp",
storeCountsPerm = FALSE,
seed = 420)
`

All the best,
Dominika

from netcomi.

Thanks @lixiaopi1985 for sending me your data. I just uploaded a new version (only at 'develop' branch for now) that should fix this error. Could you please install the latest version from develop branch (version 1.0.2.9000) and let me know whether still any errors occur?

The 'develop' branch version can be installed via:

devtools::install_github("stefpeschel/NetCoMi", ref = "develop")

from netcomi.

@stefpeschel Thank you so much for fixing this error! I did not have any problems this time on the 16S data after updating the package. I am running it on the ITS data as it showed different error messages last time, I will share the result soon.

from netcomi.

Great! Please let me know if you're still running into errors for the ITS data.

Last time, the ITS error occured when computing the group differences for the centralities. I also adapted this approach, so I hope it runs without any problems now.

from netcomi.

@stefpeschel No issue running the ITS data this time, hooray!!!

Thank you so much!

from netcomi.