# Taylor Diagram with the plotrix package

Hello,

I want to compare 6 climate models to my observations via a Taylor diagram in R using the plotrix package.

library('plotrix')

ref = c(mean(Y1\$Z1),sd(Y1\$Z1))
model1 = c(mean(Y1\$X1),sd(Y1\$X1))
model2 = c(mean(Y1\$X2),sd(Y1\$X2))
model3 = c(mean(Y1\$X3),sd(Y1\$X3))
model4 = c(mean(Y1\$X4),sd(Y1\$X4))
model5 = c(mean(Y1\$X5),sd(Y1\$X5))
model6 = c(mean(Y1\$X6),sd(Y1\$X6))

crs = c('red', 'orange', 'blue', 'red', 'orange', 'blue')

taylor.diagram(ref, model1, pch = 1, col = crs[1], normalize = F, pos.cor=TRUE)
taylor.diagram(ref, model2, pch = 2, col = crs[2], normalize = F, add = T)
taylor.diagram(ref, model3, pch = 0, col = crs[3], normalize = F, add = T)
taylor.diagram(ref, model4, pch = 16, col = crs[4], normalize = F, add = T)
taylor.diagram(ref, model5, pch = 17, col = crs[5], normalize = F, add = T)
taylor.diagram(ref, model6, pch = 15, col = crs[6], normalize = F, add = T)
[enter image description here][1]

legend('topright', legend = c('Mod1', 'Mod2', 'Mod3', 'Mod4', 'Mod5', 'Mod6'),
pch = c(1, 2, 0, 16, 17, 15,18), col = crs, bty = 'n')

and this is the plot I obtain. What did I do wrong?

This topic was automatically closed 21 days after the last reply. New replies are no longer allowed.