SpatialExtremes (version 1.8-0)

print.copula: Printing objects of class copula''

Description

A method for printing object of class ``copula''.

Usage

## S3 method for class 'copula':
print(x, digits = max(3, getOption("digits") - 3), ...)

Arguments

x

An object of class ``copula. Most often, x is the output of the fitcopula function.

digits

The number of digits to be printed.

...

Other options to be passed to the print function.

Value

• Print several information on screen.

Examples

n.site <- 30
n.obs <- 50

coord <- matrix(runif(2 * n.site, -10, 10), ncol = 2)
colnames(coord) <- c("lon", "lat")

## Generate data from a Gaussian copula model

## 1. Generate of Gaussian random field
data <- rgp(n.obs, coord, "powexp", sill = 1, range = 4, smooth = 1.2)

## 2. Transform the margins to U(0,1)
data <- pnorm(data)

## 3. Transform the margins to GEV
locs <- -5 + coord[,"lon"] / 10
scales <- 10 + coord[,"lat"] / 2
shapes <- rep(0.2, n.site)

for (i in 1:n.site)
  data[,i] <- qgev(data[,i], locs[i], scales[i], shapes[i])

## Fit a Gaussian copula model

## 1. Define trend surfaces
loc.form <- y ~ lon
scale.form <- y ~ lat
shape.form <- y ~ 1

## 2. Fit
M0 <- fitcopula(data, coord, "gaussian", "powexp", loc.form, scale.form,
                shape.form, sill = 1)
M0