Skip to main content
added 40 characters in body
Source Link
Michael Albanese
  • 19.3k
  • 9
  • 87
  • 160

I have the following problem. If iI can get some help, I would greatly appreciate it. II am trying to replicate a particular research paper and came across a problem:

Suppose X$X$ is a birth death process (represents population size) that evolves by:

$X -> X+1 $$X \to X+1$ if a birth occurs with rate $\mu$,

$X -> X-1 $$X \to X-1$ if a death occurs with rate $\theta$.

Suppose $T_A$ is first passage time of a BD process from state A$A$ to state 0$0$ and suppose $T_B$ is first passage time of another BD process from state B$B$ to state 0$0$. TheyThey are both independent.

I need to find $P(T_A < T_B)$. That is, probability that a population of size A$A$ goes to 0$0$ before population of size B$B$.

By definition:

$T_A = T_{A,A-1} + T_{A-1,A-2} + ... + T_{1,0}$$$T_A = T_{A,A-1} + T_{A-1,A-2} + \dots + T_{1,0}$$

where $T_{i,i-1}$ represents first passage time from state i,$i$ to state i-1$i-1$.

I read some articles online that mentioned that if $T = T_A - T_B$, then the CDF defined by:

$G_T(t) := P(T <=t)$$$G_T(t) := P(T \leq t)$$

is what I need. TheThe paper suggested taking inverse Laplace of a CDF to obtain the CDF and evaluate it at 0$0$. ItIt first suggested finding Laplace transform of T$T$, which is given by

$L[T] = E(e^{-ST}) = -\frac{L[T_A](s) L[T_B](-s)}{s}$$$L[T] = E(e^{-ST}) = -\frac{L[T_A](s) L[T_B](-s)}{s}.$$

Then it suggested taking Laplace of $G_T(t)$ i.e  . $L[G_T(t)]$. However,

However, $L[G_T(t)] = \frac{L[T]}{s} = \frac{L[T_A](s) L[T_B](-s)}{s}$$$L[G_T(t)] = \frac{L[T]}{s} = \frac{L[T_A](s) L[T_B](-s)}{s}.$$

Then the paper suggests taking the inverse of the above evaluated at 0$0$ to get $P(T_A < T_B)$.

QuestionQuestions:

  1. Given the Laplace transform of CDF, which is $L[G_T(t)]$, I want to use the inverse laplace to obtain $G_T(t)$ evaluated at 0. But, by definition, inverse laplace using algorithms in python are all one sided from $[0,\infty]$ My random variable T is given by difference of two first passage times, $T = T_A - T_B $. Won't this be negative?

  2. In the paper, it says they are shifting the random variable X under study by a constant c such that P(X + c > 0) is approximately 1. Then inverting the corresponding one sided Laplace transform. How would I do that in this context here?

  1. Given the Laplace transform of CDF, which is $L[G_T(t)]$, I want to use the inverse Laplace to obtain $G_T(t)$ evaluated at $0$. But, by definition, inverse Laplace using algorithms in python are all one sided from $[0,\infty]$. My random variable $T$ is given by difference of two first passage times, $T = T_A - T_B$. Won't this be negative?

Thanks

  1. In the paper, it says they are shifting the random variable $X$ under study by a constant $c$ such that $P(X + c > 0)$ is approximately $1$. Then inverting the corresponding one sided Laplace transform. How would I do that in this context here?

I have the following problem. If i can get some help, I would greatly appreciate it. I am trying to replicate a particular research paper and came across a problem:

Suppose X is a birth death process (represents population size) that evolves by:

$X -> X+1 $ if a birth occurs with rate $\mu$

$X -> X-1 $ if a death occurs with rate $\theta$

Suppose $T_A$ is first passage time of a BD process from state A to state 0 and suppose $T_B$ is first passage time of another BD process from state B to state 0. They are both independent.

I need to find $P(T_A < T_B)$. That is, probability that a population of size A goes to 0 before population of size B.

By definition:

$T_A = T_{A,A-1} + T_{A-1,A-2} + ... + T_{1,0}$

