Description Usage Format Value Source References Examples
Australian daily minimum temperature climate data for 8 different stations are provided. The data is taken from Australian Government Bureau of Meteorology.
1 |
An object of class data.table
or data.frame
.
Australian daily minimum temperature climate data for 8 different stations:
Sydney |
Sydney (Observatory Hill), taken from 1859 to 2012 |
Melbourne |
Melbourne (Regional Office), taken from 1855 to 2012 |
Boulia |
Boulia Airport, taken from 1888 to 2012 |
Cape_Otway |
Cape Otway Lighthouse, taken from 1864 to 2012 |
Gayndah |
Gayndah Post Office, taken from 1893 to 2009 |
Gunnedah |
Gunnedah Pool, taken from 1876 to 2011 |
Hobart |
Hobart (Ellerslie Road), taken from 1882 to 2012 |
Robe |
Robe Comparison, taken from 1884 to 2012 |
The daily observations are available from http://www.bom.gov.au/climate/data. Copyright Commonwealth of Australia 2010, Bureau of Meteorology. Definitions adapted from http://www.bom.gov.au/climate/dwo/IDCJDW0000.shtml
Australian Government Bureau of Meteorology
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 | library(fda)
library(reshape2)
fun_data_S = Australian_Temp$Sydney
D = 21
basis = create.fourier.basis(rangeval = c(0, 1), nbasis = D)
nas = which(is.na(fun_data_S$Days.of.accumulation.of.minimum.temperature))
fun_data_S = fun_data_S[-nas, ]
yy = unique(fun_data_S$Year)
mat.S = matrix(0, D, length(yy))
for (i in 1:length(yy)){
aa = subset(fun_data_S, Year==yy[i])
cc = aa$Minimum.temperature..Degree.C.
a = which(is.na(cc))
if (length(a)>0){
cc = cc[-which(is.na(cc))]
}else{
cc = cc
}
f_Obs = Data2fd(argvals=seq(0, 1, length = length(cc)) , cc, basisobj = basis)
mat.S[, i] = f_Obs$coefs
}
fdata = fd(mat.S, basis)
# note that the last year, has data only up to 6 months
# therefore we remove it
fdata = fdata[-length(yy)]
plot(fdata)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.