R-inla discussion group

1–30 of 3031

Welcome to this discussion group about r-inla. Please ask your questions here in case you think they will be useful for others, otherwise send them to help@r-inla.org. You are of'course free to comment on questions from others as well.

Best,
H

0 selected

9:53 AM

ar1c with replicate /group / etc

Hi all, I believe a similar question has been asked before, but wasn't answered. Is it possible,

unread,

ar1c with replicate /group / etc

Hi all, I believe a similar question has been asked before, but wasn't answered. Is it possible,

9:53 AM

Andrew, Finn Lindgren4

Mar 23

GMRFLib_2order_approx: rescue NAN/INF values in logl for idx

Hi Finn, Thanks for the quick response. I have just tried running again but with the following

unread,

GMRFLib_2order_approx: rescue NAN/INF values in logl for idx

Hi Finn, Thanks for the quick response. I have just tried running again but with the following

Mar 23

Mik, Helpdesk (Haavard Rue)5

Mar 23

Spatial INLA model computing negative predictive values

Hi Håvard, I have tried using inv.logit() in place of exp() but i am still obtaining predicted values

unread,

Spatial INLA model computing negative predictive values

Hi Håvard, I have tried using inv.logit() in place of exp() but i am still obtaining predicted values

Mar 23

Bas Matsuura, Helpdesk (Haavard Rue)2

Mar 23

How to specify control.family and control.link arguments for multiple distributions

with three likelihoods, then the control.family needs to be control.family=list( list(), list(), list

unread,

How to specify control.family and control.link arguments for multiple distributions

with three likelihoods, then the control.family needs to be control.family=list( list(), list(), list

Mar 23

Mar 22

Precision dominates individual linear.predictor standard deviation

Hello, I have noticed that my model$summary.linear.predictor$sd values are typically two orders of

unread,

Precision dominates individual linear.predictor standard deviation

Hello, I have noticed that my model$summary.linear.predictor$sd values are typically two orders of

Mar 22

TheCorinna1994, Finn Lindgren9

Mar 22

non-stationary precipitation model

Ok Ok thanks. One last question: I have my mesh and now I just imported a file containing the

unread,

non-stationary precipitation model

Ok Ok thanks. One last question: I have my mesh and now I just imported a file containing the

Mar 22

Subramanian Swaminathan

Mar 22

inla.group.cv - post-estimation

Hi prof Rue, This is in connection with ' inla.group.cv'. This has greatly simplified

unread,

inla.group.cv - post-estimation

Hi prof Rue, This is in connection with ' inla.group.cv'. This has greatly simplified

Mar 22

Evan Baker, Finn Lindgren6

Mar 22

Inlabru hierarchical model - incorrect random effect estimates

Hi, I've attached these things (summary, bruverbose = 3 output, and covnergence plots for both

unread,

Inlabru hierarchical model - incorrect random effect estimates

Hi, I've attached these things (summary, bruverbose = 3 output, and covnergence plots for both

Mar 22

Mar 21

interpretation results and minimum of data

Hi all, I would like to understand what is the interpertration of my results. I apply a bym model for

unread,

interpretation results and minimum of data

Hi all, I would like to understand what is the interpertration of my results. I apply a bym model for

Mar 21

rafa.arc...@gmail.com, Finn Lindgren4

Mar 21

Predictive scoring of Point Processes

No, CPO isn't well defined for point process likelihoods. Some node on prediction scores (that

unread,

Predictive scoring of Point Processes

No, CPO isn't well defined for point process likelihoods. Some node on prediction scores (that

Mar 21

Saran Ahluwalia, … Helpdesk (Haavard Rue)5

Mar 21

Error when running in Google's Colaboratory

looks good to me On Mon, 2023-03-20 at 08:46 -0700, Saran Ahluwalia wrote: > Incorporating both of

unread,

Error when running in Google's Colaboratory

looks good to me On Mon, 2023-03-20 at 08:46 -0700, Saran Ahluwalia wrote: > Incorporating both of

Mar 21

Lorenzo Valleggi, … Virgilio Gómez-Rubio4

Mar 20

Variance modeling

Hi, \sigma^2_{i,j}(X) = abs(X)^{2\theta_j}} where i are the observations and j are groups where

unread,

Variance modeling

Hi, \sigma^2_{i,j}(X) = abs(X)^{2\theta_j}} where i are the observations and j are groups where

Mar 20

Saran Ahluwalia

Mar 19

Running R-inla on Colaboratory Yields Error

Hi All: Here is an MVP of the issue I am running into. I am not sure what is the best way to

unread,

Running R-inla on Colaboratory Yields Error

Hi All: Here is an MVP of the issue I am running into. I am not sure what is the best way to

Mar 19

JJ Hubbard, … Finn Lindgren11

Mar 19

How can I calculate the standard deviation of the linear.predictor from the standard deviation of the model components?

sd=, not mean=. The mean should be 0. Finn On 19 Mar 2023, at 04:51, JJ Hubbard <jordanjhubbard@

unread,

How can I calculate the standard deviation of the linear.predictor from the standard deviation of the model components?

sd=, not mean=. The mean should be 0. Finn On 19 Mar 2023, at 04:51, JJ Hubbard <jordanjhubbard@

Mar 19

Prince Michael Amegbor, Helpdesk (Haavard Rue)4

Mar 17

this is old code that I do not know if it will work with inla.mode="experimental" which is

