search

khabbazian / l1ou

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

Error: object 'elist.min' not found when estimating convergent regimes using "pBIC" and "pBICess" #39

Closed Gonzobilbao closed 4 years ago

Gonzobilbao commented 6 years ago

Hi,

I am running the following piece of code

lasso1TreeB1_pBIC<-estimate_shift_configuration(AdjTree, Y[,1], criterion="pBIC")
lasso2TreeB1_pBIC <- estimate_convergent_regimes(lasso1TreeB1_pBIC, criterion="pBIC")

I get the following error

Error: object 'elist.min' not found

This seems to be also the case when I choose the criterion "pBICess" for estimating shifts. Is there any problem with my code?

Also, this is probably a silly question, but is the criterion for estimating shifts and convergence supposed to be the same? Or can it be, say pBIC for estimating shifts and BIC/AICc/ for estimating convergence (I illustrate it with an example bellow):

lasso1TreeB1_pBIC<-estimate_shift_configuration(AdjTree, Y[,1], criterion="pBIC")
lasso2TreeB1_pBIC <- estimate_convergent_regimes(lasso1TreeB1_pBIC, criterion="AICc")

Does this make any sense? Thank you!

Gonzalo

khabbazian commented 6 years ago

Hi Gonzalo, Sorry for the delay. The code seems correct. I guess the method only finds one shift (or maybe no shifts). Then when you run estimate_convergent_regimes you get the error. Please let me know if that's the case.

Gonzobilbao commented 6 years ago

Hi Mohammed, thanks for the answer. Yes, so no shifts were found by estimate_shift_configuration, that was why the error popped up. While doing this, I realized that the function summary.l1ou is not there, although it is on the manual. Maybe an issue with the manual installation?

Gonzalo