1

我需要使用该nlm函数按组插入一个大型数据框。我在单个组的 df 上使用它没有任何问题：

#example data
df <- data.frame(var= cumsum(sort(rnorm(100, mean=20, sd=4))),
                 time= seq(from=0,to=550,length.out=100))
#create function
my_function <- function(Cini, time, theta,var){
  fy <- (theta[1]-(theta[1]- Cini)*exp((-theta[2]/100000)*(time-theta[3])))
  ssq<-sum((var-fy)^2)
  return(ssq)
}
th.start <- c(77, 148, 5)   #set starting parameters

#run nlm
my_fitt <- nlm(f=my_function, Cini=400, var = df$var,
               time=df$time, p=th.start)

然后，我尝试使用该函数在具有多个组的 df 中应用该dlply函数：

#data with groups
df.2 <- data.frame(var= cumsum(sort(rnorm(300, mean=20, sd=4))),
                   time= rep(seq(from=0,to=1200,length.out=100),3),
                   groups=rep(c(1:3),each=100))
#run nlm
library(plyr)
my_fitt.2 <- dlply(df.2, .(groups),
               nlm(f=my_function, Cini=400, var  = df.2$var,time=df.2$time, p=th.start))

但是我得到消息：Error in fs[[i]](x, ...) : attempt to apply non-function。我还尝试在此示例中以及在我的原始 df （是变量之一）中删除df.2$, 获取。Error in time - theta[3] : non-numeric argument to binary operatorError in f(x, ...) : object 'time.clos' not foundtime.clos

另外，我还没有使用 dplyr 库

library(dplyr)
df.2 %>%
  group_by(groups) %>%
  nlm(f=my_function, Cini=400, v= var,
      time=time, p=th.start)

获得Error in f(x, ...) : unused argument (.). 可能是什么问题呢？

4

2 回答 2

2

考虑基本 R by（面向对象的包装器tapply），它可以按因子对数据帧进行子集化，并将子集化的数据帧传递到诸如您的nlm调用之类的方法中，所有这些都返回对象列表：

run_nlm <- function(sub_df) nlm(f=my_function, Cini=400, var=sub_df$var, 
                                time=sub_df$time, p=th.start)

# LIST OF nlm OUTPUTS (EQUAL TO NUMBER OF DISTINCT df$groups)
my_fitt_list <- by(df, df$groups, run_nlm)
于 2018-07-31T15:33:27.797 回答
1

我对环境帮不上什么忙，tidyverse因为我更像是一个基本的 R 类型的人。我认为你上次调用的问题是你正在将一个组data.frame传递给一个将function对象作为第一个参数的函数。那是行不通的。

让我向您推荐一种基本的 R 方式：

df.2 %>% 
  split(.$groups) %>% 
  lapply(function(xx) nlm(f=my_function, Cini=400, var = xx$var, time=xx$time, p=th.start))

这会产生list长度为 3 的 a（三组）和您的三个结果。

于 2018-07-31T12:56:51.870 回答