paternogbc / sensiPhy

R package to perform sensitivity analysis for comparative methods
http://onlinelibrary.wiley.com/doi/10.1111/2041-210X.12990/full
GNU General Public License v2.0
12 stars 4 forks source link

BUG: clade_phylm #192

Open MatiasIMunoz opened 1 year ago

MatiasIMunoz commented 1 year ago

Dear developers,

function

clade_phylm()

Problem

I used the clade_phylm() before, but now when trying it I come across an error indicating different number of rows in the clade to be removed. I come across the same error with my own data, and with the primates example from the vignette. The samp_phylm() function works fine with my data, but not clade_phylm().

I couldn't manage to use the reprex library, but below you can see the output from running the primate example on my laptop.

sessionInfo() R version 4.3.1 (2023-06-16) Platform: x86_64-apple-darwin20 (64-bit) Running under: macOS Big Sur 11.3.1

sensiPhy version 0.8.5

Reproducible example

data(primates)

run analysis:

clade <- clade_phylm(log(sexMaturity) ~ log(adultMass), phy = primates$phy[[1]],

  • data = primates$data, nsim = 10,clade.col = "family", track=TRUE) Used dataset has 95 species that match data and phylogeny |===================================================================| 100%

    To check summary results and most influential clades:

    summary(clade) Error in data.frame(clade removed = c, N.species = ce$N.species, estimate = numeric(length(c)), : arguments imply differing number of rows: 0, 4