R: Phosphorous Concentrations in Speed River, Monthly
GuelphP
R 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
>