
fitswavecav(cdat, cavdat, tanm = "trend1", pnames,
yrstart = 0, yrend = 0, tndbeg = 0, tndend = 0,
iwcav = c("none"), dcol = "dates",
qwcols = c("R", "P"), mclass = 1)
survreg.object
loglik numeric the log-likelihood for the model
cint numeric coefficient for model intercept
cwave numeric coefficient for the seasonal wave
ctnd numeric coefficient for the trend component of model
c[alphanumeric] numeric 0 or more coefficients for the
continuous ancillary variables
seint numeric standard error for the intercept
sewave numeric standard error for the seasonal wave
setnd numeric standard error for the trend
se[alphanumeric] numeric 0 or more standard errors for the
continuous ancillary variables
pvaltnd numeric the p-value for the trend line
} U.S. Geological Survey, 2013a, National Water Information
System: Web Interface, accessed Febaruary 26, 2013, at
U.S. Geological Survey, 2013b, Parameter code definition:
National Water Information System: Web Interface,
accessed Febaruary 26, 2013, at
Vecchia, A.V., Martin, J.D., and Gilliiom, R.J., 2008,
Modeling variability and trends in pesticide
concentrations in streams: Journal of the American Water
Resources Association, v. 44, no. 5, p. 1308-1324,
fitswavecav
calls internally:
prepData
and fitMod
.data(swData)
modMoRivOmaha<-combineData(qwdat=qwMoRivOmaha, cqwdat=cqwMoRivOmaha)
myfit1 <- fitswavecav(cdat=modMoRivOmaha, cavdat=cqwMoRivOmaha,
tanm="myfit1", pnames=c("04035", "04037", "04041"), yrstart=1995,
yrend=2003, tndbeg=1995, tndend=2003, iwcav=c("flowa30","flowa1"),
dcol="dates", qwcols=c("R","P"))
myfit2 <- fitswavecav(cdat=modMoRivOmaha, cavdat=cqwMoRivOmaha,
tanm="myfit2", pnames=c("04035", "04037", "04041"), yrstart=1995,
yrend=2003, tndbeg=1995, tndend=2003, iwcav=c("seda30","seda1"),
dcol="dates", qwcols=c("R","P"))
myfit3 <- fitswavecav(cdat=modMoRivOmaha, cavdat=cqwMoRivOmaha,
tanm="myfit3", pnames=c("04035", "04037", "04041"), yrstart=1995,
yrend=2003, tndbeg=1995, tndend=2003, iwcav=c("flowa30","flowa1",
"seda30", "seda1"), dcol="dates", qwcols=c("R","P"))
# trend model results
myfit1[[1]]
# example regression call
myfit1[[2]][[1]]
# first few lines of observed concentrations
head(myfit1[[3]])
# first few lines of predicted concentrations
head(myfit1[[4]])
# summary statistics for predicted concentrations
head(myfit1[[5]])
Run the code above in your browser using DataLab