问题描述
我正在尝试理解逻辑回归的源代码（R 中的 glm 函数）。
在summary.fit 函数的p 值计算步骤中，使用了一个名为chol2inv 的函数，该函数应该执行(X'X)^(-1)。
然而，chol2inv(X) 和 inv(X'X) 的输出不匹配。你能检查一下如何从我使用的代码中协调它们吗？
X <- matrix(c(11,12,11),2,2)
chol2inv(X)
solve(t(X) %*% X)
解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）
