Plotando uma amostra gaussiana 2D

Question

Você pode me ajudar a fazer uma amostra Guassiana 2D com médias e variações especificadas?

Eu só sei fazer uma curva guassiana :D

\documentclass{article}
\usepackage{paralist,pst-func, pst-plot, pst-math, pstricks-add,pgfplots}
\usetikzlibrary{patterns,matrix,arrows}
\begin{document}
\begin{tikzpicture}
    \begin{axis}[hide axis,clip=false,xmin=-4,xmax=4,xlabel={X},ymin=0,ymax=1] 
    \addplot[color=lime, samples=100] {1/sqrt(2)*exp(-(x+1)^2/1)} ;
    \end{axis}
\end{tikzpicture}
\end{document}

E aqui está a fórmula para isso:

E este seria um exemplo de três distribuições normais juntas:

Answer 1

Para isso, você precisa de uma forma de gerar números aleatórios normalmente distribuídos. Uma maneira de fazer isso é usar a transformada Box-Muller.

Aqui está um exemplo de uso de PGFPlots para isso (com base na minha resposta paraTikZ: Desenhando os mesmos dados com gráficos de dispersão e coordenadas paralelas). Plotei as distribuições marginais para mostrar que os números são de fato aproximadamente normais:

\documentclass{article}
\usepackage{pgfplots}
% Create a function for generating inverse normally distributed numbers using the Box–Muller transform
\pgfmathdeclarefunction{invgauss}{2}{%
  \pgfmathparse{sqrt(-2*ln(#1))*cos(deg(2*pi*#2))}%
}
\begin{document}
\begin{tikzpicture}
\begin{axis}[
    axis equal image,
    xmin=-2.5,xmax=2.5,
    ymin=-2.5, ymax=2.5,
    enlargelimits=false,
    xtick=data, 
    xticklabel=\empty,
    ytick=data,
    yticklabel=\empty,
    extra x ticks={-2,...,2},
    every extra x tick/.style={
        tick align=outside,
        xticklabel=\pgfmathprintnumber{\tick}
    },
    extra y ticks={-2,...,2},
    every extra y tick/.style={
        tick align=outside,
        yticklabel=\pgfmathprintnumber{\tick}
    }
]
\addplot [only marks, samples=100] ({invgauss(rnd,rnd)},{invgauss(rnd,rnd)});
\end{axis}
\end{tikzpicture}
\end{document}