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

1. aaa<- data.frame(group1=c(1,1,2,2,1,1),
   
2.                  group2=c("A","B","A","B","A","B"))

复制代码

写一个函数通过aaa得到dataframe bbb(4行6列)：
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


注：要求采用sum(with(aaa, group1==1 & group2=="A")) 实现。
    不能采用R中已有产生四格表的函数。

dxystata 发表于 2014-5-23 09:35
你能运行成功吗？我这里运行报错！能否把code复制到txt中，发到我的邮箱，邮箱pm你。谢谢！

嗯，我这里看到的贴在网页上的代码有问题
有些$显示不正常，应该是论坛的问题
我这回贴附件里了，你再试试吧



test.txt (889 Bytes)

2014-5-23 14:02:09 上传

代码

继续悬赏！

太高级了，没看懂，完全不知道那个4行6列是神马样子……

1. aaa<- data.frame(group1=c(1,1,2,2,1,1),group2=c("A","B","A","B","A","B"))
   
2. 
   
3. res <- data.frame(gcat=unique(aaa$group2),stringsAsFactors=FALSE)
   
4. 
   
5. for (i in 1:length(res$gcat)){
   
6.   grp <- unique(aaa$group1[which(aaa$group2==res$gcat)])
   
7.   
   
8.   for (j in 1:length(grp)){
   
9.     res[,paste("grp",j,sep="")] = grp[j]
   
10.     res[,paste("pct",j,sep="")] = sum(with(aaa, group1==j & group2==res$gcat)) / sum(with(aaa, group2==res$gcat))
    
11.   }
    
12. }
    
13. 
    
14. res$Total <- sum(with(aaa, group2==res$gcat))
    
15. 
    
16. ttl <- apply(res[,-1],2,sum)
    
17. ttl$gcat <- "Total"
    
18. 
    
19. res <- rbind(res,ttl)

复制代码

基本上算是完成了任务，只不过左下角那个Total实在是搞不定
不知道怎么让R取消认定gcat是factor，之前加了stringsAsFactors也没用

playmore 发表于 2014-5-22 10:05
aaa

谢谢！希望调用这样的函数tttab(ds=aaa,var="group1",group="group2")实现。

playmore 发表于 2014-5-22 10:05
aaa

建议代码用

1.      

复制代码

，显示出来好看一点，谢谢！

dxystata 发表于 2014-5-22 11:25
建议代码用，显示出来好看一点，谢谢！

写成函数形式了
另外还是有factor的问题，有高人看看顺便给解决了吧
我就是觉得R里的factor挺烦人的，也看不出来有什么好处

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

复制代码

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

程序运行有误，函数内部还有group1 和 group2

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

1. options(stringsAsFactors = FALSE)

复制代码