OzoneStgo(skewtools)
OzoneStgo()所属R语言包:skewtools
Ozone Time Series Data
臭氧时间序列数据
译者:生物统计家园网 机器人LoveR
描述----------Description----------
Ozone Time Series Data set of hourly observations (Santiago de Chile)
臭氧每小时的观测数据时间序列数据集(智利圣地亚哥)
用法----------Usage----------
data(OzoneStgo)
格式----------Format----------
A data frame with 744 observations on the following 9 variables.
744观察以下9个变量的数据框。
date a numeric vector of dates
date一个数字向量日期,
hour a numeric vector of hours
hour一个数字向量小时,
F a numeric vector of Independencia Station Ozone concentrations
F独立站臭氧浓度的一个数值向量
L a numeric vector of La Florida Station Ozone concentrations
L佛罗里达州的La站臭氧浓度的一个数值向量
M a numeric vector of Las Condes Station Ozone concentrations
M拉斯孔德斯站臭氧浓度的一个数值向量
N a numeric vector of Parque O'Higgins Station Ozone concentrations
N公园奥希金斯站臭氧浓度的一个数值向量
O a numeric vector of Pudahuel Station Ozone concentrations
O遭Pudahuel站臭氧浓度的一个数值向量
P a numeric vector of Cerrillos Station Ozone concentrations
P的塞里约斯站臭氧浓度的一个数值向量
Q a numeric vector of El Bosque Station Ozone concentrations
Q一个数字向量El博斯克车站的臭氧浓度,
Details
详细信息----------Details----------
This Dataset correspond to MACAM network in Seremi de Salud (March, 2006) of time series of ozone concentrations at 7 monitoring stations (F, L, M, N, O, P, Q) with n=5208 hourly observations
该数据集对应的MACAM网络的每期Seremi(2006年3月),7个监测站的臭氧浓度时间序列(F,L,M,N, O,P,Q)n=5208每小时的观测数据。
源----------Source----------
Data base available at Seremi de Salud (2006), Indices de Calidad del Aire, Santiago de Chile, http://www.seremisaludrm.cl/sitio/pag/aire/indexjs3aireindgasesdemo-prueba.asp
数据库,Seremi的每期(2006年),指数CALIDAD DEL艾尔,智利圣地亚哥,http://www.seremisaludrm.cl/sitio/pag/aire/indexjs3aireindgasesdemo-prueba.asp
参考文献----------References----------
Skew-Elliptical Distributions. Scandinavian Journal of Statistics, doi: 10.1111/j.1467-9469.2011.00774.x.
Atmospheric Environment, 37, 2337-2345.
实例----------Examples----------
data(OzoneStgo)
str(OzoneStgo)
data1=OzoneStgo[,3:9]
## Movile Average s-Function [Movile平均S-函数]
mmf <- function(X, s=1){
m=dim(X)[2]; t=dim(X)[1]
n=dim(X)[1]/s; datos.a=matrix(NA,t,m)
for(i in (s+1):t) for(j in 1:m){
aux=X[(i-s+1):i,j]
datos.a[i,j]=mean(aux[is.na(aux)==FALSE])
}
datos.a[(s+1):t,]
}
est=2 # L Station[L拨号站]
dat0=mmf(data1,s=8)[,est]
dat1=mmf(data1,s=16)[,est]
dat2=mmf(data1,s=24)[,est]
dat3=mmf(data1,s=32)[,est]
plot(data1[1 24*7),est],type="l",lty=2,
ylab=expression(paste(O[3] (mu ,g/m^3))),
xlab="time (hours)", main="Ozone - Station L")
lines(dat0[1 24*7)],type="l",lwd=2,col="yellow")
lines(dat1[1 24*7)],type="l",lwd=2,col="red")
lines(dat2[1 24*7)],type="l",lwd=2,col="green")
lines(dat3[1 24*7)],type="l",lwd=2,col="blue")
legend(122,90,c("Original",expression(MA[8]),expression(MA[16]),
expression(MA[24]),expression(MA[32])),
lty=c(2,1,1,1,1),lwd=c(1,2,2,2,2),
col=c("black","yellow","red","green","blue"))
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|