MaplePrimes - answers and comments on Question, Beta Prime Distribution in Maple

Assigning parameters

Markiyan Hirnyk — Tue, 26 Feb 2013 21:27:39 Z

This works.

with(Statistics):
U := Distribution(PDF = (proc (x) options operator, arrow; x^(alpha-1)*(1+x)^(-alpha-beta)/Beta(alpha, beta) end proc), Support = 0 .. infinity):
Z := RandomVariable(U):
alpha := 2: beta := 1.5:
DensityPlot(Z, range = 0 .. infinity);

You can also follow http://www.mapleprimes.com/questions/124046-Truncated-Normal-Distribution#comment124049

when creating U.

By Explore

Markiyan Hirnyk — Tue, 26 Feb 2013 22:24:19 Z

Also try

See ?Explore for more info.

Explore

pagan — Tue, 26 Feb 2013 23:19:47 Z

@Markiyan Hirnyk It might be called in Maple 16 like,

Explore( proc(a,b)
           Statistics:-DensityPlot(
             Statistics:-RandomVariable(
               Statistics:-Distribution(
                 PDF = (x-> x^(a-1)*(1+x)^(-a-b)/Beta(a, b)),
                 Support = 0 .. infinity)), range=0..infinity);
         end proc(alpha,beta) );

And then in the pop-up window adjust (at least some lower) end-values. But the Slider components that will appear in a new worksheet will reference names whose assignments are in a hidden document block, so it's not easy to cut and paste it as a whole back into the original open worksheet.

It might be easier to build an interactive version in the currently opened worksheet by following the few manual steps outlined here. The call to `plot` might then be replaced by your DensityPlot(Z,...) call above, for example.

Thank you for the valuable comment

Markiyan Hirnyk — Tue, 26 Feb 2013 23:23:44 Z

I am glad to see you in MaplePrimes.