2011-04-15 86 views
2

出于某种原因,我得到的结果比我预期的更高,因为升级到R-2.13.0 - 并且升级到plyr_1.5.1.tar.gz ...我在旧版本上尝试过版本plyr的(版本不确定不幸的是,因为我已经仅覆盖了吧...)ddply返回的结果太多

library(plyr) 
dd <-data.frame(matrix(rnorm(216),72,3),c(rep("A",24),rep("B",24), 
    rep("C",24)),c(rep("J",36),rep("K",36))) 
colnames(dd) <- c("v1", "v2", "v3", "dim1", "dim2") 

results1 <- ddply(dd, c("dim1","dim2"), function(df) c(m1=mean(df$v1))) 
results2 <- ddply(dd, c("dim1","dim2"), function(df) { c(m1=mean(df$v1), 
    m2=mean(df$v2)) }) 
results3 <- ddply(dd, c("dim1","dim2"), function(df) { c(m1=mean(df$v1), 
    m2=mean(df$v2), m3=mean(df$v3)) }) 

我不明白,为什么结果2有结果1行数的两倍,并有结果3的三倍多 - 其中原始结果1仅复制两次或三次。

我使用的是旧版本plyr的,我期待为结果[R版本,方便的拷贝2.11.0修补的(2010-05-01 r51907)...

> results1 
    dim1 dim2   m1 
1 A J 0.07312783 
2 B J -0.22428746 
3 B K -0.44205832 
4 C K 0.21421456 
> results2 
    dim1 dim2   m1   m2 
1 A J 0.07312783 -0.1130148 
2 B J -0.22428746 0.4394832 
3 B K -0.44205832 -0.1934018 
4 C K 0.21421456 -0.0178809 
> results3 
    dim1 dim2   m1   m2   m3 
1 A J 0.07312783 -0.1130148 -0.03175873 
2 B J -0.22428746 0.4394832 0.21581696 
3 B K -0.44205832 -0.1934018 -0.28313530 
4 C K 0.21421456 -0.0178809 -0.21948430 

结果我得到从R版本2.13.0(2011-04-13)

> results1 
    dim1 dim2   m1 
1 A J -0.2270726 
2 B J 0.5860493 
3 B K -0.5986129 
4 C K 0.3135809 
> results2 
    dim1 dim2   m1   m2 
1 A J -0.2270726 -0.19037813 
2 B J 0.5860493 -0.05385395 
3 B K -0.5986129 0.29404095 
4 C K 0.3135809 -0.26744010 
5 A J -0.2270726 -0.19037813 
6 B J 0.5860493 -0.05385395 
7 B K -0.5986129 0.29404095 
8 C K 0.3135809 -0.26744010 
> results3 
    dim1 dim2   m1   m2   m3 
1  A J -0.2270726 -0.19037813 -0.20448734 
2  B J 0.5860493 -0.05385395 -0.11190857 
3  B K -0.5986129 0.29404095 -0.27072101 
4  C K 0.3135809 -0.26744010 -0.03184949 
5  A J -0.2270726 -0.19037813 -0.20448734 
6  B J 0.5860493 -0.05385395 -0.11190857 
7  B K -0.5986129 0.29404095 -0.27072101 
8  C K 0.3135809 -0.26744010 -0.03184949 
9  A J -0.2270726 -0.19037813 -0.20448734 
10 B J 0.5860493 -0.05385395 -0.11190857 
11 B K -0.5986129 0.29404095 -0.27072101 
12 C K 0.3135809 -0.26744010 -0.03184949 

为什么有results2得到8行而不是4和results3得到12行而不是4?

感谢, 肖恩

回答

5

这将很快在plyr 1.5.2

+1

+1谢谢... – Andrie 2011-04-15 13:36:02

+0

正常工作1.4.1万一别人刚刚升级的一切太:http://cran.r-project.org/src/contrib/Archive/plyr/plyr_1.4.1.tar.gz – 2011-04-17 01:20:43

5

固定这是你的ddply内的C()()函数是造成问题。

这里有三种不同的方法,你可以写你的结果3陈述,逐渐变得简单了:

  1. 使用data.frame自己的函数中:

    ddply(dd, c("dim1","dim2"), function(df) {data.frame(m1=mean(df$v1), m2=mean(df$v2), m3=mean(df$v3)) })

  2. 使用总结:

    ddply(dd, .(dim1, dim2), summarise, m1=mean(v1), m2=mean(v2), m3=mean(v3))

  3. 使用numcolwise。

    ddply(dd, .(dim1, dim2), numcolwise(mean))

在每一种情况下,结果是你所期望的:

dim1 dim2   m1   m2   m3 
1 A J -0.04272659 -0.1468376 0.17902942 
2 B J -0.10133503 -0.1427358 -0.05241214 
3 B K 0.29698847 -0.0989732 0.14422812 
4 C K 0.04108324 0.2014864 -0.15893221