using Distributions
using PDMats
using StatPlots
ENV["PLOTS_TEST"] = "true"
"true"
gr()
plot(Gamma(10, 0.1), label="Gamma(10, 0.1)")
A = Float64[
3 -2
-2 5
]
@show eig(A)
sig = PDMat(A)
mu = Float64[1, 3]
d = MvNormal(mu, sig)
eig(A) = ([1.76393, 6.23607], [-0.850651 -0.525731; -0.525731 0.850651])
FullNormal( dim: 2 μ: [1.0, 3.0] Σ: [3.0 -2.0; -2.0 5.0] )
gr()
x = collect(linspace(mu[1]-4, mu[1]+4, 201))
y = collect(linspace(mu[2]-4, mu[2]+4, 201))
f(x, y) = pdf(d, [x, y])
contour(x, y, f, size=(450,400))
gr()
x = collect(linspace(mu[1]-4, mu[1]+4, 201))
y = collect(linspace(mu[2]-4, mu[2]+4, 201))
f(x, y) = pdf(d, [x, y])
contourf(x, y, f, size=(450,400))
pyplot()
x = collect(linspace(mu[1]-4, mu[1]+4, 21))
y = collect(linspace(mu[2]-4, mu[2]+4, 21))
f(x, y) = pdf(d, [x, y])
wireframe(x, y, f)
pyplot()
x = collect(linspace(mu[1]-4, mu[1]+4, 201))
y = collect(linspace(mu[2]-4, mu[2]+4, 201))
f(x, y) = pdf(d, [x, y])
surface(x, y, f)
plotlyjs()
x = collect(linspace(mu[1]-4, mu[1]+4, 201))
y = collect(linspace(mu[2]-4, mu[2]+4, 201))
f(x, y) = pdf(d, [x, y])
surface(x, y, f)
Plotly javascript loaded.
To load again call
init_notebook(true)
pyplot()
x = collect(linspace(0,1,21))
y = collect(linspace(0,1,21))
f(x,y) = (0< x+y < 1 ? pdf(Dirichlet([5, 5, 5]), [x, y, 1-x-y]) : -1.0)
wireframe(x, y, f, camera=(60, 30), colorbar=false)
pyplot()
x = collect(linspace(0,1,201))
y = collect(linspace(0,1,201))
f(x,y) = (0< x+y < 1 ? pdf(Dirichlet([5, 5, 5]), [x, y, 1-x-y]) : -1.0)
surface(x, y, f, camera=(60, 30), colorbar=false)
pyplot()
x = collect(linspace(0,1,201))
y = collect(linspace(0,1,201))
f(x,y) = (0< x+y < 1 ? pdf(Dirichlet([5, 5, 5]), [x, y, 1-x-y]) : -1.0)
surface(x, y, f, camera=(180,60), colorbar=false)
pyplot()
n = 2^5
x = collect(0:n)
y = collect(0:n)
f(x,y) = (x+y ≤ n::Int ? pdf(Multinomial(n::Int, [0.3, 0.5, 0.2]), [x, y, n::Int - (x+y)]) : -0.005)
surface(x, y, f, camera=(180,60), colorbar=false)
plotlyjs()
n = 2^4
x = collect(0:n)
y = collect(0:n)
f(x,y) = (x+y ≤ n::Int ? pdf(Multinomial(n::Int, [0.3, 0.5, 0.2]), [x, y, n::Int - (x+y)]) : -0.005)
surface(x, y, f, camera=(180,60), colorbar=false)