# [Rd] problem with glm (PR#452)

**Peter Dalgaard BSA
**
p.dalgaard@biostat.ku.dk

*19 Feb 2000 13:59:27 +0100*

mcneney@cs.sfu.ca writes:
>* N_100; NREPS_100; lam.e_5; logit.inv_function(x){ exp(x)/(1+exp(x)) }
*>* for(i in (1:NREPS)) {
*>* cat(paste("REP",i,"\n"))
*>* T_rexp(N,(lam.e))
*>* Delta_rbinom(N,1,logit.inv(1+T))
*>* ff_glm(Delta ~ T,family=binomial(logit))
*>* }
*>*
*>* It often completes all reps, but I occasionally get errors related to
*>* lengths of various objects, etc.:
*...
>* --please do not edit the information below--
*>*
*>* Version:
*>* platform = alpha-dec-osf5.0
*>* arch = alpha
*>* os = osf5.0
*>* system = alpha, osf5.0
*>* status = Patched
*>* major = 0
*>* minor = 99.0
*>* year = 2000
*>* month = February
*>* day = 9
*>* language = R
*>*
*>* Search Path:
*>* .GlobalEnv, Autoloads, package:base
*
I just did the above ~10 times on i386 and nothing happened. We've had
previous reports on trouble in the lm/glm/model.matrix etc. region on
Alpha.
If you compile with all warnings on, are you getting warnings on
pointer/integer comparisons or assignments?
Can you create a reproducible case (with set.seed to keep the random
generator sequence constant) and would you know how run a debugger on
your version so that we can figure out what is going wrong?
--
O__ ---- Peter Dalgaard Blegdamsvej 3
c/ /'_ --- Dept. of Biostatistics 2200 Cph. N
(*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard@biostat.ku.dk) FAX: (+45) 35327907
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._