landerlabcode / BigSurR

MIT License
3 stars 0 forks source link

Missing Rmpfr? #2

Open DanielWendisch opened 3 weeks ago

DanielWendisch commented 3 weeks ago

Hi and thank you for this project! Just a note:

I got..

Error in mpfr(-(x[3]^2)/2, precBits = 128) : could not find function "mpfr"

when running..

output <- BigSur(example.seurat.subset, variable.features = F, correlations=T, depthlist = depths)

It disappeared after installing Rmpfr and executing library(Rmpfr).

I had installed all dependencies but still got the error on my windows machine.

Hope this is helpful. Best wishes, Daniel

[1] "Pipeline started execution." [1] "Modified corrected Pearson residuals calculated." [1] "Modified corrected Fano factors calculated." [1] "Beginning correlation calculation." [1] "Modified-corrected Pearson Correlation Coefficients calculated." [1] "Inverse sqrt moments calculated." [1] "PCC cumulants calculated." [1] "PCC Cornish Fisher coefficients calculated." [1] "Beginning root finding process for Cornish Fisher." [1] "First pruning complete. Removed 0 insignificant correlations." [1] "Second pruning complete. 91 correlations remain." [1] "Beginning root finding." [1] "Root finding complete." [1] "Estimating p-values." Error in mpfr(-(x[3]^2)/2, precBits = 128) : could not find function "mpfr" [1] "Error in mpfr(-(x[3]^2)/2, precBits = 128) : \n could not find function \"mpfr\"\n" [1] "Traceback:" [1] "6: ifelse(abs(x[3]) >= 38.4, as.double(-log10((0.5 exp(mpfr(-(x[3]^2)/2, "
[2] " precBits = 128))))), -log10(-x[4]/log(10)))"
[3] "5: ifelse(abs(x[3]) < 8.2, -log10(1 - exp(x[4])), ifelse(abs(x[3]) >= "
[4] " 38.4, as.double(-log10((0.5
exp(mpfr(-(x[3]^2)/2, precBits = 128))))), "
[5] " -log10(-x[4]/log(10))))"
[6] "4: FUN(newX[, i], ...)"
[7] "3: apply(p.matrix, 1, function(x) {"
[8] " ifelse(abs(x[3]) < 8.2, -log10(1 - exp(x[4])), ifelse(abs(x[3]) >= "
[9] " 38.4, as.double(-log10((0.5 * exp(mpfr(-(x[3]^2)/2, precBits = 128))))), " [10] " -log10(-x[4]/log(10))))"
[11] " })"
[12] "2: CF.PCC.pval(cor.roots)"
[13] "1: BigSur(example.seurat.subset, variable.features = F, correlations = T, "
[14] " depthlist = depths)"

silksmoothie commented 2 weeks ago

Thanks for letting us know! It looks like I forgot to include rmpfr as an import. This should be fixed. I will close the issue when I test that it works on Windows.