group by in R, ddply con weighted.mean

Question

Sto provando a fare una media ponderata in stile "group by" in R. Con alcuni media di base il codice seguente (usando il pacchetto plyr di Hadley) ha funzionato bene.

ddply(mydf,.(period),mean) 

Se uso lo stesso approccio con weighted.mean ottengo il seguente errore " 'x' e 'w' deve avere la stessa lunghezza", che non capisco perché la parte weighted.mean lavora fuori ddply.

weighted.mean(mydf$mycol,mydf$myweight) # works just fine 
ddply(mydf,.(period),weighted.mean,mydf$mycol,mydf$myweight) # returns the erros described above 
ddply(mydf,.(period),weighted.mean(mydf$mycol,mydf$myweight)) # different code same story 

ho pensato di scrivere una funzione personalizzata invece di utilizzare weighted.mean e poi passarlo a ddply o anche scrivere qualcosa di nuovo da zero con sottoinsieme. Nel mio caso sarebbe troppo lavoro spero, ma ci dovrebbe essere una soluzione più intelligente con quello che c'è già.

thx per eventuali suggerimenti in anticipo!

Answer 1

Utilizzare una funzione anonima:

> ddply(iris,"Species",function(X) data.frame(wmn=weighted.mean(X$Sepal.Length, 
+                X$Petal.Length), 
+            mn=mean(X$Sepal.Length))) 
    Species  wmn mn 
1  setosa 5.016963 5.006 
2 versicolor 5.978075 5.936 
3 virginica 6.641535 6.588 
> 

Questo calcola una media pesata di Sepal.Length (ponderato per Petal.Length) così come media non ponderata e restituisce entrambi.

Answer 2

Usa riassumere (o riassumere):

ddply(iris, "Species", summarise, 
    wmn = weighted.mean(Sepal.Length, Petal.Length), 
    mn = mean(Sepal.Length))