Runtime error with spatio-temporal model

[Adam Green]

I am attempting to fit a poisson spatial model and I keep getting the "This application has requested the runtime to terminate in an unusual way: please contact the application's support team" error that it seems some others on the forum have encountered before. I have tried some of the other suggestions from previous posts like including control.fixed = list(prec.intercept = 0.01) and control.inla(diagonal=100) in the inla call. I also tried to get starting values by fitting the model with control.inla=list(int.strategy='eb',strategy='gaussian') but that doesn't run either. Here is my code, and I can provide the data and graph files if necessary. Thanks for your help.

formula.tmp=max.male~s_6400+f(year,model='ar1',replicate=FID2)+f(FID2,model='besag',graph='W.dat.txt',replicate=year)
formula.1=as.formula(formula.tmp)

inla.pois=inla(formula.1,family='poisson',data=data,verbose=T,control.inla=list(int.strategy='eb',strategy='gaussian'))

Adam Green

[INLA help]

On Fri, 2014-07-25 at 12:30 -0700, Adam Green wrote:
> I am attempting to fit a poisson spatial model and I keep getting the
> "This application has requested the runtime to terminate in an unusual
> way: please contact the application's support team" error that it
> seems some others on the forum have encountered before. I have tried
> some of the other suggestions from previous posts like including
> control.fixed = list(prec.intercept = 0.01) and
> control.inla(diagonal=100) in the inla call. I also tried to get
> starting values by fitting the model with
> control.inla=list(int.strategy='eb',strategy='gaussian') but that
> doesn't run either. Here is my code, and I can provide the data and
> graph files if necessary. Thanks for your help.
>
> formula.tmp=max.male~s_6400
> +f(year,model='ar1',replicate=FID2)+f(FID2,model='besag',graph='W.dat.txt',replicate=year)
> formula.1=as.formula(formula.tmp)

this is a 'shot in the dark', but the argument of 'replicate' is
integers. if max(year)=2014, say, then ALL replications between 1...2014
are generated (internally) as well. so if you have years 2010,2011,2012,
say, then code them as 1,2 and 3 and pass these into replicate. this
might be the reason is 'crashing' due to lack of memory.

H


--
Håvard Rue
he...@r-inla.org

[Adam Green]

Thanks for your help. That seemed to fix the problem. I'm still getting an error (status 253) when I fit the model with my entire data set. It works with a subset of the data though, so I think it's just a memory issue.

Adam

[Adam Green]

Ok. Sorry to keep bothering you, but a computer with more memory didn't fix the problem. I have 2176 sites, and inla runs on a subset of at least 2000 of those sites. Maybe it's not a memory issue after all, especially since I only have about 23,500 observations. I have pasted the output and error messages below. Any ideas you have would be greatly appreciated. Thanks.

Adam

        hgid: 85351f348950  date: Fri Jul 04 08:57:14 2014 +0100
