通过aaa得到dataframe bbb(4行6列)

dxystata 发表于 2014-5-22 13:49
程序运行有误，函数内部还有group1 和 group2

嗯，粗心了
你照着原样改下就好了

另外，我不知道你为什么要自己做个函数，用现成的不好吗？
如果想要写个更出色的
可以先看下table()函数的源代码，然后在它的基础上进行改进比较好

yywan0913 发表于 2014-5-22 14:03

多谢，这个设置我都用了，没有效果啊
另外我函数里还加了stringsAsFactors=FALSE，也没用，搞不明白了

playmore 发表于 2014-5-22 15:43
多谢，这个设置我都用了，没有效果啊
另外我函数里还加了stringsAsFactors=FALSE，也没用，搞不明白了

放在所有data.frame的前面设置。
之后再aaa=date.frame();   res=data.frame()

这是题目么？
从数据管理的角度，
gcat    grp1  pct1  grp2  pct2   Total (列名)
  A        2   66.67    1   33.33    3
  B        2   66.67    1   33.33    3
Total     4               2              6
这个最后一行的total是毒药。一个干净的数据表，每一行都应该是一个数据，这里A B 和total的性质完全不一样，直接导致对这个数据表进行处理的时候得不到应该有的值。
如果不用得到最后一行的话，直接用 reshape2库中的 dcast函数直接完成你要的转换。

直接用 reshape2 库中的 dcast 函数就可以。

另外，
gcat    grp1  pct1  grp2  pct2   Total (列名)
  A        2   66.67    1   33.33    3
  B        2   66.67    1   33.33    3
Total     4               2              6

这个表的最后一行非常没有必要。从一个干净数据库的角度，total这一行破坏了数据表的完整性，使得这个表不能在此使用。 比如你summary 这个新的table，得到的 grp1有8个。

playmore 发表于 2014-5-22 15:42
嗯，粗心了
你照着原样改下就好了

如何看table()函数的源代码？谢谢！

playmore 发表于 2014-5-22 13:24
写成函数形式了
另外还是有factor的问题，有高人看看顺便给解决了吧
我就是觉得R里的factor挺烦人的，也 ...

能否把你的程序修改一下能运行，谢谢！

playmore 发表于 2014-5-22 13:24
写成函数形式了
另外还是有factor的问题，有高人看看顺便给解决了吧
我就是觉得R里的factor挺烦人的，也 ...

能否把你的程序修改一下能运行，谢谢！

dxystata 发表于 2014-5-22 22:36
能否把你的程序修改一下能运行，谢谢！

你直接运行下面的代码，然后看bbb是什么

1. options(stringsAsFactors = FALSE)
   
2. 
   
3. tttab <- function(ds,var="group1",group="group2"){
   
4.   
   
5.   stopifnot(is.data.frame(ds))
   
6.   
   
7.   res <- data.frame(gcat=unique(ds[,group]),stringsAsFactors=FALSE)
   
8.   
   
9.   for (i in 1:length(res$gcat)){
   
10.     grp <- unique(ds[,var][which(ds[,group]==res$gcat[i])])
    
11.    
    
12.     for (j in 1:length(grp)){
    
13.       res[,paste("grp",j,sep="")] = grp[j]
    
14.       res[,paste("pct",j,sep="")] = eval(parse(text=paste("sum(with(ds,",var,"==j & ",group,"==res$gcat[i])) / sum(with(ds,",group,"==res$gcat[i]))",sep="")))
    
15.     }
    
16.   }
    
17.   
    
18.   res$Total <- eval(parse(text=paste("sum(with(ds,",group,"==res$gcat[i]))",sep="")))
    
19.   
    
20.   ttl <- apply(res[,-1],2,sum)
    
21.   ttl$gcat <- "Total"
    
22.   
    
23.   res <- rbind(res,ttl)
    
24.   
    
25.   return(res)
    
26. }
    
27.   
    
28. aaa <- data.frame(group1=c(1,1,2,2,1,1),group2=c("A","B","A","B","A","B"))
    
29. 
    
30. bbb <- tttab(ds=aaa,var="group1",group="group2")

复制代码

dxystata 发表于 2014-5-22 20:52
如何看table()函数的源代码？谢谢！

在R里直接table回车