khabbazian / l1ou

Detection of evolutionary shifts in Ornstein-Uhlenbeck models

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

problem when runing estimate_convergent_regimes with pBIC option

Gonzobilbao opened this issue · comments

Hi,
I am running the following piece of code

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

[1] 4

lasso2TreeB1_pBIC <- estimate_convergent_regimes(lasso1TreeB1_pBIC, criterion="pBIC")

Error: $ operator is invalid for atomic vectors

However, when I use AICc or BIC instead of pBIC, it seems to work and $nShifts gives:

[1] 4

I have checked the manual for the function but not sure what is not working. Is there any problem with pBIC when using estimate_convergent_regimes? Something I am overlooking? Thank you!

close it because of not having enough info to fix to reproduce ...