bomsoi                 package:DAAG                 R Documentation

_S_o_u_t_h_e_r_n _O_s_c_i_l_l_a_t_i_o_n _I_n_d_e_x _D_a_t_a

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

     The Southern Oscillation Index (SOI) is the difference in
     barometric   pressure at sea level between Tahiti and Darwin. 
     Annual SOI and Australian rainfall data, for the years 1900-2001,
     are given. Australia's annual mean rainfall is an area-weighted
     average of the total  annual precipitation at approximately 370
     rainfall stations  around the country.

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

     bomsoi

_F_o_r_m_a_t:

     This data frame contains the following columns:

     _Y_e_a_r a numeric vector

     _J_a_n average January SOI values for each year

     _F_e_b average February SOI values for each year

     _M_a_r average March SOI values for each year

     _A_p_r average April SOI values for each year

     _M_a_y average May SOI values for each year

     _J_u_n average June SOI values for each year

     _J_u_l average July SOI values for each year

     _A_u_g average August SOI values for each year

     _S_e_p average September SOI values for each year

     _O_c_t average October SOI values for each year

     _N_o_v average November SOI values for each year

     _D_e_c average December SOI values for each year

     _S_O_I a numeric vector consisting of average annual SOI values

     _a_v_r_a_i_n a numeric vector consisting of a weighted average annual
          rainfall at a large number of Australian sites

_S_o_u_r_c_e:

     Australian Bureau of Meteorology web pages:

     http://www.bom.gov.au/climate/change/rain02.txt and
     http://www.bom.gov.au/climate/current/soihtm1.shtml

_R_e_f_e_r_e_n_c_e_s:

     Nicholls, N., Lavery, B., Frederiksen, C. and Drosdowsky, W. 1996.
     Recent apparent changes in relationships between the El Nino -
     southern oscillation and Australian rainfall and temperature.
     Geophysical Research Letters 23: 3357-3360.

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

      
     plot(ts(bomsoi[, 15:14], start=1900),
          panel=function(y,...)panel.smooth(1900:2001, y,...))

     # Check for skewness by comparing the normal probability plots for 
     # different a, e.g.
     par(mfrow = c(2,3))
     for (a in c(50, 100, 150, 200, 250, 300))
     qqnorm(log(bomsoi[, "avrain"] - a))
       # a = 250 leads to a nearly linear plot
     par(mfrow = c(1,1))

     plot(bomsoi$SOI, log(bomsoi$avrain - 250), xlab = "SOI",
          ylab = "log(avrain = 250)")
     lines(lowess(bomsoi$SOI)$y, lowess(log(bomsoi$avrain - 250))$y, lwd=2)
       # NB: separate lowess fits against time
     lines(lowess(bomsoi$SOI, log(bomsoi$avrain - 250)))

     detsoi <- data.frame(
       detSOI = bomsoi[, "SOI"] - lowess(bomsoi[, "SOI"])$y,
       detrain = log(bomsoi$avrain - 250) - lowess(log(bomsoi$avrain - 250))$y)
     row.names(detsoi) <- paste(1900:2001)

     par(mfrow = c(1,2))  
     plot(log(avrain-250) ~ SOI, data = bomsoi, ylab = 
      "log(Average rainfall - 250)")
     lines(lowess(bomsoi$SOI, log(bomsoi$avrain-250)))
     plot(detrain ~ detSOI, data = detsoi,
       xlab="Detrended SOI", ylab = "Detrended log(Rainfall-250)")
     lines(lowess(detsoi$detrain ~ detsoi$detSOI)) 
     par(mfrow = c(1,1))

