For starters, use 'pmin' (parallel min) instead of 'min'. substitute(MIN(psi/K14,EXP(((ABS(H14)/peak)^omega)*LN(psi/K14))), list(MIN=quote(pmin), K14=quote(VaR), ABS=quote(abs), EXP=quote(exp), LN=quote(log), H14=quote(Bmax), omega=quote(w2))) # pmin(psi/VaR, exp(((abs(Bmax)/peak)^w2) * log(psi/VaR)))
Bill Dunlap TIBCO Software wdunlap tibco.com On Mon, May 2, 2016 at 11:31 AM, T.Riedle <tr...@kent.ac.uk> wrote: > Dear R users, > > I am trying to generate a 3D plot using the wireframe() function in the > lattice package. > > The corresponding formula in Excel looks as follows and is applied to the > wireframe() function: > > MIN(psi/K14,EXP(((ABS(H14)/peak)^omega)*LN(psi/K14))) > > I tried to "translate" this formula in R and the code looks as follows > > > min(psi/VaR,exp(((abs(B)/Bmax)^w2)*log((psi/VaR),2.718182))) > > > > In this case I get only one value which is NAN but I should get a time > series with 13000 values. Hence, I deleted the min() function and tried the > wireframe() function: > > > > wireframe(inflator~exp((abs(B)/Bmax)^w2)*log((psi/VaR),base=2.718182),data=data_3Dplot) > > However, it doesn't work as the formula for the inflator is incorrect. It > produces NANs although it shouldn't do that and actually doesn't do that in > Excel. Furthermore, the results are totally different to those calculated > by Excel. > > Can anyone help me with the formula and translate it correctly in terms of > R? How do I get the time series? How do I get the values I get in Excel? > > I have attached the data as csv file to this email. > > Thank you very much for your help. > > ______________________________________________ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > 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. > [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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.