pmf

Author: bloodberrys

binom_dist.m

@@ -0,0 +1,7 @@
+function pmf = binom_dist(N,p,k)
+  nValues = numel(k);
+  pmf = zeros(1,nValues);
+  for i = 1:nValues
+    pmf(i) = nchoosek(N,k(i))*p^k(i)*(1-p)^(N-k(i));
+  end
+end

newpmf.m

@@ -0,0 +1,16 @@
+
+format longG
+
+data = readtable('/Users/mymac/Documents/SCRIPTSHEET/SKRIPSI/data_jmeter/4/LatenciesOverTime.csv', 'ReadVariableNames', false, 'HeaderLines', 1);
+% x = data{:,2};
+ x = -1000:0.5:1000;
+ B = data{:,2}';
+ B = B(:)';
+ %x = B;
+h=x/sum(x);
+ %y = randn(10000,1);
+y = h(:);
+disp(x)
+n = hist(y);
+pmf = n/sum(n);
+plot(pmf,'o');

pmf.m

@@ -1,47 +1,27 @@
-format longG
-
-data = readtable('/Users/mymac/Documents/SCRIPTSHEET/SKRIPSI/data_jmeter/4/LatenciesOverTime.csv', 'ReadVariableNames', false, 'HeaderLines', 1);
-
-% ini untuk format date seperti halnya di jmeter ya
-% x = seconds((datenum(datestr(data.ElapsedTime, 'yyyy-mm-dd hh:MM:ss.fff')) - datenum(datestr(data{1,1}, 'yyyy-mm-dd hh:MM:ss.fff'))) * 100000);
-% x.Format = 'hh:mm:ss';
-
-markers = {'+','*','.','o','x','v','d','^','s','>','<','v','p','h','p','v','<','>','s','^','d','v','x','o','.','*'};
-colors = {'r','b','m','k','y','c','g','c','y','k','m','b'};
-lines = {'-','--',':','-.',':','--'};
-line_width = 0.9;
-marker_size = 5;
-marker_counter = 1;
-color_counter = 1;
-line_counter = 1;
-ylabels='Response Time (ms)';
-ylabels1='Number of Threads';
-ylabels2='Mega Bytes (MB)';
-ylabels3='TCP Connections';
-legend_base_name = 'data-';
-
-% kalo ini format number dalam menit elapsed time
-% x = (datenum(datestr(data{:,1}, 'yyyy-mm-dd hh:MM:ss.fff')) - datenum(datestr(data{1,1}, 'yyyy-mm-dd hh:MM:ss.fff'))) * 100000/60;
-
-% ini untuk data x yang pure number aja, bukan elapsed time
-x = data{:,2};
-
-hold on
-title('Any Plot');
-xlabel('Elapsed time (minutes), Granulation: 500 ms','FontSize',12,'FontWeight','bold');
-ylabel(ylabels3,'FontSize',12,'FontWeight','bold');
-disp(data)
-pmfs=zeros(size(x(:))); 
-    for i=1:length(x) 
-        pmfs(i)= sum(px(find(sx==x(i)))); 
+function [vals, freqs] = pmf(X)
+%#PMF Return the probability mass function for a vector/matrix.
+%#
+%#INPUTS:
+%#   X       Input matrix
+%#
+%#OUTPUTS:
+%#   VALS    Vector of unique values
+%#   FREQS   Vector of frequencies of occurence of each value.
+%#
+
+    [vals,junk,idx] = unique(X);
+
+    vals   = vals(:);
+    frequs = NaN(length(vals),1);
+
+    for i = 1:length(vals)
+        freqs(i) = mean(idx == i);
     end
-% pmfs = histcounts(data{:,2},[unique(data{:,2}) Inf],'Normalization','probability');
-histogram(pmfs);
-
 
+    %# If 0 or 1 output is requested, put the values and counts in two columns
+    %# of a matrix.
+    if nargout < 2
+        vals = [vals freqs];
+    end
 
-box on;
-ax = gca;
-ax.YAxis.Exponent = 0;
-legend('show');
-hold off
+end

pmf1.m