Report bugs to <he...@r-inla.org>
Processing file [C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/Model.ini] max_threads=[1]
inla_build...
        number of sections=[8]
        parse section=[0] name=[INLA.Model] type=[PROBLEM]
        inla_parse_problem...
                name=[INLA.Model]
                openmp.strategy=[default]
        store results in directory=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/results.files]
                output:
                        cpo=[0]
                        po=[0]
                        dic=[0]
                        kld=[1]
                        mlik=[1]
                        q=[0]
                        graph=[0]
                        gdensity=[0]
                        hyperparameters=[1]
                        summary=[1]
                        return.marginals=[1]
                        nquantiles=[3]  [ 0.025 0.5 0.975 ]
                        ncdf=[0]  [ ]
        parse section=[2] name=[Predictor] type=[PREDICTOR]
        inla_parse_predictor ...
                section=[Predictor]
                dir=[predictor]
                PRIOR->name=[loggamma]
                PRIOR->from_theta=[function (x) <<NEWLINE>>exp(x)]
                PRIOR->to_theta = [function (x) <<NEWLINE>>log(x)]
                PRIOR->PARAMETERS=[1, 1e-005]
                initialise log_precision[11]
                fixed=[1]
                user.scale=[1]
                n=[23383]
                m=[0]
                ndata=[23383]
                compute=[0]
                read offsets from file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f04f907b6a]
                read n=[46766] entries from file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f04f907b6a]
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f04f907b6a] 0/23383  (idx,y) = (0, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f04f907b6a] 1/23383  (idx,y) = (1, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f04f907b6a] 2/23383  (idx,y) = (2, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f04f907b6a] 3/23383  (idx,y) = (3, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f04f907b6a] 4/23383  (idx,y) = (4, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f04f907b6a] 5/23383  (idx,y) = (5, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f04f907b6a] 6/23383  (idx,y) = (6, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f04f907b6a] 7/23383  (idx,y) = (7, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f04f907b6a] 8/23383  (idx,y) = (8, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f04f907b6a] 9/23383  (idx,y) = (9, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f04f907b6a] 10/23383  (idx,y) = (10, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f04f907b6a] 11/23383  (idx,y) = (11, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f04f907b6a] 12/23383  (idx,y) = (12, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f04f907b6a] 13/23383  (idx,y) = (13, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f04f907b6a] 14/23383  (idx,y) = (14, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f04f907b6a] 15/23383  (idx,y) = (15, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f04f907b6a] 16/23383  (idx,y) = (16, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f04f907b6a] 17/23383  (idx,y) = (17, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f04f907b6a] 18/23383  (idx,y) = (18, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f04f907b6a] 19/23383  (idx,y) = (19, 0)
                Aext=[(null)]
                AextPrecision=[1e+008]
                output:
                        summary=[1]
                        return.marginals=[1]
                        nquantiles=[3]  [ 0.025 0.5 0.975 ]
                        ncdf=[0]  [ ]
        parse section=[1] name=[INLA.Data1] type=[DATA]
        inla_parse_data [section 1]...
                tag=[INLA.Data1]
                family=[POISSON]
                likelihood=[POISSON]
                file->name=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f01ca468e6]
                file->name=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f0733f345]
                read n=[70149] entries from file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f01ca468e6]
                        0/23383  (idx,a,y,d) = (0, 1, 18, 1)
                        1/23383  (idx,a,y,d) = (1, 1, 17, 1)
                        2/23383  (idx,a,y,d) = (2, 1, 19, 1)
                        3/23383  (idx,a,y,d) = (3, 1, 14, 1)
                        4/23383  (idx,a,y,d) = (4, 1, 30, 1)
                        5/23383  (idx,a,y,d) = (5, 1, 18, 1)
                        6/23383  (idx,a,y,d) = (6, 1, 23, 1)
                        7/23383  (idx,a,y,d) = (7, 1, 52, 1)
                        8/23383  (idx,a,y,d) = (8, 1, 25, 1)
                        9/23383  (idx,a,y,d) = (9, 1, 23, 1)
                        10/23383  (idx,a,y,d) = (10, 1, 37, 1)
                        11/23383  (idx,a,y,d) = (11, 1, 26, 1)
                        12/23383  (idx,a,y,d) = (12, 1, 17, 1)
                        13/23383  (idx,a,y,d) = (13, 1, 25, 1)
                        14/23383  (idx,a,y,d) = (14, 1, 63, 1)
                        15/23383  (idx,a,y,d) = (15, 1, 18, 1)
                        16/23383  (idx,a,y,d) = (16, 1, 17, 1)
                        17/23383  (idx,a,y,d) = (17, 1, 15, 1)
                        18/23383  (idx,a,y,d) = (18, 1, 25, 1)
                        19/23383  (idx,a,y,d) = (19, 1, 15, 1)
                use variant [0]
                        bit 0 is off
                        bit 1 is off
                        bit 2 is off
                        bit 3 is off
                Link model  [LOG]
                Link order  [-1]
                Link ntheta [0]
                mix.use[0]
        parse section=[5] name=[year2] type=[FFIELD]
        inla_parse_ffield...
                section=[year2]
                dir=[random.effect00000001]
                model=[ar1]
                PRIOR0->name=[loggamma]
                PRIOR0->from_theta=[function (x) <<NEWLINE>>exp(x)]
                PRIOR0->to_theta = [function (x) <<NEWLINE>>log(x)]
                PRIOR0->PARAMETERS0=[1, 5e-005]
                PRIOR1->name=[normal]
                PRIOR1->from_theta=[function (x) <<NEWLINE>>2 * exp(x)/(1 + exp(x)) - 1]
                PRIOR1->to_theta = [function (x) <<NEWLINE>>log((1 + x)/(1 - x))]
                PRIOR1->PARAMETERS1=[0, 0.15]
                constr=[0]
                diagonal=[0]
                id.names=<not present>
                compute=[1]
                nrep=[2176]
                ngroup=[1]
                read covariates from file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f02ac31763]
                read n=[46766] entries from file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f02ac31763]
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f02ac31763] 0/23383  (idx,y) = (0, 43312)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f02ac31763] 1/23383  (idx,y) = (1, 43316)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f02ac31763] 2/23383  (idx,y) = (2, 43319)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f02ac31763] 3/23383  (idx,y) = (3, 43315)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f02ac31763] 4/23383  (idx,y) = (4, 43306)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f02ac31763] 5/23383  (idx,y) = (5, 43317)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f02ac31763] 6/23383  (idx,y) = (6, 43324)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f02ac31763] 7/23383  (idx,y) = (7, 43307)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f02ac31763] 8/23383  (idx,y) = (8, 43303)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f02ac31763] 9/23383  (idx,y) = (9, 43323)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f02ac31763] 10/23383  (idx,y) = (10, 43310)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f02ac31763] 11/23383  (idx,y) = (11, 43305)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f02ac31763] 12/23383  (idx,y) = (12, 43322)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f02ac31763] 13/23383  (idx,y) = (13, 43320)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f02ac31763] 14/23383  (idx,y) = (14, 43308)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f02ac31763] 15/23383  (idx,y) = (15, 43313)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f02ac31763] 16/23383  (idx,y) = (16, 43321)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f02ac31763] 17/23383  (idx,y) = (17, 43314)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f02ac31763] 18/23383  (idx,y) = (18, 43304)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f02ac31763] 19/23383  (idx,y) = (19, 43318)
                n=[29]
                cyclic=[0]
                file for locations=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f06fa676f0]
                        nlocations=[29]
                        locations[0]=[1]
                        locations[1]=[2]
                        locations[2]=[3]
                        locations[3]=[4]
                        locations[4]=[5]
                        locations[5]=[6]
                        locations[6]=[7]
                        locations[7]=[8]
                        locations[8]=[9]
                        locations[9]=[10]
                        locations[10]=[11]
                        locations[11]=[12]
                        locations[12]=[13]
                        locations[13]=[14]
                        locations[14]=[15]
                        locations[15]=[16]
                        locations[16]=[17]
                        locations[17]=[18]
                        locations[18]=[19]
                        locations[19]=[20]
                initialise log_precision[4]
                fixed=[0]
                initialise phi_intern[2]
                fixed=[0]
                computed/guessed rank-deficiency = [0]
                output:
                        summary=[1]
                        return.marginals=[1]
                        nquantiles=[3]  [ 0.025 0.5 0.975 ]
                        ncdf=[0]  [ ]
        parse section=[6] name=[FID2] type=[FFIELD]
        inla_parse_ffield...
                section=[FID2]
                dir=[random.effect00000002]
                model=[besag]
                PRIOR->name=[loggamma]
                PRIOR->from_theta=[function (x) <<NEWLINE>>exp(x)]
                PRIOR->to_theta = [function (x) <<NEWLINE>>log(x)]
                PRIOR->PARAMETERS=[1, 5e-005]
                constr=[0]
                diagonal=[9.99978e-007]
                id.names=<not present>
                compute=[1]
                nrep=[29]
                ngroup=[1]
                read covariates from file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f06b9bd5]
                read n=[46766] entries from file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f06b9bd5]
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f06b9bd5] 0/23383  (idx,y) = (0, 34133)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f06b9bd5] 1/23383  (idx,y) = (1, 42837)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f06b9bd5] 2/23383  (idx,y) = (2, 49365)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f06b9bd5] 3/23383  (idx,y) = (3, 40661)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f06b9bd5] 4/23383  (idx,y) = (4, 21077)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f06b9bd5] 5/23383  (idx,y) = (5, 45013)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f06b9bd5] 6/23383  (idx,y) = (6, 60245)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f06b9bd5] 7/23383  (idx,y) = (7, 23253)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f06b9bd5] 8/23383  (idx,y) = (8, 14549)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f06b9bd5] 9/23383  (idx,y) = (9, 58069)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f06b9bd5] 10/23383  (idx,y) = (10, 29781)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f06b9bd5] 11/23383  (idx,y) = (11, 18901)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f06b9bd5] 12/23383  (idx,y) = (12, 55893)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f06b9bd5] 13/23383  (idx,y) = (13, 51541)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f06b9bd5] 14/23383  (idx,y) = (14, 25429)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f06b9bd5] 15/23383  (idx,y) = (15, 36309)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f06b9bd5] 16/23383  (idx,y) = (16, 53717)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f06b9bd5] 17/23383  (idx,y) = (17, 38485)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f06b9bd5] 18/23383  (idx,y) = (18, 16725)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f06b9bd5] 19/23383  (idx,y) = (19, 47189)
                read graph from file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f08619e2]
                file for locations=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f070e3d05]
                        nlocations=[2176]
                        locations[0]=[1]
                        locations[1]=[2]
                        locations[2]=[3]
                        locations[3]=[4]
                        locations[4]=[5]
                        locations[5]=[6]
                        locations[6]=[7]
                        locations[7]=[8]
                        locations[8]=[9]
                        locations[9]=[10]
                        locations[10]=[11]
                        locations[11]=[12]
                        locations[12]=[13]
                        locations[13]=[14]
                        locations[14]=[15]
                        locations[15]=[16]
                        locations[16]=[17]
                        locations[17]=[18]
                        locations[18]=[19]
                        locations[19]=[20]
                initialise log_precision[4]
                fixed=[0]
                adjust.for.con.comp[1]
                scale.model[0]
                read extra constraint from file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f0479144a3]
                Constraint[0]
                        A[0] = 1.000000
                        A[1] = 1.000000
                        A[2] = 1.000000
                        A[3] = 1.000000
                        A[4] = 1.000000
                        A[5] = 1.000000
                        A[6] = 1.000000
                        A[7] = 1.000000
                        A[8] = 1.000000
                        A[9] = 1.000000
                        A[10] = 1.000000
                        A[11] = 1.000000
                        A[12] = 1.000000
                        A[13] = 1.000000
                        A[14] = 1.000000
                        A[15] = 1.000000
                        A[16] = 1.000000
                        A[17] = 1.000000
                        A[18] = 1.000000
                        A[19] = 1.000000
                        A[20] = 1.000000
                        e[0] = 0.000000
                Constraint[1]
                        A[264] = 1.000000
                        A[265] = 1.000000
                        A[266] = 1.000000
                        e[1] = 0.000000
                Constraint[2]
                        A[688] = 1.000000
                        A[714] = 1.000000
                        e[2] = 0.000000
                Constraint[3]
                        A[1118] = 1.000000
                        A[1120] = 1.000000
                        A[1123] = 1.000000
                        A[1124] = 1.000000
                        A[1131] = 1.000000
                        e[3] = 0.000000
                Constraint[4]
                        A[1493] = 1.000000
                        A[1494] = 1.000000
                        A[1495] = 1.000000
                        A[1496] = 1.000000
                        A[1497] = 1.000000
                        A[1498] = 1.000000
                        A[1501] = 1.000000
                        A[1502] = 1.000000
                        A[1503] = 1.000000
                        A[1504] = 1.000000
                        A[1505] = 1.000000
                        A[1506] = 1.000000
                        A[1507] = 1.000000
                        e[4] = 0.000000
                Constraint[5]
                        A[1499] = 1.000000
                        A[1500] = 1.000000
                        e[5] = 0.000000
                rank-deficiency is *defined* [10]
                output:
                        summary=[1]
                        return.marginals=[1]
                        nquantiles=[3]  [ 0.025 0.5 0.975 ]
                        ncdf=[0]  [ ]
        section=[3] name=[(Intercept)] type=[LINEAR]
        inla_parse_linear...
                section[(Intercept)]
                dir=[fixed.effect00000001]
                file for covariates=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f05cb36e5c]
                read n=[46766] entries from file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f05cb36e5c]
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f05cb36e5c] 0/23383  (idx,y) = (0, 1)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f05cb36e5c] 1/23383  (idx,y) = (1, 1)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f05cb36e5c] 2/23383  (idx,y) = (2, 1)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f05cb36e5c] 3/23383  (idx,y) = (3, 1)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f05cb36e5c] 4/23383  (idx,y) = (4, 1)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f05cb36e5c] 5/23383  (idx,y) = (5, 1)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f05cb36e5c] 6/23383  (idx,y) = (6, 1)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f05cb36e5c] 7/23383  (idx,y) = (7, 1)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f05cb36e5c] 8/23383  (idx,y) = (8, 1)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f05cb36e5c] 9/23383  (idx,y) = (9, 1)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f05cb36e5c] 10/23383  (idx,y) = (10, 1)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f05cb36e5c] 11/23383  (idx,y) = (11, 1)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f05cb36e5c] 12/23383  (idx,y) = (12, 1)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f05cb36e5c] 13/23383  (idx,y) = (13, 1)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f05cb36e5c] 14/23383  (idx,y) = (14, 1)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f05cb36e5c] 15/23383  (idx,y) = (15, 1)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f05cb36e5c] 16/23383  (idx,y) = (16, 1)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f05cb36e5c] 17/23383  (idx,y) = (17, 1)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f05cb36e5c] 18/23383  (idx,y) = (18, 1)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f05cb36e5c] 19/23383  (idx,y) = (19, 1)
                prior mean=[0]
                prior precision=[0]
                compute=[1]
                output:
                        summary=[1]
                        return.marginals=[1]
                        nquantiles=[3]  [ 0.025 0.5 0.975 ]
                        ncdf=[0]  [ ]
        section=[4] name=[da_800] type=[LINEAR]
        inla_parse_linear...
                section[da_800]
                dir=[fixed.effect00000002]
                file for covariates=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f03e361c2b]
                read n=[46766] entries from file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f03e361c2b]
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f03e361c2b] 0/23383  (idx,y) = (0, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f03e361c2b] 1/23383  (idx,y) = (1, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f03e361c2b] 2/23383  (idx,y) = (2, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f03e361c2b] 3/23383  (idx,y) = (3, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f03e361c2b] 4/23383  (idx,y) = (4, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f03e361c2b] 5/23383  (idx,y) = (5, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f03e361c2b] 6/23383  (idx,y) = (6, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f03e361c2b] 7/23383  (idx,y) = (7, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f03e361c2b] 8/23383  (idx,y) = (8, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f03e361c2b] 9/23383  (idx,y) = (9, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f03e361c2b] 10/23383  (idx,y) = (10, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f03e361c2b] 11/23383  (idx,y) = (11, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/data.files/file19f03e361c2b] 12/23383  (idx,y) = (12, 0)
                file=[C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19fError 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>.
In addition: Warning messages:
1: In INLA::f(FID2, model = "besag", graph = "W.dat.txt", replicate = year2) :
  The graph for the model besag has 8 connected components!!! Model is revised accordingly.
2: running command '"C:/Program Files/R/R-3.0.0/library/INLA/bin/windows/inla64.exe"  -b -s -v  -t 1  "C:/Users/agreen/AppData/Local/Temp/1/RtmpiQpyEv/file19f011397219/Model.ini"' had status 253

[Gustavo Machado]

I do have the same problem, can someone help?

> model.par <- inla(formula.par,family="poisson",data=data,E=exp, + control.predictor=list(compute=TRUE), + control.compute=list(dic=TRUE,cpo=TRUE),control.fixed = list(prec.intercept = 0.01),verbose=T,control.inla=list(int.strategy='eb',strategy='gaussian')) hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 Report bugs to <he...@r-inla.org> Processing file [C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/Model.ini] max_threads=[8] inla_build... number of sections=[9] parse section=[8] name=[INLA.Expert] type=[EXPERT] inla_parse_expert... section[INLA.Expert] disable.gaussian.check=[0] cpo.manual=[0] jp.R_HOME=[(null)] jp.Rfile=[(null)] jp.RData=NULL jp.func=[(null)] parse section=[0] name=[INLA.Model] type=[PROBLEM] inla_parse_problem... name=[INLA.Model] openmp.strategy=[default] store results in directory=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/results.files] output: cpo=[1] po=[0] dic=[1] kld=[1] mlik=[1] q=[0] graph=[0] gdensity=[0] hyperparameters=[1] summary=[1] return.marginals=[1] nquantiles=[3] [ 0.025 0.5 0.975 ] ncdf=[0] [ ] parse section=[2] name=[Predictor] type=[PREDICTOR] inla_parse_predictor ... section=[Predictor] dir=[predictor] PRIOR->name=[loggamma] hyperid=[53001|Predictor] PRIOR->from_theta=[function (x) <<NEWLINE>>exp(x)] PRIOR->to_theta = [function (x) <<NEWLINE>>log(x)] PRIOR->PARAMETERS=[1, 1e-005] initialise log_precision[11] fixed=[1] user.scale=[1] n=[44512] m=[0] ndata=[44512] compute=[1] read offsets from file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf485d222219] read n=[89024] entries from file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf485d222219] file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf485d222219] 0/44512 (idx,y) = (0, 0) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf485d222219] 1/44512 (idx,y) = (1, 0) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf485d222219] 2/44512 (idx,y) = (2, 0) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf485d222219] 3/44512 (idx,y) = (3, 0) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf485d222219] 4/44512 (idx,y) = (4, 0) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf485d222219] 5/44512 (idx,y) = (5, 0) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf485d222219] 6/44512 (idx,y) = (6, 0) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf485d222219] 7/44512 (idx,y) = (7, 0) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf485d222219] 8/44512 (idx,y) = (8, 0) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf485d222219] 9/44512 (idx,y) = (9, 0) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf485d222219] 10/44512 (idx,y) = (10, 0) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf485d222219] 11/44512 (idx,y) = (11, 0) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf485d222219] 12/44512 (idx,y) = (12, 0) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf485d222219] 13/44512 (idx,y) = (13, 0) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf485d222219] 14/44512 (idx,y) = (14, 0) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf485d222219] 15/44512 (idx,y) = (15, 0) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf485d222219] 16/44512 (idx,y) = (16, 0) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf485d222219] 17/44512 (idx,y) = (17, 0) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf485d222219] 18/44512 (idx,y) = (18, 0) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf485d222219] 19/44512 (idx,y) = (19, 0) Aext=[(null)] AextPrecision=[1e+008] output: summary=[1] return.marginals=[1] nquantiles=[3] [ 0.025 0.5 0.975 ] ncdf=[0] [ ] parse section=[1] name=[INLA.Data1] type=[DATA] inla_parse_data [section 1]... tag=[INLA.Data1] family=[POISSON] likelihood=[POISSON] file->name=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48c31fbd] file->name=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf486512653d] read n=[133536] entries from file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48c31fbd] 0/44512 (idx,a,y,d) = (0, 0.186571, 0, 1) 1/44512 (idx,a,y,d) = (1, 0.691096, 0, 1) 2/44512 (idx,a,y,d) = (2, 0.0482872, 0, 1) 3/44512 (idx,a,y,d) = (3, 0.601, 0, 1) 4/44512 (idx,a,y,d) = (4, 0.130252, 0, 1) 5/44512 (idx,a,y,d) = (5, 0.1422, 0, 1) 6/44512 (idx,a,y,d) = (6, 0.0671798, 0, 1) 7/44512 (idx,a,y,d) = (7, 0.104621, 0, 1) 8/44512 (idx,a,y,d) = (8, 0.219744, 0, 1) 9/44512 (idx,a,y,d) = (9, 0.31862, 0, 1) 10/44512 (idx,a,y,d) = (10, 0.397778, 0, 1) 11/44512 (idx,a,y,d) = (11, 0.891932, 0, 1) 12/44512 (idx,a,y,d) = (12, 0.238147, 0, 1) 13/44512 (idx,a,y,d) = (13, 0.152013, 0, 1) 14/44512 (idx,a,y,d) = (14, 0.290105, 0, 1) 15/44512 (idx,a,y,d) = (15, 0.258699, 0, 1) 16/44512 (idx,a,y,d) = (16, 1.#QNAN, 0, 1) 17/44512 (idx,a,y,d) = (17, 0.170715, 0, 1) 18/44512 (idx,a,y,d) = (18, 0.0253636, 0, 1) 19/44512 (idx,a,y,d) = (19, 0.387399, 0, 1) likelihood.variant=[0] Link model [LOG] Link order [-1] Link variant [-1] Link ntheta [0] mix.use[0] parse section=[5] name=[county] type=[FFIELD] inla_parse_ffield... section=[county] dir=[random.effect00000001] model=[bym] PRIOR0->name=[loggamma] hyperid=[10001|county] PRIOR0->from_theta=[function (x) <<NEWLINE>>exp(x)] PRIOR0->to_theta = [function (x) <<NEWLINE>>log(x)] PRIOR0->PARAMETERS0=[1, 0.0005] PRIOR1->name=[loggamma] hyperid=[10002|county] PRIOR1->from_theta=[function (x) <<NEWLINE>>exp(x)] PRIOR1->to_theta = [function (x) <<NEWLINE>>log(x)] PRIOR1->PARAMETERS1=[1, 0.0005] correct=[-1] constr=[0] diagonal=[9.99978e-007] id.names=<not present> compute=[1] nrep=[1] ngroup=[1] read covariates from file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4810b0311f] read n=[89024] entries from file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4810b0311f] file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4810b0311f] 0/44512 (idx,y) = (0, 0) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4810b0311f] 1/44512 (idx,y) = (1, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4810b0311f] 2/44512 (idx,y) = (2, 2) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4810b0311f] 3/44512 (idx,y) = (3, 3) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4810b0311f] 4/44512 (idx,y) = (4, 4) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4810b0311f] 5/44512 (idx,y) = (5, 5) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4810b0311f] 6/44512 (idx,y) = (6, 6) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4810b0311f] 7/44512 (idx,y) = (7, 7) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4810b0311f] 8/44512 (idx,y) = (8, 8) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4810b0311f] 9/44512 (idx,y) = (9, 9) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4810b0311f] 10/44512 (idx,y) = (10, 10) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4810b0311f] 11/44512 (idx,y) = (11, 11) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4810b0311f] 12/44512 (idx,y) = (12, 12) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4810b0311f] 13/44512 (idx,y) = (13, 13) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4810b0311f] 14/44512 (idx,y) = (14, 14) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4810b0311f] 15/44512 (idx,y) = (15, 15) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4810b0311f] 16/44512 (idx,y) = (16, 16) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4810b0311f] 17/44512 (idx,y) = (17, 17) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4810b0311f] 18/44512 (idx,y) = (18, 18) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4810b0311f] 19/44512 (idx,y) = (19, 19) read graph from file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4872fd57f0] file for locations=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48648c29fb] nlocations=[5564] locations[0]=[1] locations[1]=[2] locations[2]=[3] locations[3]=[4] locations[4]=[5] locations[5]=[6] locations[6]=[7] locations[7]=[8] locations[8]=[9] locations[9]=[10] locations[10]=[11] locations[11]=[12] locations[12]=[13] locations[13]=[14] locations[14]=[15] locations[15]=[16] locations[16]=[17] locations[17]=[18] locations[18]=[19] locations[19]=[20] initialise log_precision (iid component)[4] fixed=[0] initialise log_precision (spatial component)[4] fixed=[0] adjust.for.con.comp[0] scale.model[0] read extra constraint from file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4833c74030] Constraint[0] A[5564] = 1.000000 A[5565] = 1.000000 A[5566] = 1.000000 A[5567] = 1.000000 A[5568] = 1.000000 A[5569] = 1.000000 A[5570] = 1.000000 A[5571] = 1.000000 A[5572] = 1.000000 A[5573] = 1.000000 A[5574] = 1.000000 A[5575] = 1.000000 A[5576] = 1.000000 A[5577] = 1.000000 A[5578] = 1.000000 A[5579] = 1.000000 A[5580] = 1.000000 A[5581] = 1.000000 A[5582] = 1.000000 A[5583] = 1.000000 A[5584] = 1.000000 e[0] = 0.000000 computed/guessed rank-deficiency = [1] output: summary=[1] return.marginals=[1] nquantiles=[3] [ 0.025 0.5 0.975 ] ncdf=[0] [ ] parse section=[6] name=[county1] type=[FFIELD] inla_parse_ffield... section=[county1] dir=[random.effect00000002] model=[iid] PRIOR->name=[loggamma] hyperid=[1001|county1] PRIOR->from_theta=[function (x) <<NEWLINE>>exp(x)] PRIOR->to_theta = [function (x) <<NEWLINE>>log(x)] PRIOR->PARAMETERS=[1, 5e-005] correct=[-1] constr=[1] diagonal=[9.99978e-007] id.names=<not present> compute=[1] nrep=[1] ngroup=[1] read covariates from file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4851287760] read n=[89024] entries from file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4851287760] file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4851287760] 0/44512 (idx,y) = (0, 0) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4851287760] 1/44512 (idx,y) = (1, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4851287760] 2/44512 (idx,y) = (2, 2) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4851287760] 3/44512 (idx,y) = (3, 3) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4851287760] 4/44512 (idx,y) = (4, 4) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4851287760] 5/44512 (idx,y) = (5, 5) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4851287760] 6/44512 (idx,y) = (6, 6) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4851287760] 7/44512 (idx,y) = (7, 7) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4851287760] 8/44512 (idx,y) = (8, 8) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4851287760] 9/44512 (idx,y) = (9, 9) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4851287760] 10/44512 (idx,y) = (10, 10) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4851287760] 11/44512 (idx,y) = (11, 11) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4851287760] 12/44512 (idx,y) = (12, 12) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4851287760] 13/44512 (idx,y) = (13, 13) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4851287760] 14/44512 (idx,y) = (14, 14) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4851287760] 15/44512 (idx,y) = (15, 15) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4851287760] 16/44512 (idx,y) = (16, 16) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4851287760] 17/44512 (idx,y) = (17, 17) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4851287760] 18/44512 (idx,y) = (18, 18) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4851287760] 19/44512 (idx,y) = (19, 19) read weights from file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48594b27db] read n=[89024] entries from file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48594b27db] file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48594b27db] 0/44512 (idx,y) = (0, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48594b27db] 1/44512 (idx,y) = (1, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48594b27db] 2/44512 (idx,y) = (2, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48594b27db] 3/44512 (idx,y) = (3, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48594b27db] 4/44512 (idx,y) = (4, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48594b27db] 5/44512 (idx,y) = (5, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48594b27db] 6/44512 (idx,y) = (6, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48594b27db] 7/44512 (idx,y) = (7, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48594b27db] 8/44512 (idx,y) = (8, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48594b27db] 9/44512 (idx,y) = (9, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48594b27db] 10/44512 (idx,y) = (10, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48594b27db] 11/44512 (idx,y) = (11, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48594b27db] 12/44512 (idx,y) = (12, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48594b27db] 13/44512 (idx,y) = (13, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48594b27db] 14/44512 (idx,y) = (14, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48594b27db] 15/44512 (idx,y) = (15, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48594b27db] 16/44512 (idx,y) = (16, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48594b27db] 17/44512 (idx,y) = (17, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48594b27db] 18/44512 (idx,y) = (18, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48594b27db] 19/44512 (idx,y) = (19, 1) file for locations=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48702dfe7] nlocations=[5564] locations[0]=[1] locations[1]=[2] locations[2]=[3] locations[3]=[4] locations[4]=[5] locations[5]=[6] locations[6]=[7] locations[7]=[8] locations[8]=[9] locations[9]=[10] locations[10]=[11] locations[11]=[12] locations[12]=[13] locations[13]=[14] locations[14]=[15] locations[15]=[16] locations[16]=[17] locations[17]=[18] locations[18]=[19] locations[19]=[20] cyclic=[0] initialise log_precision[4] fixed=[0] computed/guessed rank-deficiency = [1] output: summary=[1] return.marginals=[1] nquantiles=[3] [ 0.025 0.5 0.975 ] ncdf=[0] [ ] section=[3] name=[(Intercept)] type=[LINEAR] inla_parse_linear... section[(Intercept)] dir=[fixed.effect00000001] file for covariates=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4867971e9d] read n=[89024] entries from file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4867971e9d] file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4867971e9d] 0/44512 (idx,y) = (0, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4867971e9d] 1/44512 (idx,y) = (1, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4867971e9d] 2/44512 (idx,y) = (2, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4867971e9d] 3/44512 (idx,y) = (3, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4867971e9d] 4/44512 (idx,y) = (4, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4867971e9d] 5/44512 (idx,y) = (5, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4867971e9d] 6/44512 (idx,y) = (6, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4867971e9d] 7/44512 (idx,y) = (7, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4867971e9d] 8/44512 (idx,y) = (8, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4867971e9d] 9/44512 (idx,y) = (9, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4867971e9d] 10/44512 (idx,y) = (10, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4867971e9d] 11/44512 (idx,y) = (11, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4867971e9d] 12/44512 (idx,y) = (12, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4867971e9d] 13/44512 (idx,y) = (13, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4867971e9d] 14/44512 (idx,y) = (14, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4867971e9d] 15/44512 (idx,y) = (15, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4867971e9d] 16/44512 (idx,y) = (16, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4867971e9d] 17/44512 (idx,y) = (17, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4867971e9d] 18/44512 (idx,y) = (18, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf4867971e9d] 19/44512 (idx,y) = (19, 1) prior mean=[0] prior precision=[0.01] compute=[1] output: summary=[1] return.marginals=[1] nquantiles=[3] [ 0.025 0.5 0.975 ] ncdf=[0] [ ] section=[4] name=[years] type=[LINEAR] inla_parse_linear... section[years] dir=[fixed.effect00000002] file for covariates=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48232474db] read n=[89024] entries from file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48232474db] file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48232474db] 0/44512 (idx,y) = (0, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48232474db] 1/44512 (idx,y) = (1, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48232474db] 2/44512 (idx,y) = (2, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48232474db] 3/44512 (idx,y) = (3, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48232474db] 4/44512 (idx,y) = (4, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48232474db] 5/44512 (idx,y) = (5, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48232474db] 6/44512 (idx,y) = (6, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48232474db] 7/44512 (idx,y) = (7, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48232474db] 8/44512 (idx,y) = (8, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48232474db] 9/44512 (idx,y) = (9, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48232474db] 10/44512 (idx,y) = (10, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48232474db] 11/44512 (idx,y) = (11, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48232474db] 12/44512 (idx,y) = (12, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48232474db] 13/44512 (idx,y) = (13, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48232474db] 14/44512 (idx,y) = (14, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48232474db] 15/44512 (idx,y) = (15, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48232474db] 16/44512 (idx,y) = (16, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48232474db] 17/44512 (idx,y) = (17, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48232474db] 18/44512 (idx,y) = (18, 1) file=[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/data.files/file3cf48232474db] 19/44512 (idx,y) = (19, 1) prior mean=[0] prior precision=[0.001] compute=[1] output: summary=[1] return.marginals=[1] nquantiles=[3] [ 0.025 0.5 0.975 ] ncdf=[0] [ ] Index table: number of entries[5], total length[61206] tag start-index length Predictor 0 44512 county 44512 11128 county1 55640 5564 (Intercept) 61204 1 years 61205 1 parse section=[7] name=[INLA.Parameters] type=[INLA] inla_parse_INLA... section[INLA.Parameters] lincomb.derived.only = [Yes] lincomb.derived.correlation.matrix = [No] global_node.factor = 2.000 global_node.degree = 2147483647 reordering = -1 Contents of ai_param 0302EEC8 Optimiser: DEFAULT METHOD Option for domin-BFGS: epsx = 0.005 Option for domin-BFGS: epsf = 1e-005 (rounding error) Option for domin-BFGS: epsg = 0.005 Option for GSL-BFGS2: tol = 0.1 Option for GSL-BFGS2: step_size = 1 Option for GSL-BFGS2: epsx = 0.005 Option for GSL-BFGS2: epsf = 0.000353553 Option for GSL-BFGS2: epsg = 0.005 Restart: 0 Mode known: No Gaussian approximation: abserr_func = 0.0005 abserr_step = 0.0005 optpar_fp = 0 optpar_nr_step_factor = -0.1 Gaussian data: No Strategy: Use the Gaussian approximation Fast mode: On Use linear approximation to log(|Q +c|)? Yes Method: Compute the derivative exact Parameters for improved approximations Number of points evaluate: 9 Step length to compute derivatives numerically: 0.000100002 Stencil to compute derivatives numerically: 5 Cutoff value to construct local neigborhood: 0.0001 Log calculations: On Log calculated marginal for the hyperparameters: On Integration strategy: Use only the modal configuration (EMPIRICAL_BAYES) f0 (CCD only): 1.100000 dz (GRID only): 0.750000 Adjust weights (GRID only): On Difference in log-density limit (GRID only): 4.000000 Skip configurations with (presumed) small density (GRID only): On Gradient is computed using Central difference with step-length 0.010000 Hessian is computed using Central difference with step-length 0.100000 Hessian matrix is forced to be a diagonal matrix? [No] Compute effective number of parameters? [Yes] Perform a Monte Carlo error-test? [No] Interpolator [Auto] CPO required diff in log-density [3] Stupid search mode: Status [On] Max iter [1000] Factor [1.05] Numerical integration of hyperparameters: Maximum number of function evaluations [100000] Relative error ....................... [1e-005] Absolute error ....................... [1e-006] To stabilise the numerical optimisation: Minimum value of the -Hessian [-1.#INF] CPO manual calculation[No] Laplace-correction is Disabled. inla_build: check for unused entries in[C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/Model.ini] inla_INLA... Strategy = [DEFAULT] Size is [61206] Chose OpenMP-strategy [HUGE] Chose density-strategy [LOW] Size of graph=[61206] constraints=[2] Found optimal reordering=[amdc] nnz(L)=[330139] and use_global_nodes(user)=[no] List of hyperparameters: theta[0] = [Log precision for county (idd component)] theta[1] = [Log precision for county (spatial component)] theta[2] = [Log precision for county1] Optimise using DEFAULT METHOD file: smtp-taucs.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712, Thread: 0 Variable evaluates to NAN or INF. idx=(16,16). I will try to fix it... file: smtp-taucs.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712, Thread: 0 Variable evaluates to NAN or INF. idx=(16,16). I will try to fix it... file: smtp-taucs.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712, Thread: 1 Variable evaluates to NAN or INF. idx=(15302,15302). I will try to fix it... file: smtp-taucs.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712, Thread: 1 Variable evaluates to NAN or INF. idx=(15302,15302). I will try to fix it... file: smtp-taucs.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712, Thread: 0 Variable evaluates to NAN or INF. idx=(16,16). I will try to fix it... file: smtp-taucs.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712, Thread: 1 Variable evaluates to NAN or INF. idx=(15302,15302). I will try to fix it... file: smtp-taucs.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712, Thread: 0 Variable evaluates to NAN or INF. idx=(16,16). I will try to fix it... file: smtp-taucs.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712, Thread: 2 Variable evaluates to NAN or INF. idx=(30666,30666). I will try to fix it... file: smtp-taucs.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712, Thread: 1 Variable evaluates to NAN or INF. idx=(15302,15302). I will try to fix it... file: smtp-taucs.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712, Thread: 0 Variable evaluates to NAN or INF. idx=(16,16). I will try to fix it... file: smtp-taucs.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712, Thread: 0 Variable evaluates to NAN or INF. idx=(16,16). I will try to fix it... file: smtp-taucs.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712, Thread: 1 Variable evaluates to NAN or INF. idx=(15302,15302). I will try to fix it... file: smtp-taucs.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712, Thread: 1 Variable evaluates to NAN or INF. idx=(15302,15302). I will try to fix it... file: smtp-taucs.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712, Thread: 0 Variable evaluates to NAN or INF. idx=(16,16). I will try to fix it... file: smtp-taucs.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712, Thread: 0 Variable evaluates to NAN or INF. idx=(16,16). I will try to fix it... file: smtp-taucs.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712, Thread: 1 Variable evaluates to NAN or INF. idx=(15302,15302). I will try to fix it... file: smtp-taucs.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712, Thread: 0 Variable evaluates to NAN or INF. idx=(16,16). I will try to fix it... file: smtp-taucs.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712, Thread: 0 Variable evaluates to NAN or INF. idx=(16,16). I will try to fix it... file: smtp-taucs.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712, Thread: 1 Variable evaluates to NAN or INF. idx=(15302,15302). I will try to fix it... file: smtp-taucs.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712, Thread: 0 Variable evaluates to NAN or INF. idx=(16,16). I will try to fix it... file: smtp-taucs.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712, Thread: 0 Variable evaluates to NAN or INF. idx=(16,16). I will try to fix it... file: smtp-taucs.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712, Thread: 1 Variable evaluates to NAN or INF. idx=(15302,15302). I will try to fix it... file: smtp-taucs.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712, Thread: 0 Variable evaluates to NAN or INF. idx=(16,16). I will try to fix it... file: smtp-taucs.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712, Thread: 1 Variable evaluates to NAN or INF. idx=(15302,15302). I will try to fix it... GMRFLib version 3.0-0-snapshot, has recived error no [12] Reason : The Newton-Raphson optimizer did not converge Message : Condition `lambda < 1.0 / lambda_lim' is not TRUE Function : GMRFLib_init_GMRF_approximation_store__intern File : approx-inference.c Line : 2856 RCSId : file: approx-inference.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 GMRFLib version 3.0-0-snapshot, has recived error no [12] Reason : The Newton-Raphson optimizer did not converge Message : Condition `lambda < 1.0 / lambda_lim' is not TRUE Function : GMRFLib_init_GMRF_approximation_store__intern File : approx-inference.c Line : 2856 RCSId : file: approx-inference.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 GMRFLib version 3.0-0-snapshot, has recived error no [12] Reason : The Newton-Raphson optimizer did not converge Message : Condition `lambda < 1.0 / lambda_lim' is not TRUE Function : GMRFLib_init_GMRF_approximation_store__intern File : approx-inference.c Line : 2856 RCSId : file: approx-inference.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200 GMRFLib version 3.0-0-snapshot, has recived error no [12] Reason : The Newton-Raphson optimizer did not converge Message : Condition `lambda < 1.0 / lambda_lim' is not TRUE Function : GMRFLib_init_GMRF_approximation_store__intern File : approx-inference.c Line : 2856 RCSId : file: approx-inference.c hgid: e6dd688ed77e date: Sat Jun 11 22:34:15 2016 +0200

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>. In addition: Warning messages:

1: In writeBin(as.double(as.vector(A)), fp) : NAs introduced by coercion 2: running command '"C:/Users/macha058/Documents/R/R-3.0.2/library/INLA/bin/windows/32bit/inla.exe" -b -s -v "C:/Users/macha058/AppData/Local/Temp/RtmpYLMJa9/file3cf4835b85698/Model.ini"' had status 3

[Håvard Rue]

Hi.

this is saying that something in the latent field goes to +- infinity.
could it be that the graph is not connected and you have islands where
all the observations are zero? 

H

On Wed, 2016-06-15 at 21:45 -0700, Gustavo Machado wrote:
> Optimise using DEFAULT METHOD
>
> file: smtp-taucs.c  hgid: e6dd688ed77e  date: Sat Jun 11
> 22:34:15 2016 +0200
> Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712,
> Thread: 0
> Variable evaluates to NAN or INF. idx=(16,16). I will try to
> fix it...
> file: smtp-taucs.c  hgid: e6dd688ed77e  date: Sat Jun 11
> 22:34:15 2016 +0200
> Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712,
> Thread: 0
> Variable evaluates to NAN or INF. idx=(16,16). I will try to
> fix it...
> file: smtp-taucs.c  hgid: e6dd688ed77e  date: Sat Jun 11
> 22:34:15 2016 +0200
> Function: GMRFLib_build_sparse_matrix_TAUCS(), Line: 712,
> Thread: 1

--

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

[Haakon Bakka]

It looks like the starting values for the hyper-parameters represent a model that does not make sense (the computed precision matrix fails).

If that is the case, you can try fixing the hyperparameters at values you know are reasonable. And, if that works, set the initial values for the optimizer to those reasonable values, through:

control.mode=list(restart=T, theta=c(4,1,-1))

where the theta vector is the internal hyper-parameters.

Regards,

Haakon Bakka

--
You received this message because you are subscribed to the Google Groups "R-inla discussion group" group.
To unsubscribe from this group and stop receiving emails from it, send an email to r-inla-discussion...@googlegroups.com.
To post to this group, send an email to r-inla-disc...@googlegroups.com.

Visit this group at https://groups.google.com/group/r-inla-discussion-group.
For more options, visit https://groups.google.com/d/optout.