jtrecenti

require(dplyr)
require(ggplot2)

# http://atlasbrasil.org.br/2013/dadosbrutos/atlas2013_dados_brutos.xls
# Salve em csv ou use o pacote gdata
dados <- read.csv2('atlas2013_dados_brutos.csv', as.is=T)

bd <- select(dados, ANO, pesotot, Município, IDHM, UFN, RDPC) %.% 
  filter(ANO==2010) %.%
  arrange(desc(pesotot)) %.%

jtrecenti

require(stringr)

pega_do <- function(yyyy,mm,dd) {
  wget <- sprintf('wget --no-check-certificate \'https://www.dje.tjsp.jus.br/cdje/downloadCaderno.do?dtDiario=%02d/%02d/%d&cdCaderno=14\' -O pdf_do/%d_%02d_%02d.pdf', dd, mm, yyyy, yyyy, mm, dd)
  if(!file.exists(sprintf('pdf_do/%d_%02d_%02d.pdf',yyyy,mm,dd))) {
    system(wget, ignore.stdout=T, ignore.stderr=T, wait=T)  
  }
}
pdf2txt_do <- function(arq, verbose=F) {
  if(file.info(arq)$size > 20000 & !file.exists(gsub('pdf','txt',arq))) {

jtrecenti

\documentclass{article}

\usepackage[brazil]{babel}
\usepackage[utf8]{inputenc}


\begin{document}
\SweaveOpts{concordance=TRUE}

\section{Simulação Priest and Klein}

jtrecenti

# Calcula a previsão de um modelo. O mesmo que "fitted", 
# mas aceita um modelo com coeficientes modificados
pred <- function(model) {
  pred <- model.matrix(model) %*% model$coefficients
  pred <- exp(pred) / (1+exp(pred))
  return(pred)
}

# Calcula o valor Z do teste de hipóteses para uma variável no modelo glm
zval <- function(model, v) {

jtrecenti

# Exemplo de amostrador de gibbs

# Ronaldo, quando era magro, marcava um número de gols por mês que seguia distribuição Poisson(2).
# Ao passar do tempo, teve um momento em que o número de gols de Ronaldo mudou de distribuição.
# Essa distribuição continuou sendo poisson, mas mudou o parâmetro...
# Eu acho que essa mudança pode ter acontecido em qualquer mês entre os anos de 1998 e 2010, de forma
# igual, e também ouso dizer que a média dele agora não deve ser muito diferente de uma exponencial(1).
# pergunto: quando foi que o infeliz (mas rico) mudou a forma de jogar, e pra quanto foi a sua média de gols?
# minha função de perda para acertar o mês é 0-1 e minha função de perda para acertar a média é quadrática.

jtrecenti

# Exemplo de amostrador de gibbs

Ronaldo, quando era magro, marcava um número de gols por mês que seguia distribuição Poisson(2). Ao passar do tempo, teve um momento em que o número de gols de Ronaldo mudou de distribuição. Essa distribuição continuou sendo poisson, mas mudou o parâmetro... Eu acho que essa mudança pode ter acontecido em qualquer mês entre os anos de 1998 e 2010, de forma igual, e também ouso dizer que a média dele agora não deve ser muito diferente de uma exponencial(1). Pergunto: quando foi que o infeliz (mas rico) mudou a forma de jogar, e pra quanto foi a sua média de gols? Minha função de perda para acertar o mês é 0-1 e minha função de perda para acertar a média é quadrática.

```{r}
require(ggplot2)
require(dplyr)
```

```{r}

jtrecenti

x <- rnorm(1000)
z <- runif(1000)
d <- data.frame(x = x, z = z, y= x + z/10 + rnorm(1000, sd=.001))
cor(d)
summary(lm(y~x+z, data=d))
summary(lm(y~x, data=d))
summary(lm(y~z, data=d))

jtrecenti

require(abjutils)
require(ggmap)
# isso aqui é um bd que nós já preparamos há algum tempo e deixamos pronto caso alguém precise
data(br_uf_map)

# isso aqui é um bd auxiliar com os estados nas linhas e onde 
# cada coluna é uma variável que vc quer mostrar no mapa
estados <- data.frame(est=c('AC','AL','AM','AP','BA','CE','DF','ES','GO',
                            'MA','MG','MS','MT','PA','PB','PE','PI','PR',
                            'RJ','RN','RO','RR','RS','SC','SE','SP','TO'),

jtrecenti

library(jurimetria)
library(abjutils)
library(ggplot2)

data(br_uf_map)
data(dados_aj)

d_aj <- dados_aj %>% 
  filter(tipo=='Justiça Estadual e DF') %>%
  mutate(nome_vara2=rm_accent(tolower(nome_vara))) %>%

jtrecenti

#' faz o download de uma tabela do Anuário Estatístico da USP
#' secao: numero da secao (1 a 11)
#' id: id da tabela (depende da secao)
#' retorna um data.frame com os dados.
rusp_tab <- function(secao, id) {
  link <- sprintf("https://uspdigital.usp.br/anuario/br/tabelas/XLS/2013/T%d.%02d.xls", secao, id)
  cat('Acessando ', link, '...\n', sep='')
  tmp <- tempfile()
  r <- RCurl::getBinaryURL(link, ssl.verifypeer=FALSE, encoding='latin1')
  writeBin(as.vector(r), tmp)