- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我正在尝试使用数据掩码评估最大似然表达式。这个想法是允许在函数内按名称调用参数和变量,同时避免多次调用attach()
。和 detach()
.这是一个非常简单的小例子,实际功能非常庞大和复杂。
set.seed(1)
# Data
db <- data.frame(
x = runif(10),
y = runif(10),
z = sample(c(0, 1), 10, replace = TRUE)
)
# Log likelihood function
ll_lik <- function(param) {
pr_1 <- 1 / (1 + exp(-(param[1]*x - param[2]*y)))
pr_2 <- 1 - pr_1
lik <- z * pr_1 + (1 - z) * pr_2
log(lik)
}
# Parameters
param <- c(p1 = 0.1, p2 = 0.2)
# Run the model with attach()/detach()
attach(db)
model <- maxLik::maxLik(ll_lik, start = param)
detach(db)
summary(model)
attach()
和
detach()
.要按名称访问参数,我需要转换
param
到对数似然函数内的一个列表,然后调用
attach()
/
detach()
.这不仅是困惑的,而且对于大型函数和数据,它会产生不必要的开销。我一直在研究的一种可能性是使用
rlang
package 和 wrapper 函数主要是为表达式的整洁评估而编写的。
mask <- as_data_mask(db)
eval_tidy(quo(maxLik::maxLik(ll_lik, start = param)), mask)
Error in fnOrig(theta, ...) : object 'x' not found
) 中的对象。也许问题出在
maxLik
,但我什至无法评估
ll_lik()
这给出了同样的错误:
eval_tidy(quo(ll_lik(param)), mask)
eval_tidy(quo(x*3), mask)
ll_lik()
有“错误”的父级,这就是为什么我的数据掩码可能不在函数的搜索路径中,因此它找不到变量。现在,
as_data_mask()
的帮助部分确实提供了一些示例,说明如何通过创建顶层、中层和底层环境来“嵌套”环境。好的,让我们看看我是否可以将我的函数创建为数据掩码结构的一部分:
call_stack <- function() {lobstr::cst()}
# Create a new environment (child of empty) that takes a list of objects to populate it
top <- new_environment(list(ll_lik = ll_lik, call_stack = call_stack))
# Create a child of the "top" environment"
middle <- env(top)
# Create a child of the "middle environment and add the data object to it
bottom <- env(middle, db=db)
# Create a data_mask where the bottom contains the masking elements and the top
# the last element of the data_mask.
new_mask <- new_data_mask(bottom, top = top)
x
.我什至没有尝试
maxLik
在这里发挥作用。所以为了尝试更深入地挖掘,我开始搞乱调用堆栈。
eval_tidy(call_stack(), data = new_mask)
█
1. ├─rlang::eval_tidy(call_stack(), data = new_mask)
2. └─global::call_stack()
3. └─lobstr::cst()
maxLik
中按名称调用参数无需调用
attach()
/
detach()
, 那将是真棒。
最佳答案
一种选择是创建一个包装器来评估 ll_lik
的主体。作为表达式,使用 db
作为上下文:
llwrap <- function(param) {
eval( body(ll_lik), db )
}
model <- maxLik::maxLik(llwrap, start=param) # Works
body()
返回一个表达式,因此您可以在该表达式中使用您想要的任何名称,只要您在评估时提供适当的上下文。但是,如果您将函数的主体与其参数列表完全分离,为什么不从一开始就将其定义为表达式呢?
ll_expr <- rlang::expr({ # An expression, not a function
pr_1 <- 1 / (1 + exp(-(p1*x - p2*y))) # <-- now using p1, p2
pr_2 <- 1 - pr_1
lik <- z * pr_1 + (1 - z) * pr_2
log(lik)
})
llwrap2 <- function(param) {
ctx <- c( as.list(db), as.list(param) ) # Combine param and db into one context
eval( ll_expr, ctx ) # No longer need body()
}
model <- maxLik::maxLik(llwrap2, start=param) # Works
关于r - 使用 R 中的数据掩码评估最大似然表达式,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/61834186/
初学者 android 问题。好的,我已经成功写入文件。例如。 //获取文件名 String filename = getResources().getString(R.string.filename
我已经将相同的图像保存到/data/data/mypackage/img/中,现在我想显示这个全屏,我曾尝试使用 ACTION_VIEW 来显示 android 标准程序,但它不是从/data/dat
我正在使用Xcode 9,Swift 4。 我正在尝试使用以下代码从URL在ImageView中显示图像: func getImageFromUrl(sourceUrl: String) -> UII
我的 Ubuntu 安装 genymotion 有问题。主要是我无法调试我的数据库,因为通过 eclipse 中的 DBMS 和 shell 中的 adb 我无法查看/data/文件夹的内容。没有显示
我正在尝试用 PHP 发布一些 JSON 数据。但是出了点问题。 这是我的 html -- {% for x in sets %}
我观察到两种方法的结果不同。为什么是这样?我知道 lm 上发生了什么,但无法弄清楚 tslm 上发生了什么。 > library(forecast) > set.seed(2) > tts lm(t
我不确定为什么会这样!我有一个由 spring data elasticsearch 和 spring data jpa 使用的类,但是当我尝试运行我的应用程序时出现错误。 Error creatin
在 this vega 图表,如果我下载并转换 flare-dependencies.json使用以下 jq 到 csv命令, jq -r '(map(keys) | add | unique) as
我正在提交一个项目,我必须在其中创建一个带有表的 mysql 数据库。一切都在我这边进行,所以我只想检查如何将我所有的压缩文件发送给使用不同计算机的人。基本上,我如何为另一台计算机创建我的数据库文件,
我有一个应用程序可以将文本文件写入内部存储。我想仔细看看我的电脑。 我运行了 Toast.makeText 来显示路径,它说:/数据/数据/我的包 但是当我转到 Android Studio 的 An
我喜欢使用 Genymotion 模拟器以如此出色的速度加载 Android。它有非常好的速度,但仍然有一些不稳定的性能。 如何从 Eclipse 中的文件资源管理器访问 Genymotion 模拟器
我需要更改 Silverlight 中文本框的格式。数据通过 MVVM 绑定(bind)。 例如,有一个 int 属性,我将 1 添加到 setter 中的值并调用 OnPropertyChanged
我想向 Youtube Data API 提出请求,但我不需要访问任何用户信息。我只想浏览公共(public)视频并根据搜索词显示视频。 我可以在未经授权的情况下这样做吗? 最佳答案 YouTube
我已经设置了一个 Twilio 应用程序,我想向人们发送更新,但我不想回复单个文本。我只是想让他们在有问题时打电话。我一切正常,但我想在发送文本时显示传入文本,以确保我不会错过任何问题。我正在使用 p
我有一个带有表单的网站(目前它是纯 HTML,但我们正在切换到 JQuery)。流程是这样的: 接受用户的输入 --- 5 个整数 通过 REST 调用网络服务 在服务器端运行一些计算...并生成一个
假设我们有一个名为 configuration.js 的文件,当我们查看内部时,我们会看到: 'use strict'; var profile = { "project": "%Projec
这部分是对 Previous Question 的扩展我的: 我现在可以从我的 CI Controller 成功返回 JSON 数据,它返回: {"results":[{"id":"1","Sourc
有什么有效的方法可以删除 ios 中 CBL 的所有文档存储?我对此有疑问,或者,如果有人知道如何从本质上使该应用程序像刚刚安装一样,那也会非常有帮助。我们正在努力确保我们的注销实际上将应用程序设置为
我有一个 Rails 应用程序,它与其他 Rails 应用程序通信以进行数据插入。我使用 jQuery $.post 方法进行数据插入。对于插入,我的其他 Rails 应用程序显示 200 OK。但在
我正在为服务于发布请求的 API 调用运行单元测试。我正在传递请求正文,并且必须将响应作为帐户数据返回。但我只收到断言错误 注意:数据是从 Azure 中获取的 spec.js const accou
我是一名优秀的程序员,十分优秀!