where $T_{i,i-1}$ represents first passage time from state i, to state i-1.

I read some articles online that mentioned that if $T = T_A - T_B$, then the CDF defined by:

$G_T(t) := P(T <=t)$

is what I need. The paper suggested taking inverse Laplace of a CDF to obtain the CDF and evaluate it at 0. It first suggested finding Laplace transform of T, which is given by

$L[T] = E(e^{-ST}) = -\frac{L[T_A](s) L[T_B](-s)}{s}$

Then it suggested taking Laplace of $G_T(t)$ i.e  $L[G_T(t)]$

However, $L[G_T(t)] = \frac{L[T]}{s} = \frac{L[T_A](s) L[T_B](-s)}{s}$

Then the paper suggests taking the inverse of above evaluated at 0 to get $P(T_A < T_B)$.

Question:

  1. Given the Laplace transform of CDF, which is $L[G_T(t)]$, I want to use the inverse laplace to obtain $G_T(t)$ evaluated at 0. But, by definition, inverse laplace using algorithms in python are all one sided from $[0,\infty]$ My random variable T is given by difference of two first passage times, $T = T_A - T_B $. Won't this be negative?

  2. In the paper, it says they are shifting the random variable X under study by a constant c such that P(X + c > 0) is approximately 1. Then inverting the corresponding one sided Laplace transform. How would I do that in this context here?

Thanks

I have the following problem. If I can get some help, I would greatly appreciate it. I am trying to replicate a particular research paper and came across a problem:

Suppose $X$ is a birth death process (represents population size) that evolves by:

$X \to X+1$ if a birth occurs with rate $\mu$,

$X \to X-1$ if a death occurs with rate $\theta$.

Suppose $T_A$ is first passage time of a BD process from state $A$ to state $0$ and suppose $T_B$ is first passage time of another BD process from state $B$ to state $0$. They are both independent.

I need to find $P(T_A < T_B)$. That is, probability that a population of size $A$ goes to $0$ before population of size $B$.

By definition:

$$T_A = T_{A,A-1} + T_{A-1,A-2} + \dots + T_{1,0}$$

where $T_{i,i-1}$ represents first passage time from state $i$ to state $i-1$.

I read some articles online that mentioned that if $T = T_A - T_B$, then the CDF defined by:

$$G_T(t) := P(T \leq t)$$

is what I need. The paper suggested taking inverse Laplace of a CDF to obtain the CDF and evaluate it at $0$. It first suggested finding Laplace transform of $T$, which is given by

$$L[T] = E(e^{-ST}) = -\frac{L[T_A](s) L[T_B](-s)}{s}.$$

Then it suggested taking Laplace of $G_T(t)$ i.e. $L[G_T(t)]$. However,

$$L[G_T(t)] = \frac{L[T]}{s} = \frac{L[T_A](s) L[T_B](-s)}{s}.$$

Then the paper suggests taking the inverse of the above evaluated at $0$ to get $P(T_A < T_B)$.

Questions:

  1. Given the Laplace transform of CDF, which is $L[G_T(t)]$, I want to use the inverse Laplace to obtain $G_T(t)$ evaluated at $0$. But, by definition, inverse Laplace using algorithms in python are all one sided from $[0,\infty]$. My random variable $T$ is given by difference of two first passage times, $T = T_A - T_B$. Won't this be negative?
  1. In the paper, it says they are shifting the random variable $X$ under study by a constant $c$ such that $P(X + c > 0)$ is approximately $1$. Then inverting the corresponding one sided Laplace transform. How would I do that in this context here?
Removed the (tag-removed) tag (The question has been bumped anyway - by a new answer.)
Source Link
Martin Sleziak
  • 4.7k
  • 4
  • 35
  • 40

I have the following problem. If i can get some help, I would greatly appreciate it. I am trying to replicate a particular research paper and came across a problem:

Suppose X is a birth death process (represents population size) that evolves by:

$X -> X+1 $ if a birth occurs with rate $\mu$

$X -> X-1 $ if a death occurs with rate $\theta$

