How do I simulate Bernoulli observations?

4 views (last 30 days)
Chris Chen
Chris Chen on 15 Oct 2021
Answered: Prachi Kulkarni on 18 Oct 2021
I want to simulate 1000 Bernoulli observations with p = 0.2 and compute an MLE estimate of p along with standard errors and t-statistics, what's the easiest way to do so?

Answers (1)

Prachi Kulkarni
Prachi Kulkarni on 18 Oct 2021
You can generate the Bernoulli observation and get the statistics as follows.
n = 1000;
p = 0.2;
data = binornd(1,p,n,1); % Generate 1000 observatiosn of Bernoulli Trials with p = 0.2
phat = mle(data,Distribution="Bernoulli"); % Maximum Likelihood Estimate of p
se = std(data)/sqrt(n); % Standard Error of Mean
[hypothesis,pvalue,ci,stats] = ttest(data); % t-statistics
For more information, see the following documentation.


Find more on Statistics and Machine Learning Toolbox in Help Center and File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by