A dataset containing linear trend from an ensemble of observed or reconstructed temperatures
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
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)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.