<div dir="ltr">Nilson,<div><br></div><div>O que você deseja fazer tem um nome que nossos antepassados que já enfrentaram esse problema deram: "construção da ROC pelo método bi-normal".</div><div><br></div><div>Por transitividade matemática, você pode descobrir a curva logística que identifica o valor da sua variável que discrimina melhor as duas condições.</div><div><br></div><div>Note que embora do ponto de vista <u>matemático</u> esses procedimentos todos podem ser interpretados e compreensíveis (<i>in extremis</i> você poderia fazer uma busca do valor da variável que maximize ou um alfa ou um beta de um teste de hipótese versus as duas subdistribuições), usar isso como busca <i>post hoc</i> nos seus dados para provar uma hipótese pode ser considerado inválido.</div><div><br></div><div>Se seu problema for de classificação, você precisa se proteger usando as técnicas desse campo de estudo, como separação de parte dos dados para treino e outra para teste, etc.</div><div><br></div><div>HTH</div><div>--</div><div>Cesar Rabak</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Nov 1, 2019 at 8:50 PM Nilson Guiotoku por (R-br) <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Pessoal, <div><br></div><div>eu tenho uma variável que segue uma distribuição bimodal, e eu preciso encontrar o melhor ponto que separa essas bimodais, de forma lógica matematicamente. </div><div><br></div><div>Então eu pensei em gerar a curva de densidade e encontrar a derivada dessa curva. O local em que a derivada for zero eu teria o valor critico e também o ponto de mínimo, portanto, o melhor valor para determinar o corte e encontrar os dois grupos distintos. </div><div><br></div><div>Como encontrar a derivada primeira e segunda da densidade?</div><div><br></div><div># Código mínimo executável:</div><div><br></div><div>X = c(rnorm(100, 10, 3), rnorm(50, 40, 5))  # Variável bimodal<br></div><div>hist(X, freq = F)<br>denfuncao = density(X)<br>lines(denfuncao, lwd = 2)</div><div>  # como encontrar a derivada de uma curva se não conheço a função de </div><div>  # densidade?<br>deriv(expr = denfuncao$x, namevec = denfuncao$x)  <br></div></div>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br>
</blockquote></div>