Suppose $T_A$ is first passage time of a BD process from state A to state 0 and suppose $T_B$ is first passage time of another BD process from state B to state 0. They are both independent.

I need to find $P(T_A < T_B)$. That is, probability that a population of size A goes to 0 before population of size B.

By definition:

$T_A = T_{A,A-1} + T_{A-1,A-2} + ... + T_{1,0}$

where $T_{i,i-1}$ represents first passage time from state i, to state i-1.

I read some articles online that mentioned that if $T = T_A - T_B$, then the CDF defined by:

$G_T(t) := P(T <=t)$

is what iI need. The paper suggested taking inverse laplaceLaplace of a CDF to obtain the CDF and evaluate it at 0. It first suggested finding laplaceLaplace transform of T, which is given by

$L[T] = E(e^{-ST}) = -\frac{L[T_A](s) L[T_B](-s)}{s}$

Then it suggested taking laplaceLaplace of $G_T(t)$ i.e $L[G_T(t)]$

However, $L[G_T(t)] = \frac{L[T]}{s} = \frac{L[T_A](s) L[T_B](-s)}{s}$

Then the paper suggests taking the inverse of above evaluated at 0 to get $P(T_A < T_B)$.

Question:

  1. Given the Laplace transform of CDF, which is $L[G_T(t)]$, I want to use the inverse laplace to obtain $G_T(t)$ evaluated at 0. But, by definition, inverse laplace using algorithms in python are all one sided from $[0,\infty]$ My random variable T is given by difference of two first passage times, $T = T_A - T_B $. Won't this be negative?

  2. In the paper, it says they are shifting the random variable X under study by a constant c such that P(X + c > 0) is approximately 1. Then inverting the corresponding one sided Laplace transform. How would iI do that in this context here?

Thanks

I have the following problem. If i can get some help, I would greatly appreciate it. I am trying to replicate a particular research paper and came across a problem:

Suppose X is a birth death process (represents population size) that evolves by:

$X -> X+1 $ if a birth occurs with rate $\mu$

$X -> X-1 $ if a death occurs with rate $\theta$

Suppose $T_A$ is first passage time of a BD process from state A to state 0 and suppose $T_B$ is first passage time of another BD process from state B to state 0. They are both independent.

I need to find $P(T_A < T_B)$. That is, probability that a population of size A goes to 0 before population of size B.

By definition:

$T_A = T_{A,A-1} + T_{A-1,A-2} + ... + T_{1,0}$

where $T_{i,i-1}$ represents first passage time from state i, to state i-1.

I read some articles online that mentioned that if $T = T_A - T_B$, then the CDF defined by:

$G_T(t) := P(T <=t)$

is what i need. The paper suggested taking inverse laplace of a CDF to obtain the CDF and evaluate it at 0. It first suggested finding laplace transform of T, which is given by

$L[T] = E(e^{-ST}) = -\frac{L[T_A](s) L[T_B](-s)}{s}$

Then it suggested taking laplace of $G_T(t)$ i.e $L[G_T(t)]$

However, $L[G_T(t)] = \frac{L[T]}{s} = \frac{L[T_A](s) L[T_B](-s)}{s}$

Then the paper suggests taking the inverse of above evaluated at 0 to get $P(T_A < T_B)$.

Question:

  1. Given the Laplace transform of CDF, which is $L[G_T(t)]$, I want to use the inverse laplace to obtain $G_T(t)$ evaluated at 0. But, by definition, inverse laplace using algorithms in python are all one sided from $[0,\infty]$ My random variable T is given by difference of two first passage times, $T = T_A - T_B $. Won't this be negative?

  2. In the paper, it says they are shifting the random variable X under study by a constant c such that P(X + c > 0) is approximately 1. Then inverting the corresponding one sided Laplace transform. How would i do that in this context here?

Thanks

I have the following problem. If i can get some help, I would greatly appreciate it. I am trying to replicate a particular research paper and came across a problem:

Suppose X is a birth death process (represents population size) that evolves by:

$X -> X+1 $ if a birth occurs with rate $\mu$

$X -> X-1 $ if a death occurs with rate $\theta$