unread,

this is old code that I do not know if it will work with inla.mode="experimental" which is

Mar 17

Natoya j, Helpdesk (Haavard Rue)2

Mar 17

approximate variance of data-level prediction

the easy way would be to include year 2021 in the model, but with data equal to NA In this way you

unread,

approximate variance of data-level prediction

the easy way would be to include year 2021 in the model, but with data equal to NA In this way you

Mar 17

Pooja B S, Helpdesk (Haavard Rue)2

Mar 17

unable to fit the shared multiscale model for the count data using two datas sets(one is of district level and another state level)

This seems straight forward. and its look like its something minor only. btw, you use 'g1'

unread,

unable to fit the shared multiscale model for the count data using two datas sets(one is of district level and another state level)

This seems straight forward. and its look like its something minor only. btw, you use 'g1'

Mar 17

Daniel Johnson, … Johnson, Daniel P9

Mar 15

Error Running Blangiardo Camelleti Piemonte SPDE S-T Example

Havard, Insang, and Elias- This last bit from Havard did the trick on the HPC system. I will load

unread,

Error Running Blangiardo Camelleti Piemonte SPDE S-T Example

Havard, Insang, and Elias- This last bit from Havard did the trick on the HPC system. I will load

Mar 15

Sylvan Benaksas, INLA help3

Mar 15

model doesn't run with larger mesh

Hi Haavard, Following another group post https://groups.google.com/g/r-inla-discussion-group/c/

unread,

model doesn't run with larger mesh

Hi Haavard, Following another group post https://groups.google.com/g/r-inla-discussion-group/c/

Mar 15

aynaz lotfata, Elias T. Krainski2

Mar 15

Hi, The 'Ntrials' argument is for 'binomial' while the equivalent for 'nbinomial

unread,

Hi, The 'Ntrials' argument is for 'binomial' while the equivalent for 'nbinomial

Mar 15

Insang Song, Helpdesk (Haavard Rue)4

Mar 13

Successful inla.cpo results still include failures

Dear Dr. Rue: Thank you for the prompt response. I will try inla.group.cv. Best, Insang Song On

unread,

Successful inla.cpo results still include failures

Dear Dr. Rue: Thank you for the prompt response. I will try inla.group.cv. Best, Insang Song On

Mar 13

JJ Hubbard, … Håvard Rue3

Mar 13

Is it possible to continue training a fitted INLA model?

with more data, then likely the latent dimension change, hence you can only carry over the

unread,

Is it possible to continue training a fitted INLA model?

with more data, then likely the latent dimension change, hence you can only carry over the

Mar 13

Josh Cullen, Finn Lindgren3

Mar 13

Hierarchical Gaussian Process Regression

Hi Finn, Thanks for the quick response! I'll take a look at trying to fit this model via inlabru.

unread,

Hierarchical Gaussian Process Regression

Hi Finn, Thanks for the quick response! I'll take a look at trying to fit this model via inlabru.

Mar 13

Vera Tibot, Finn Lindgren2

Mar 13

Larger values of the posterior standard deviation at the mesh nodes?

Hi, this is due to the conditionally deterministic behaviour inside each triangle. For visualisation,

unread,

Larger values of the posterior standard deviation at the mesh nodes?

Hi, this is due to the conditionally deterministic behaviour inside each triangle. For visualisation,

Mar 13

Moisés Augusto, Helpdesk (Haavard Rue)3

Mar 13

Another distribution

Hi, to be short, I want to incorporate the Bell distribution and use it to fit a spatio-temporal

unread,

Another distribution

Hi, to be short, I want to incorporate the Bell distribution and use it to fit a spatio-temporal

Mar 13

Geir Storvik, Håvard Rue2

Mar 6

Marginal likelihoods for the generic0 model

On Mon, 2023-03-06 at 01:22 -0800, Geir Storvik wrote: > As stated in the documentation of the

unread,

Marginal likelihoods for the generic0 model

On Mon, 2023-03-06 at 01:22 -0800, Geir Storvik wrote: > As stated in the documentation of the

Mar 6

Mar 4

Compare classical and Bayesian linear regression

Hi, I have doubt that whether we can compare the Ordinary linear regression and Bayesian regression

unread,

Compare classical and Bayesian linear regression

Hi, I have doubt that whether we can compare the Ordinary linear regression and Bayesian regression

Mar 4

Geir Storvik, … FredrikLAa3

Mar 3

Reparametrization from x to eta

Hi, In the classical way of estimating models with INLA, in the Gaussian approx of pi(Latent field|

unread,

Reparametrization from x to eta

Hi, In the classical way of estimating models with INLA, in the Gaussian approx of pi(Latent field|

Mar 3

Christos Kokaliaris, Finn Lindgren6

Mar 2

spatially varying effect in a coregionalisation model (inlabru)

Yes, in inlabru, the "copy" feature is nearly only needed when estimating a single scaling

unread,

spatially varying effect in a coregionalisation model (inlabru)

Yes, in inlabru, the "copy" feature is nearly only needed when estimating a single scaling

Mar 2

Feb 28

INLA Keeps crashing

Below is the code I run: ## Without interactions count_total_gv ~ gini + urban + minority +

unread,

INLA Keeps crashing

Below is the code I run: ## Without interactions count_total_gv ~ gini + urban + minority +

Feb 28

Search

Clear search

Close search

Google apps

Main menu