- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我已经为一些数据拟合了一个多项式模型,我想从该模型中提取公式以找到它的最大值。我能够从 lm 对象中提取一个公式作为字符串,但是我无法从该字符串创建一个将与 optimize 函数一起使用的新函数。
## function for generating data
f1 = function(x) 1 + x^2 - x^3
## random variable from normal distribution
yran = rnorm(500, 1, .025)
## create data by fitting function f to x points times the random variable
dt = data.frame(x = seq(.01, 1, .01) * yran,
y = f1(seq(.01, 1, .01))*yran)
## sort data frame by x
dt = dt[order(dt$x, decreasing = FALSE), ]
## plot the generated data
plot(dt, ylim = c(.9, 1.24))
## create a polynomial model
fit3 = lm(y ~ poly(x, 3), dt)
## plot the models over the data
lines(x = dt$x, predict(fit3, data.frame(x = dt$x)), col = "red", lwd = 2)
## fit the original model for comparison
lines(x = dt$x, f1(dt$x), lwd = 2, lty = 2)
我创建的函数只是提取系数并将公式粘贴在一起。我的挑战是能够从模型字符串创建一个函数。
ext.mdl = function(lm) {
int = paste(lm$coefficients[[1]][[1]])
coef = paste(lm$coefficients[2:length(lm[[1]])])
out = as.character()
for (i in 1:length(coef)) {
out = paste(out, coef[i], "*x^", i, " + ", sep = "")
}
out = gsub('.{3}$', '', out)
out = paste(int, '+', out)
return(out)
}
> ext.mdl(fit3)
[1] "1.08475891509144 + 0.599668223720749*x^1
> + -0.822484955777266*x^2 + -0.377150292824362*x^3"
理想情况下,我希望能够为从 ext.mdl() 中提取的任何内容分配一个新函数,这样我就可以使用 optimize() 找到最大值在函数中。最终,我需要能够通过“function(x) [model string]”进行优化。
> optimize(function(x) 1.08475891509144 + 0.599668223720749*x^1
+ + -0.822484955777266*x^2 + -0.377150292824362*x^3,
+ interval = c(0,1), maximum = TRUE)
$maximum
[1] 0.3018792
$objective
[1] 1.180457
有什么想法吗?
最佳答案
#parse the string and then evaluate the expression
optimize(function(x) eval(parse(text=ext.mdl(fit3))),
interval = c(0,1), maximum = TRUE)
$maximum
[1] 0.3007581
$objective
[1] 1.179404
关于R将字符串公式转换为函数,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/28908494/
如何使用 SPListCollection.Add(String, String, String, String, Int32, String, SPListTemplate.QuickLaunchO
我刚刚开始使用 C++ 并且对 C# 有一些经验,所以我有一些一般的编程经验。然而,似乎我马上就被击落了。我试过在谷歌上寻找,以免浪费任何人的时间,但没有结果。 int main(int argc,
这个问题已经有答案了: In Java 8 how do I transform a Map to another Map using a lambda? (8 个回答) Convert a Map>
我正在使用 node + typescript 和集成的 swagger 进行 API 调用。我 Swagger 提出以下要求 http://localhost:3033/employees/sear
我是 C++ 容器模板的新手。我收集了一些记录。每条记录都有一个唯一的名称,以及一个字段/值对列表。将按名称访问记录。字段/值对的顺序很重要。因此我设计如下: typedef string
我需要这两种方法,但j2me没有,我找到了一个replaceall();但这是 replaceall(string,string,string); 第二个方法是SringBuffer但在j2me中它没
If string is an alias of String in the .net framework为什么会发生这种情况,我应该如何解释它: type JustAString = string
我有两个列表(或字符串):一个大,另一个小。 我想检查较大的(A)是否包含小的(B)。 我的期望如下: 案例 1. B 是 A 的子集 A = [1,2,3] B = [1,2] contains(A
我有一个似乎无法解决的小问题。 这里...我有一个像这样创建的输入... var input = $(''); 如果我这样做......一切都很好 $(this).append(input); 如果我
我有以下代码片段 string[] lines = objects.Split(new string[] { "\r\n", "\n" }, StringSplitOptions.No
这可能真的很简单,但我已经坚持了一段时间了。 我正在尝试输出一个字符串,然后输出一个带有两位小数的 double ,后跟另一个字符串,这是我的代码。 System.out.printf("成本:%.2
以下是 Cloud Firestore 列表查询中的示例之一 citiesRef.where("state", ">=", "CA").where("state", "= 字符串,我们在Stack O
我正在尝试检查一个字符串是否包含在另一个字符串中。后面的代码非常简单。我怎样才能在 jquery 中做到这一点? function deleteRow(locName, locID) { if
这个问题在这里已经有了答案: How to implement big int in C++ (14 个答案) 关闭 9 年前。 我有 2 个字符串,都只包含数字。这些数字大于 uint64_t 的
我有一个带有自定义转换器的 Dozer 映射: com.xyz.Customer com.xyz.CustomerDAO customerName
这个问题在这里已经有了答案: How do I compare strings in Java? (23 个回答) 关闭 6 年前。 我想了解字符串池的工作原理以及一个字符串等于另一个字符串的规则是
我已阅读 this问题和其他一些问题。但它们与我的问题有些无关 对于 UILabel 如果你不指定 ? 或 ! 你会得到这样的错误: @IBOutlet property has non-option
这两种方法中哪一种在理论上更快,为什么? (指向字符串的指针必须是常量。) destination[count] 和 *destination++ 之间的确切区别是什么? destination[co
This question already has answers here: Closed 11 years ago. Possible Duplicates: Is String.Format a
我有一个Stream一个文件的,现在我想将相同的单词组合成 Map这很重要,这个词在 Stream 中出现的频率. 我知道我必须使用 collect(Collectors.groupingBy(..)
我是一名优秀的程序员,十分优秀!