[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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._