2012-11-02 55 views
0

我有几分定义的数据帧如下(有真200+变量)熔融:熔体/ dcast /与边缘不工作

class_level var1 var2 
      1 4 7 
      1 6 7 
      1 7 6 
      4 3 1 
      4 1 3 

其中class_level是1或4,和VAR1和VAR2已值为1-7,这是来自调查的likert-like响应级别。

我想,其中包括数以供类级别每个变种水平的数据帧,以总计行每个变量,像:

class_level variable Levels Students 
      1  var1  1  10 
      1  var1  2  7 
      1  var1  3  28 
      1  var1  4  15 
      1  var1  5  54 
      1  var1  6  38 
      1  var1  7  16 
      1  var1 Total  168 
      4  var1  1  58 
      .  .  .  . 
      .  .  .  . 
      .  .  .  . 
      4  var1  7  33 
      4  var1 Total  294 

我曾尝试以下:

df.m <- melt(df, id.vars=c("class_level"), na.rm=TRUE) 
head(df.m) 
# class_level variable value 
#1   4  var1  4 
#2   4  var1  6 
#3   1  var1  7 
#4   4  var1  3 
#5   1  var1  5 
#6   4  var1  6 

df.c <- dcast(df.m, class_level+variable ~ value, 
       fun.aggregate=length, 
       subset=.(variable %in% c("var1","var2")), 
       margins=TRUE 
      ) 
head(df.c) 
# class_level variable 1 2 3 4 5 6 7 (all) 
#1   1  var1 1 1 8 24 56 101 32 223 
#2   1  var2 2 4 4 22 49 79 56 216 
#3   4  var1 4 5 11 38 91 114 76 339 
#4   4  var2 2 6 11 35 73 98 106 331 

df.o <- melt(df.c, id.vars=c("class_level","variable"), 
       variable.name="Levels", value.name="Students" 
      ) 
head(df.o) 
# class_level variable Levels Students 
#1   1  var1  1  1 
#2   1  var2  1  2 
#3   4  var1  1  4 
#4   4  var2  1  2 
#5   1  var1  2  1 
#6   1  var2  2  4 

正如您所看到的,这会按每个变量的级别生成计数,但不会生成总计行。如何获取最终数据集(df.o)中的总计行数。任何帮助将不胜感激。

大卫

+0

请提供[重复性好示例](http://stackoverflow.com/questions/5963269/how-to-制作一个伟大的可再现的例子)以及理想的结果。我现在还不清楚你想要做什么。 – Dason

+1

可能你在'dcast'调用中只是缺少一个')'? – joran

+0

@joran - 是的,更正。 – dmonder

回答

1

我会倾向于使用plyr的帮助:

df.m <- melt(df, id.vars=c("class_level"), na.rm=TRUE) 
df.m$value <- factor(df.m$value, levels=1:7) # To ensure 0 counts as well 
df.c <- ddply(df.m, .(class_level, variable), 
      function(x) c(table(df.m$value), Total=length(df$m.value))) 
df.o <- melt(df.c, id.vars=c("class_level", "variable")) 
+0

刚刚看到@ joran的评论:看起来像OP的代码毕竟有效! – seancarmody

+0

啊,但OP会通知...? :) – joran

+0

在你的帮助下,让我们希望如此! – seancarmody