Displaying 1 result from an estimated 1 matches for "t_rexp".
2000 Feb 19
1
problem with glm (PR#452)
...th gcc. In both cases I've built with
g77.
I want to do some simulations with a call to glm in each rep. A simplified
version would be something like:
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.:
REP 1
REP 2
<snip>
REP 9
Error in model.matrix(t, data) : variable lengths differ
REP 1...