我想知道是否有人知道在R中加载按计数加权的摘要数据的方法,如SAS中的PROC FREQ。
例如在SAS中加载数据:
data q1;
input resid $ stress $ out $ count @@; datalines;
urban low fav 55 urban low unfav 34
urban med fav 66 urban med unfav 64
urban high fav 107 urban high unfav 89
rural low fav 40 rural low unfav 25
rural med fav 57 rural med unfav 47
rural high fav 88 rural high unfav 60
;
run;
允许您通过proc freq:
获得加权数据表PROC FREQ data=q1;
weight count;
tables fav*urban;
run;
或运行逻辑回归,如:
Proc logistic data=q1;
freq count;
class stress(ref='low') resid(ref='urban')/ param=ref;
model out(ref='unfav')=stress resid;
run;
任何指导都是非常感谢的。
使用tidyr::uncount
,您可以执行以下操作来回答上述问题
q1<-tibble(resid=c(rep(c("urban"),6),rep("rural",6)),
stress=c(rep(c("low","low","med","med","high","high"),2)),
out=c(rep(c("fav","unfav"),6)),
n=c(55,34,66,64,107,89,40,25,57,47,88,60))
sum(q1$n)
q1_expanded<-tidyr::uncount(q1,n)
,你最终从汇总数据制作一个可用的数据框架。