加權表格數據框與plyr

Question

我正在處理由多個問題（y1，y2，y3，...）的整數值響應和分配給每個應答者的加權計數組成的調查數據，如下所示：

foo <- data.frame(wcount = c(10, 1, 2, 3),  # weighted counts 
        y1 = sample(1:5, 4, replace=T), # numeric responses 
        y2 = sample(1:5, 4, replace=T), # 
        y3 = sample(1:5, 4, replace=T)) # 
>foo 
    wcount y1 y2 y3 
1  10 5 5 5 
2  1 1 4 4 
3  2 1 2 5 
4  3 2 5 3 

我想將其轉換爲加權表的統一數據框版本，第一列代表響應值，接下來的3列代表加權計數。

library(Hmisc) 
ty1 <- wtd.table(foo$y1, foo$wcount) 
ty2 <- wtd.table(foo$y2, foo$wcount) 
ty3 <- wtd.table(foo$y3, foo$wcount) 

bar <- merge(ty1, ty2, all=T, by="x") 
bar <- merge(bar, ty3, all=T, by="x") 

names(bar) <- c("x", "ty1", "ty2", "ty3") 
bar[is.na(bar)]<-0 
>bar 
    x ty1 ty2 ty3 
1 1 3 0 0 
2 2 3 2 0 
3 3 0 0 3 
4 4 0 1 1 
5 5 10 13 12 

我懷疑有與plyr和numcolwise或ddply這種自動化的方法：可以明確地列使用來完成。例如，以下內容即將結束，但我不確定還需要完成該任務：

library(plyr) 
bar2 <- numcolwise(wtd.table)(foo[c("y1","y2","y3")], foo$wcount) 
>bar2 
     y1  y2  y3 
1 1, 2, 5 2, 4, 5 3, 4, 5 
2 3, 3, 10 2, 1, 13 3, 1, 12 

有什麼想法？

Answer 1

不是plyr的答案，但是這給我的印象重塑/聚合，可以直截了當地用從包reshape2功能來解決問題。

首先，melt數據集，使響應值的一列可以命名爲x（y1 - y3中的唯一值）。

library(reshape2) 
dat2 = melt(foo, id.var = "wcount", value.name = "x") 

現在，這可以與廣投dcast回，使用sum作爲聚合功能。這將y1 - y3作爲列，總計爲wcount，每個值爲x。

# Cast back wide using the values within y1-y3 as response values 
    # and filling with the sum of "wcount" 
dcast(dat2, x ~ variable, value.var = "wcount", fun = sum) 

給予

x y1 y2 y3 
1 1 3 0 0 
2 2 3 2 0 
3 3 0 0 3 
4 4 0 1 1 
5 5 10 13 12 

Answer 2

你所描述的是使用重複權的一項調查數據集。看到http://asdfree.com/許多很多的例子，但對於recs，做這樣的事情：

library(survey) 

x <- read.csv("http://www.eia.gov/consumption/residential/data/2009/csv/recs2009_public.csv") 
rw <- read.csv("http://www.eia.gov/consumption/residential/data/2009/csv/recs2009_public_repweights.csv") 
y <- merge(x , rw) 

# create a replicate-weighted survey design object 
z <- svrepdesign(data = y , weights = ~NWEIGHT , repweights = "brr_weight_[0-9]") 

# now run all of your analyses on the object `z` .. 
# see the `survey` package homepage for details 

# distribution 
svymean(~ factor(BASEHEAT) , z) 

# mean 
svymean(~ TOTHSQFT , z)