3

apply最近、2つのリストで関数を使用することについて質問しました。各リストは、大きなデータフレームを分割して作成されたデータフレームのリストです。関数が実行されるたびに、の最初の要素（データフレーム）からベクトルを取得し、の最初の要素（データフレーム）mylist1からいくつかのベクトルを取得してmylist2、それらを相互に回帰します。次に、次のmylist1要素とmylist2要素に移動します。事実上、この関数は同じ数の要素を持つ2つのリストを取得し、ペア（各リストから1つ）を取得してそれらを操作します。

次のことを試しましたが、得られる結果は私が望むものではありません。

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 に答える