plotfields          package:ProbForecastGOP          R Documentation

_P_l_o_t _o_f _w_e_a_t_h_e_r _f_i_e_l_d_s

_D_e_s_c_r_i_p_t_i_o_n:

     Plot weather fields.

_U_s_a_g_e:

     plotfields(field, x.lim, y.lim, title)

_A_r_g_u_m_e_n_t_s:

   field: numeric square matrix containing the values of the weather
          field on a grid.

   x.lim: numeric vector giving the smallest and the largest longitude
          to be displayed.

   y.lim: numeric vector giving the smallest and the largest latitude
          to be displayed.

   title: character string with the title for the plot.

_V_a_l_u_e:

     The function returns a graphical display of the weather field on a
     region delimited by the lower and upper bound for the longitude
     and the latitude.

_N_o_t_e:

     This function uses the package 'fields'.

_A_u_t_h_o_r(_s):

     Gel, Y., Raftery, A. E., Gneiting, T., Berrocal, V. J.
     veronica@stat.washington.edu.

_S_e_e _A_l_s_o:

     The package 'fields' for display of spatial data, and 'US' and
     'world' for a map of the US and the world.

_E_x_a_m_p_l_e_s:

     ## Loading data
     library(fields)
     library(RandomFields)
     data(slp)
     data(gridlong)
     data(gridlat)
     data(forecast.grid)
     day <- slp$date.obs
     id <- slp$id.stat
     coord1 <- slp$lon.stat
     coord2 <- slp$lat.stat
     obs <- slp$obs
     forecast <- slp$forecast
     coord1.grid <- gridlong$gridded.long
     coord2.grid <- gridlat$gridded.lat
     forecast.grid <- forecast.grid$gridded.forecast

     ## Computing the empirical variogram
     variogram <- Emp.variog(day,obs,forecast,id,coord1,coord2,cut.points=NULL,
     max.dist=NULL,nbins=NULL)

     ## Estimating parameters
     param.est <- Variog.fit(variogram,"exponential",max.dist.fit=NULL,
     init.val=NULL,fix.nugget=FALSE)

     ## Simulating realizations of the weather random field
     simul <- Field.sim(obs, forecast, coord1.grid, coord2.grid, forecast.grid,
     variog.model="exponential", param.est=c(param.est$nugget,param.est$variance,
     param.est$range), n.sim=4, n.displ=0, qt.displ=c(10,50,90))

     ##Plotting one of the simulated weather random fields
     par(mfrow=c(1,1))
     plotfields(simul$sim.fields[,,1],x.lim=c(min(coord1.grid),max(coord1.grid)),
     y.lim=c(min(coord2.grid),max(coord2.grid)),title="Simulated weather field")

     ## Plotting one of the percentiles of the weather field
     par(mfrow=c(1,1))
     plotfields(simul$pct.fields[,,1],x.lim=c(min(coord1.grid),max(coord1.grid)),
     y.lim=c(min(coord2.grid),max(coord2.grid)),title="10th percentile")

