khabbazian / l1ou

Detection of evolutionary shifts in Ornstein-Uhlenbeck models
GNU General Public License v2.0
19 stars 7 forks source link

error related with R version #53

Open nhazzi opened 1 year ago

nhazzi commented 1 year ago

Dear developers of l1ou package

I have use your package l1ou before with previous data, but now I am getting this error: Error in cmp_sqrt_OU_covariance(my.edge.list, length(tre$tip.label), tre$root.edge) : internal error!@sqrt_OU_covariance.cpp.

I am almost certain that the error has to be related with R version, because after I got the error, I tried to run the analysis in a older R version and it work. But now week after, I am getting the error again, so I think that is some kind of bug error between the package and the R version. I would really appreciatte any help that you can provided me, because this analysis is a very important part of my PhD.

thanks in advance

cecileane commented 1 year ago

Sorry to hear that. issue #50 was also caused by a change in the version of R.

@nhazzi : it would help if you could document under which version of R you got or didn't get the error --for anyone who may run into the same error, or anyone with the time to troubleshoot the error. Also, it would help to post a minimal working example, to help someone else troubleshoot the problem, so they can reproduce the error first.

nhazzi commented 1 year ago

Hi. thanks for your fast response.

I am using R-4.2.2.

and this is the code that I am using:

nstall_github("glmgen/genlasso") install_github("khabbazian/l1ou") install.packages('l1ou') library(l1ou) library(genlasso) library(ape) library(phytools) library(phangorn)

dat<- read.csv("pca.csv",row.names = 1) tree<-read.tree("tree_pruned_final") is.ultrametric(tree) nnls<-nnls.tree(cophenetic(tree),tree,rooted=TRUE) is.ultrametric(nnls) ltree<-adjust_data(nnls,dat[,]) shifts<-estimate_shift_configuration(ltree$tree, ltree$Y)

thanks in advance

Nicolas