Error in my.RSS.table[as.character(i), 3:4] <- c(pot.index[opt], break.RSS[opt]) : replacement has length zero
Nobody has claimed this yet.
Assessment
- Difficulty
- 4/5
- Estimated time
- 3-5 days
- Newbie friendliness
- 35/100
Research direction
Start by reproducing the reported breakpoints() call with the attached CSV, using bfastpp() and bfastts() as shown. Trace how recresid(), X1, break.RSS, and opt produce the NaN and empty replacement; done means the example no longer fails with a replacement-length error and the behavior for this edge case is covered by a regression test.
Written by the indexing model from the issue text.
Description
This is a very odd edge case. With a specific h value and specifically when using both trend and harmon components of bfastpp, for specific data, I get the error:
Error in my.RSS.table[as.character(i), 3:4] <- c(pot.index[opt], break.RSS[opt]) :
replacement has length zero
Also a lot of warnings about sqrt() producing NaN values. Looks like that line fails because opt is an integer(0), and that is because break.RSS is a single NaN (which probably comes from the sqrt problem, which in turn comes from fr in recresid being a negative value, which comes from X1 having negative values, and I don't know what X1 stands for).
To reproduce, see attached CSV, and then run the code:
breakpoints(formula(response~trend+harmon), data=
bfastpp(bfastts(CSV$evi, CSV$time), "irregular"), order=3), h=46)
- Dominant language
- R
- Stars
- 6
- Forks
- 5
- PR merge metrics
- No merged PRs in 30d
Contributor guide
No contributing guide indexed for this repository
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
More from bfast2/strucchangeRcpp
-
Tweak at the end Openquestion
Difficulty 4/5 3-5 days Newbie friendliness 20/100
bfast2/strucchangeRcpp#38 · 1 comment ·
-
enhancement
bfast2/strucchangeRcpp#32 · 2 comments · 1 assignee ·
-
bug help wanted
Difficulty 5/5 Over a week Newbie friendliness 28/100
bfast2/strucchangeRcpp#22 · 2 comments ·
-
help wanted question
Difficulty 5/5 Over a week Newbie friendliness 25/100
bfast2/strucchangeRcpp#20 · 1 comment ·
-
enhancement
Difficulty 4/5 3-5 days Newbie friendliness 35/100
bfast2/strucchangeRcpp#16 ·
All issues in bfast2/strucchangeRcpp
Similar issues
-
documentation pkg infrastructure
Difficulty 2/5 1-3 hours Newbie friendliness 76/100
epiverse-trace/epiparameter#511 ·
-
function:write_dwc
Difficulty 2/5 1-3 hours Newbie friendliness 72/100
-
Urgent request: Due to vulnerabilities move to API version 12.6.1 (12.6.2 eventually) or 13.1.1 Open
Difficulty 2/5 1-3 hours Newbie friendliness 65/100
jbkunst/highcharter#849 ·
-
Difficulty 2/5 1-3 hours Newbie friendliness 82/100
r-lib/pkgdepends#485 · 3 comments ·
-
Difficulty 1/5 Under an hour Newbie friendliness 92/100