<div>
<div>Obrigado colega!<br /><br />E aproveitando o insejo, encontrei um artigo que aborta sobre Modelos de Mistura, o qual disponibiliza a rotina R utilizada no mesmo, entretanto, como ainda não estou bem familiarizado com o assunto, não entendi qual o objetivo da rotina. Mas de qualquer forma, a rotina é a seguinte:<br /><br />
<div>library(mixtools) # mixtools package must be installed first</div>
<div>set.seed(123) # Ensure that results are exactly reproducible</div>
<div>#Generate data:</div>
<div>mu <- matrix(c(0, 15), 2, 3)</div>
<div>sigma <- matrix(c(1, 5), 2, 3)</div>
<div>x <- rmvnormmix(300, lambda = c(.4,.6), mu = mu, sigma = sigma)</div>
<div># npEM algorithm results:</div>
<div>a <- npEM(x, mu0 = 2, blockid = rep(1,3), samebw = TRUE)</div>
<div>b <- npEM(x, mu0 = 2, blockid = rep(1,3), samebw = FALSE)</div>
<div># Produce plots like in Figure 1:</div>
<div>s <- seq(-10, 40, len = 200)</div>
<div>plot(a, xlim=c(-10, 40), ylim = c(0, .16), xlab = "", breaks = 30)</div>
<div>lines(s, dnorm(s)*.4 + dnorm((s-15)/5)/5*.6, lwd = 2, lty = 2)</div>
<div>plot(b, xlim=c(-10, 40), ylim = c(0, .16), xlab = "", breaks = 30)</div>
<div>lines(s, dnorm(s)*.4 + dnorm((s-15)/5)/5*.6, lwd = 2, lty = 2)</div>
<div># Display npEM bandwidths, minimum values of max_j p_{ij}</div>
<div>a$bandwidth</div>
<div>b$bandwidth</div>
<div>min(apply(a$posteriors, 1, max))</div>
<div>min(apply(b$posteriors, 1, max))</div>
<div># spEM algorithm results:</div>
<div>d <- spEM(x, mu0 = 2, blockid = rep(1,3), constbw = FALSE)</div>
<div>d2 <- spEM(x, mu0 = 2, blockid = rep(1,3), constbw = TRUE)</div>
<div># Produce plot like in Figure 2:</div>
<div>plot(d, xlim=c(-10, 40), ylim = c(0, .16), xlab = "", breaks = 30,</div>
<div>addlegend=FALSE)</div>
<div>plot(d2, newplot=FALSE, addlegend=FALSE, lty=2, dens.col=1)</div>
<div># Display spEM bandwidths</div>
<div>d$bandwidth</div>
<div>d2$bandwidth<br /><br />Caso se interesse pelo artigo, posso enviá-lo nesse mesmo e-mail.<br /><br /><em>Att.</em><br /><em>André</em> </div>
</div>
</div>
<hr style="border-top: 1px solid #ccc;" />
<div><br /><strong>De:</strong> felipe.e.barletta@gmail.com<br /><strong>Enviada:</strong> Sábado, 23 de Agosto de 2014 19:58<br /><strong>Para:</strong> r-br@listas.c3sl.ufpr.br<br /><strong>Assunto:</strong> [R-br] Dúvida Simples - Histograma<br /><br /> Prezado Andre,<br /> <br /> Pode usar a estimação Kernel que é um método não-paramétrico para estimação de densidades que pode ser visto como um aperfeiçoamento do histograma pois gera uma função suavizada.<br /> Precisa defir a função Kernel a ser utilizada e o tamanho da janela(h) ou tamanho da largura da banda(binwidth).<br /> O fator preponderante é o valor de h, pois este define o formato da curva estimada, e sua escolha é arbitrária.<br /> As funções Kernel mais conhecidas são:<br /> Uniforme<br /> Tricúbica<br /> Gaussiana<br /> Epanechnikov<br /> <br /> #############Exemplos no R#####################################################<br /> <br /> ##Exemplo 1<br /> data(precip)<br /> hist(precip, prob=T)<br /> precip.d <- density(precip)<br /> lines(precip.d)<br /> <br /> <br /> ##Exemplo2<br /> plot(density(precip,n=1000))<br /> rug(precip)<br /> lines(density(precip,bw='nrd'),col='red')<br /> lines(density(precip,bw='ucv'),col='green')<br /> lines(density(precip,bw='SJ-ste'),col='yellow')<br /> lines(density(precip,bw='SJ-dpi'),col='blue')<br /> <br /> <br /> plot(density(precip,n=1000))<br /> rug(precip)<br /> lines(density(precip,bw=4),col='red')<br /> lines(density(precip,bw=5),col='green')<br /> lines(density(precip,bw=6),col='yellow')<br /> lines(density(precip,bw=7),col='blue')<br /> <br /> ##Exemplo 3<br /> #Estimação de densidades<br /> #bandwidth=largura da janela<br /> data(faithful)<br /> attach(faithful)<br /> hist(eruptions,breaks="scott",prob=T)<br /> den<-density(eruptions)<br /> lines(den)<br /> den # x valor do tempo de erupção, y valor de densidade estimado, bw ótimo, não precisa se preocupar para achá-lo, o próprio pacote o encontra.<br /> names(den)<br /> summary(eruptions)<br /> summary(den$x)<br /> <br /> <br /> ##Exemplo 4<br /> require(MASS)<br /> data(Boston)<br /> attach(Boston)<br /> hist(rm, prob=T)<br /> rm.d <- density(rm)###### Kernel gaussiano, o defaut do R <br /> rm.d$bw###### tamanho da largura da banda<br /> rug(rm)<br /> lines(rm.d)<br /> curve(dnorm(x,mean(rm),sd(rm)),0,10,add=T,col="blue")##comando para jogar a curva da normal<br /> <br /> </div>
<div class="moz-cite-prefix">On 23-08-2014 19:21, Andre Oliveira wrote:</div>
<blockquote cite="mid:1408832494.5254.YahooMailNeo@web125901.mail.ne1.yahoo.com">
<div style="color: #000; background-color: #fff; font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial,
Lucida Grande, sans-serif; font-size: 10pt;">
<div><span>Walmes,</span></div>
<div style="color: #000000; font-size: 13.3333px; font-family: HelveticaNeue,Helvetica
Neue,Helvetica,Arial,Lucida Grande,sans-serif; background-color: transparent; font-style: normal;"><span>sei que tem várias funções para </span>construir os gráficos de densidade empírica, mas qual função você tem usado com esta finalidade no momento e que tem dado melhores resultados?</div>
<div style="color: #000000; font-size: 13.3333px; font-family: HelveticaNeue,Helvetica
Neue,Helvetica,Arial,Lucida Grande,sans-serif; background-color: transparent; font-style: normal;">obrigado </div>
<div> </div>
<div> </div>
<div><br /> </div>
<div><br /> André Oliveira Souza</div>
<div class="qtdSeparateBR"><br /> </div>
<div class="yahoo_quoted" style="display: block;">
<div style="font-family: HelveticaNeue, Helvetica
Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 10pt;">
<div style="font-family: HelveticaNeue, Helvetica
Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 12pt;">
<div dir="ltr"><span style="font-family: Arial; font-size: small;"> Em Sexta-feira, 22 de Agosto de 2014 17:42, walmes . <a class="moz-txt-link-rfc2396E" href="../../../undefined//compose?to=walmeszeviani@gmail.com" target="_blank"><walmeszeviani@gmail.com></a> escreveu:<br /> </span></div>
<br /> <br />
<div>
<div id="yiv8663867924">
<div dir="ltr">
<div style="font-family: trebuchet ms,
sans-serif;">Para ser bem sincero, o histograma, seja com os melhores elementos/recursos, de um pacote ou outro, está começando a ficar ultrapassado e deixando o lugar para os gráficos de densidade empírica. Nas minhas aulas eu tenho dado mais enfase ao gráfico de densidade que o histograma. E a principal razão é bem simples, é quando desejamos comparar a distribuição a partir de amostras de diferentes populações. Com histogramas há uma certa dificuldade de sobrepor na mesma janela os resultados para mais de duas amostras porque os retângulos, sem o devido cuidado, podem esconder outros de menor altura atrás deles, comprometendo a visualização. Claro que é possível usar retângulo com semitransparência mas não se pode assumir esse conhecimento de todo usuário. Com os gráficos de densidade, que traçam as linhas da função de densidade empírica, permitem sobreposição mais facilmente. Outras coisa interessante é a facilidade na identificação da moda da distribuição como sendo o pico da função traçada. No histograma é fácil achar a classe modal, mas um valor para a moda é menos óbvio, apesar de ser possível interpolar considerando a frequência das classes vizinhas. O bom do histograma é que se pode cobrar em provas, basta o aluno saber classificar dados e obter frequências ao passo que a densidade empírica tem aplicabilidade limitada sem computador.</div>
<div style="font-family: trebuchet ms,
sans-serif;"> </div>
<div style="font-family: trebuchet ms,
sans-serif;">À disposição.</div>
<div style="font-family: trebuchet ms,
sans-serif;">Walmes.</div>
<div style="display: inline;"></div>
</div>
</div>
<br /> _______________________________________________<br /> R-br mailing list<br /> <a href="../../../undefined//compose?to=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" 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" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br /> </div>
</div>
</div>
</div>
</div>
<br /><fieldset class="mimeAttachmentHeader"></fieldset><br />
<pre>_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="../../../undefined//compose?to=R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<div><br /> _______________________________________________<br />R-br mailing list<br />R-br@listas.c3sl.ufpr.br<br /><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" 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" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</div>