Description Major functions Toy example functions References Examples
The package cope computes and plots CoPE sets defined in Sommerfeld, Sain and Schwartzman (2015) for 2D functions. CoPE sets for a real-valued target function μ(s) on a two-dimensional domain are designed to bound the excursion set {μ(s) >= c} from above and below with a predefined probability. The target function can be a parameter in spatially indexed linear regression Y(s)=X*b(s)+ ε(s), where s is the spatial location, X is a known n by p design matrix, ε(s) is an error field and Y(s) is the observed data.
ComputeCope
Computes CoPE sets for 2D data.
PlotCope
Plots CoPE sets.
ToySignal
Generates a toy signal.
ToyNoise1
, ToyNoise2
,
ToyNoise3
Generates realizations of toy noise fields.
M. Sommerfeld, S. Sain and A. Schwartzman. Confidence regions for excursion sets in asymptotically Gaussian random fields, with an application to climate. Preprint, 2015.
1 2 3 4 5 6 7 8 9 10 11 12 | # An example using the ToySignal and the Toy Noise1 of this package.
# Set sample size.
n = 30
# Generate n realizations of the toy noise field.
Data = ToyNoise1(n = n)
# Add the toy signal to the noise.
Data$z = Data$z + rep(ToySignal()$z, n)
# Compute the CoPE sets.
CopeSet = ComputeCope(Data,level=4/3, mu=ToySignal()$z)
# Plot the result.
PlotCope(CopeSet)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.