inst/doc/Introduction.R

## ------------------------------------------------------------------------
library(systemicrisk)

## ------------------------------------------------------------------------
l <- c(1,2.5,3)
a <- c(0.7,2.7,3.1)

## ------------------------------------------------------------------------
L <- sample_ERE(l,a,p=0.5,lambda=0.25,nsamples=200,thin=20,burnin=10)

## ------------------------------------------------------------------------
L[[1]]
L[[2]]
L[[3]]

## ----fig.width=7,fig.height=4--------------------------------------------
plot(sapply(L,function(x)x[1,2]),type="b")

## ----fig.width=7,fig.height=4--------------------------------------------
acf(sapply(L,function(x)x[1,2]))

## ------------------------------------------------------------------------
ea <- c(1,1,1)
el <- c(1,1,1)

## ------------------------------------------------------------------------
default(L[[1]],ea=ea,el=el)$defaultind
default(L[[2]],ea=ea,el=el)$defaultind
default(L[[3]],ea=ea,el=el)$defaultind

## ------------------------------------------------------------------------
rowMeans(sapply(L, function(Lakt) default(Lakt,ea=ea,el=el)$defaultind))

## ------------------------------------------------------------------------
rowMeans(sapply(L, function(Lakt) default(Lakt,ea=ea,el=el,alpha=0.98,beta=0.98)$defaultind))

Try the systemicrisk package in your browser

Any scripts or data that you put into this service are public.

systemicrisk documentation built on Nov. 17, 2017, 6:13 a.m.