Bonjour,
Soit une distribution de loi F doublement non centrée définit avec 4 paramètres, les degrés de liberté n1 et n2 et les paramètres de non centralité mu1 et mu2. La loi F doublement non centrée est en fait le ratio de 2 loi chi2 non centrée divisée par le nombre de ddl :
F(n1,n1,mu1,mu2) = (Chi2(n1,mu1)/n1) / (Chi2(n2,mu2)/n2)
L'approximation de Patnaik permet d'exprimer une loi CHi2 non centrée en loi CHi2 centrée : Chi2(n1,mu1) = c Chi2(n3) avec c=(n1+2mu1)/(n1+mu1) et n3=n1 + (mu1²)/(n1+2mu1)
Mon problème est que je n'arrive pas à implémenter sous Matlab, l'approximation d'une loi F doublement non centrée en loi F simplement non centrée :
F(n1,n2,mu1,mu2) = (1/(1+mu2/n2))F(n1,n3,mu1)
Quelqu'un a t-il déja essayé d'implémenter cela sous Matlab ? L'approximation de chi2 non centrée en chi2, j'y arrive mais pas pour la loi F.
Merci d'avance pour vos éclairements
-----