You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Hello!
This is not a real issue but just a suggestion.
I have noticided that EnvCpt is not robust to outliers (see code below). It might be useful to have an additional argument in envcpt() to choose if to perform a robust changepoint analysis or not (See 'robust' argument in the stl() function of the stats R package).
Thank you!
without outlier
library(EnvCpt)
y = c(rnorm(50,100,sd = 5),rnorm(50,50,sd = 5))
out = envcpt(y)
which.min(BIC(out))
plot(out$meancpt,type="p")
with outlier
y[20]=200
out = envcpt(y)
which.min(BIC(out))
plot(out$meancpt,type="p")
The text was updated successfully, but these errors were encountered:
Hello!
This is not a real issue but just a suggestion.
I have noticided that EnvCpt is not robust to outliers (see code below). It might be useful to have an additional argument in envcpt() to choose if to perform a robust changepoint analysis or not (See 'robust' argument in the stl() function of the stats R package).
Thank you!
without outlier
library(EnvCpt)
y = c(rnorm(50,100,sd = 5),rnorm(50,50,sd = 5))
out = envcpt(y)
which.min(BIC(out))
plot(out$meancpt,type="p")
with outlier
y[20]=200
out = envcpt(y)
which.min(BIC(out))
plot(out$meancpt,type="p")
The text was updated successfully, but these errors were encountered: