Is it possible to create a graph file for INLA containing region names

[James Cranston]

I've been trying to link a column in my data containing these regions as a factor, to an INLA model with a graph describing their spatial arrangement.

#something like this f(regions, model = "bym2", graph = gb_regions)

Is there a way to get the region names into the graph file?

Spdep ignores region.id in class nb  when using spdep::nb2INLA?

[INLA help]

You have to map names to integers 1:n manually…

Get Outlook for iOS

---

From: r-inla-disc...@googlegroups.com <r-inla-disc...@googlegroups.com> on behalf of James Cranston <james.c...@ons.gov.uk>
Sent: Tuesday, August 9, 2022 11:42:20 AM
To: R-inla discussion group <r-inla-disc...@googlegroups.com>
Subject: [r-inla] Is it possible to create a graph file for INLA containing region names

 

--
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 view this discussion on the web, visit https://groups.google.com/d/msgid/r-inla-discussion-group/3e134dcc-a0dd-4246-b8c5-adc4b376a8dbn%40googlegroups.com.

[James Cranston]

Thanks, I've created a graph like this manually, I've included a small example here, should this run or am I missing something else?

graph

$n 

3

$nnbs

1 2 1

$nbs

nbs$north

south

nbs$south

north west

nbs$west

south

$cc

$cc$id

1 1 1

$cc$n 

1 

$cc$nodes

north south west

A balanced dataset with all regions presents

data.frame(

y = 1:30, 

x = 1:10,

region =  c("north", "south", "west")

)

y ~ x +

f(region, model = "bym2", graph = graph)

[INLA help]

You need to use 1,2,3 as region , not north, south,west 

Get Outlook for iOS

---

From: r-inla-disc...@googlegroups.com <r-inla-disc...@googlegroups.com> on behalf of James Cranston <james.c...@ons.gov.uk>

Sent: Wednesday, August 10, 2022 4:15:33 PM

To: R-inla discussion group <r-inla-disc...@googlegroups.com>

Subject: Re: [r-inla] Is it possible to create a graph file for INLA containing region names

 

To view this discussion on the web, visit https://groups.google.com/d/msgid/r-inla-discussion-group/3454acd8-88ab-4c6a-a205-c7498efe725fn%40googlegroups.com.

[James Cranston]

So the answer to the original question is no, you can't have a graph file with region names in it? Sorry just checking I understand.

[INLA help]

Yes.  No names, only integers 1:n 

Get Outlook for iOS

---

From: r-inla-disc...@googlegroups.com <r-inla-disc...@googlegroups.com> on behalf of James Cranston <james.c...@ons.gov.uk>

Sent: Friday, August 12, 2022 9:56:15 AM

To: R-inla discussion group <r-inla-disc...@googlegroups.com>

To view this discussion on the web, visit https://groups.google.com/d/msgid/r-inla-discussion-group/01b15b23-9311-49fa-8751-89633cc6331dn%40googlegroups.com.