<html>
<head>
<meta content="text/html; charset=UTF-8" http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
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>
<br>
<div class="moz-cite-prefix">On 23-08-2014 19:21, Andre Oliveira
wrote:<br>
</div>
<blockquote
cite="mid:1408832494.5254.YahooMailNeo@web125901.mail.ne1.yahoo.com"
type="cite">
<div style="color:#000; background-color:#fff;
font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial,
Lucida Grande, sans-serif;font-size:10pt">
<div style="" class=""><span style="" class="">Walmes,</span></div>
<div class="" style="color: rgb(0, 0, 0); font-size: 13.3333px;
font-family: HelveticaNeue,Helvetica
Neue,Helvetica,Arial,Lucida Grande,sans-serif;
background-color: transparent; font-style: normal;"><span
style="" class="">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 class="" style="color: rgb(0, 0, 0); font-size: 13.3333px;
font-family: HelveticaNeue,Helvetica
Neue,Helvetica,Arial,Lucida Grande,sans-serif;
background-color: transparent; font-style: normal;">obrigado <br>
</div>
<div style="" class=""> </div>
<div style="" class="">
</div>
<div style="" class=""><br style="" class="">
<br style="" class="">
</div>
<div style="" class="">
<br style="" class="">
André Oliveira Souza</div>
<div class="qtdSeparateBR"><br>
<br>
</div>
<div style="display: block;" class="yahoo_quoted">
<div class="" style="font-family: HelveticaNeue, Helvetica
Neue, Helvetica, Arial, Lucida Grande, sans-serif;
font-size: 10pt;">
<div class="" style="font-family: HelveticaNeue, Helvetica
Neue, Helvetica, Arial, Lucida Grande, sans-serif;
font-size: 12pt;">
<div style="" class="" dir="ltr"> <font style="" class=""
face="Arial" size="2"> Em Sexta-feira, 22 de Agosto de
2014 17:42, walmes . <a class="moz-txt-link-rfc2396E" href="mailto:walmeszeviani@gmail.com"><walmeszeviani@gmail.com></a>
escreveu:<br style="" class="">
</font> </div>
<br style="" class="">
<br style="" class="">
<div style="" class="">
<div style="" class="" id="yiv8663867924">
<div style="" class="" dir="ltr">
<div class="" 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 class="" style="font-family:trebuchet ms,
sans-serif;"><br style="" class="">
</div>
<div class="" style="font-family:trebuchet ms,
sans-serif;">À disposição.</div>
<div class="" style="font-family:trebuchet ms,
sans-serif;">
Walmes.</div>
<div class="" style="display:inline;"></div>
</div>
</div>
<br style="" class="">
_______________________________________________<br
style="" class="">
R-br mailing list<br style="" class="">
<a moz-do-not-send="true" style="" class=""
ymailto="mailto:R-br@listas.c3sl.ufpr.br"
href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br
style="" class="">
<a moz-do-not-send="true" style="" class=""
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
style="" class="">
Leia o guia de postagem (<a moz-do-not-send="true"
style="" class=""
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 style=""
class="">
<br style="" class="">
</div>
</div>
</div>
</div>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">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">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">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
</body>
</html>