Python中文
首页
教程
问答
标签
搜索
登录
注册
极值的logit和逆logit函数
回答此问题可获得
20
贡献值,回答如果被采纳可获得
50
分。
<p>我需要logit和逆logit函数,以便<code>logit(inv_logit(n)) == n</code>。我用的是numpy,这是我所拥有的:</p> <pre><code>import numpy as np def logit(p): return np.log(p) - np.log(1 - p) def inv_logit(p): return np.exp(p) / (1 + np.exp(p)) </code></pre> <p>以下是价值观:</p> <pre><code>print logit(inv_logit(2)) 2.0 print logit(inv_logit(10)) 10.0 print logit(inv_logit(20)) 20.000000018 #well, pretty close print logit(inv_logit(50)) Warning: divide by zero encountered in log inf </code></pre> <p>现在让我们测试负数</p> <pre><code>print logit(inv_logit(-10)) -10.0 print logit(inv_logit(-20)) -20.0 print logit(inv_logit(-200)) -200.0 print logit(inv_logit(-500)) -500.0 print logit(inv_logit(-2000)) Warning: divide by zero encountered in log -inf </code></pre> <p>所以我的问题是:实现这些功能的正确方法是什么,以便要求<code>logit(inv_logit(n)) == n</code>在尽可能宽的范围内(至少是[-1e4;1e4)对任何<code>n</code>都有效?</p> <p>另外(我肯定这和第一个有关),为什么我的函数在负值的情况下比正值更稳定?</p>
0 条评论
分类:
Python问答
请先
登录
后评论
默认排序
时间排序
1 个回答
匿名
1天前
擅长:python、mysql、java
<p>你正在接近IEEE 754双精度浮点的精度限制。如果需要更大的范围和更精确的域,则需要使用更高精度的数字和操作。</p> <pre><code>>>> 1 + np.exp(-37) 1.0 >>> 1 + decimal.Decimal(-37).exp() Decimal('1.000000000000000085330476257') </code></pre>
请先
登录
后评论
针对此问题:
更多的回答
关注
89
关注
收藏
1
收藏,
216
浏览
网友 提问于 2天前
相关Python问题
VirtualEnvRapper错误:路径python2(来自python=python2)不存在
1 回答
virtualenvs上的pyinstaller,没有名为导入错误的模块
2 回答
virtualenvs是否可以退回到用户包而不是系统包?
8 回答
virtualenvwrapper CentOS7
5 回答
virtualenvwrapper IOError:[Errno 13]权限被拒绝
8 回答
virtualenvwrapper mkproject和shell在windows中的启动问题?
7 回答
virtualenvwrapper mkvirtualenv不工作但没有错误
5 回答
Virtualenvwrapper python bash
7 回答
virtualenvwrapper:“workon”何时更改到项目目录?
9 回答
virtualenvwrapper:mkvirtualenv可以工作,但是rmvirtualenv返回bash:没有这样的文件或目录
1 回答
virtualenvwrapper:virtualenv信息存储在哪里?
3 回答
virtualenvwrapper:命令“python设置.pyegg_info“失败,错误代码为1
1 回答
virtualenvwrapper:如何将mkvirtualenv的默认Python版本/路径更改为ins
2 回答
Virtualenvwrapper:模块“pkg_resources”没有属性“iter_entry_points”
6 回答
Virtualenvwrapper:没有名为virtualenvwrapp的模块
2 回答
Virtualenvwrapper.bash_profi的正确设置
1 回答
Virtualenvwrapper.hook:权限被拒绝
9 回答
virtualenvwrapper.sh:fork:资源暂时不可用Python/Djang
1 回答
Virtualenvwrapper.shlssitepackages命令不工作
5 回答
Virtualenvwrapper.sh函数在bash sh中不可用
6 回答