data(mscstn)
data(mscdata)
head(mscdata)
stnids <- levels(mscdata$id)
# show stations and station names available in this data frame
data.frame(stnids, name=getstnname(stnids))
dat <- mksub(mscdata, id=1108447)
dat$fact <- mkfact(dat, "mon")
plot(mean_t ~ date, data=dat, type="l")
plot(mean_t ~ date, data=dat, subset=fact=="Dec", type="l")
plot.seas.temp(dat)
plot.year(dat)
Run the code above in your browser using DataLab