Suppose $T_A$ is first passage time of a BD process from state A to state 0 and suppose $T_B$ is first passage time of another BD process from state B to state 0. They are both independent.

I need to find $P(T_A < T_B)$. That is, probability that a population of size A goes to 0 before population of size B.

By definition:

$T_A = T_{A,A-1} + T_{A-1,A-2} + ... + T_{1,0}$

where $T_{i,i-1}$ represents first passage time from state i, to state i-1.

I read some articles online that mentioned that if $T = T_A - T_B$, then the CDF defined by:

$G_T(t) := P(T <=t)$

is what I need. The paper suggested taking inverse Laplace of a CDF to obtain the CDF and evaluate it at 0. It first suggested finding Laplace transform of T, which is given by

$L[T] = E(e^{-ST}) = -\frac{L[T_A](s) L[T_B](-s)}{s}$

Then it suggested taking Laplace of $G_T(t)$ i.e $L[G_T(t)]$

However, $L[G_T(t)] = \frac{L[T]}{s} = \frac{L[T_A](s) L[T_B](-s)}{s}$

Then the paper suggests taking the inverse of above evaluated at 0 to get $P(T_A < T_B)$.

Question:

  1. Given the Laplace transform of CDF, which is $L[G_T(t)]$, I want to use the inverse laplace to obtain $G_T(t)$ evaluated at 0. But, by definition, inverse laplace using algorithms in python are all one sided from $[0,\infty]$ My random variable T is given by difference of two first passage times, $T = T_A - T_B $. Won't this be negative?

  2. In the paper, it says they are shifting the random variable X under study by a constant c such that P(X + c > 0) is approximately 1. Then inverting the corresponding one sided Laplace transform. How would I do that in this context here?

Thanks

deleted 23 characters in body
Source Link

I have the following problem. If i can get some help, I would greatly appreciate it. I am trying to replicate a particular research paper and came across a problem:

Suppose X is a birth death process (represents population size) that evolves by:

$X -> X+1 $ if a birth occurs with rate $\mu$

$X -> X-1 $ if a death occurs with rate $\theta$

Suppose $T_A$ is first passage time of a BD process from state A to state 0 and suppose $T_B$ is first passage time of another BD process from state B to state 0. They are both independent.

I need to find $P(T_A < T_B)$. That is, probability that a population of size A goes to 0 before population of size B.

By definition:

$T_A = T_{A,A-1} + T_{A-1,A-2} + ... + T_{1,0}$

where $T_{i,i-1}$ represents first passage time from state i, to state i-1.

I read some articles online that mentioned that if $T = T_A - T_B$, then the CDF defined by:

$G_T(t) := P(T <=t)$

is what i need. The paper suggested taking inverse laplace of a CDF to obtain the CDF and evaluate it at 0. It first suggested finding laplace transform of T, which is given by

$L[T] = E(e^{-ST}) = -\frac{L[T_A](s) L[T_B](-s)}{s}$

Then it suggested taking laplace of $G_T(t)$ i.e $L[G_T(t)]$

However, $L[G_T(t)] = \frac{L[T]}{s} = \frac{L[T_A](s) L[T_B](-s)}{s}$

Then the paper suggests taking the inverse of above evaluated at 0 to get $P(T_A < T_B)$.

Question:

  1. Given the Laplace transform of CDF, which is $L[G_T(t)]$, I want to use the inverse laplace to obtain $G_T(t)$ evaluated at 0. But, by definition, inverse laplace using algorithms in python are all one sided from $[0,\infty]$ My random variable T is given by difference of two first passage times, $T = T_A - T_B $. Won't this be negative?

  2. In the paper, it says they are shifting the random variable X under study by a constant c such that P(X + c > 0) is approximately 1. Then inverting the corresponding one sided Laplace transform. They say choose c = 0.5 How How would i do that in this context here?

Thanks

I have the following problem. If i can get some help, I would greatly appreciate it. I am trying to replicate a particular research paper and came across a problem:

Suppose X is a birth death process (represents population size) that evolves by:

