function [P] = gmpe_fcc_FN(iwyn)
    m = numel(iwyn(1).b) - 1 ;
    n = size(iwyn(1).x,2) ;
    M = numel(iwyn) ;
    P = zeros(M,1) ;
    df = ones(M,1)*n ;
    df0 = n ;
    for j = 1:M ;
        df0 = df0 - m - j ;
        df(j) = df0 ;
    end
    for k2 = 2:M
        for k1 = 1:k2
            F = ((iwyn(k1).SSE - iwyn(k2).SSE) / (df(k1) - df(k2))) / (iwyn(k2).SSE / df(k2)) ;
            P(k2)= max(P(k2),1-fcdf(F,(df(k1) - df(k2)),df(k2)));
        end
    end
    for j = 1:M
        mp = min(P(j+1:end)) ;
        if P(j) > mp
            P(j) = mp ;
        end
    end
end