Inla Program crashing

[K.e]

Hi All,

The following code is causing my INLA program to crash.

ormula=(PfPR ~ -1 + Intercept + Long+ Lat + prec+ f(field, model = spde, group = field.group, control.group = list(model="ar",order=1)))

Ntrials=c(mau.data$Ex, rep(NA,10857))


result_mauritania_final2=inla(formula, data = inla.stack.data(stack, spde=spde),
                       family = "binomial", Ntrials=Ntrials,
                       control.predictor = list(A=inla.stack.A(stack),link = 1,
                                                compute=TRUE),control.compute=list(cpo=T,dic=TRUE), 
                       keep=FALSE,verbose=TRUE)

I am getting the following error. What could be the cause

Segmentation fault (core dumped)
Error in inla.inlaprogram.has.crashed() : 
  The inla-program exited with an error. Unless you interupted it yourself, please rerun with verbose=TRUE and check the output carefully.
  If this does help; please contact the developers at <he...@r-inla.org>.

[INLA help]

On Tue, 2014-01-21 at 23:33 -0800, K.e wrote:
> Hi All,
>
> The following code is causing my INLA program to crash.

Hi,

we need more info, like if you run with option inla(...,verbose=TRUE),
and if we need to rerun it here, we'll need code/data so we can do
that.

thanks

H

http://www.r-inla.org/contact-us#TOC-Report-an-error
--
Håvard Rue
he...@r-inla.org

[K.e]

Thanks for reply,

I am running with verbose=TRUE.  The model is running well in windows server but crashes when run in Linux server.
I have realized  the  problem may have arisen  because  I terminated a  running code which   I had specified  control.group as  "ar2" instead of "ar1".
How can I correct this problem.

Many Thanks

[Håvard Rue]

On Wed, 2014-01-22 at 01:24 -0800, K.e wrote:
> Thanks for reply,
>
> I am running with verbose=TRUE. The model is running well in windows
> server but crashes when run in Linux server.
> I have realized the problem may have arisen because I terminated a
> running code which I had specified control.group as "ar2" instead
> of "ar1".
> How can I correct this problem.

if the (group) model was set to "ar2" which does not exists, you should
have got an error...

H

Read 3921 items
> r = inla(y~f(idx, group=g, control.group = list(model="ar2")), data =
data.frame(y=1:3,g=c(1,1,2), idx=c(1,2,1)))
Error in inla.model.properties.generic(inla.trim.family(model),
(mm[names(mm) == :
Name 'ar2' is unknown. Valid choices are: exchangeable ar1 ar rw1 rw2
besag.


--
Håvard Rue
Department of Mathematical Sciences
Norwegian University of Science and Technology
N-7491 Trondheim, Norway
Voice: +47-7359-3533 URL : http://www.math.ntnu.no/~hrue
Mobile: +47-9260-0021 Email: havar...@math.ntnu.no

R-INLA: www.r-inla.org

[K.e]

Thanks very much Rue. I had specified control.group  wrongly. INLA is not crashing after implementing changes you have recommended