OZrain(VLMC)
OZrain()所属R语言包：VLMC
Daily Rainfall in Melbourne, Australia, 1981-1990
日雨量在澳大利亚墨尔本，1981年至1990年
译者：生物统计家园网 机器人LoveR
描述----------Description----------
Amount of daily rainfall in Melbourne, Australia, 1981-1990, measured in millimeters. The amounts are integers with many zeros and three days of more than 500mm rain.
日降雨量在澳大利亚墨尔本，1981-1990年，金额单位为毫米。的金额都是整数，有许多零三天超过500毫米雨。
用法----------Usage----------
data(OZrain)
格式----------Format----------
A time-series of length 3653 with the amount of daily rainfall in mm. Because of the two leap years 1984 and '88, we have constructed it with ts(*, start=1981, frequency=365.25, end = 1981+ (3653 - 1)/365.25).
一个时间序列的长度3653毫米的日雨量。由于两个飞跃，1984年和88年的，我们已经建立了与ts(*, start=1981, frequency=365.25, end = 1981+ (3653 - 1)/365.25)。
注意----------Note----------
There must be one extra observation since for the ten years with two leap years, there are only 3652 days. In 61 out of 100 days, there's no rain.
必须有一个额外的观测以来的10年，两个闰年，有3652天的。在61个100天，有没有下雨。
源----------Source----------
"rainfall.dat" in Rob J. Hyndman's Time Series Data Library, http://www-personal.buseco.monash.edu.au/~hyndman/TSDL/
“”rainfall.dat罗伯J.海德门的时间序列数据库，http://www-personal.buseco.monash.edu.au/海德门/ TSDL /
originally, Australian Bureau of Meteorology, http://www.abs.gov.au.
最初，澳大利亚气象http://www.abs.gov.au局。
实例----------Examples----------
data(OZrain)
(n <- length(OZrain)) ## should be 1 more than[＃应该是1年多]
ISOdate(1990,12,31) - ISOdate(1981, 1,1)## but it's 2 ..[＃但它2。]
has.rain <- OZrain > 0
summary(OZrain[has.rain])# Median = 18, Q3 = 50[中位数= 18，Q3 = 50]
table(rain01 <- as.integer(has.rain))
table(rain4c <- cut(OZrain, c(-.1, 0.5, 18.5, 50.1, 1000)))
AIC(v1 <- vlmc(rain01))# cutoff = 1.92[截止= 1.92]
AIC(v00 <- vlmc(rain01, cut = 1.4))
AIC(v0 <- vlmc(rain01, cut = 1.5))
hist(OZrain)
hist(OZrain, breaks = c(0,1,5,10,50,1000), xlim = c(0,100))
plot(OZrain, main = "Rainfall 1981-1990 in Melbourne")
## work around bug in plot.ts() for R <= 1.3.1 :[＃解决错误plot.ts（），R <= 1.3.1：]
newer1.4 <- {v <- R.version; v$major > 1 || as.numeric(v$minor) >= 4.2}
if(newer1.4) {
plot(OZrain, log="y", main = "Non-0 Rainfall [LOG scale]")
} else
plot.ts(time(OZrain)[has.rain],
OZrain[has.rain], log="y", main = "Non-0 Rainfall [LOG scale]")
lOZ <- lowess(log10(OZrain[has.rain]), f= .05)
lines(time(OZrain)[has.rain], 10^lOZ$y, col = 2, lwd = 2)
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注：
注1：为了方便大家学习，本文档为生物统计家园网机器人LoveR翻译而成，仅供个人R语言学习参考使用，生物统计家园保留版权。
注2：由于是机器人自动翻译，难免有不准确之处，使用时仔细对照中、英文内容进行反复理解，可以帮助R语言的学习。
注3：如遇到不准确之处，请在本贴的后面进行回帖，我们会逐渐进行修订。
|
发表于 2012-10-1 16:24:05
