Hi meshesurf
As per my understanding, you want to find the plot for the joint probability density functions.
Here is an example code that you can modify according to your need:
X = [1;2;1;1;2;1;1];
Y = [1;1;2;2;1;1;2];
%First find all list of values in X and Y seperately
% In this case I have taken only 1,2 as the possible values,
x_vals = [1,2];
y_vals = [1,2];
values={x_vals y_vals};
count=hist3([X,Y],'ctrs',values);%first count the values
p=count/sum(sum(count)); % then divide by the sum on all matrix value to get probability
%%
imagesc(p);
