[S] error message with plot.gam

Renaud LANCELOT (lancelot@telecomplus.sn)
Fri, 24 Jul 1998 09:23:45 +0100


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