@@ -0,0 +1,51 @@
+format longG
+
+data = readtable('/Users/mymac/Documents/SCRIPTSHEET/SKRIPSI/data_jmeter/4/LatenciesOverTime.csv', 'ReadVariableNames', false, 'HeaderLines', 1);
+
+% ini untuk format date seperti halnya di jmeter ya
+% x = seconds((datenum(datestr(data.ElapsedTime, 'yyyy-mm-dd hh:MM:ss.fff')) - datenum(datestr(data{1,1}, 'yyyy-mm-dd hh:MM:ss.fff'))) * 100000);
+% x.Format = 'hh:mm:ss';
+
+markers = {'+','*','.','o','x','v','d','^','s','>','<','v','p','h','p','v','<','>','s','^','d','v','x','o','.','*'};
+colors = {'r','b','m','k','y','c','g','c','y','k','m','b'};
+lines = {'-','--',':','-.',':','--'};
+line_width = 0.9;
+marker_size = 5;
+marker_counter = 1;
+color_counter = 1;
+line_counter = 1;
+ylabels='Response Time (ms)';
+ylabels1='Number of Threads';
+ylabels2='Mega Bytes (MB)';
+ylabels3='Probability';
+legend_base_name = 'data-';
+
+% kalo ini format number dalam menit elapsed time
+% x = (datenum(datestr(data{:,1}, 'yyyy-mm-dd hh:MM:ss.fff')) - datenum(datestr(data{1,1}, 'yyyy-mm-dd hh:MM:ss.fff'))) * 100000/60;
+
+% ini untuk data x yang pure number aja, bukan elapsed time
+x = fix(data{:,2});
+
+hold on
+title('Any Plot');
+xlabel('Delay','FontSize',12,'FontWeight','bold');
+ylabel(ylabels3,'FontSize',12,'FontWeight','bold');
+
+ %code
+  b=x/sum(x);
+%  b = b';
+%  b = b(:)';
+    
+
+B = x';
+B = B(:)';
+
+disp(sum(b))
+
+% stem(B,b)
+
+box on;
+ax = gca;
+ax.YAxis.Exponent = 0;
+legend('show');
+hold off

pmf1.m~

@@ -0,0 +1,63 @@
+format longG
+
+data = readtable('/Users/mymac/Documents/SCRIPTSHEET/SKRIPSI/data_jmeter/4/LatenciesOverTime.csv', 'ReadVariableNames', false, 'HeaderLines', 1);
+
+% ini untuk format date seperti halnya di jmeter ya
+% x = seconds((datenum(datestr(data.ElapsedTime, 'yyyy-mm-dd hh:MM:ss.fff')) - datenum(datestr(data{1,1}, 'yyyy-mm-dd hh:MM:ss.fff'))) * 100000);
+% x.Format = 'hh:mm:ss';
+
+markers = {'+','*','.','o','x','v','d','^','s','>','<','v','p','h','p','v','<','>','s','^','d','v','x','o','.','*'};
+colors = {'r','b','m','k','y','c','g','c','y','k','m','b'};
+lines = {'-','--',':','-.',':','--'};
+line_width = 0.9;
+marker_size = 5;
+marker_counter = 1;
+color_counter = 1;
+line_counter = 1;
+ylabels='Response Time (ms)';
+ylabels1='Number of Threads';
+ylabels2='Mega Bytes (MB)';
+ylabels3='Probability';
+legend_base_name = 'data-';
+
+% kalo ini format number dalam menit elapsed time
+% x = (datenum(datestr(data{:,1}, 'yyyy-mm-dd hh:MM:ss.fff')) - datenum(datestr(data{1,1}, 'yyyy-mm-dd hh:MM:ss.fff'))) * 100000/60;
+
+% ini untuk data x yang pure number aja, bukan elapsed time
+x = data{:,2};
+
+hold on
+title('Any Plot');
+xlabel('Delay','FontSize',12,'FontWeight','bold');
+ylabel(ylabels3,'FontSize',12,'FontWeight','bold');
+
+ %code
+  b=x/sum(x);
+  c=[x b];
+  
+
+ 
+B = data{:,2}';
+B = B(:)';
+
+ pmfs = accumarray(x(:),1)./numel(B);
+ disp(pmfs)
+
+% plot(B,b);
+  
+ 
+
+% [res,a] = pmf(x);
+% pmfa = accumarray(x(:),1)./numel(x);
+% pmfs = pmfa./sum(pmfa);
+% disp(pmfs)
+% pmfs = histcounts(data{:,2},[unique(data{:,2}) Inf],'Normalization','probability');
+% stem(tes,y);
+% size(pmfs);
+
+
+box on;
+ax = gca;
+ax.YAxis.Exponent = 0;
+legend('show');
+hold off