Hi, how are you?
I have an attached file (StationIDs) that lists 74 observation sites
that I need to process using the EGRET/WRTDS R-package. For each site, I
need to load information into the Sample, INFO, and Daily fields.
The "sitenumber" in the following code refers to one of the 74
observation sites.
Is there anyway to process the following commands through a function
that exports the set of 74 AnnualResults (each of the 4 season results
in order from Fall to Summer) to a .csv file?
I would like to thank you in advance.
Irucka Embry
file <- "nwis_wqdata_sc_n.csv"
Sample <- getSampleDataFromFilecustom(file, hasHeader = TRUE, separator
=",")
Sample <- subset(Sample,station_id==sitenumber)
Sample <- subset(Sample, select=c(date, remark, value))
write.csv(Sample, file = "sitenumber_nitrogen.csv")
file <- "sitenumber_nitrogen.csv"
Sample <- getCustomSampleDataFromFilet(file, hasHeader = TRUE, separator
=",")
Sample <- subset(Sample,Julian>=55060)
Sample <- subset(Sample,Julian<=56520)
Sample<-removeDuplicates(Sample)
summary(Sample)
length(Sample$Date)
file <-"info_n.csv"
INFO<- getMetaDataFromFile(file)
INFO<- subset(INFO,site.no=="sitenumber")
file <- "sitenumber.rdb.gz"
Daily <- getDailyDataFromFileCustom(file)
Daily <- subset(Daily,Julian>=55060)
Daily <- subset(Daily,Julian<=56520)
summary(Daily)
Sample<-mergeReport()
modelEstimation(localDaily = Daily,localSample = Sample, localINFO =
INFO, windowY=10, windowQ=2, windowS=0.5,minNumObs=10,minNumUncen=10,
env=parent.frame())
INFO<-setPA(window = 4)
annualSeries<-makeAnnualSeries()
plotFlowSingle(istat=5)
plotFlowSingle(istat=2)
printSeries(istat=5)
printSeries(istat=2)
plotConcTime()
AnnualResults<-setupYears()
AnnualResults
tableResults(localAnnualResults = AnnualResults, localINFO = INFO, qUnit
= 1, fluxUnit = 13)
yearPoints<-c(2001,2002,2003,2004)
tableChange(localAnnualResults = AnnualResults, localINFO = INFO,
fluxUnit = 13, yearPoints)
AnnualResults<- setupYears(paLong = 3, paStart = 10, localDaily = Daily)
AnnualResults
AnnualResults<- setupYears(paLong = 3, paStart = 1, localDaily = Daily)
AnnualResults
AnnualResults<- setupYears(paLong = 3, paStart = 4, localDaily = Daily)
AnnualResults
AnnualResults<- setupYears(paLong = 3, paStart = 7, localDaily = Daily)
AnnualResults
<span id=m2wTl><p><font face="Arial, Helvetica, sans-serif" size="2"
style="font-size:13.5px">_______________________________________________________________<BR>Get
the Free email that has everyone talking at <a href=http://www.mail2world.com
target=new>http://www.mail2world.com</a><br> <font color=#999999>Unlimited
Email Storage – POP3 – Calendar – SMS – Translator –
Much More!</font></font></span>
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.