[datatable-help] stratified weighted means

Matthew Dowle mdowle at mdowle.plus.com
Wed Dec 19 18:38:59 CET 2012


Deja vu?!

http://stackoverflow.com/a/13454651/403310

https://r-forge.r-project.org/tracker/index.php?func=detail&aid=2381&group_id=240&atid=975


On 19.12.2012 17:31, David Freedman wrote:
> I apologize for the simple question, but something that worked in v 
> 1.8.4
> fails in v 1.8.6 - it involves getting stratified weighted means from 
> a data
> table.  The code is
>
> library(data.table)
> (dt <- data.table(id     = rep(c(1,2), each=9),
>                  month  = rep(100L:108L, times=2),
>                  value=rnorm(18),
>                  swt=sample(1:18,18) ))
>
> dt[,lapply(.SD,mean),keyby=list(id)] #fine - stratified means
> dt[,lapply(.SD,function(x)weighted.mean(x,swt))] #fine - overall 
> weighted
> mean
> dt[,lapply(.SD,function(x)weighted.mean(x,swt)),keyby=list(id)] 
> #fails
> 'object "swt" not found'
>
> I'm using v 1.8.6 of data.table, and I am sure that this statement 
> for
> calculating stratified weight means worked in v 1.8.4.
>
> Output of sessionInfo is:
> R version 2.15.2 (2012-10-26)
> Platform: i386-w64-mingw32/i386 (32-bit)
>
> locale:
> [1] LC_COLLATE=English_United States.1252  LC_CTYPE=English_United
> States.1252
> [3] LC_MONETARY=English_United States.1252 LC_NUMERIC=C
>
> [5] LC_TIME=English_United States.1252
>
> attached base packages:
> [1] stats     graphics  grDevices utils     datasets  methods   base
>
> other attached packages:
> [1] data.table_1.8.6
>
> loaded via a namespace (and not attached):
> [1] tools_2.15.2
>
> Any help would be greatly appreciated as this is something I do 
> almost
> daily.
>
> And I realize that I could always use something like
> dt[,list(wm(month,swt),wm(value,swt)),keyby=list(id)]
> But this becomes difficult if I want the stratified weighted means 
> for a
> large number of variables
>
> thanks, david freedman
>
>
>
>
>
>
> --
> View this message in context:
> 
> http://r.789695.n4.nabble.com/stratified-weighted-means-tp4653556.html
> Sent from the datatable-help mailing list archive at Nabble.com.
> _______________________________________________
> datatable-help mailing list
> datatable-help at lists.r-forge.r-project.org
> 
> https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/datatable-help


More information about the datatable-help mailing list