3

我最近问了一个关于apply在两个列表上使用函数的问题。每个列表都是通过拆分大型数据框创建的数据框列表。每次函数运行时，我都想从第一个元素（数据帧）中获取向量，并从第一个元素（数据帧）中mylist1获取一些向量，mylist2并将它们相互回归。然后移动到下一个mylist1元素和mylist2元素。实际上，该函数采用两个具有相同数量元素的列表并采用一对（每个列表中的一个）并使用它们。

我尝试了以下方法，但得到的结果不是我想要的：

a1<-c(1:5,rep(0,5))
a2<-c(1:5,10:6)
b2<-c(rep(100,5),rep(50,5))
z<-c(rep("part1",5),rep("part2",5))
df1<-data.frame(a1,z)
df2<-data.frame(a2,b2,z)

mylist1<-split(df1,z)
mylist2<-split(df2,z)


myfunction<-function(x,y) 
{

meana <- mean(x$a)
meanb <- mean(y$b)
model<-lm((x$a)~(y$a))
return(c(model$coefficients[2],meana=meana,meanb=meanb))
}

result <- mapply(myfunction,x=mylist,y=mylist2)

 #result
#        x   y
#y$a     1  -1
#meana   3   8
#meanb 100  50

我想要的是：

#y$a     1   0   
#meana   3   0
#meanb   100 50


#e.g. the results in the first row are from lm((mylist1[[1]][,1])~(mylist2[[1]][,1]))  and lm((mylist1[[2]][,1])~(mylist2[[2]][,1]))
4

1 回答 1

2

我运行了你的代码并得到了

> result <- mapply(myfunction,x=mylist1,y=mylist2)
> result
      part1 part2
y$a       1     0
meana     3     0
meanb   100    50

你有错字

result <- mapply(myfunction,x=mylist,y=mylist2)

我改成

result <- mapply(myfunction,x=mylist1,y=mylist2)

也许这就是问题

于 2012-06-16T10:47:42.157 回答