Hi Arun, thank-you very much! The 2nd option worked perfectly. That was what I wanted.
Now, I have another question. I am using the R packages dataRetrieval and EGRET from https://github.com/USGS-CIDA/WRTDS. I have 2 objects Daily and Sample that have the naming convention (Names = "21NC02WQ.C1000000" or whatevver the list of site names happens to be) that I need to have after running modelEstimation, but the names are stripped after running modelEstimation (see below). The new names after running modelEstimation for both Sample and Daily are the names of the columns rather than the site name. ## modelEstimation -- this is the code that I use to do the loop through modelEstimation for(i in Daily) { modelEstimation(localDaily = Daily[[i]], localSample = Sample[[i]], localINFO = INFO[[i]], windowY=10, windowQ=2, windowS=0.5,minNumObs=10,minNumUncen=6, env=parent.frame()) } # shows the correct naming convention before being ran through modelEstimation > dput(Daily) structure(list(`21NC02WQ.C1000000` = structure(list(Date = structure(c(11231, 11232, 11233, 11234, 11235, 11236), class = "Date"), Q = c(0.886534950681092, 0.857937049046218, 0.829339147411344, 0.80074124577647, 0.772143344141596, 0.743545442506722), Julian = c(55060, 55061, 55062, 55063, 55064, 55065), Month = c(10, 10, 10, 10, 10, 10), Day = c(275, 276, 277, 278, 279, 280), DecYear = c(2000.75, 2000.75273224044, 2000.75546448087, 2000.75819672131, 2000.76092896175, 2000.76366120219), MonthSeq = c(1810, 1810, 1810, 1810, 1810, 1810), i = 1:6, LogQ = c(-0.120434728772853, -0.153224551595844, -0.187126103271526, -0.222217423082796, -0.258585067253671, -0.296325395236517)), .Names = c("Date", "Q", "Julian", "Month", "Day", "DecYear", "MonthSeq", "i", "LogQ"), row.names = c(NA, -6L), class = "data.frame")), .Names = "21NC02WQ.C1000000") # shows the correct naming convention before being ran through modelEstimation > dput(Sample) structure(list(`21NC02WQ.C1000000` = structure(list(Date = structure(numeric(0), class = "Date"), ConcLow = numeric(0), ConcHigh = numeric(0), Uncen = numeric(0), ConcAve = numeric(0), Julian = numeric(0), Month = numeric(0), Day = numeric(0), DecYear = numeric(0), MonthSeq = numeric(0), SinDY = numeric(0), CosDY = numeric(0)), .Names = c("Date", "ConcLow", "ConcHigh", "Uncen", "ConcAve", "Julian", "Month", "Day", "DecYear", "MonthSeq", "SinDY", "CosDY"), row.names = integer(0), class = "data.frame")), .Names = "21NC02WQ.C1000000") # shows the Daily incorrect naming convention after being ran through modelEstimation )))), .Names = c("Date", "Q", "Julian", "Month", "Day", "DecYear", "MonthSeq", "i", "LogQ", "Q7", "Q30", "yHat", "SE", "ConcDay", "FluxDay", "FNConc", "FNFlux"), row.names = c(NA, -1461L), class = "data.frame") # shows the Sample incorrect naming convention after being ran through modelEstimation , .Names = c("Date", "ConcLow", "ConcHigh", "Uncen", "ConcAve", "Julian", "Month", "Day", "DecYear", "MonthSeq", "SinDY", "CosDY", "Q", "LogQ", "yHat", "SE", "ConcHat" ), row.names = 176:198, class = "data.frame") Is there anyway to make the Name of the object Daily and Sample the site name after being processed through modelEstimation, i.e. making the current names the column names which they should be? I also have a question about the INFO object that I will address in the next e-mail. Thank you for all of your continued great assistance Arun! Irucka ## the modelEstimation function and its major components are below modelEstimation function (localDaily = Daily, localSample = Sample, localINFO = INFO, windowY = 10, windowQ = 2, windowS = 0.5, minNumObs = 100, minNumUncen = 50, env = parent.frame()) { cat("\n first step running estCrossVal may take about 1 minute") Sample1 <- estCrossVal(SampleCrossV = localSample, windowY, windowQ, windowS, minNumObs, minNumUncen) surfaceIndexParameters <- surfaceIndex(localDaily = localDaily) localINFO$bottomLogQ <- surfaceIndexParameters[1] localINFO$stepLogQ <- surfaceIndexParameters[2] localINFO$nVectorLogQ <- surfaceIndexParameters[3] localINFO$bottomYear <- surfaceIndexParameters[4] localINFO$stepYear <- surfaceIndexParameters[5] localINFO$nVectorYear <- surfaceIndexParameters[6] localINFO$windowY <- windowY localINFO$windowQ <- windowQ localINFO$windowS <- windowS localINFO$minNumObs <- minNumObs localINFO$minNumUncen <- minNumUncen cat("\nNext step running estSurfaces with survival regression:\n") surfaces1 <- estSurfaces(localDaily = localDaily, localSample = localSample, windowY, windowQ, windowS, minNumObs, minNumUncen) Daily1 <- estDailyFromSurfaces(localDaily = localDaily, localINFO = localINFO, localsurfaces = surfaces1) env$Daily <- Daily1 env$INFO <- localINFO env$Sample <- Sample1 env$surfaces <- surfaces1 cat("\nDone with modelEstimation,\nnow do AnnualResults<-setupYears()\nor if using a period of analysis other than Water Year specify the arguments paStart and paLong in call to setupYears ") } <environment: namespace:EGRET> > estSurfaces function (localDaily = Daily, localSample = Sample, windowY = 10, windowQ = 2, windowS = 0.5, minNumObs = 100, minNumUncen = 50) { bottomLogQ <- min(localDaily$LogQ) - 0.05 topLogQ <- max(localDaily$LogQ) + 0.05 stepLogQ <- (topLogQ - bottomLogQ)/13 vectorLogQ <- seq(bottomLogQ, topLogQ, stepLogQ) stepYear <- 1/16 bottomYear <- floor(min(localDaily$DecYear)) topYear <- ceiling(max(localDaily$DecYear)) vectorYear <- seq(bottomYear, topYear, stepYear) nVectorYear <- length(vectorYear) estPtLogQ <- rep(vectorLogQ, nVectorYear) estPtYear <- rep(vectorYear, each = 14) resultSurvReg <- runSurvReg(estPtYear, estPtLogQ, localSample, windowY, windowQ, windowS, minNumObs, minNumUncen) surfaces <- array(0, dim = c(14, nVectorYear, 3)) for (iQ in 1:14) { for (iY in 1:nVectorYear) { k <- (iY - 1) * 14 + iQ surfaces[iQ, iY, ] <- resultSurvReg[k, ] } } return(surfaces) } <environment: namespace:EGRET> > estDailyFromSurfaces function (localDaily = Daily, localINFO = INFO, localsurfaces = surfaces) { library("fields") LogQ <- seq(localINFO$bottomLogQ, by = localINFO$stepLogQ, length.out = localINFO$nVectorLogQ) Year <- seq(localINFO$bottomYear, by = localINFO$stepYear, length.out = localINFO$nVectorYear) localDaily$yHat <- interp.surface(obj = list(x = LogQ, y = Year, z = localsurfaces[, , 1]), loc = data.frame(localDaily$LogQ, localDaily$DecYear)) localDaily$SE <- interp.surface(obj = list(x = LogQ, y = Year, z = localsurfaces[, , 2]), loc = data.frame(localDaily$LogQ, localDaily$DecYear)) localDaily$ConcDay <- interp.surface(obj = list(x = LogQ, y = Year, z = localsurfaces[, , 3]), loc = data.frame(localDaily$LogQ, localDaily$DecYear)) localDaily$FluxDay <- localDaily$ConcDay * localDaily$Q * 86.4 allLogQsByDayOfYear <- split(localDaily$LogQ, localDaily$Day) allLogQsByDayOfYear[["365"]] <- c(unlist(allLogQsByDayOfYear["365"]), unlist(allLogQsByDayOfYear["366"])) allLogQsByDayOfYear["366"] <- allLogQsByDayOfYear["365"] allLogQsReplicated <- allLogQsByDayOfYear[localDaily$Day] allDatesReplicated <- rep(localDaily$DecYear, lapply(allLogQsReplicated, length)) allConcReplicated <- interp.surface(obj = list(x = LogQ, y = Year, z = localsurfaces[, , 3]), loc = data.frame(unlist(x = allLogQsReplicated), y = allDatesReplicated)) allFluxReplicated <- allConcReplicated * exp(unlist(allLogQsReplicated)) * 86.4 localDaily$FNConc <- tapply(allConcReplicated, allDatesReplicated, "mean") localDaily$FNFlux <- tapply(allFluxReplicated, allDatesReplicated, "mean") return(localDaily) } <environment: namespace:EGRET> <-----Original Message-----> >From: arun kirshna [via R] [ml-node+s789695n4661331...@n4.nabble.com] >Sent: 3/14/2013 7:30:54 AM >To: iruc...@mail2world.com >Subject: Re: "[[i]]$" <- "" indexing and lapply > >Hi Irucka, > >Not sure how `Q` should be attached to Data object. >Try this: >const<- list(1.01, 1.00, 0.99) >DailyModeled1<-mapply(cbind,DailyModeled,const=const,SIMPLIFY=FALSE) > >DailyModeled$Q<-lapply(seq_along(Daily),function(i) Daily[[i]]$Q*DailyModeled1[[i]]$const) > str(DailyModeled) #different than what you showed. Also, I noticed that 'Q' is the same in your 'str' >#List of 4 ># $ :'data.frame': 4 obs. of 1 variable: ># ..$ Date: Date[1:4], format: "2000-10-01" "2000-10-02" ... ># $ :'data.frame': 4 obs. of 1 variable: ># ..$ Date: Date[1:4], format: "2000-10-01" "2000-10-02" ... ># $ :'data.frame': 4 obs. of 1 variable: ># ..$ Date: Date[1:4], format: "2000-10-01" "2000-10-02" ... ># $ Q:List of 3 ># ..$ : num [1:4] 1.316 1.03 0.944 0.858 ># ..$ : num [1:4] 0.397 0.369 0.312 0.281 ># ..$ : num [1:4] 76.3 78.8 78.5 75.7 > >#Or >#Maybe you want something like this: >DailyModeled[]<-lapply(seq_along(DailyModeled1),function(i) {DailyModeled[[i]]$Q<- >Daily[[i]]$Q*DailyModeled1[[i]]$const;DailyModeled[[i]]} ) > DailyModeled >#[[1]] ># Date Q >#1 2000-10-01 1.3156007 >#2 2000-10-02 1.0296005 >#3 2000-10-03 0.9438005 >#4 2000-10-04 0.8580004 > >#[[2]] ># Date Q >#1 2000-10-01 0.3969416 >#2 2000-10-02 0.3686247 >#3 2000-10-03 0.3119910 >#4 2000-10-04 0.2808425 > >#[[3]] ># Date Q >#1 2000-10-01 76.25160 >#2 2000-10-02 78.77463 >#3 2000-10-03 78.49430 >#4 2000-10-04 75.69093 > > str(DailyModeled) >#List of 3 ># $ :'data.frame': 4 obs. of 2 variables: ># ..$ Date: Date[1:4], format: "2000-10-01" "2000-10-02" ... ># ..$ Q : num [1:4] 1.316 1.03 0.944 0.858 ># $ :'data.frame': 4 obs. of 2 variables: ># ..$ Date: Date[1:4], format: "2000-10-01" "2000-10-02" ... ># ..$ Q : num [1:4] 0.397 0.369 0.312 0.281 ># $ :'data.frame': 4 obs. of 2 variables: ># ..$ Date: Date[1:4], format: "2000-10-01" "2000-10-02" ... ># ..$ Q : num [1:4] 76.3 78.8 78.5 75.7 >A.K. > > > > > >If you reply to this email, your message will be added to the discussion below: >http://r.789695.n4.nabble.com/i-indexing-and-lapply-tp4661286p4661331.h tml >To unsubscribe from "[[i]]$" <- "" indexing and lapply, click here. >NAML <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> -- View this message in context: http://r.789695.n4.nabble.com/i-indexing-and-lapply-tp4661286p4662103.html Sent from the R help mailing list archive at Nabble.com. [[alternative HTML version deleted]] ______________________________________________ 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.
