Имитация случайного поля с помощью пакета RandomFields в r - PullRequest
Новая
       17

Имитация случайного поля с помощью пакета RandomFields в r

0 голосов
/ 04 декабря 2018

У меня положительно определенная матрица B, и мне нужно сгенерировать случайное поле X, используя матрицу B в качестве ковариационной матрицы поля.У меня есть этот код для генерации 10000 реализаций поля 

library(RandomFields)
RFoptions(seed=NA) 
library(Matrix)
A<-matrix(c(1,2,7,2,5,6,7,6,9),3,3)
B<-nearPD(A,keepDiag =T,maxit = 200)
B<-matrix(B$mat,3,3)
model<-RMfixcov(B)
X<-RFsimulate(model,x=seq(0, 1, length=100),y=seq(0, 1, length=100))

Я получил ошибку:

Ошибка в rfInit (model = list ("Simulate", setseed = eval (parse (text = "quote (set.seed (seed = seed))")),: 'RFsimulate': количество данных (9), не кратное количеству местоположений (10000 x 0) ^ 2

Что я делаю не так?

...