使用来自不同数据帧的值作为嵌套函数的输入

Question

我如何运行以下代码，使用来自数据框I 的值的函数ftn4 的输出被newtonraphson 函数与来自数据框R0 的值结合使用。我想将输出存储在与I 和R0 具有相同维度（ncol = 5，nrow = 5）的新数据框（R）中。

我被卡住了，有什么想法吗？任何提示都非常感谢

c           = c(0.3)      
Ccf         = c(0.3)      
Acf         = c(3*10^-6)  
E           = c(0.00006) 
Ke          = c(1000)    

I           = as.data.frame(matrix(sample(seq(150,300,0.1),25), ncol = 5, nrow = 5))
R0          = as.data.frame(matrix(sample(seq(0.01,0.2,0.001),25), ncol = 5, nrow = 5))

library(spuRs) 

ftn4 <- function(x) {
  f = { x^2*exp(c*x)- (Ccf*Acf*E*(I/Ke+I))}
  f1 = { 2 * x * exp(c * x) + x^2 * (exp(c * x) * c)}
        return(c(f,f1))
}

R = as.data.frame(newtonraphson(ftn4, R0, tol = 1e-09 ,max.iter = 1000))

Answer 1

按照 cmets 的规定，您应该将参数传递给函数 ftn4，该函数产生一个值。所以ftn4，应该改成：

ftn4 <- function(x,I) {
  f = { x^2*exp(c*x)- (Ccf*Acf*E*(I/Ke+I))}
  f1 = { 2 * x * exp(c * x) + x^2 * (exp(c * x) * c)}
  return(c(f,f1))
}

现在遍历您的矩阵 I 和 R0。

R <- data.frame(matrix(NA,nrow=nrow(I), ncol(I)))
for(i in 1:nrow(I)){
  for(j in 1:ncol(I)){
    R[i,j] <- newtonraphson(ftn  = function(x) ftn4(x, I[i,j]), 
                            x0  = R0[i,j], 
                            tol = 1e-09, 
                            max.iter = 1000)
  }
}