Simulate manifold indexed Gaussian field by the Fieldsim method

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

fieldsim(process,Ne,nbNeighbor)

Arguments

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

# 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.

Results

R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(FieldSim)
Loading required package: rgl
Loading required package: RColorBrewer
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/FieldSim/fieldsim.Rd_%03d_medium.png", width=480, height=480)
> ### Name: fieldsim
> ### Title: Simulate manifold indexed Gaussian field by the Fieldsim method
> ### Aliases: fieldsim
> 
> ### ** Examples
> 
> # Load FieldSim library
> library(FieldSim)
> 
> # Fractional Brownian field on [0,1]^2
> plane.fBm<-setProcess("fBm-plane",0.7)
> str(plane.fBm)
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(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.
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>