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
- Codice:
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)