Changes between Version 14 and Version 15 of udg/ecoms/RPackage/examples/continentalSelection

Timestamp:

Sep 2, 2014 6:08:45 PM (7 years ago)

Author:

juaco

Comment:

--

udg/ecoms/RPackage/examples/continentalSelection

@@ -81 +81 @@
 }}}
 
-[[Image(image-20140902-172548.png)]]
+
+[[Image(image-20140902-172529.png)]]
 
 
@@ -88 +89 @@
 {{{
 #!text/R
-ref <- apply(obs.regridded$Data, MARGIN = c(3,2), mean, na.rm = TRUE)
+> ref <- apply(obs.regridded$Data, MARGIN = c(3,2), mean, na.rm = TRUE)
 }}}
 
+The following lines of code compute the bias of each member w.r.t. the reference and plot them:
 
+{{{
+#!text/R
+# Now we compute the difference agains each of the multimember spatial means:
+> require(fields) 
+> n.members <- dim(ex2$Data)[1]
+> par(mfrow = c(1,2))
+> for (i in 1:n.members) {
++       member <- apply(ex2$Data[i, , , ], MARGIN = c(3,2), mean, na.rm = TRUE)
++       bias <- member - ref      
++       image.plot(ex2$xyCoords$x, ex2$xyCoords$y, bias, xlab = "lon", ylab = "lat", asp = 1)
++       title(paste("Bias member", i))
++       world(add = TRUE)
++ }
+> par(mfrow = c(1,1)) # To reset the graphical window
+}}}
+
+[[Image(image-20140902-180232.png)]]
+
+
+