search

statdivlab / rigr

Regression, Inference, and General Data Analysis Tools for R
Other
10 stars 3 forks source link

null predictor space: "Error in zzs$coefficients[u, 1] : subscript out of bounds" #149

Open adw96 opened 1 year ago

adw96 commented 1 year ago

The cryptic error message Error in zzs$coefficients[u, 1] : subscript out of bounds arises when the single predictor variable is all zeroes. (This sounds trivial that may happen in automated large scale data analysis.) 

### ideal behaviour
lm(rnorm(10) ~ rep(0, 10))
glm(rpois(10, 10) ~ rep(0, 10), family="poisson")
rigr::regress("rate",
              formula = rpois(10, 10) ~ rep(0, 10) + rnorm(10),
              data=df)
rigr::regress("mean",
              formula = rpois(10, 10) ~ rep(0, 10) + rnorm(10),
              data=df)

### error with unhelpful error message
rigr::regress("rate",
              formula = rpois(10, 10) ~ rep(0, 10),
              data=df)
rigr::regress("mean",
              formula = rpois(10, 10) ~ rep(0, 10),
              data=df)