Spikegraph:=proc(distributionlist,xmin,xmax,show,c)
local num,linelist,j,i,templist:
templist:=[]:
if ((distributionlist[1][1]<xmin) or
(distributionlist[nops(distributionlist)][1]>xmax)
) then
lprint(`Note: some data values lie outside the user-defined interval.`)
fi:
for i from 1 to nops(distributionlist) do
if not((distributionlist[i][1]<xmin) or (distributionlist[i][1]>xmax)) then
templist:=[op(templist),
[distributionlist[i][1],
distributionlist[i][2]
]
]:
fi:
od:
num:=nops(templist):
linelist:=[]:
for j from 1 to num do
linelist:=[op(linelist),[templist[j][1],0]]:
linelist:=[op(linelist),[templist[j][1],templist[j][2]]]:
linelist:=[op(linelist),[templist[j][1],0]]:
od:
if show then plot(linelist,style=LINE,color=c) fi;
end:
Bin:=proc(n,p,k)
(binomial(n,k))*(p^k)*((1-p)^(n-k));
end:
with(plots):
"This programs computes for us the probabilty that our result is due to chance error givne the null hypothesis is true. Its
ofeten called a significane measure. ":
BinomialPlot:=proc(n,m,p)
local length,i,j1,j2,j3,j4,j5,spikelist1,spikelist2,spikelist3,
spikelist4,spikelist5,A1,A2,A3,A4,A5:
spikelist1:=[]:
spikelist2:=[]:
spikelist3:=[]:
spikelist4:=[]:
spikelist5:=[]:
if m <= p*n then
for j2 from 0 to n
do
spikelist2:=[op(spikelist2),Bin(n,p,j2)] ;
od:
for j1 from 0 to m
do
spikelist1:=[op(spikelist1),[j1,Bin(n,p,j1)]] ;
od:
for j5 from m+1 to n do
spikelist5:=[op(spikelist5),[j5,Bin(n,p,j5)]] ;
od:
A1:= Spikegraph(spikelist1,0,m,true,red);
A5:=Spikegraph(spikelist5,m+1,n,true,black);
A2:=Spikegraph([[m,max(op(spikelist2))]],m,m,true,blue);
fi;
if m > p*n then
for j2 from 0 to n
do
spikelist2:=[op(spikelist2),Bin(n,p,j2)] ;
od:
for j1 from 0 to m-1
do
spikelist1:=[op(spikelist1),[j1,Bin(n,p,j1)]] ;
od:
for j5 from m to n do
spikelist5:=[op(spikelist5),[j5,Bin(n,p,j5)]] ;
od:
A1:= Spikegraph(spikelist1,0,m-1,true,black);
A5:=Spikegraph(spikelist5,m,n,true,red);
A2:=Spikegraph([[m,max(op(spikelist2))]],m,m,true,blue);
fi;
display(A2,A1,A5);
end:
Hyp:=proc(n,m,p)
local r,s;
r:=sum(
binomial(n,k)*(p^k)*((1-p)^(n-k)),k=m...n
);
s:=sum(
binomial(n,k)*(p^k)*((1-p)^(n-k)),k=0...m
);
if m > p*n then
lprint("their is a",100*r,"percent chance that we would see this or something even more dramatic due to chance alone.");
fi;
if m <= p*n then
lprint("their is a",100*s,"percent chance that we would see this or something even more dramatic due to chance alone.");
fi;
lprint("In the following picture the blue line represents the cut off value and the red region indicates is the more dramatic possibilities.");
BinomialPlot(n,m,p);
end:Hyp(58,28,.375);