ensemble_temperature: Observed or reconstructed temperatures (e.g., from CRU or...

Description Format Examples

Description

A dataset containing linear trend from an ensemble of observed or reconstructed temperatures

Format

A n by N matrix, where n=2107 is the number of grid cells over globe, and N=394 is the number of ensemble memembers

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
# estimate linear trend
## Not run: 
data("ensemble_temperature")
data("GCM_runs")
y <- ensemble_temperature
x <- GCM_runs
N <- dim(y)[2]
avg <- function(x) {return(apply(x, 1, mean))}
xbar <- sapply(x, avg)
Y <- c(y)
X <- kronecker(diag(rep(1,N),N,1), xbar)

# OLS estimates
fit <- lm(Y~X)
summary(fit)

## End(Not run)

# plot the first 10 ensemble memebers
boxplot(ensemble_temperature[ ,1:10], xlab="ensemble member", 
        ylab="temperature increase", cex=0.5, cex.axis=0.6, cex.lab=0.8, 
        pars=list(boxwex=0.5), medcol="red", medlwd=1)

mapn/DAbayes documentation built on May 21, 2017, 5:47 p.m.