從R中的lm函數記錄對數似然

Question

我想使用最大對數似然函數做一個概率模型。 我已經在R中擁有一個多項式lm函數，並希望從那里開始。

這是lm函數，它是數字和因子變量的數字（0,1）y變量：

ee2002ext $ participe2〜enf18.f + enf3.f + ee2002ext $ agd + ee2002ext $ expot1，data = ee2002ext，with（ee2002ext，s == 1））

我該怎么做呢？

我暫時這樣做：

    function(beta0, beta1, beta2, beta3, beta4, mu, sig) {
R=participe - enf3.f *beta1 - enf18.f*beta2 - age*beta3 - expot*beta4 - beta0 
R= suppressWarnings(dnorm(R, mu, sig))
-sum(R)
}

謝謝，

亞歷山大

Answer 1

您可以使用glm（）運行概率模型。 就你而言

reg <- glm(participe2 ~ enf18.f + enf3.f + agd + expot1 , data = ee2002ext, family = binomial(link = "probit"))

應該管用。 要獲得良好的輸出，請使用summary（reg），但您可能已經知道了。