fieldsim: Simulate manifold indexed Gaussian field by the Fieldsim...

Description Usage Arguments Value Author(s) References See Also Examples

View source: R/fieldsim.R

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.

See Also

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
Loading required package: RColorBrewer
Warning messages:
1: In rgl.init(initValue, onlyNULL) : RGL: unable to open X11 display
2: 'rgl_init' failed, running with rgl.useNULL = TRUE 
3: .onUnload failed in unloadNamespace() for 'rgl', details:
  call: fun(...)
  error: object 'rgl_quit' not found 
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.