Записать вероятность из функции lm в R

Question

Я хотел бы сделать пробную модель, используя функцию максимального логарифмического правдоподобия.У меня уже есть полиноминальная функция lm в R, и я хотел бы начать с нее.

Вот функция lm, представляющая собой числовую (0,1) y переменную для числовых и факторных переменных:

ee2002ext $partipe2 ~ enf18.f + enf3.f + ee2002ext $ agd + ee2002ext$ expot1, data = ee2002ext, с (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"))

должно работать.Чтобы получить хороший результат, используйте сводку (рег), но вы, вероятно, уже это знали.