> On May 3, 2016, at 3:54 AM, T.Riedle <tr...@kent.ac.uk> wrote: > > Something is wrong here. The formula > > pmin(psi/VaR,exp(((abs(B)/Bmax)^w2)*log((psi/VaR),2.718182))) > > > provides a time series. Nevertheless, the returned values are incorrect and > it produces NANs.
The only thing we can conclude at the moment is that you have not provided sufficient information. The standard mailing list "sig" is there in an effort to prevent just this situation. > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. -- David. > > > > ________________________________ > From: William Dunlap <wdun...@tibco.com> > Sent: 02 May 2016 20:04 > To: T.Riedle > Cc: r-help@r-project.org > Subject: Re: [R] Generating 3Dplot in lattice package > > 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<http://tibco.com> > > On Mon, May 2, 2016 at 11:31 AM, T.Riedle > <tr...@kent.ac.uk<mailto: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. > > > https://stat.ethz.ch/mailman/listinfo/r-help David Winsemius Alameda, CA, USA ______________________________________________ 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.
