二変量3パラメータ関数のMLEの求め方R

Question

以下の関数で3つのパラメータのMLEを求める必要があります（X & Yのランダム2次元点の集合）。私はRの中に複雑な関数を入力する方法を

？

ご協力いただければ幸いです。

Answer 1

私はこれがあなたが探していることを望みます。 dnormは、normalの確率密度関数です（meanとsdが指定されていない場合は、デフォルトは標準のnormalになります）。 pnormは（平均場合は、標準正規にもデフォルトと供給され、私たちが望むものである下テール領域を与えないSD）、通常の累積分布関数である：

# Construct function 
MLE= function(x, y, mu1, mu2, theta){ 
    numerator = dnorm(x-mu1)*dnorm(y-mu2)*theta*(1-exp(-theta))*exp(-theta*(pnorm(x-mu1)+pnorm(y-mu2))) 
    denominator = (1-exp(-theta)-(1-exp(-theta*pnorm(x-mu1)))*(1-exp(-theta*pnorm(y-mu2))))^2 
    result = numerator/denominator 
    return(result) 
}