[MatLab] media, deviazione standard ecc di variabili di poisson

[lalla8998]

Salve a tutti!
Potete dirmi se questo pezzo di codice è stato scritto bene?

Generare N=100 numeri casuali da una distribuzione di Poisson con parametro 10, ripetere usando N=1000, 10000, 100000; generare N=100 coppie di numeri casuali da due Poisson con parametri 10 e 20; per tutte le variabili calcolare media, deviazione standard, varianza, mediana, decili, percentili; fare boxplot

Lambda=10; %imponiamo il valore del parametro lambda

R=poissrnd(Lambda, 100, 1); %genero un vettore colonna di 100 numeri casuali da una distribuzione di Poisson di parametro 10

R=poissrnd(Lambda, 1000, 1); 

R=poissrnd(Lambda, 10000, 1);

R=poissrnd(Lambda, 100000, 1);

Lambda1=10;
Lambda2=20;

R1=poissrnd(Lambda1, 100, 1);

R2=poissrnd(Lambda2, 100, 1);

media=mean(R) %media

media1=mean(R1)

media2=mean(R2)

%%
devS=std(R) %deviazione standard

devS1=std(R1)

devS2=std(R2)

%%
varianza=(devS)^2 %varianza

varianza1=(devS1)^2

varianza2=(devS2)^2

%%
mediana=median(R)

mediana1=median(R1)

mediana2=median(R2)

%%
quantile(R,[0.1 0.5]) % vediamo il primo e il quinto decile (il decile divide la distribuzione in 10 parti)

quantile(R,[0.01 0.1 0.5]) %vediamo il primo, il decimo e il cinquantesimo percentile (il percentile divide la distribuzione in 100 parti)
 
%%
quantile(R, 9) %vediamo tutti i decili

quantile(R,99) %vediamo tutti i percentili

%%
boxplot(R)

boxplot(R1)

boxplot(R2)