R help - May 2016 - Generating 3Dplot in lattice package

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.

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 <tr206 at 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 at 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]]