2015-05-15 2 views
1

Раньше я использовал INLA с помощью SpatialPolygonsDataFrame, но я никогда не использовал его с растром. Я много часов прочесываю в интернете, читая виньетки и все, что могу найти, но все еще не могу понять, как использовать INLA с растром. Ниже всего лишь простой пример того, что мои данные, как:Использование INLA с растром в R

library("raster") 
library("INLA") 

lat = c(10, 10, 10, 10, 14, 14, 14, 14, 10, 10, 10, 10, 14, 14, 14, 14) 
lon = c(20 ,24, 28, 32, 20, 24, 28, 32, 20, 24, 28, 32, 20, 24, 28, 32) 
response = c(17.3, 17.4, 17.6, 17.9, 17.1, 17.0, 17.4, 17.5, 17.4, 17.4, 
17.7, 17.8, 17.2, 17.1, 17.7, 17.9) 
explan1 = c(31, 35, 33, 36, 32, 37, 36, 40, 32, 34, 33, 37, 35, 35, 39, 38) 
explan2 = c(112, 116, 111, 114, 115, 117, 117, 112, 110, 114, 113, 117, 112, 
113, 115, 116) 

data.df = data.frame(lat = lat, lon = lon, response = response, explan1 = 
explan1, explan2 = explan2) 

myGrid = raster(ncol = 5, nrow = 5) 

Так что я брошусь перед всеми вами, умоляя о пощаде и помощи таким образом. Заранее благодарю вас за то, что вышли из моих страданий!

+0

Что вы хотите сделать? И, пожалуйста, покажите нам пример кода, который работает. 'a = c (8 случайных значений)' не является правильным кодом 'R' – RobertH

+0

Я улучшил свой пример. Я пытаюсь сделать пространственно-временную модель, которая будет пытаться моделировать будущее событие с помощью INLA. Без использования растровых данных я бы просто создал SpatialPolygonsDataFrame, а затем использовал poly2nb и построил таблицу пространства-времени, прежде чем выполнять спецификации формул и запустить модель INLA. Любая помощь в том, как это сделать с использованием растровых данных, будет с благодарностью оценена. – user2221184

ответ

1

Для установки данных используйте свои данные.df (и вашу модель времени пространства).

Для прогнозирования/проецирования на растре используйте inla.mesh.project. Вы можете проецировать результат $ summary.linear.predictor, result $ summary.fitted.values ​​или result $ summary.random, в зависимости от вашего кода и того, что вы хотите получить.

Если вам нужны образцы вместо резюме, используйте inla.posterior.sample и проецируйте правильную часть (в зависимости от вашей модели).

Пространственно-временные модели сложны, а рабочие примеры очень длинные. См. Учебное пособие spde http://www.r-inla.org/examples/tutorials/spde-tutorial

Смежные вопросы