使用 pgfplots 绘制理论二项分布

Question 1

将@John Kormylos 和 @cgnieders 的评论放在一起，我得到以下解决方案，ybar interval而不是hist：

\documentclass{article}
\usepackage{pgfplots}
\usepackage{python}
\begin{document}


\begin{python}%
import numpy as np
from scipy import special
N,p=150,.5
def f(n):
    return special.binom(N,n)*p**n*(1-p)**(N-n)
s = zip(range(1,N),map(f,range(1,N)))
np.savetxt("bernoulli.dat",s,fmt='%0.5f')
\end{python}%


\begin{tikzpicture}
  \begin{axis}[width=10cm]
    \addplot[%
    ybar interval,
    fill=blue!30,
    draw opacity=0.5,
    line width=0.4pt,
    ]%
    table {bernoulli.dat};
  \end{axis}
\end{tikzpicture}

\end{document}

输出：

Answer

将@John Kormylos 和 @cgnieders 的评论放在一起，我得到以下解决方案，ybar interval而不是hist：

\documentclass{article}
\usepackage{pgfplots}
\usepackage{python}
\begin{document}


\begin{python}%
import numpy as np
from scipy import special
N,p=150,.5
def f(n):
    return special.binom(N,n)*p**n*(1-p)**(N-n)
s = zip(range(1,N),map(f,range(1,N)))
np.savetxt("bernoulli.dat",s,fmt='%0.5f')
\end{python}%


\begin{tikzpicture}
  \begin{axis}[width=10cm]
    \addplot[%
    ybar interval,
    fill=blue!30,
    draw opacity=0.5,
    line width=0.4pt,
    ]%
    table {bernoulli.dat};
  \end{axis}
\end{tikzpicture}

\end{document}

输出：

Question 2

问题已由 Anna 解决，但您可能有兴趣知道，使用sagetex包可以得到相同的结果，而无需在此过程中创建数据文件。Sage 为您提供 Python 以及对数学有用的额外命令。

\documentclass{article}
\usepackage{sagetex}
\usepackage{pgfplots}
\pagestyle{empty}
\begin{document}
\begin{sagesilent}
p = .5
N = 140

output =""
output += r"\begin{tikzpicture}[scale=1.0]"
output += r"\begin{axis}[width=10cm]"
output += r"\addplot[ybar interval,fill=blue!30,draw opacity=0.5,line   width=0.4pt] coordinates{"
for i in range(0,N):
    output += r"(%d, %f) "%(i, binomial(N,i)*p^i*(1-p)^(N-i))

output += r"};"
output += r"\end{axis}"
output += r"\end{tikzpicture}"
\end{sagesilent}
\sagestr{output}
\end{document}

无需数据文件即可运行此输出萨基马云：在此处输入图片描述

Answer

问题已由 Anna 解决，但您可能有兴趣知道，使用sagetex包可以得到相同的结果，而无需在此过程中创建数据文件。Sage 为您提供 Python 以及对数学有用的额外命令。

\documentclass{article}
\usepackage{sagetex}
\usepackage{pgfplots}
\pagestyle{empty}
\begin{document}
\begin{sagesilent}
p = .5
N = 140

output =""
output += r"\begin{tikzpicture}[scale=1.0]"
output += r"\begin{axis}[width=10cm]"
output += r"\addplot[ybar interval,fill=blue!30,draw opacity=0.5,line   width=0.4pt] coordinates{"
for i in range(0,N):
    output += r"(%d, %f) "%(i, binomial(N,i)*p^i*(1-p)^(N-i))

output += r"};"
output += r"\end{axis}"
output += r"\end{tikzpicture}"
\end{sagesilent}
\sagestr{output}
\end{document}

无需数据文件即可运行此输出萨基马云：在此处输入图片描述

使用 pgfplots 绘制理论二项分布

答案1

答案2

相关内容