Não foi possível enviar o arquivo. Será algum problema com as permissões?
Diferenças

Diferenças

Aqui você vê as diferenças entre duas revisões dessa página.

Link para esta página de comparações

Ambos lados da revisão anterior Revisão anterior
Próxima revisão
Revisão anterior
Próxima revisão Ambos lados da revisão seguinte
pessoais:eder [2011/06/05 19:56]
eder [section 4]
pessoais:eder [2011/06/05 21:13]
eder [section 5]
Linha 12: Linha 12:
   *  Estatística Espacial   *  Estatística Espacial
   * [[http://​www.leg.ufpr.br/​doku.php/​projetos:​gem2|GEM²]] Grupo de estudos em modelos mistos   * [[http://​www.leg.ufpr.br/​doku.php/​projetos:​gem2|GEM²]] Grupo de estudos em modelos mistos
 +===== Disciplinas 2011/1 ===== 
 +  * [[http://​www.leg.ufpr.br/​doku.php/​disciplinas:​ce210-2010-02|CE-210:​ Inferência estatística II]]
 +  * [[http://​www.leg.ufpr.br/​doku.php/​disciplinas:​ce718|CE-718:​ Métodos Computacionalmente Intensivos]]
 +
  
 ===== Minicursos =====  ===== Minicursos ===== 
Linha 74: Linha 78:
 MCcirculo(1,​seq(5,​5000,​by=1000),​plotS=FALSE) MCcirculo(1,​seq(5,​5000,​by=1000),​plotS=FALSE)
 ###​-----------------------------------------------------------------###​ ###​-----------------------------------------------------------------###​
-### inversão de p 
 ### Inversão de Probabilidade ### Inversão de Probabilidade
 +### OBJ: gerar x~exp transformando de uma uniforme
 NS <- 10000 NS <- 10000
 +lam <- 0.5
 +#​f(x)=exp(lam) F(x)=1-exp(-lam*x),​ logo: F^-1(x)= -lam^-1*log(1-x)
 +Gexp <- function(x,​lam){-(log(1-U))/​lam}
 +
 U <- runif(NS) U <- runif(NS)
-X <- - log(U) +X <- Gexp(U,lam
-Y <- rexp(NS)+Y <- rexp(NS,lam) 
 par(mfrow=c(1,​3)) par(mfrow=c(1,​3))
 hist(U,​freq=FALSE,​main='​Uniforme',​col='​lightblue'​) hist(U,​freq=FALSE,​main='​Uniforme',​col='​lightblue'​)
 lines(density(U),​col='​red',​lwd=2) lines(density(U),​col='​red',​lwd=2)
 +
 hist(X,​freq=FALSE,​main='​Expoencial via uniforme',​col='​lightblue'​) hist(X,​freq=FALSE,​main='​Expoencial via uniforme',​col='​lightblue'​)
 lines(density(X),​col='​red',​lwd=2) lines(density(X),​col='​red',​lwd=2)
-lines(curve(dexp(x,​1),​min(X),​max(X),​add=TRUE),​col='​blue',​lwd=2)+lines(curve(dexp(x,​lam),​min(X),​max(X),​add=TRUE),​col='​blue',​lwd=2) 
 hist(Y,​freq=FALSE,​main='​Expoencial do R',​col='​lightblue'​) hist(Y,​freq=FALSE,​main='​Expoencial do R',​col='​lightblue'​)
 lines(density(Y),​col='​red',​lwd=2) lines(density(Y),​col='​red',​lwd=2)
-lines(curve(dexp(x,​1),​min(Y),​max(Y),​add=TRUE),​col='​blue',​lwd=2)+lines(curve(dexp(x,​lam),​min(Y),​max(Y),​add=TRUE),​col='​blue',​lwd=2)
 ###​-----------------------------------------------------------------###​ ###​-----------------------------------------------------------------###​
 ### Regressão Beta ### Regressão Beta

QR Code
QR Code pessoais:eder (generated for current page)