<font size="3" face="Times New Roman">

</font><p style="margin: 0cm 0cm 10pt;" class="MsoNormal"><font size="3"><font face="Calibri">Benilton,</font></font></p><div style="margin: 0cm 0cm 10pt;" class="MsoNormal"><font size="3"><font face="Calibri">vou tentar
explicar melhor o que eu preciso.</font></font></div><div style="margin: 0cm 0cm 10pt;" class="MsoNormal"><font size="3"><font face="Calibri">Eu tenho um arquivo com as variáveis codigo, peso, val,
pond, mes e ano. Esse arquivo de +/- 1.500.000 observações, tem <strong>n</strong> 
observações para cada código por mês e ano. Na minha rotina eu esqueço a
informação do mês, pois eu quero essa rotina por ano. Eu preciso calcular os
quartis (1º e 3º) por código da variável pond ponderados pela variável val. Em
seguida preciso calcular a diferença interquartílica e os limites (superior e
inferior) com os resultados dos quartis. Depois, junto esse arquivo com os
limites com o arquivo com todas as informações e retiro dele as observações
cujo pond não esteja entre os limites. Executo essa<span>  </span>rotina 2 vezes.</font></font></div><div style="margin: 0cm 0cm 10pt;" class="MsoNormal"><font size="3"><font face="Calibri">Estou mandando os novos links com o código e os dados.</font></font></div>
<div style="margin: 0cm 0cm 10pt;" class="MsoNormal"><a href="http://www.datafilehost.com/download-7fe849b1.html"><font color="#0000ff" size="3" face="Calibri">http://www.datafilehost.com/download-7fe849b1.html</font></a></div>
<div style="margin: 0cm 0cm 10pt;" class="MsoNormal"><font size="3" face="Times New Roman">

</font><p style="margin: 0cm 0cm 10pt;" class="MsoNormal"><a href="http://www.datafilehost.com/download-0437cad6.html"><font color="#0000ff" size="3" face="Calibri">http://www.datafilehost.com/download-0437cad6.html</font></a></p>
<font size="3" face="Times New Roman">

</font></div><div style="margin: 0cm 0cm 10pt;" class="MsoNormal"><font face="Times New Roman"></font> </div><div style="margin: 0cm 0cm 10pt;" class="MsoNormal">Em 22 de setembro de 2011 10:51, Carlos Mendonça <span dir="ltr"><<a href="mailto:csaeslpv@centroin.com.br" target="_blank">csaeslpv@centroin.com.br</a>></span> escreveu:<br>
</div><div class="gmail_quote">
<blockquote style="margin: 0px 0px 0px 0.8ex; padding-left: 1ex; border-left-color: rgb(204, 204, 204); border-left-width: 1px; border-left-style: solid;" class="gmail_quote">
<br><br><div class="gmail_quote"><div>Em 31 de agosto de 2011 09:30, Carlos Mendonça <span dir="ltr"><<a href="mailto:csaeslpv@centroin.com.br" target="_blank">csaeslpv@centroin.com.br</a>></span> escreveu:<br>
</div><blockquote style="margin: 0px 0px 0px 0.8ex; padding-left: 1ex; border-left-color: rgb(204, 204, 204); border-left-width: 1px; border-left-style: solid;" class="gmail_quote">
<div>Caros, estou utilizando a função abaixo que calcula os quartis ponderados que funcionou muito bem até agora. </div><div><div></div><div><div> </div><div>tmp = subset(arq1, select=c(codigo, qtd, val))</div>
<div>equantileByCounts <- function(x, counts, qs){<br>

    tot <- sum(counts)<br>    i <- order(x)<br>    x <- x[i]<br>    counts <- counts[i]<br>    partial <- cumsum(counts)<br>    qsObs <- qs*tot<br>    start <- floor(qsObs)<br>    end <- ceiling(qsObs)<br>




    i <- length(qs)<br>    aux1 <- sapply(1:i,<br>                  function(.x, start, end, x){<br>                      coord <- c(start[.x], end[.x])<br>                      bins <- c(1, partial)<br>                      grps <- cut(coord, bins, labels=FALSE, include.lowest=T)<br>




                      mean(x[grps])<br>                  }, start=start, end=end, x=x)<br>    names(aux1) <- names(qs)<br>    return(aux1)<br>}</div><div>aux1 <- do.call(rbind,<br>               lapply(split(tmp, tmp$codigo),<br>




                      function(y)<br>                      equantileByCounts(y$qtd, y$val, c(.25, .75))))</div><div> </div><div> </div><div>O meu processamento é por trimestres e em apenas um deles deu a seguinte mensagem de erro:</div>




<div> </div><div> "Erro em cut.default(coord, bins, labels = FALSE, include.lowest = T) : <br>  'breaks' não são únicos"<br></div><div>Alguém saberia me dizer o que posso fazer?<br clear="all"><br></div>




<div>                Um abraço,</div>
<div> </div><font color="#888888">
<div>                Carlos Mendonça.</div><br>
</font></div></div></blockquote></div><font color="#888888"><br><br clear="all"><br>-- <br><div>                Um abraço,</div>
<div> </div>
<div>                Mendonça</div><br>
</font></blockquote></div><br><br clear="all"><br>-- <br><div>                Um abraço,</div>
<div> </div>
<div>                Mendonça</div><br>