Estou tentando criar um ECDF (e um limite de confiança) a partir de dados em Python. Posso gerar o ECDF facilmente com a numpy
classificação e o uso linspace
. No entanto, não estou totalmente certo de quais são os limites de confiança apropriados e não parece haver nenhuma biblioteca statsmodels
interna que calcule os limites ( parece apenas fornecer o ECDF).
Se eu quero uma confiança pontual vinculada aé apropriado usar a desigualdade DKW para calcular minha região com
Onde é o número de observações na minha amostra? Assim, se é meu ECDF, meus limites superior e inferior seriam
O MATLAB possui uma função interna ECDF , mas não tive muita sorte em entender como aplicar a fórmula de Greenwood (mencionada na parte inferior) para gerar os limites.
Respostas:
No console do Matlab, digite:
Abre o código fonte no editor.
Vá para a linha 194:
Este é o início do bloco de código que calcula os mais baixos - limites (de confiança) e superior:
[Flo, Fup]
. O bloco de código tem 30 linhas e é bastante direto. Postado abaixo para sua conveniência:A raiz quadrada da fórmula de Greenwood, ou seja,
é implementado na linha 203 como:
Você pode pegar daqui? Avise-se me.
fonte