- #1
jcates7
- 3
- 0
Hello,
I'm trying to calculate Fisher information (and eventually the Cramer-Rao lower bound) for this particular pdf with Mathematica:
\text{pte}[t,\Theta ] = \frac{P_{\text{ec}}}{\tau _d-\tau _r}\left[e^{\frac{-(t-\Theta )}{\tau _d}}-e^{\frac{-(t-\Theta )}{\tau _r}}\right]; \text{domain}[\text{pte}] = \{t,-\infty,\infty\}\&\&\{\Theta >0\}
So I want to find I(t|\Theta )
This seems relatively straightforward with:
\text{Integrate}\left[D[\text{Log}[\text{pte}],\Theta ]^2,\{t,-\infty ,\infty \}\right]
However, Mathematica doesn't want to compute the integral. It just returns the integral itself:
\int_{-\infty}^{\infty} \frac{\left(\frac{e^{-\frac{t-\Theta }{\tau _d}}}{\tau _d}-\frac{e^{-\frac{t-\Theta }{\tau _r}}}{\tau _r}\right){}^2 \left(\frac{P_{\text{ec}}}{\tau _d-\tau _r}\right)'\left[e^{-\frac{t-\Theta }{\tau _d}}-e^{-\frac{t-\Theta }{\tau _r}}\right]{}^2}{\frac{P_{\text{ec}}}{\tau _d-\tau _r}\left[e^{-\frac{t-\Theta }{\tau _d}}-e^{-\frac{t-\Theta }{\tau _r}}\right]{}^2} \, dt
Initially my thought is that there isn't a closed-form solution, but this is something I have seen calculated in journal papers with the same pdf. I'm not an experience Mathematica user. Is there something I have missed in the input (syntax or additional options) or are there any general simplifications that Mathematica would need?
Thanks!
I'm trying to calculate Fisher information (and eventually the Cramer-Rao lower bound) for this particular pdf with Mathematica:
\text{pte}[t,\Theta ] = \frac{P_{\text{ec}}}{\tau _d-\tau _r}\left[e^{\frac{-(t-\Theta )}{\tau _d}}-e^{\frac{-(t-\Theta )}{\tau _r}}\right]; \text{domain}[\text{pte}] = \{t,-\infty,\infty\}\&\&\{\Theta >0\}
So I want to find I(t|\Theta )
This seems relatively straightforward with:
\text{Integrate}\left[D[\text{Log}[\text{pte}],\Theta ]^2,\{t,-\infty ,\infty \}\right]
However, Mathematica doesn't want to compute the integral. It just returns the integral itself:
\int_{-\infty}^{\infty} \frac{\left(\frac{e^{-\frac{t-\Theta }{\tau _d}}}{\tau _d}-\frac{e^{-\frac{t-\Theta }{\tau _r}}}{\tau _r}\right){}^2 \left(\frac{P_{\text{ec}}}{\tau _d-\tau _r}\right)'\left[e^{-\frac{t-\Theta }{\tau _d}}-e^{-\frac{t-\Theta }{\tau _r}}\right]{}^2}{\frac{P_{\text{ec}}}{\tau _d-\tau _r}\left[e^{-\frac{t-\Theta }{\tau _d}}-e^{-\frac{t-\Theta }{\tau _r}}\right]{}^2} \, dt
Initially my thought is that there isn't a closed-form solution, but this is something I have seen calculated in journal papers with the same pdf. I'm not an experience Mathematica user. Is there something I have missed in the input (syntax or additional options) or are there any general simplifications that Mathematica would need?
Thanks!
