fieldsim: Simulate manifold indexed Gaussian field by the Fieldsim... In FieldSim: Random Fields (and Bridges) Simulations

Description

The function `fieldsim` yields simulation of sample path of a manifold indexed Gaussian field (or bridge) following the procedure described in Brouste et al. (2007, 2010, 2014).

Usage

 `1` ```fieldsim(process,Ne,nbNeighbor) ```

Arguments

 `process` an S4 object process `Ne` a positive integer corresponding to the number of points to simulate with the accurate simulation step `nbNeighbor` a positive integer (between 2 and 32) corresponding to the number of neighbors to use in the second refined step of the algorithm.

Value

The function returns in the slot `values` of the object process the values of the process on the manifold atlas

Author(s)

Alexandre Brouste (http://perso.univ-lemans.fr/~abrouste/) and Sophie Lambert-Lacroix (http://membres-timc.imag.fr/Sophie.Lambert/).

References

A. Brouste, J. Istas and S. Lambert-Lacroix (2007). On Gaussian random fields simulations. A. Brouste, J. Istas and S. Lambert-Lacroix (2010) On simulation of manifold indexed fractional Gaussian fields. A. Brouste, J. Istas and S. Lambert-Lacroix (2014) Fractional Gaussian bridges with the package `FieldSim`.

`process-class`, `setProcess`.

Examples

 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21``` ```# Load FieldSim library library(FieldSim) # Fractional Brownian field on [0,1]^2 plane.fBm<-setProcess("fBm-plane",0.7) str(plane.fBm) fieldsim(plane.fBm) plot(plane.fBm) # Sphere indexed fractional Brownian field #sphere.fBm<-setProcess("fBm-sphere",0.3) #fieldsim(sphere.fBm) #plot(sphere.fBm) # Bridge associated to the Fractional Brownian field on [0,1]^2 #Gamma<-matrix(c(1,0,0,0,1,1,1,1,1,1/2,1/2,0.5),3,4) #bridge.plane.fBm<-setProcess("bridge-fBm-plane",list(Gamma=Gamma,par=0.9)) #fieldsim(bridge.plane.fBm) #plot(bridge.plane.fBm) # Other examples can be found in the setProcess documentation. ```

Example output

```Loading required package: rgl
Warning messages:
1: In rgl.init(initValue, onlyNULL) : RGL: unable to open X11 display
2: 'rgl_init' failed, running with rgl.useNULL = TRUE
call: fun(...)
Formal class 'process' [package "FieldSim"] with 5 slots
..@ name     : chr "fBm"
..@ values   : num 0
..@ manifold :Formal class 'manifold' [package "FieldSim"] with 5 slots
.. .. ..@ name    : chr "plane"
.. .. ..@ atlas   : num [1:2, 1:289] 0 0 0 1 1 0 1 1 0.5 0.5 ...
.. .. ..@ gridtype: chr "visualization"
.. .. ..@ distance:function (xi, xj)
.. .. ..@ origin  : num [1:2, 1] 0 0
..@ covf     :function (xi, xj)
..@ parameter: num 0.7
```

FieldSim documentation built on May 29, 2017, 2:10 p.m.