Ho scritto il seguente codice per tracciare la funzione Logit dove beta = 1.
beta = 1
e(:,1) = 0:0.01:1;
e(:,2) = log(e(:,1)./(1-e(:,1)))+beta;
plot(e(:,2),e(:,1))
Per quanto posso dire che sembra a posto, ma non lo facciosapere come aggiungere la funzione Probit. Capisco come mettere in grafici in una figura, ma non capisco come generare dati da questa distribuzione. C'è un modo conveniente per farlo in MATLAB?
In definitiva, sto cercando di creare qualcosa come il seguente, in cui si può vedere che Probit e Logit sono molto simili.
risposte:
3 per risposta № 1beta = 1
e(:,1) = 0:0.01:1;
e(:,2) = log(e(:,1)./(1-e(:,1)))+beta; % logit fn + beta
e(:,3) = sqrt(2) * erfinv(2*e(:,1) - 1) ; % probit fn
plot(e(:,1),e(:,2:3))
Come suggerito da @patrick se si dispone della "tool machine and statistic toolbox", sarebbe più pulito utilizzare il file norminv
funzione, cioè sostituire la quarta riga con:
e(:,3) = norminv(e(:,1)) ; % probit fn