function gpdisc(call)
% This is the machine-generated representation of a MATLAB object
% and its children. Note that handle values may change when these
% objects are re-created. This may cause problems with some callbacks.
% The command syntax may be supported in the future, but is currently
% incomplete and subject to change.
%
% To re-open this system, just type the name of the m-file at the MATLAB
% prompt. The M-file and its associtated MAT-file must be on your path.
if nargin==0, call='start'; end
switch call
case 'start'
load gpdisc
a = figure('Color',[0.8 0.8 0.8], ...
'Colormap',mat0, ...
'Name','Learning about p using discrete models', ...
'PointerShapeCData',mat1, ...
'Position',[-2 30 1028 702], ...
'Tag','Fig1');
b = axes('Parent',a, ...
'Box','on', ...
'ButtonDownFcn','gpdisc(''adjust_bars'')', ...
'CameraUpVector',[0 1 0], ...
'CameraUpVectorMode','manual', ...
'Color',[1 1 1], ...
'ColorOrder',mat2, ...
'NextPlot','add', ...
'Position',[0.116803 0.697039 0.727459 0.15262], ...
'Tag','prior plot', ...
'UserData',mat3, ...
'XColor',[0 0 0], ...
'XLim',[-0.1 1.1], ...
'XLimMode','manual', ...
'YColor',[0 0 0], ...
'YLim',[0 0.181818], ...
'YLimMode','manual', ...
'ZColor',[0 0 0]);
c = text('Parent',b, ...
'Color',[0 0 0], ...
'HandleVisibility','callback', ...
'HorizontalAlignment','center', ...
'Position',[0.499196 -0.0454545 0], ...
'String','p', ...
'Tag','Text1', ...
'VerticalAlignment','cap');
set(get(c,'Parent'),'XLabel',c);
c = text('Parent',b, ...
'Color',[0 0 0], ...
'HandleVisibility','callback', ...
'HorizontalAlignment','center', ...
'Position',[-0.169973 0.0917666 0], ...
'Rotation',90, ...
'String','prob', ...
'Tag','Text2', ...
'VerticalAlignment','baseline');
set(get(c,'Parent'),'YLabel',c);
c = text('Parent',b, ...
'Color',[0 0 0], ...
'HandleVisibility','callback', ...
'HorizontalAlignment','right', ...
'Position',[-0.293834 0.362778 0], ...
'Tag','Text3', ...
'Visible','off');
set(get(c,'Parent'),'ZLabel',c);
c = text('Parent',b, ...
'Color',[0 0 0], ...
'HandleVisibility','callback', ...
'HorizontalAlignment','center', ...
'Position',[0.499196 0.198113 0], ...
'String','PRIOR DISTRIBUTION', ...
'Tag','Text4', ...
'VerticalAlignment','bottom');
set(get(c,'Parent'),'Title',c);
c = line('Parent',b, ...
'Color',[0 0 1], ...
'Tag','pbars', ...
'UserData',mat4, ...
'XData',mat5, ...
'YData',mat6);
b = axes('Parent',a, ...
'Units','points', ...
'Box','on', ...
'CameraUpVector',[0 1 0], ...
'CameraUpVectorMode','manual', ...
'Color',[1 1 1], ...
'ColorOrder',mat7, ...
'Position',[70.7586 129.724 466.759 75.7241], ...
'Tag','post plot', ...
'UserData',mat8, ...
'XColor',[0 0 0], ...
'XLim',[-0.1 1.1], ...
'XLimMode','manual', ...
'YColor',[0 0 0], ...
'YLim',[0 0.255409], ...
'YLimMode','manual', ...
'ZColor',[0 0 0]);
c = line('Parent',b, ...
'Color',[0 0 1], ...
'Tag','post bars', ...
'UserData',0, ...
'XData',mat9, ...
'YData',mat10);
c = text('Parent',b, ...
'Color',[0 0 0], ...
'HandleVisibility','callback', ...
'HorizontalAlignment','center', ...
'Position',[0.5 -0.0559367 0], ...
'String','p', ...
'Tag','Text5', ...
'VerticalAlignment','cap');
set(get(c,'Parent'),'XLabel',c);
c = text('Parent',b, ...
'Color',[0 0 0], ...
'HandleVisibility','callback', ...
'HorizontalAlignment','center', ...
'Position',[-0.155126 0.127704 0], ...
'Rotation',90, ...
'String','prob', ...
'Tag','Text6', ...
'VerticalAlignment','baseline');
set(get(c,'Parent'),'YLabel',c);
c = text('Parent',b, ...
'Color',[0 0 0], ...
'HandleVisibility','callback', ...
'HorizontalAlignment','right', ...
'Position',[-0.281358 1.03958 0], ...
'Tag','Text7', ...
'Visible','off');
set(get(c,'Parent'),'ZLabel',c);
c = text('Parent',b, ...
'Color',[0 0 0], ...
'HandleVisibility','callback', ...
'HorizontalAlignment','center', ...
'Position',[0.5 0.275462 0], ...
'String','POSTERIOR DISTRIBUTION', ...
'Tag','Text8', ...
'VerticalAlignment','bottom');
set(get(c,'Parent'),'Title',c);
b = uicontrol('Parent',a, ...
'Units','points', ...
'Callback','gpdisc(''get_data'')', ...
'Position',[127.241 240.828 49.0345 19.2414], ...
'String','Get Data', ...
'Tag','get data', ...
'UserData',[0 0]);
b = uicontrol('Parent',a, ...
'Units','points', ...
'BackgroundColor',[0.752941 0.752941 0.752941], ...
'Position',[266.276 243.931 54 15.5172], ...
'String','s=0 f=0', ...
'Style','text', ...
'Tag','text_data', ...
'UserData',[1 1]);
b = uicontrol('Parent',a, ...
'Units','points', ...
'Callback','gpdisc(''update'')', ...
'Position',[385.448 241.448 50.2759 19.8621], ...
'String','Update', ...
'Tag','Pushbutton1', ...
'UserData',[23 100]);
b = uicontrol('Parent',a, ...
'Units','points', ...
'Callback','gpdisc(''define_prior'')', ...
'Position',[248.897 399.724 107.379 16.1379], ...
'String','Define Prior Grid', ...
'Tag','define prior', ...
'UserData',[0 1 11]);
b = uicontrol('Parent',a, ...
'Units','points', ...
'BackgroundColor',[0.752941 0.752941 0.752941], ...
'Position',[369.931 36 89.3793 50.8966], ...
'Style','text', ...
'Tag','summaries');
b = uicontrol('Parent',a, ...
'Units','points', ...
'BackgroundColor',[1 1 1], ...
'Position',[102.414 33.5172 137.793 54], ...
'Style','listbox', ...
'Tag','list probs', ...
'Value',1);
b = uicontrol('Parent',a, ...
'Units','points', ...
'BackgroundColor',[0.752941 0.752941 0.752941], ...
'Position',[119.793 92.4828 94.9655 12.4138], ...
'String','PROBABILITY TABLE', ...
'Style','text', ...
'Tag','StaticText1');
b = uicontrol('Parent',a, ...
'Units','points', ...
'Position',[372.414 93.7241 81.3103 11.1724], ...
'String','SUMMARIES', ...
'Style','text', ...
'Tag','StaticText1');
case 'define_prior'
prompt={'Enter smallest value for p::','Enter largest value for p:',...
'Enter number of values:'};
def={'','',''};
title='Define prior grid';
lineNo=1;
answer=inputdlg(prompt,title,lineNo,def);
lo=eval(answer{1}); hi=eval(answer{2}); n=eval(answer{3}); w=hi-lo;
set(gco,'userdata',[lo hi n]);
p=linspace(lo,hi,n);
prior=1/n+0*p;
t=findobj('tag','pbars');
X=[p;p;p]; Y=[0*prior;prior;0*prior+NaN];
set(t,'ydata',Y(:));
set(t,'xdata',X(:));
t=findobj('tag','prior plot');
set(t,'xlim',[lo-w/(n-1) hi+w/(n-1)]);
set(t,'ylim',[0 2*prior(1)]);
set(t,'userdata',[p;prior]);
t=findobj('tag','post bars');
t2=findobj('tag','post plot');
pdist=get(t2,'userdata');
p=pdist(1,:);
y=[0*p;0*p;NaN+0*p];
set(t,'ydata',y(:))
set(t2,'userdata',[p;0*p]);
set(findobj('tag','summaries'),'string','')
set(findobj('tag','list probs'),'string','')
case 'get_data'
prompt={'Enter the number of successes:','Enter the number of failures:'};
def={'',''};
title='Input Binomial Data';
lineNo=1;
answer=inputdlg(prompt,title,lineNo,def);
if prod(size(answer))==2
if (length(answer{1})>0)&(length(answer{2})>0)==1
set(gco,'userdata',[eval(answer{1}),eval(answer{2})])
t=findobj('tag','text_data');
set(t,'string',['s=' answer{1} ' f=' answer{2}])
t=findobj('tag','post bars');
flag=get(t,'userdata');
if flag==1
set(t,'userdata',0)
t2=findobj('tag','post plot');
pdist=get(t2,'userdata');
p=pdist(1,:);
y=[0*p;0*p;NaN+0*p];
set(t,'ydata',y(:))
set(t2,'userdata',[p;0*p]);
set(findobj('tag','summaries'),'string','')
set(findobj('tag','list probs'),'string','')
end
end
end
case 'update'
t=findobj('tag','prior plot');
pdist=get(t,'userdata');
t=findobj('tag','get data');
data=get(t,'userdata');
post=p_disc(pdist(1,:),pdist(2,:),data);
x=ones(3,1)*pdist(1,:);
y=[0*post;post;NaN+0*post;];
t=findobj('tag','post bars');
set(t,'xdata',x(:))
set(t,'ydata',y(:))
set(t,'userdata',1)
t=findobj('tag','post plot');
t2=findobj('tag','define prior');
lohi=get(t2,'userdata');
lo=lohi(1); hi=lohi(2); w=hi-lo; n=lohi(3);
set(t,'xlim',[lo-w/(n-1) hi+w/(n-1)])
set(t,'ylim',[0 max(post)*1.1])
postdist=[pdist(1,:);post];
set(t,'userdata',postdist);
clear
t=findobj('tag','post plot');
S=c_summ(get(t,'userdata'));
set(findobj('tag','summaries'),'string',S)
set(findobj('tag','list probs'),'string', ...
num2str(get(t,'userdata')','%5.3f %5.3f'));
case 'adjust_bars'
t=findobj('tag','define prior');
lohi=get(t,'userdata');
lo=lohi(1); hi=lohi(2); n=lohi(3);
p=get(gca,'currentpoint'); % reads point at mouseclick
dx=(hi-lo)/(n-1);
x=round((p(1,1)-lo)/dx+1); % finds the closest value of p
y=p(1,2);
t=findobj('tag','prior plot');
pdist=get(t,'userdata');
pr=pdist(2,:);
pr(x)=y; % put new value of p in probability vector
pr=pr/sum(pr); % normalizes
Y=[0*pr;pr;0*pr+NaN];
prior_handle=findobj('tag','pbars');
set(t,'ylim',[0 max(pr)*1.2])
set(prior_handle,'ydata',Y(:)); % resets current ydata of line
pdist=get(t,'userdata');
p=pdist(1,:); pdist=[p;pr];
set(gca,'userdata',pdist);
t=findobj('tag','post bars');
flag=get(t,'userdata');
if flag==1
set(t,'userdata',0)
x=[p;p;NaN+p];
y=[0*p;0*p;NaN+0*p];
set(t,'xdata',x(:))
set(t,'ydata',y(:))
t=findobj('tag','post plot');
set(t,'userdata',[pdist(1,:);0*p]);
set(findobj('tag','summaries'),'string','')
set(findobj('tag','list probs'),'string','')
end
end
function S=c_summ(postdist)
mn=sum(postdist(1,:).*postdist(2,:));
sd=sqrt(sum(postdist(1,:).^2.*postdist(2,:))-mn^2);
set=disc_int(postdist',.9); lo=min(set); hi=max(set);
t=findobj('tag','summaries');
S=strvcat(['Mean = ' num2str(mn)],...
['St dev = ' num2str(sd)], ...
['90% interval: (' num2str(lo) ',' num2str(hi) ')']);
