Advertisement
Guest User

Untitled

a guest
Jul 21st, 2017
62
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.37 KB | None | 0 0
  1. t := 1
  2. logQlb[a_] := -a^2/2 + Log[a] - Log[1 + a^2] - Log[Sqrt[2*[Pi]]]
  3. logQub[a_] := -a^2/2 - Log[a] - Log[Sqrt[2*[Pi]]]
  4. ratefn[t_] := t^2/2
  5. Plot[{-n*ratefn[t], logQlb[t*Sqrt[n]], logQub[t*Sqrt[n]]}, {n, 1,
  6. 100}, PlotLegends -> {"Large Deviations (Cramer) estimate",
  7. "Exact (lower bound) via Q fn", "Exact (upper bound) via Q fn"},
  8. AxesLabel -> {"n", "Log[P[Sn>=t]]"}]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement