Dear S+ users,
Yesterday I sent the following message:
*********************
I met for the first time a strange problem with plot.gam() which I often
use. I want to plot a simple gam (logistic) model with the following data:
> test <- gam(mort ~ lo(tmin, 1/3), data=ovlouga, family=binomial(logit))
> plot.gam(test)
Error: Expressions nested beyond limit (256) -- increase limit with
options(expressions=...)
only 30 of 213 frames dumped
Dumped
There is no problem with the fitted model, and the increase of
options(expression=n) up to 4096 does not change anything.
Did anyone meet the same problem ? How to fix it ?
********************
In fact, it seems there is a conflict between one of the loaded library and
the plot.gam() function:
> m1 <- gam(mort ~ lo(pds, 1/4), data=ovlouga, family=binomial(logit))
> plot.gam(m1) # no problem
> library(yags)
> plot.gam(m1) # error message
Error: Expressions nested beyond limit (256) -- increase limit with
options(expressions=...)
only 30 of 214 frames dumped
Dumped
The yags library is a gee solver. The problem does not seem to occur with
other libraries.
Could someone explain me why this happens, and how to prevent it ?
Thanks in advance,
Renaud
Renaud Lancelot
Programme PPR
ISRA-URA II
BP 2057
Dakar-Hann (Senegal)
Tel (221) 832 49 02
Fax (221) 821 18 79 / (221) 832 33 34
-----------------------------------------------------------------------
This message was distributed by s-news@wubios.wustl.edu. To unsubscribe
send e-mail to s-news-request@wubios.wustl.edu with the BODY of the
message: unsubscribe s-news
|