Register to reply

Delta method fails. Any suggestions how to calculate E(Y)?

by Hejdun
Tags: delta, fails, method, suggestions
Share this thread:
Hejdun
#1
Mar31-12, 08:54 PM
P: 25
Hi,

Let Z be a r.v. from N(0, σ^2) and let Y = Z^2*exp(Z)/(1 + exp(Z))^2.
What is the exptected value of Y, E(Y)?

The delta method (Taylor expansion) doesn't work since
the errors seem to accumulate.

Any suggestions? This is a non-standard problem related to my research and I am totally stuck with this one...

// H
Phys.Org News Partner Science news on Phys.org
Security CTO to detail Android Fake ID flaw at Black Hat
Huge waves measured for first time in Arctic Ocean
Mysterious molecules in space
chiro
#2
Mar31-12, 10:21 PM
P: 4,572
Quote Quote by Hejdun View Post
Hi,

Let Z be a r.v. from N(0, σ^2) and let Y = Z^2*exp(Z)/(1 + exp(Z))^2.
What is the exptected value of Y, E(Y)?

The delta method (Taylor expansion) doesn't work since
the errors seem to accumulate.

Any suggestions? This is a non-standard problem related to my research and I am totally stuck with this one...

// H
Hey Hejdun and welcome to the forums.

If it's hard to get an analytic answer for the distribution, you should do a simulation and then create a simulated distribution for Y in which you can can easily calculate E[Y].

I would use R to do this: it's free and it would take you probably ten minutes to get this result up on screen.

The commands you need to use

z = rnorm(10000,0,sd=whatever) and then y = z^2*exp(z)/(1 + exp(z))^2 and then mean(y). It should take only ten seconds or less to execute.
Hejdun
#3
Apr1-12, 05:31 AM
P: 25
Quote Quote by chiro View Post
Hey Hejdun and welcome to the forums.

If it's hard to get an analytic answer for the distribution, you should do a simulation and then create a simulated distribution for Y in which you can can easily calculate E[Y].

I would use R to do this: it's free and it would take you probably ten minutes to get this result up on screen.

The commands you need to use

z = rnorm(10000,0,sd=whatever) and then y = z^2*exp(z)/(1 + exp(z))^2 and then mean(y). It should take only ten seconds or less to execute.

Thanks for your help. However, I am not looking for a numerical solution. I am aware of that no exact analytical solution exists, but maybe there is some approximation that is "good enough".

If you plot the function in R, then it is clear that Y=f(Z) is a symmetric bimodal distribution. Perhaps it simplifies if we put a restriction and look at Y when Z>0.

But anyway, I am stuck right now. Perhaps not solution exist.

chiro
#4
Apr1-12, 05:48 AM
P: 4,572
Delta method fails. Any suggestions how to calculate E(Y)?

Quote Quote by Hejdun View Post
Thanks for your help. However, I am not looking for a numerical solution. I am aware of that no exact analytical solution exists, but maybe there is some approximation that is "good enough".

If you plot the function in R, then it is clear that Y=f(Z) is a symmetric bimodal distribution. Perhaps it simplifies if we put a restriction and look at Y when Z>0.

But anyway, I am stuck right now. Perhaps not solution exist.
Instead of using a normal, why don't you use the square root of a chi-square. This distribution will act like a positive normal distribution and is only defined from 0 onwards.
Hejdun
#5
Apr1-12, 06:02 AM
P: 25
Quote Quote by chiro View Post
Instead of using a normal, why don't you use the square root of a chi-square. This distribution will act like a positive normal distribution and is only defined from 0 onwards.
Thanks again for you interest in my problem and your suggestions.
However, I am not sure how your suggestion in this case would help. What chi-square distribution and how many degrees of freedom?

BTW, I tried doing a Padé Approximation instead of Taylor expansion, but that didn't help either.

/H
awkward
#6
Apr1-12, 08:30 AM
P: 327
Since Y(Z) = Y(-Z) and Z's distribution is symmetric about 0, doesn't this imply that E(Y) = 0?
Hejdun
#7
Apr1-12, 08:39 AM
P: 25
Quote Quote by awkward View Post
Since Y(Z) = Y(-Z) and Z's distribution is symmetric about 0, doesn't this imply that E(Y) = 0?
Thanks for your answer.

That is true. I think the assumption Y(Z)=Y(-Z) could be wrong however (even though they are very similar), since when I run simulations to check E(Y) then E(Y) ≠ 0.
Hurkyl
#8
Apr1-12, 09:09 AM
Emeritus
Sci Advisor
PF Gold
Hurkyl's Avatar
P: 16,092
Quote Quote by Hejdun View Post
Thanks for your help. However, I am not looking for a numerical solution. I am aware of that no exact analytical solution exists, but maybe there is some approximation that is "good enough".
I don't understand this. Why can't a numerical solution be "good enough"?
Hejdun
#9
Apr1-12, 09:13 AM
P: 25
Quote Quote by Hurkyl View Post
I don't understand this. Why can't a numerical solution be "good enough"?
Because this expectation is a part of a larger derivation of a proof and needs to be more general.

But thanks for you comment.
awkward
#10
Apr1-12, 09:56 AM
P: 327
Quote Quote by Hejdun View Post
Thanks for your answer.

That is true. I think the assumption Y(Z)=Y(-Z) could be wrong however (even though they are very similar), since when I run simulations to check E(Y) then E(Y) ≠ 0.
It's not an assumption, it's an identity. See if you can prove it.

[edit] Oops... Although it's true that Y(Z) = Y(-Z), unfortunately that's not what you would need to get a mean of zero. For that you would need Y(Z) = -Y(Z). So please ignore what I wrote. [/edit]


Register to reply

Related Discussions
Calculate delta H for the combustion of 3-methylpentan-2-one. Biology, Chemistry & Other Homework 0
I need a method to calculate ln(x) for small x, other than Taylor series method General Math 12
Delta-epsilon method Calculus 9
Need help designing a lab to calculate delta H? Biology, Chemistry & Other Homework 1
How do you calculate delta G Chemistry 2