[英]Log likelihood from a lm function in R
我想使用最大对数似然函数做一个概率模型。 我已经在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)
}
谢谢,
亚历山大
您可以使用glm()运行概率模型。 就你而言
reg <- glm(participe2 ~ enf18.f + enf3.f + agd + expot1 , data = ee2002ext, family = binomial(link = "probit"))
应该管用。 要获得良好的输出,请使用summary(reg),但您可能已经知道了。
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.