$X -> X+1 $ if a birth occurs with rate $\mu$

$X -> X-1 $ if a death occurs with rate $\theta$

Suppose $T_A$ is first passage time of a BD process from state A to state 0 and suppose $T_B$ is first passage time of another BD process from state B to state 0. They are both independent.

I need to find $P(T_A < T_B)$. That is, probability that a population of size A goes to 0 before population of size B.

By definition:

$T_A = T_{A,A-1} + T_{A-1,A-2} + ... + T_{1,0}$

where $T_{i,i-1}$ represents first passage time from state i, to state i-1.

I read some articles online that mentioned that if $T = T_A - T_B$, then the CDF defined by:

$G_T(t) := P(T <=t)$

is what i need. The paper suggested taking inverse laplace of a CDF to obtain the CDF and evaluate it at 0. It first suggested finding laplace transform of T, which is given by

$L[T] = E(e^{-ST}) = -\frac{L[T_A](s) L[T_B](-s)}{s}$

Then it suggested taking laplace of $G_T(t)$ i.e $L[G_T(t)]$

However, $L[G_T(t)] = \frac{L[T]}{s} = \frac{L[T_A](s) L[T_B](-s)}{s}$

Then the paper suggests taking the inverse of above evaluated at 0 to get $P(T_A < T_B)$.

Question:

  1. Given the Laplace transform of CDF, which is $L[G_T(t)]$, I want to use the inverse laplace to obtain $G_T(t)$ evaluated at 0. But, by definition, inverse laplace using algorithms in python are all one sided from $[0,\infty]$ My random variable T is given by difference of two first passage times, $T = T_A - T_B $. Won't this be negative?

  2. In the paper, it says they are shifting the random variable X under study by a constant c such that P(X + c > 0) is approximately 1. Then inverting the corresponding one sided Laplace transform. They say choose c = 0.5 How would i do that in this context here?

Thanks

I have the following problem. If i can get some help, I would greatly appreciate it. I am trying to replicate a particular research paper and came across a problem:

Suppose X is a birth death process (represents population size) that evolves by:

$X -> X+1 $ if a birth occurs with rate $\mu$

$X -> X-1 $ if a death occurs with rate $\theta$

Suppose $T_A$ is first passage time of a BD process from state A to state 0 and suppose $T_B$ is first passage time of another BD process from state B to state 0. They are both independent.

I need to find $P(T_A < T_B)$. That is, probability that a population of size A goes to 0 before population of size B.

By definition:

$T_A = T_{A,A-1} + T_{A-1,A-2} + ... + T_{1,0}$

where $T_{i,i-1}$ represents first passage time from state i, to state i-1.

I read some articles online that mentioned that if $T = T_A - T_B$, then the CDF defined by:

$G_T(t) := P(T <=t)$

is what i need. The paper suggested taking inverse laplace of a CDF to obtain the CDF and evaluate it at 0. It first suggested finding laplace transform of T, which is given by

$L[T] = E(e^{-ST}) = -\frac{L[T_A](s) L[T_B](-s)}{s}$

Then it suggested taking laplace of $G_T(t)$ i.e $L[G_T(t)]$

However, $L[G_T(t)] = \frac{L[T]}{s} = \frac{L[T_A](s) L[T_B](-s)}{s}$

Then the paper suggests taking the inverse of above evaluated at 0 to get $P(T_A < T_B)$.

Question:

  1. Given the Laplace transform of CDF, which is $L[G_T(t)]$, I want to use the inverse laplace to obtain $G_T(t)$ evaluated at 0. But, by definition, inverse laplace using algorithms in python are all one sided from $[0,\infty]$ My random variable T is given by difference of two first passage times, $T = T_A - T_B $. Won't this be negative?

  2. In the paper, it says they are shifting the random variable X under study by a constant c such that P(X + c > 0) is approximately 1. Then inverting the corresponding one sided Laplace transform. How would i do that in this context here?

Thanks

added 149 characters in body
Source Link
Loading
added 876 characters in body
Source Link
Loading
added 4 characters in body
Source Link
Loading
Source Link
Loading