作者celestialgod (天)
看板R_Language
标题Re: [问题] 利用newton-raphson method 求回归系数
时间Sat Dec 5 19:10:06 2015
我懒得推文了...
问题如下：
1. 你知道logistic的MLE怎麽算吗？
2. 你有了牛顿法的程式，但是你知道f3放什麽？x0放什麽吗？
3. 我看不懂你创的矩阵，[Y-pi1]是从哪里来的？
4. 从你的叙述里面，似乎都对newton这个函数不熟悉，你要怎麽跑这个程式？
5. 你知道R可以直接用glm求logistic regression的回归系数吗？
※ 引述《ritajen (asdfge)》之铭言：
: 如题，利用牛顿法求logistic regression回归系数的最大概似估计值，
: 目前我有牛顿法的雏型，
: rm(list = ls())
: newton <- function(f3, x0, tol = 1e-9, n.max = 100) {
: x <- x0
: f3.x <- f3(x)
: n <- 0
: while ((max(abs(f3.x[[1]])) > tol) & (n < n.max)) {
: x <- x - solve(f3.x[[2]]) %*% f3.x[[1]]
: f3.x <- f3(x)
: n <- n + 1
: }
: if (n == n.max) {
: cat('newton failed to converge\n')
: } else {
: return(x)
: }
: }
: 以下是我的data，我目前了解需先创立一个新的矩阵
: X=[1 gpa gre] (共有n行是一个n*3的矩阵)
: 再创造一个
: [Y-pi1](共有n行是一个n*1的矩阵，其中pi是机率)
: 但这个地方我不太明白pi要是多少?
: 是利用上面所创作出来的矩阵相乘得到f3，在代回上述求解吗?
: 这样解会是一个矩阵? (回归系数)
: 麻烦厉害的人求解 感谢!
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.109.73.236
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1449313809.A.BBC.html
1F：推 obarisk: 别这麽凶嘛，感觉像是数值课的作业 12/05 19:22
我只是平铺直叙的问问题啦QQ
2F：→ obarisk: 文字比较容易有这种感觉Orz 12/06 09:09
文字生硬阿XD
※ 编辑: celestialgod (140.109.73.234), 12/06/2015 13:21:47