如何将多个列作为参数传递给uniroot

Question

我现在的任务是解决一个方程，如eqn(a,b,c) ，从中，a，b，c分别是数据中的一列。 如何将这些uniroot以同时求解根？ 我想也许在这里sapply很有用。

正如listista所建议的那样，我将详细信息粘贴如下：

等式如下：

     eqnt <- function(T,l,fn,m,EI) { T - fn^2*m/l^2 - EI }

现在将所有参数（l，fn，m，EI）加载到一个数据帧中，例如df 。 我打算用df的参数求解所有方程式。

    L        m       EI       fn
 1 10 6.190004 9988.997  6.59710
 2 10 6.190004 9988.997  8.01847
 3 10 6.190004 9988.997  9.21858
 4 15 6.190004 9988.997 10.27676
 5 15 6.190004 9988.997 11.23391
 6 15 6.190004 9988.997 12.11441

我该如何？

当我尝试遵循代码时，定义一个函数并使用mapply ：

rootfun <- function(l,fn,m,EI){
     uniroot.all(eqnt,l,fn,m,EI,lower=0.8e5, upper=3.5e5,
                 interval=c(0.8e5,3.8e5))
 }

 res <- mapply(rootfun,l=data$L,fn=data$f.4,m=data$m,EI=data$EI)

错误出现：

 Error in f(xseq, ...) (from #8) : argument "mm" is missing, with no default

Answer 1

您可以使用pmap但没有一些数据来解决它的问题。 参见例如：

library(purrr)

df <- data.frame(
  x = c("apple", "banana", "cherry"),
  pattern = c("p", "n", "h"),
  replacement = c("x", "f", "q"),
  stringsAsFactors = FALSE
)
pmap(df, gsub)

说明-为您的数据命名，需要与函数的输入进行内联循环（在本例中为uniroot ）。 将数据框与函数一起传递给pmap。