| mergeStation {clim.pact} | R Documentation |
Merges two series from different sources, eg from NACD
(getnacd) and DNMI (getdnmi). The code is
useful for updating long climate series with new observations from a
different database. The routine cokpares data for overlapping times
and prints out diagnostics about the two data sets.
mergeStation(x.1,x.2,plot=FALSE,print=TRUE,rescale=TRUE)
x.1 |
1st series. |
x.2 |
2nd series. |
plot |
'TRUE' plots the overlap. |
print |
'TRUE' prints diagnostics. |
rescale |
'TRUE' uses lm to determine amplitude of second sequence otherwise the amplitude of second sequence is not changed. |
A climate station series object.
R.E. Benestad
## Not run:
oslo.1 <- getnacd("OSLO-BLINDERN")
oslo.2 <- getdnmi("oslo")
print(range(oslo.1$yy))
#[1] 1890 1990
print(range(oslo.2$yy))
#[1] 1937 2002
oslo <- mergeStation(oslo.1,oslo.2)
#[1] "Time intervals:"
#[1] 1890 1990
#[1] 1937 2002
#[1] 1937.042 1990.958
#[1] "RMSE: 0.04"
#
#Call:
#lm(formula = y ~ 1 + x, data = ovrlp)
#
#Residuals:
# Min 1Q Median 3Q Max
#-7.24005 -0.03271 0.01161 0.06006 7.61593
#Coefficients:
# Estimate Std. Error t value Pr(>|t|)
#(Intercept) 0.029044 0.047482 0.612 0.541
#x 0.993886 0.004866 204.231 <2e-16 ***
#---
#Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1
#
#Residual standard error: 0.9738 on 644 degrees of freedom
#Multiple R-Squared: 0.9848, Adjusted R-squared: 0.9848
#F-statistic: 4.171e+04 on 1 and 644 DF, p-value: < 2.2e-16
print(range(oslo$yy))
#[1] 1890 2002
## End(Not run)