#C7.txt: ExpMap class Help:=proc(): print(` PostD(Pri,t,a,b,Lik,k,Data1)`): print(`PostD(t^(3-1)*(1-t)^(3-1),t,0,1,t^k*(1-t)^(1-k),k,[0,1,0,1,0]);`): ptint(`dgamma(t,a,b)`): end: #PostD(Pri,t,a,b,Lik,k,Data1): inputs a CONT. pdf Pri in the variable t, #for the Prior distribution, whose domain is t=a to t=b, a #DISCRETE r.v. Lik in terms of k, and a LIST of the data, Data1, uses #the Bayesian methodology to output the PDF of the POSTERIOR distribution #For example, #PostD(t^(3-1)*(1-t)^(3-1),t,0,1,t^k*(1-t)^(1-k),k,[0,1,0,1,0]); PostD:=proc(Pri,t,a,b,Lik,k,Data1) local Pri1,Post1,i: Pri1:=Pri/int(Pri,t=a..b): Post1:=Pri1*mul( subs(k=Data1[i],Lik),i=1..nops(Data1)): Post1:=simplify(Post1/int(Post1,t=a..b)): Post1: end: #dgamma(t,a,b): The Gamma distribution dgamma:=proc(a,b,t) local A: A:=b^a*t^(a-1)*exp(-b*t)/GAMMA(a): end: