Last data update: 2014.03.03

R: Phosphorous Concentrations in Speed River, Monthly
GuelphPR Documentation

Phosphorous Concentrations in Speed River, Monthly

Description

Monthly time series of phosphorous (P) concentrations in mg/l, Speed River, Guelph, Ontario, 1972.1-1977.1."

Usage

data(GuelphP)

Format

The format is: Time-Series [1:72] from 1972 to 1978: 0.47 0.51 0.35 0.19 0.33 NA 0.365 0.65 0.825 1 ... - attr(*, "title")= chr "Phosphorous Data,Speed River,Guelph,1972.1-1977.1"

Details

See Hipel and McLeod (2005, pp.707-709 and p.986). Hipel and McLeod used intervention analysis to estimate the effect of a phosphorous removal scheme. Intervention analysis was also used to estimate the missing values.

Source

Hipel, K.W. and McLeod, A.I., (2005). Time Series Modelling of Water Resources and Environmental Systems. Electronic reprint of our book orginally published in 1994. http://www.stats.uwo.ca/faculty/aim/1994Book/.

Examples

data(GuelphP)
plot(GuelphP)
#replace missing values in the series with estimates obtained from Hipel and McLeod
# (2005, p. 986)
missingEst<-c(0.1524, 0.2144, 0.3064, 0.1342)
GuelphP2<-GuelphP
GuelphP2[is.na(GuelphP)]<-missingEst

Results


R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(Kendall)
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/Kendall/GuelphP.Rd_%03d_medium.png", width=480, height=480)
> ### Name: GuelphP
> ### Title: Phosphorous Concentrations in Speed River, Monthly
> ### Aliases: GuelphP
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(GuelphP)
> plot(GuelphP)
> #replace missing values in the series with estimates obtained from Hipel and McLeod
> # (2005, p. 986)
> missingEst<-c(0.1524, 0.2144, 0.3064, 0.1342)
> GuelphP2<-GuelphP
> GuelphP2[is.na(GuelphP)]<-missingEst
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>