[R] doubt with function on R software

[Pete Brecknock]

monicamir88 wrote
> Hello!
> I have a doubt with the R software. I have this function:
> 
> results <- function(bCODBOD, BOD, VSS, COD, sBOD, sCOD, TSS, TKNpa, T,
> NH3Ne, DO, Q, TKN, MLSS, NO3Ne, RAS, tanoxic1, tanoxic2, rbCOD, SDNR1,
> SDNR2, tdanoxic, tanaerobic, IRp, P) {
> bCOD <- bCODBOD*BOD
> nbCOD <- COD-bCOD
> nbVSS <- VSS*(1-((bCODBOD)*(BOD-sBOD)/(COD-sCOD)))
> iTSS <- TSS-VSS
> SF <- TKNpa
> knt <- 0.74*1.053^(T-20)
> kdnt <- 0.08*1.04^(T-20)
> unmt <- 0.75*1.07^(T-20)
> un <- (unmt*NH3Ne)/(knt+NH3Ne)*(DO/(0.5+DO))-kdnt
> kst <- 20*1^(T-20)
> kdt <- 0.12*1.04^(T-20)
> umt <- 6*1.07^(T-20)
> SRTtheo <- 1/un
> SRT <- SRTtheo*SF
> S <- kst*(1+kdt*SRT)/(SRT*(umt-kdt)-1)
> Pxbh <- Q*0.4*(bCOD-S)/(1+kdt*SRT)
> Pxdead <-  0.15*kdt*Q*0.4*(bCOD-S)*SRT/(1+kdt*SRT)
> Nox1 <- 0.8*TKN
> Pxbnit1 <- Q*0.12*Nox1/(1+kdnt*SRT)
> Pxbio1 <- Pxbh+Pxdead+Pxbnit1
> Nox2 <- TKN-NH3Ne-0.12*Pxbio1/Q
> Pxbnit2 <- Q*0.12*Nox2/(1+kdnt*SRT)
> Pxbio2 <- Pxbh+Pxdead+Pxbnit2
> Nox3 <- TKN-NH3Ne-0.12*Pxbio2/Q
> Pxbnit3<- Q*0.12*Nox3/(1+kdnt*SRT)
> Pxbio3 <- Pxbh+Pxdead+Pxbnit3
> Pxvss <- Pxbio3+Q*nbVSS
> Pxtss <- Pxbio3/0.88+Q*nbVSS+Q*(TSS-VSS)
> MassMLVSS <- Pxvss*SRT
> MassMLSS <- Pxtss*SRT
> Vaerobic <- MassMLSS/MLSS
> taerobic <- Vaerobic/Q
> thaerobic <- taerobic*24
> VSSfraction <- MassMLVSS/MassMLSS
> MLVSS <- MLSS*VSSfraction
> observedyield <- Pxtss/(Q*(bCOD-S))*bCODBOD
> Xb <- SRT*Q/Vaerobic*(0.4*bCOD/(1+kdt*SRT))
> IR <- Nox3/NO3Ne-1-RAS
> Noxfeed <- (IR*Q+RAS*Q)*NO3Ne
> Vanoxic1 <- tanoxic1*Q
> Vanoxic2 <- tanoxic2*Q
> FMb1 <- Q*BOD/(Vanoxic1*Xb)
> FMb2 <- Q*BOD/(Vanoxic2*Xb)
> rbCODbCOD <- rbCOD/bCOD
> SDNRt1 <- SDNR1*1.026^(T-20)
> SDNRt2 <- SDNR2*1.026^(T-20)
> NOr1 <- Vanoxic1*SDNRt1*Xb
> NOr2 <- Vanoxic2*SDNRt2*Xb
> Excesscap1 <- NOr1/Noxfeed
> Excesscap2 <- NOr2/Noxfeed
> Vanoxic <- tdanoxic*Q
> Vanaerobic <- tanaerobic*Q
> NO3reac <- NO3Ne*IRp/(1+IRp)
> rbCODremov <- NO3reac*6.6
> rbCODava <- rbCOD-rbCODremov
> biorem <- rbCODava/10
> Pbiomassg <- (Pxbh+Pxbnit3)*0.015/Q
> Premov <- Pbiomassg+biorem
> Peff <- P-Premov
>   
> res <- c(bCOD, nbCOD, nbVSS, iTSS, SF, knt, kdnt, unmt, un, kst, kdt, umt,
> SRTtheo, SRT, S, Pxbh, Pxdead, Nox1, Pxbnit1, Pxbio1, Nox2, Pxbnit2,
> Pxbio2, Nox3, Pxbnit3, Pxbio3 , Pxvss, Pxtss, MassMLVSS, MassMLSS,
> Vaerobic, taerobic, thaerobic, VSSfraction, MLVSS, observedyield, Xb, IR,
> Noxfeed, Vanoxic1, Vanoxic2, FMb1, FMb2, rbCODbCOD, SDNRt1, SDNRt2, NOr1,
> NOr2, Excesscap1, Excesscap2, Vanoxic, Vanaerobic, NO3reac, rbCODremov,
> rbCODava, biorem, Pbiomassg, Premov, Peff)
>   
> names(res) <- c("bCOD", "nbCOD", "nbVSS", "iTSS", "SF", "knt", "kdnt",
> "unmt", "un", "kst", "kdt", "umt", "SRTtheo", "SRT", "S", "Pxbh",
> "Pxdead", "Nox1"," Pxbnit1", "Pxbio1", "Nox2", "Pxbnit2", "Pxbio2",
> "Nox3", "Pxbnit3", "Pxbio3" , "Pxvss", "Pxtss", "MassMLVSS", "MassMLSS",
> "Vaerobic", "taerobic", "thaerobic", "VSSfraction", "MLVSS",
> "observedyield", "Xb", "IR", "Noxfeed", "Vanoxic1", "Vanoxic2", "FMb1",
> "FMb2", "rbCODbCOD", "SDNRt1", "SDNRt2", "NOr1", "NOr2", "Excesscap1",
> "Excesscap2", "Vanoxic", "Vanaerobic", "NO3reac", "rbCODremov",
> "rbCODava", "biorem", "Pbiomassg", "Premov", "Peff")
>   
> return(res) }
> 
> 
>  
> Where, for example, the values of the variables are:
> results(1.6, 140, 60, 300, 70, 132, 70, 1.5, 12, 0.5, 2, 22464, 35, 3000,
> 6, 0.6, 0.104, 0.0625, 80, 0.22, 0.31, 0.0625, 0.0625, 2, 6)
> 
> I want to take the values, for example, from a data.frame. Is it possible?
> 
> What I mean is, if I have a table with the values:
> 
> Variable	Value
>      Q	       22464
>      T                	12
>   BOD	        140
>      ...	           ...
> 
> 
> When I want to call the function, I indicate where to go for each value of
> the variable in the table. Thus, when changing the value of the variables
> in the table would be easier than doing in the function.
> 
> Thanks!

Hi

I took the liberty of simplifying things ..... 

# Define Your Function
myfun <- function(a,b,c){
 A = a + 1
 B = b + 2
 C = c + 3
 output <- data.frame(A,B,C)
 return(output)
}

# Test The Function
myfun(0,0,0)

# Input Data
inputData <- data.frame(a=c(1,2,3,4), b=c(10,20,30,40),
c=c(100,200,300,400))

# Ouput Data
outputData  <- do.call(myfun, inputData)

print(outputData)

HTH

Pete



--
View this message in context: http://r.789695.n4.nabble.com/doubt-with-function-on-R-software-tp4658973p4658990.html
Sent from the R help mailing list archive at Nabble.com.