spatial_data.R
In RTMB: 'R' Bindings for 'TMB'

n <-
100
y <-
c(2, 5, 0, 2, 4, 2, 8, 6, 8, 5, 0, 2, 2, 4, 0, 3, 9, 8, 4, 4, 
1, 3, 3, 3, 1, 2, 4, 5, 6, 2, 3, 3, 4, 4, 1, 5, 6, 4, 1, 1, 2, 
1, 3, 11, 2, 4, 13, 4, 4, 10, 3, 6, 1, 1, 0, 8, 4, 4, 5, 3, 1, 
1, 1, 1, 4, 3, 1, 4, 4, 4, 0, 0, 3, 7, 5, 4, 2, 2, 5, 2, 3, 1, 
6, 6, 4, 3, 0, 5, 4, 3, 3, 6, 4, 1, 5, 0, 5, 2, 5, 2)
X <-
structure(c(0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 
0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 
0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 
0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 
0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 
0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 
0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 
0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, -0.0784305278075794, 
0.0431035142746222, -0.104877716200653, -0.0653639052797231, 
-0.00207364333418504, 0.14174297334029, 0.149126846472008, 0.0248120288058295, 
0.023958580335832, 0.0441664660258962, 0.084214520441791, 0.0781646215611708, 
-0.0140286791875, 0.108077655265921, 0.0693553828017935, -0.0249721650120785, 
-0.0854416251152327, -0.104265396974621, -0.11326280842651, -0.148281223149694, 
-0.127827104476326, -0.0413755411977594, -0.0977355708057351, 
0.116678837176254, 0.145064258925477, 0.0846687181295023, -0.0373373593063127, 
0.00561201526140945, -0.0561050454241364, 0.0518086254198018, 
0.120955495574931, 0.0979538109685988, -0.116962483427178, -0.0423237799587788, 
0.0890394355516125, 0.135984843443705, -0.144150565765741, -0.0940597800527178, 
0.124292453948915, 0.105890162989044, 0.0981768409483373, -0.141881203948403, 
0.12264808661279, 0.0461216263368792, 0.0248917491993786, -0.0259374726724203, 
-0.0119806957542054, -0.15595680733957, -0.136177087359117, -0.06191070786313, 
-0.0772899201946506, 0.0306790923750974, 0.0520776107478291, 
-0.142080712145883, -0.0839762733374298, -0.0964707543100864, 
0.125455320902941, -0.133587684630007, 0.0885559762013782, 0.133777728683172, 
-0.136156120950409, -0.13656354375424, 0.0888431758998778, 0.0417775756918825, 
-0.150056308164483, -0.147051661221829, 0.0402184600367766, 0.0535249328171864, 
0.0356723244992564, -0.0302477301665515, 0.141994732988922, 0.0593378578340184, 
0.00262609269040793, 0.105238187499615, -0.159267782491596, 0.102440857341435, 
-0.143790649635233, 0.125100749161356, 0.102454647129236, 0.0576697080864678, 
-0.134575349069181, 0.0667664474308766, -0.153575801186939, 0.138065395831201, 
0.137650408221086, -0.00160093829184136, -0.15425002103278, 0.14877912862145, 
-0.158567579796551, 0.0202554752585219, 0.0721763435862204, 0.0814037847299997, 
-0.0209023918240075, -0.111590523129262, -0.125927907569247, 
0.11884670667127, 0.0962652563280377, 0.00503878019688148, 0.0192770700134143, 
-0.0982308285460897), .Dim = c(100L, 2L), .Dimnames = list(NULL, 
    c("V1", "V2")))
Z <-
structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 
2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 
3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 5L, 
5L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 7L, 
7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 
8L, 8L, 8L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 10L, 10L, 
10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 1L, 2L, 3L, 4L, 5L, 6L, 
7L, 8L, 9L, 10L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 1L, 
2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 1L, 2L, 3L, 4L, 5L, 6L, 
7L, 8L, 9L, 10L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 1L, 
2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 1L, 2L, 3L, 4L, 5L, 6L, 
7L, 8L, 9L, 10L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 1L, 
2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 1L, 2L, 3L, 4L, 5L, 6L, 
7L, 8L, 9L, 10L), .Dim = c(100L, 2L), .Dimnames = list(NULL, 
    c("V1", "V2")))
tab.expected <-
structure(c(11.4081202944103, -0.139646809704954, 0.0054467612225716, 
-0.803098852187008, 1.00546162179489, 0.447938715274447, 1.16934361844857, 
0.676632782685668, 0.489633064450651, 0.183943489005993, 0.492307255066955, 
0.0823954101484441), dim = c(6L, 2L), dimnames = list(c("b", 
"b", "loga", "logsigma", "a", "sigma"), c("Estimate", "Std. Error"
)))