Learn R Programming

spphpr (version 1.0.0)

BJDAT: Daily Air Temperature Data of Beijing from 1952 to 2012.

Description

The data include the mean, minimum, and maximum daily temperatures (in \({}^{\circ}\)C) of Beijing between 1952 and 2012. Data source: China Meteorological Data Service Centre (https://data.cma.cn/en).

Usage

data(BJDAT)

Arguments

Details

In the data set, there are seven columns of vectors: Year, Month, Day, DOY, MDT, MinDT, and MaxDT. Year saves the recording years; Month saves the recording months; Day saves the recording days; DOY saves the dates in day-of-year; MDT saves the mean daily temperatures (in \({}^{\circ}\)C) corresponding to DOY; MinDT saves the minimum daily temperatures (in \({}^{\circ}\)C) corresponding to DOY; MaxDT saves the maximum daily temperatures (in \({}^{\circ}\)C) corresponding to DOY.

References

Guo, L., Xu, J., Dai, J., Cheng, J., Wu, H., Luedeling, E. (2015) Statistical identification of chilling and heat requirements for apricotflower buds in Beijing, China. Scientia Horticulturae 195, 138\(-\)144. tools:::Rd_expr_doi("10.1016/j.scienta.2015.09.006")

Examples

Run this code
data(BJDAT)
attach(BJDAT)

x    <- as.numeric( tapply(DOY, DOY, mean) )
y    <- as.numeric( tapply(MDT, DOY, mean) )
y.sd <- as.numeric( tapply(MDT, DOY, sd) )

dev.new()
par1 <- par(family="serif")
par2 <- par(mar=c(5, 5, 2, 2))
par3 <- par(mgp=c(3, 1, 0))
plot( x, y, cex = 1.5, xlim = c(0, 367), ylim = c(-10, 30), 
      cex.lab = 1.5, cex.axis = 1.5, type = "n", xlab = "Day-of-year", 
      ylab = expression(paste("Mean daily temperature (", degree, "C)", sep="")) )
for(i in 1:length(x)){
  lines(c(x[i], x[i]), c(y[i]-y.sd[i], y[i]+y.sd[i]), col=4)
}
points(x, y, cex = 1.5)
par(par1)
par(par2)
par(par3)

# graphics.off()

Run the code above in your browser using DataLab