Re: [問題] R跑probit模型
※ 引述《k84526991 (哇哈哈)》之銘言：
: 要如何用r跑probit模型呢??
: 看計量理論的書 必須要用maxliklihood 去找估計系數值
: 請問這要麼用r跑呢??
: 我不要用glm(.....)
: 有沒有人可以教我詳細的語法
我沒有實際用過R來解 但是有用過Matlab來解
下面可能可以提供給你一些概念
以probit model為例
1. 先建立probit model的log-likelihood function
2. 然後以R的package來求解 剛去查好像蠻多可以解ex :maxlik
maxLik(logLik, grad = NULL, hess = NULL, start, method,
logLik 是log-likelihood function
grad,hess 是gradient of log-likelihood及hessian of log-likelihood
start 是起始值
method 是最大化的方法 ex: "Newton-Raphson", "BFGS", "BHHH",
3. 求解完成後 通常你會需要的有 function value at maximum
估計係數值(estimated parameters)
Hessian at the maximum
4. 根據求解得到的值 你可以算出
Esti. Coef. Std. Err. t-value p-value
一般統計軟體還會進一步求出marginal effect及其std. err.
不過你可以先試看看做到第四步
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.86.204
討論串 (同標題文章)