If it helps... R is not Matlab. There is no '.*' operator in R.
/H
On Thu, Apr 26, 2012 at 2:37 PM, Sarah Goslee <sarah.gos...@gmail.com> wrote:
> On Thu, Apr 26, 2012 at 5:34 PM, Guaramy _ <guar...@hotmail.com> wrote:
>> Thanks for your answer, but how can i correct that, the mathematical
>> expression is correct that way with that g. Do you know how can i program it
>> ?
>
> Then you need to define g. and how it differs from g alone.
>
> g. and g are two separate objects in R.
>
> Sarah
>
> PS You also need to convert your email program to plain text, or
> whatever caused the horrid formatting below.
>
>>
>>> Date: Thu, 26 Apr 2012 13:57:36 -0400
>>> Subject: Re: [R] ErrError in f(x, ...) : object 'g.' not found
>>> From: sarah.gos...@gmail.com
>>> To: guar...@hotmail.com
>>> CC: r-help@r-project.org
>>>
>>> Just what it says:
>>>
>>> You define g but refer to a variable g. in the next line.
>>>
>>> Just get rid of the typo.
>>>
>>> Sarah
>>>
>>> On Thu, Apr 26, 2012 at 1:43 PM, Guaramy <guar...@hotmail.com> wrote:
>>> > Hi , R is a new language for me so sorry in advance if this error is to
>>> > basic
>>> > for posting. I have tried the R manual and search online for quite a
>>> > few, if
>>> > anyone could help i would be very thankful.
>>> > Here is ámy code.
>>> >
>>> > kappa = 1.1
>>> > theta = 0.1
>>> > sigma = 0.4
>>> > rho = -0.6
>>> > v0 = 0.2
>>> > r = 0.05
>>> > T = 0.5
>>> > s0 = 1
>>> > K = 0.5
>>> > type = 1
>>> > Hestoncall = function(kappa,theta,sigma,rho,v0,r,T,s0,K,type)
>>> > {
>>> >
>>> > á á á á á á á á á á á á á á á á u = 0.5
>>> > á á á á á á á á á á á á á á á á b = kappa-rho*sigma
>>> > á á á á á á á á á á á á á á á á a = kappa*theta
>>> > á á á á á á á á á á á á á á á á x = log(s0)
>>> >
>>> > á á á á á á á áHestf = function(phi)
>>> > á á á á á á á á{
>>> >
>>> > á á á á á á á á á á á á á á á á d =
>>> > sqrt((b-rho*sigma*phi*complex(1,0,1)-b )^2-sigma^2*(2*u*phi^2))
>>> > á á á á á á á á á á á á á á á á g =
>>> > (b-rho*sigma*phi*complex(1,0,1)+d)/(b-rho*sigma*phi*complex(1,0,1)-d)
>>> > á á á á á á á á á á á á á á á á C = r*phi*complex(1,0,1)*T +
>>> >
>>> > a/sigma^2*((b-rho*sigma*phi*complex(1,0,1)+d)*T-2*log((1-g.*exp(d*T))/(1-g)))
>>> > á á á á á á á á á á á á á á á á D =
>>> > (b-rho*sigma*phi*complex(1,0,1)+d)/sigma^2*((1-exp(d*T))/(1-g*exp(d*T)))
>>> >
>>> > á á á á á á á á á á á á á á á á f = exp(C+D*v0 +complex(1,0,1)*phi*x);
>>> >
>>> >
>>> > á á á á á á á á á á á á áHestonPintegrand á= real
>>> > (exp(-1i*phi*log(K))*f/(1i*phi))
>>> >
>>> > á á á á á á á á á á á á áHestonPintegrand
>>> >
>>> > á á á á á á á á á á á á}
>>> >
>>> > á á á á #int.fn = function(t){sapply(t,FUN=Hestf)}
>>> > á á á á IH =
>>> > integrate(Hestf,lower=0,upper=Inf,rel.tol=1e-10,subdivisions=1000000)
>>> > á á á á ret = 0.5 + 1/pi*IH
>>> >
>>> >
>>> >
>>> > á á á ácall = s0*ret - K*exp(r-T)*ret
>>> > á á á ácall
>>> >
>>> > }
>>> >
>>> > thanks in advance
>>> >
>
> --
> Sarah Goslee
> http://www.functionaldiversity.org
>
> ______________________________________________
> 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.
______________________________________________
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.
