Question

想象一下这两个列表如下：

seq<-seq(1,16)
L1<-list(A=seq, B=seq, C=seq, D=seq)

v1=c(11,15,17,19,21,22,24,25,26,27,28,28,29,29,30,30)
v2=c(14,19,23,24,26,27,28,29,30,30,31,32,32,32,32,33)
v3=c(2,4,5,6,6,7,8,8,9,9,9,10,10,10,11,11)
v4=c(8,13,17,20,22,24,26,27,28,29,30,31,32,33,33,34)
L2<-list(A=v1,B=v2,C=v3,D=v4)

现在我要做的是使这些列表适合使用nls公式的特定模型，L1 $ A使用L2 $ A，L1 $ B使用L2 $ B等等。我可以单独使用每个向量，例如：

nls.A<-nls(L2.A~(a/b)*(1-exp(-b*L1.A)),
    start=list(a=1, b=0.1),
    trace= TRUE, data=data.frame(L1$A, L2$A))

但我想对每个列表的四个四个元素进行检索，并检索一个对象，该对象也是一个包含所有四个（$ A，$ B，$ C，$ D）拟合结果的列表。

我尝试使用for循环虽然不是很推荐，但是这样

for (i in L1) {
  nls(L2.i~(a/b)*(1-exp(-b*L1.i)),
      start=list(a=1, b=0.1),trace= TRUE,
  data=data.frame(L1$i, L2$i))
}

但那也不能回答我的问题（我真的需要这个该死的名单）我认为拉普利的东西是理想的，但不知道如何使用lapply做这样的事情。有什么想法吗？

Answer 1

将mapply与匿名函数一起使用

res <- mapply(function(x,y){
  nls(y~(a/b)*(1-exp(-b*x)),
      start=list(a=1, b=0.1),
      trace= TRUE, data=data.frame(x, y))
},L1,L2, SIMPLIFY=FALSE)

将nls应用于列表中的元素

1 个答案: