[R] GLMMGibbs crashes on seeds data

[Peter Baker]

Hi all

I know GLMMGibbs is still in beta but has anyone experienced (and
solved ;-) this problem?

I decided to look at the seeds example but I get a core dump on two
intel linux boxes and also a sun workstation. All are running R1.3.0
but different hardware/OS's so I think I've done something wrong

> library(GLMMGibbs)
> data(seeds)
> seeds$plate <- as.factor(1:21)
> seeds$plate <- Ra(data=seeds$plate, shape=0.001, scale=0.001,
+                   type="identity", contrast="sum", of.interest=F)
> gseeds.1 <- glmm(formula=cbind(r,n-r) ~ x1 + x2, data=seeds, family=binomial,
+                store.results=T, icm=500, burnin=1000, keep=20000,
+                progress.info=500, thin=5)
Segmentation fault (core dumped)

> version
         _
platform i686-pc-linux-gnu
arch     i686
os       linux-gnu
system   i686, linux-gnu
status
major    1
minor    3.0
year     2001
month    06
day      22
language R

$ cat /etc/issue
Red Hat Linux release 6.2 (Zoot)
Kernel 2.2.16-3 on an i686

Suggestions welcome

Thanks
Peter

-- 
Dr Peter Baker, Statistician (Bioinformatics/Genetics),
CSIRO Mathematical & Information Sciences, 233 Middle St/
PO Box 120, CLEVELAND, QLD 4163 Australia.
Email: Peter.Baker at cmis.csiro.au
Phone:+61 7 3826 7263 Fax:+61 7 3826 7304



-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help 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-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._