Como posso fazer para achar os quartis e a diferença interquartílica, por uma determinada variável e ponderados por uma outra.

Por exemplo, nos dados a seguir, eu queria os quartis da variável valor, por código e ponderados pela variável quantidade.

Data frame de entrada
    codigo        valor       quantidade
    1.01     35867,14             125
    1.01      5867,05            1045
    1.01     62100,14               91
    1.01     35000,99               21
    1.02         401,87           2650
    1.02         815,7           10345

Não sei se ajuda, mas no SAS eu utilizava o seguinte código

PROC UNIVARIATE NOPRINT DATA = arq1 ;

BY codigo ;

VAR valor ;

WEIGHT quantidade ;

OUTPUT OUT = arq2 Q3=Q3 Q1=Q1 ;

Obrigado,

Carlos Mendonça.