- 积分
- 3601
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-9-8
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 hillside 于 2013-4-19 23:27 编辑
包的简介如下:
Package ‘hydroTSM’
February 15, 2013
Type Package
Title Time series management, analysis and interpolation for hydrological modelling
Version 0.3-6
Date 2011-10-18
Author Mauricio Zambrano-Bigiarini [aut, cre]
下面摘录《日降水数据分析使用手册》(见于附件)片断
Tutorial for Introductory Analysis of Daily Precipitation Data with hydroTSM
Mauricio Zambrano-Bigiarini
Oct 2010, Apr 2011, Aug 2011
1 Installation
Installing hydroTSM:
> install.packages("hydroTSM")
2 Setting Up the Environment
1. Loading the hydroTSM library, which contains data and functions used in
this analysis.
> library(hydroTSM)
2. Loading daily precipitation data at the station San Martino di Castrozza,
Trento Province, Italy, with data from 01/Jan/1921 to 31/Dec/1990.
> data(SanMartinoPPts)
3. Selecting only a 6-years time-slice for the analysis
> x <- window(SanMartinoPPts, start=as.Date("1985-01-01"))
4. Monthly values of precipitation
> ( m <- daily2monthly(x, FUN=sum) )
1985-01-01 1985-02-01 1985-03-01 1985-04-01 1985-05-01 1985-06-01 1985-07-01
141.2 7.0 140.6 72.0 175.6 131.4 85.4
1985-08-01 1985-09-01 1985-10-01 1985-11-01 1985-12-01 1986-01-01 1986-02-01
159.4 27.2 58.4 101.8 54.8 75.8 131.6
1986-03-01 1986-04-01 1986-05-01 1986-06-01 1986-07-01 1986-08-01 1986-09-01
59.6 237.8 108.2 144.8 81.2 141.0 69.8
1986-10-01 1986-11-01 1986-12-01 1987-01-01 1987-02-01 1987-03-01 1987-04-01
38.2 44.4 20.4 46.8 111.0 45.6 98.4
1987-05-01 1987-06-01 1987-07-01 1987-08-01 1987-09-01 1987-10-01 1987-11-01
212.0 153.8 221.8 175.0 90.6 278.8 164.8
1987-12-01 1988-01-01 1988-02-01 1988-03-01 1988-04-01 1988-05-01 1988-06-01
29.8 118.0 49.8 22.4 100.6 187.4 193.0
1
1988-07-01 1988-08-01 1988-09-01 1988-10-01 1988-11-01 1988-12-01 1989-01-01
120.4 149.2 61.2 136.4 10.0 59.4 0.0
1989-02-01 1989-03-01 1989-04-01 1989-05-01 1989-06-01 1989-07-01 1989-08-01
152.6 46.2 365.4 77.4 241.6 302.8 114.4
1989-09-01 1989-10-01 1989-11-01 1989-12-01 1990-01-01 1990-02-01 1990-03-01
65.4 12.8 145.0 110.6 51.6 12.4 65.8
1990-04-01 1990-05-01 1990-06-01 1990-07-01 1990-08-01 1990-09-01 1990-10-01
127.0 74.4 175.0 143.8 90.8 106.0 153.0
1990-11-01 1990-12-01
326.6 106.0
5. Dates of the daily values of 'x'
> dates <- time(x)
6. Amount of years in 'x' (needed for computations)
> ( nyears <- yip(from=start(x), to=end(x), out.type="nmbr" ) )
[1] 6
3 Basic Exploratory Data Analysis
1. Summary statistics
> smry(x)
Index x
Min. 1985-01-01 0.0000
1st Qu. 1986-07-02 0.0000
Median 1988-01-01 0.0000
Mean 1988-01-01 3.7470
3rd Qu. 1989-07-01 2.6000
Max. 1990-12-31 122.0000
IQR <NA> 2.6000
sd <NA> 10.0428
cv <NA> 2.6800
Skewness <NA> 5.3512
Kurtosis <NA> 39.1619
NA's <NA> 0.0000
n <NA> 2191.0000
2. Using the hydroplot function, which (by default) plots 9 di erent graphs: 3
ts plots, 3 boxplots and 3 histograms summarizing 'x'. For this example,
only daily and monthly plots are produced, and only data starting on
01-Jan-1987 are plotted.
此 |
评分
-
查看全部评分
|