Displaying 2 results from an estimated 2 matches for "jnext".
Did you mean:
next
2017 Oct 09
1
Using response variable in interaction as explanatory variable in glm crashes R
...fine in the integer case.
Part of the culprit code may be this (from line 717),
with the isLogical(.) which in our case, shifts the pointer by
1 in the call to firstfactor() :
int adj = isLogical(var_i)?1:0;
// avoid overflow of jstart * nn PR#15578
firstfactor(&rx[jstart * nn], n, jnext - jstart,
REAL(contrast), nrows(contrast),
ncols(contrast), INTEGER(var_i)+adj);
then in firstfactor(), we see the segfault (when running R with
'-d gdb') :
> model.matrix(dob_mon ~ dob_day*dob_mon, data = tab)
Program received signal SIGSEGV, Segmentation fault....
2017 Oct 06
2
Using response variable in interaction as explanatory variable in glm crashes R
The following code crashes R (I know I shouldn't try to estimate such a
model; this was a bug in some code of mine). I also tried with R-devel;
same result.
tab <- structure(list(dob_day = c(FALSE, FALSE, FALSE, FALSE, TRUE,
TRUE, TRUE, TRUE), dob_mon = c(FALSE, FALSE, TRUE, TRUE, FALSE,
FALSE, TRUE, TRUE), dob_year = c(FALSE, TRUE, FALSE, TRUE, FALSE,
TRUE, FALSE, TRUE), n =