Search Comments and Ratings

View All Comments and Ratings by:
Most Recent
Rank

Comments and Ratings

E-mail this View
			1 - 4 of 4
Date	File	Comment by	Comment	Rating
11 Sep 2009	RANDP Random integers with given probabilities ( (v2.0) Author: Jos	Pozzi, Francesco	I've tried this small experiment: N = 100000; % sample size w = 1:100; % quite steep distribution w = w / sum(w); % distribution (pdf) Y_std = sqrt(w .* (1 - w)); % Standard deviations associated to each probability %%%%%%%%%%%%%%%%% Y1 = randp(w, N, 1); % Try Jos function Y1_w = histc(Y1, [1:100]); % absolute frequencies table Y1_w = Y1_w / N; % transform into relative frequencies Y1_w_std = sqrt(Y1_w .* (1 - Y1_w)); % Empirical standard deviations associated to each probability %%%%%%%%%%%%%%%%% Y2 = gDiscrPdfRnd(w, N, 1); % Try Gianluca Dorini's function Y2_w = histc(Y2, [1:100]); % absolute frequencies table Y2_w = Y2_w / N; % transform into relative frequencies Y2_w_std = sqrt(Y2_w .* (1 - Y2_w)); % Empirical standard deviations associated to each probability %%%%%%%%%%%%%%%%% % Check if empirical standard deviations are coherent with theretical ones plot(Y_std - Y1_w_std', '.b') % plot empirical differences for Jos function hold on plot(Y_std - Y2_w_std', '.r') % plot empirical differences for Gianluca Dorini's function %%%%%%%%%%%%%%%%% It looks as though both functions statistically behave very well. They are both excellent functions. Thank you.
11 Sep 2009	Performing random numbers generator from a generic discrete distribution This function extracts random numbers distributed over a discrete set; the PDF is user-defined Author: Gianluca Dorini	Pozzi, Francesco	I've tried this small experiment: N = 100000; % sample size w = 1:100; % quite steep distribution w = w / sum(w); % distribution (pdf) Y_std = sqrt(w .* (1 - w)); % Standard deviations associated to each probability %%%%%%%%%%%%%%%%% Y1 = randp(w, N, 1); % Try Jos function Y1_w = histc(Y1, [1:100]); % absolute frequencies table Y1_w = Y1_w / N; % transform into relative frequencies Y1_w_std = sqrt(Y1_w .* (1 - Y1_w)); % Empirical standard deviations associated to each probability %%%%%%%%%%%%%%%%% Y2 = gDiscrPdfRnd(w, N, 1); % Try Gianluca Dorini's function Y2_w = histc(Y2, [1:100]); % absolute frequencies table Y2_w = Y2_w / N; % transform into relative frequencies Y2_w_std = sqrt(Y2_w .* (1 - Y2_w)); % Empirical standard deviations associated to each probability %%%%%%%%%%%%%%%%% % Check if empirical standard deviations are coherent with theretical ones plot(Y_std - Y1_w_std', '.b') % plot empirical differences for Jos function hold on plot(Y_std - Y2_w_std', '.r') % plot empirical differences for Gianluca Dorini's function %%%%%%%%%%%%%%%%% It looks as though both functions statistically behave very well. They are both excellent functions. Thank you.
04 Sep 2009	RANDP Random integers with given probabilities ( (v2.0) Author: Jos	Pozzi, Francesco	Wow, it really is a nice function, very fast and efficient. It is exactly what I was looking for. Thank you.
17 Jan 2009	Correlation Percentiles CORRPERC estimates percentiles & standard deviations of a correlation matrix, performing a bootstrap Author: Francesco Pozzi	Pozzi, Francesco	Sorry about that. Yes, you need the Statistics Toolbox, I forgot to mention it. Thank you for reminding. I ignore if percentiles can be computed efficiently otherwise.
1 - 4 of 4
E-mail this View

MATLAB Central Terms of Use

NOTICE: Any content you submit to MATLAB Central, including personal information, is not subject to the protections which may be afforded information collected under other sections of The MathWorks, Inc. Web site. You are entirely responsible for all content that you upload, post, e-mail, transmit or otherwise make available via MATLAB Central. The MathWorks does not control the content posted by visitors to MATLAB Central and, does not guarantee the accuracy, integrity, or quality of such content. Under no circumstances will The MathWorks be liable in any way for any content not authored by The MathWorks, or any loss or damage of any kind incurred as a result of the use of any content posted, e-mailed, transmitted or otherwise made available via MATLAB Central. Read the complete Terms prior to use.