
23 Ago
2011
23 Ago
'11
03:30
Fiz esta rotina para o domínio 1,infinito. Porém como eu faço para a curva da densidade começar no 1 assim como no histograma? d2=runif(100)^-1/alfa d3=subset(d2,d2>=1) hist(d3,probability=T,nclass=60,col='green') dens=density(d3,n=200,na.rm=T) points(dens,type='l',col=2) abline(v=1,lty=4)