O máximo de dois normais não idênticos pode ser expresso como uma distribuição normal de distorção de Azzalini. Ver, por exemplo, um documento / apresentação de trabalho de 2007 de Balakrishnan
Um olhar distorcido das estatísticas de ordens bivariadas e multivariadas
Prof. N. Balakrishnan
Working paper / presentation (2007)
Um artigo recente de ( Nadarajah e Kotz - visível aqui ) fornece algumas propriedades de max :( X, Y)
Nadarajah, S. e Kotz, S. (2008), "Distribuição exata das máximas / mínimas de duas variáveis aleatórias gaussianas", TRANSAÇÕES IEEE EM SISTEMAS DE INTEGRAÇÃO DE ESCALA MUITO GRANDE (VLSI), VOL. 16, NO. 2, fevereiro de 2008
Para trabalhos anteriores, consulte:
AP Basu e JK Ghosh, "Identificabilidade das distribuições multinormais e outras sob o modelo de riscos concorrentes", J. Multivariate Anal., Vol. 8, pp. 413–429, 1978
HN Nagaraja e NR Mohan, "Sobre a independência da distribuição da vida do sistema e a causa da falha", Scandinavian Actuarial J., pp. 188-198, 1982.
YL Tong, a distribuição normal multivariada. Nova York: Springer-Verlag, 1990.
Pode-se também usar um sistema de álgebra computacional para automatizar o cálculo. Por exemplo, dado com pdf e com pdf :X∼ N( μ1 1, σ21 1)f( X )Y∼ N( μ2, σ22)g( y)

... o pdf de Z= m a x ( X, Y) é:

onde estou usando a Maximum
função do pacote mathStatica do Mathematica e Erf
denota a função de erro.