ref = ifelse(VALOR > top30, 1, 0)
) %>%
dplyr::select(TÉCNICO, franqueado, safra, nome, MUNICÍPIO, top30, ref) -> ref
banco$bd %>%
left_join(
banco$bd %>%
count(franqueado, nome, FRENTE, INDICADOR), by = c('franqueado', 'nome', 'FRENTE', 'INDICADOR')
) %>%
filter(n > 1) %>%
select(nome, safra, franqueado,id_franquia) %>%
distinct() %>%
group_by(nome,franqueado,id_franquia) %>%
summarise(n = n()) %>%
ungroup %>%
mutate(faz_id = 1:nrow(.)) -> banco[['faz_id_data']]  ## return
banco$bd %<>%
left_join(
banco$faz_id_data %>% select(nome, franqueado, faz_id), by = c('nome', 'franqueado')
)
banco$ref <- banco$bd
return(banco)
}
b <- carregarDados()
r
head(r)
bancos$bd
dados$VALOR
dados %>% mutate(VALOR = as.numeric(VALOR))
dados %>% filter(is.na(VALOR))
dados %>% filter(is.na(VALOR)) %>% select(VALOR)
dados %>% filter(!is.na(VALOR)) %>% mutate(VALOR = as.numeric(VALOR)) %>% select(VALOR)
banco$bd %<>%
filter(!is.na(VALOR)) %>%
mutate(VALOR = gsub("\\D+",  NA, VALOR)) %>%
mutate(VALOR =  as.numeric(VALOR))
dados %<>%
filter(!is.na(VALOR)) %>%
mutate(VALOR = gsub("\\D+",  NA, VALOR)) %>%
mutate(VALOR =  as.numeric(VALOR))
dados
dados$VALOR
dados %<>%
filter(!is.na(VALOR)) %>%
mutate(VALOR = gsub("\\D+",  NA, VALOR)) %>%
mutate(VALOR =  as.numeric(VALOR))
dados %<>%
filter(!is.na(VALOR)) %>%
#mutate(VALOR = gsub("\\D+",  NA, VALOR)) %>%
mutate(VALOR =  as.numeric(VALOR))
dados
dados$VALOR
carregarDados <- function(){
banco <- list()
banco[['ind_top']] <- read_xlsx('bd/ind_top.xlsx')
#load('bd/banco.rData')
banco[['bd']] <- dados %>% ## return
dplyr::rename(sistema = `SISTEMA PECUÁRIA`
,safra = SAFRA
,franqueado = FRANQUEADO
,estado = ESTADO
,nome = `NOME DA FAZENDA`
,pastoconf = `PASTO/CONFINAMENTO`
,agricultura = `AGRICULTURA PROPRIA  (COM/SEM)`
) %>%
mutate(franqueado = if_else(franqueado == 1, 'TERRA', franqueado)) %>%
left_join(
readxl::read_xlsx(
'bd/franqueados.xlsx',1
)
)
banco$bd %<>%
filter(!is.na(VALOR)) %>%
#mutate(VALOR = gsub("\\D+",  NA, VALOR)) %>%
mutate(VALOR =  as.numeric(VALOR))
#banco$ref %>%
banco$bd %>%
filter(INDICADOR == "RESULTADO OPERAÇÃO PECUÁRIA (R$/HA) SE") %>%
group_by(safra) %>%
mutate(
top30 = quantile(VALOR,probs = c(.7)),
ref = ifelse(VALOR > top30, 1, 0)
) %>%
dplyr::select(TÉCNICO, franqueado, safra, nome, MUNICÍPIO, top30, ref) -> ref
banco$bd %>%
left_join(
banco$bd %>%
count(franqueado, nome, FRENTE, INDICADOR), by = c('franqueado', 'nome', 'FRENTE', 'INDICADOR')
) %>%
filter(n > 1) %>%
select(nome, safra, franqueado,id_franquia) %>%
distinct() %>%
group_by(nome,franqueado,id_franquia) %>%
summarise(n = n()) %>%
ungroup %>%
mutate(faz_id = 1:nrow(.)) -> banco[['faz_id_data']]  ## return
banco$bd %<>%
left_join(
banco$faz_id_data %>% select(nome, franqueado, faz_id), by = c('nome', 'franqueado')
)
banco$ref <- banco$bd
return(banco)
}
load("~/Área de Trabalho/banco.rData")
b <- carregarDados()
r$hpd
reguaTodas$hpd
debugSource('~/projects/safra-inttegra-19-20/apps/farmview/R/utils.R', echo=TRUE)
r <- processarReguas(b)
x
hpd
hpd
calc_max_hpd <- function(x) {
tryCatch({
#hpd<-suppressWarnings(HPDinterval(obj = as.mcmc(na.omit(x)),prob = .9))
hpd<-HPDinterval(obj = as.mcmc(na.omit(x)),prob = .9)
return(hpd[2])
}, error = function(e){
return(NA)
})
}
r <- processarReguas(b)
r$hpd
reguaTodas$hpd
HPDInterval(obj = as.mcmc(na.omit(x)), prob = .9)
HPDinterval(obj = as.mcmc(na.omit(x)), prob = .9)
library(HPDinterval)
install.packages(HPDinterval)
library(coda)
HPDinterval(obj = as.mcmc(na.omit(x)), prob = .9)
HPDinterval(obj = as.mcmc(na.omit(c(1,2,3,4,5,6,7,8,9))), prob = .9)
r <- processarReguas(b)
r$hpd
r$hpd - reguaTodas$hpd
View(r)
r$q25 - reguaTodas$q25
source('~/projects/safra-inttegra-19-20/apps/farmview/R/preProcessamento.R', echo=TRUE)
r$hpd
r$hpd - reguaTodas$hpd
load("~/Área de Trabalho/banco.rData")
source('~/projects/safra-inttegra-19-20/apps/farmview/R/preProcessamento.R', echo=TRUE)
r$hpd - reguaTodas$hpd
load("~/projects/safra-inttegra-19-20/apps/farmview/bd/reguaTodas.rData")
r$hpd - reguaTodas$hpd
r$max - reguaTodas$max
source('~/projects/safra-inttegra-19-20/apps/farmview/R/preProcessamento.R', echo=TRUE)
source('~/projects/safra-inttegra-19-20/apps/farmview/R/preProcessamento.R', echo=TRUE)
r$max - reguaTodas$max
load("~/Área de Trabalho/banco.rData")
source('~/projects/safra-inttegra-19-20/apps/farmview/R/preProcessamento.R', echo=TRUE)
load("~/projects/safra-inttegra-19-20/apps/farmview/bd/reguaTodas.rData")
r$max - reguaTodas$max
source('~/projects/safra-inttegra-19-20/apps/farmview/R/preProcessamento.R', echo=TRUE)
count(r$max - reguaTodas$max > 0)
count((r$max - reguaTodas$max) > 0)
table(r$max - reguaTodas$max)
z <- r$max - reguaTodas$max
z > 0
z == 0
count(z == 0)
table(z == 0)
source('~/projects/safra-inttegra-19-20/apps/farmview/R/preProcessamento.R', echo=TRUE)
table(z == 0)
count((r$hpd - reguaTodas$hpd) == 0)
table((r$hpd - reguaTodas$hpd) == 0)
source('~/projects/safra-inttegra-19-20/apps/farmview/R/preProcessamento.R', echo=TRUE)
source('~/projects/safra-inttegra-19-20/apps/farmview/R/preProcessamento.R', echo=TRUE)
table((r$hpd - reguaTodas$hpd) == 0)
source('~/projects/safra-inttegra-19-20/apps/farmview/R/preProcessamento.R', echo=TRUE)
dados$hpd
reguaTodas$hpd
data.frame(reguaTodas$hpd, rTodasOld$hpd)
data.frame(hpd_one = reguaTodas$hpd, hpd_two = rTodasOld$hpd) %>% filter(hpd_one != hpd_two)
data.frame(hpd_one = reguaTodas$max, hpd_two = rTodasOld$max) %>% filter(hpd_one != hpd_two)
data.frame(hpd_one = reguaTodas$max, hpd_two = rTodasOld$max) %>% filter(hpd_one != hpd_two) %>% mutate(diferenca = hpd_two - hpd_one)
data.frame(hpd_one = reguaTodas$q25, hpd_two = rTodasOld$q25) %>% filter(hpd_one != hpd_two) %>% mutate(diferenca = hpd_two - hpd_one)
source('~/apps/farmview/R/preProcessamento.R')
source('~/apps/farmview/R/preProcessamento.R')
load("~/apps/farmview/bd/reguaTodas.rData")
reguaTodas$hpd = rTodasOld$hpd
#todos
rTodasOld <- processarReguas(a)
reguaTodas$hpd == rTodasOld$hpd
reguaTodas$max == rTodasOld$max
dataa.frame(one = reguaTodas$max, two = rTodasOld$max) %>% filter(reguaTodas$max != rTodasOld$max)
data.frame(one = reguaTodas$max, two = rTodasOld$max) %>% filter(reguaTodas$max != rTodasOld$max)
data.frame(one = reguaTodas$max, two = rTodasOld$max) %>% filter(reguaTodas$max != rTodasOld$max) %>% mutate(dife = one-two, percentual = dife/one * 100)
data.frame(one = reguaTodas$hpd, two = rTodasOld$hpd) %>% filter(reguaTodas$max != rTodasOld$max) %>% mutate(dife = one-two, percentual = dife/one * 100)
data.frame(one = reguaTodas$hpd, two = rTodasOld$hpd) %>% filter(one != two) %>% mutate(dife = one-two, percentual = dife/one * 100)
data.frame(one = reguaTodas$max, two = rTodasOld$max) %>% filter(one != two) %>% mutate(dife = one-two, percentual = dife/one * 100)
data.frame(one = reguaTodas$q25, two = rTodasOld$q25) %>% filter(one != two) %>% mutate(dife = one-two, percentual = dife/one * 100)
data.frame(one = reguaTodas$q50, two = rTodasOld$q50) %>% filter(one != two) %>% mutate(dife = one-two, percentual = dife/one * 100)
rTodasOld$TOP == reguaTodas$TOP
table(rTodasOld$TOP == reguaTodas$TOP)
table(rTodasOld$INDICADOR == reguaTodas$INDICADOR)
table(rTodasOld$FRENTE == reguaTodas$FRENTE)
table(rTodasOld$hpd == reguaTodas$hpd)
table(rTodasOld$max == reguaTodas$max)
load("~/apps/farmview/bd/reguaCria.rData")
str(reguaCria)
reguaCria
rCOld
head(rCOld)
head(reguaCria)
dplyr::all_equal(reguaCria, rCOld)
dplyr::all_equal(select(reguaCria, min:ref), select(rCOld, min:ref))
dplyr::all_equal(select(reguaCria, data:ref), select(rCOld, data:ref))
select(rCOld, data:ref)
shiny::runApp()
install.packages('broom')
runApp()
load("~/apps/farmview/bd/reguaTodas.rData")
View(reguaTodas)
load("~/apps/farmview/bd/banco.rData")
View(dados)
load("~/apps/farmview/bd/banco.rData")
View(dados)
load("~/apps/farmview/bd/clusters_all.rData")
View(clusters_all)
load("~/apps/farmview/bd/reguaTodas.rData")
View(reguaTodas)
load("~/apps/farmview/bd/reguaTodas.rData")
View(reguaTodas)
load("~/apps/farmview/bd/reguaRecriaEngorda.rData")
View(reguaRecriaEngorda)
# Requisitos
# Arquivos necessários:
# + d/banco.rData
# +
#
library(readxl)
library(dplyr)
library(magrittr)
library(tidyr)
library(coda)
source('R/calc_regua.R')
carregarDados <- function(bd = NULL, indtop = NULL, franqueados = NULL, cluster = NULL){
banco <- list()
banco[['ind_top']] <- read_xlsx(indtop) #'../../../ind_top_old.xlsx')
load(bd) #'bd/banco.rData')
banco[['bd']] <- dados %>% ## return
dplyr::rename(sistema = `SISTEMA PECUÁRIA`
,safra = SAFRA
,franqueado = FRANQUEADO
,estado = ESTADO
,nome = `NOME DA FAZENDA`
,pastoconf = `PASTO/CONFINAMENTO`
,agricultura = `AGRICULTURA PROPRIA  (COM/SEM)`
) %>%
mutate(franqueado = if_else(franqueado == 1, 'TERRA', franqueado)) %>%
left_join(
readxl::read_xlsx(
franqueados, 1 #'../../../franqueados_old.xlsx'
)
)
banco$bd %<>%
filter(!is.na(VALOR)) %>%
#mutate(VALOR = gsub("\\D+",  NA, VALOR)) %>%
mutate(VALOR =  as.numeric(VALOR))
#banco$ref %>%
banco$bd %>%
filter(INDICADOR == "RESULTADO OPERAÇÃO PECUÁRIA (R$/HA) SE") %>%
group_by(safra) %>%
mutate(
top30 = quantile(VALOR,probs = c(.7)),
ref = ifelse(VALOR > top30, 1, 0)
) %>%
dplyr::select(TÉCNICO, franqueado, safra, nome, MUNICÍPIO, top30, ref) -> ref
banco$bd %>%
left_join(
banco$bd %>%
count(franqueado, nome, FRENTE, INDICADOR), by = c('franqueado', 'nome', 'FRENTE', 'INDICADOR')
) %>%
filter(n > 1) %>%
select(nome, safra, franqueado,id_franquia) %>%
distinct() %>%
group_by(nome,franqueado,id_franquia) %>%
summarise(n = n()) %>%
ungroup %>%
mutate(faz_id = 1:nrow(.)) -> banco[['faz_id_data']]  ## return
banco$bd %<>%
left_join(
banco$faz_id_data %>% select(nome, franqueado, faz_id), by = c('nome', 'franqueado')
)
banco$ref <- banco$bd
return(banco)
}
processarReguas <- function(b, sis = NULL){
if(!is.null(sis)){
b$ref <- b$bd %>%  filter(sistema == sis)
}
b$ref %>%
filter(INDICADOR == "RESULTADO OPERAÇÃO PECUÁRIA (R$/HA) SE") %>%
group_by(safra) %>%
mutate(
top30 = quantile(VALOR,probs = c(.7)),
ref = ifelse(VALOR > top30, 1, 0)
) %>%
dplyr::select(TÉCNICO, franqueado, safra, nome, MUNICÍPIO, top30, ref) -> ref
b$ref %>%
filter(FRENTE != "REPRODUÇÃO") %>%
bind_rows(
b$ref %>%
filter(FRENTE == "REPRODUÇÃO") %>%
mutate(safra_ref = ifelse(FRENTE == "REPRODUÇÃO", stringr::str_extract(INDICADOR,"\\d\\d.\\d\\d"), NA)) %>%
mutate(ano1 = ifelse(!is.na(safra_ref), as.numeric(stringr::str_extract(safra_ref, "\\d\\d")), NA)) %>%
mutate(INDICADOR = ifelse(FRENTE == "REPRODUÇÃO", gsub("\\(SF([^\\)]+)\\)","", INDICADOR), INDICADOR)) %>%
mutate(INDICADOR = ifelse(FRENTE == "REPRODUÇÃO", gsub("SF(.*)","", INDICADOR), INDICADOR)) %>%
mutate(INDICADOR = trimws(INDICADOR, which  = "both")) %>%
arrange(INDICADOR, ano1) %>%
group_by(INDICADOR) %>%
filter(ano1 == last(ano1) | is.na(ano1))
) %>%
left_join(
ref, by = c("TÉCNICO", "franqueado", "safra", "nome")
) %>%
right_join(
b$ind_top, by = c("FRENTE", "INDICADOR")
) %>%
group_by(INDICADOR, safra, FRENTE, TOP) %>%
nest %>%
mutate(regua = purrr::pmap(.l=list(data, top = TOP), .f = calc_regua)) %>%
unnest(regua) %>%
arrange(safra, FRENTE, INDICADOR) -> regua
regua %>% ungroup
}
# a <- carregarDados( bd = 'bd/banco.rData', indtop = '../../ind_top_old.xlsx',
#                     franqueados = '../../franqueados_old.xlsx')
#
#
# #todos
# rTodasOld <- processarReguas(a)
# #CRIA
# rCOld <- processarReguas(a, sis = "CRIA")
# #CRIA-RECRIA-ENGORDA
# rCREOld <- processarReguas(a, sis = "CRIA-RECRIA-ENGORDA")
# #RECRIA-ENGORDA
# rREOld <- processarReguas(a, sis = "RECRIA-ENGORDA")
#
# load('bd/reguaTodas.rData')
#
b <- carregarDados( bd = 'bd/banco.rData', indtop = 'bd/ind_top.xlsx',
franqueados = 'bd/franqueados.xlsx')
#todos
reguaTodas <- processarReguas(b)
View(reguaTodas)
# Requisitos
# Arquivos necessários:
# + d/banco.rData
# +
#
library(readxl)
library(dplyr)
library(magrittr)
library(tidyr)
library(coda)
source('R/utils.R')
source('R/calc_regua.R')
carregarDados <- function(bd = NULL, indtop = NULL, franqueados = NULL, cluster = NULL){
banco <- list()
banco[['ind_top']] <- read_xlsx(indtop) #'../../../ind_top_old.xlsx')
load(bd) #'bd/banco.rData')
banco[['bd']] <- dados %>% ## return
dplyr::rename(sistema = `SISTEMA PECUÁRIA`
,safra = SAFRA
,franqueado = FRANQUEADO
,estado = ESTADO
,nome = `NOME DA FAZENDA`
,pastoconf = `PASTO/CONFINAMENTO`
,agricultura = `AGRICULTURA PROPRIA  (COM/SEM)`
) %>%
mutate(franqueado = if_else(franqueado == 1, 'TERRA', franqueado)) %>%
left_join(
readxl::read_xlsx(
franqueados, 1 #'../../../franqueados_old.xlsx'
)
)
banco$bd %<>%
filter(!is.na(VALOR)) %>%
#mutate(VALOR = gsub("\\D+",  NA, VALOR)) %>%
mutate(VALOR =  as.numeric(VALOR))
#banco$ref %>%
banco$bd %>%
filter(INDICADOR == "RESULTADO OPERAÇÃO PECUÁRIA (R$/HA) SE") %>%
group_by(safra) %>%
mutate(
top30 = quantile(VALOR,probs = c(.7)),
ref = ifelse(VALOR > top30, 1, 0)
) %>%
dplyr::select(TÉCNICO, franqueado, safra, nome, MUNICÍPIO, top30, ref) -> ref
banco$bd %>%
left_join(
banco$bd %>%
count(franqueado, nome, FRENTE, INDICADOR), by = c('franqueado', 'nome', 'FRENTE', 'INDICADOR')
) %>%
filter(n > 1) %>%
select(nome, safra, franqueado,id_franquia) %>%
distinct() %>%
group_by(nome,franqueado,id_franquia) %>%
summarise(n = n()) %>%
ungroup %>%
mutate(faz_id = 1:nrow(.)) -> banco[['faz_id_data']]  ## return
banco$bd %<>%
left_join(
banco$faz_id_data %>% select(nome, franqueado, faz_id), by = c('nome', 'franqueado')
)
banco$ref <- banco$bd
return(banco)
}
processarReguas <- function(b, sis = NULL){
if(!is.null(sis)){
b$ref <- b$bd %>%  filter(sistema == sis)
}
b$ref %>%
filter(INDICADOR == "RESULTADO OPERAÇÃO PECUÁRIA (R$/HA) SE") %>%
group_by(safra) %>%
mutate(
top30 = quantile(VALOR,probs = c(.7)),
ref = ifelse(VALOR > top30, 1, 0)
) %>%
dplyr::select(TÉCNICO, franqueado, safra, nome, MUNICÍPIO, top30, ref) -> ref
b$ref %>%
filter(FRENTE != "REPRODUÇÃO") %>%
bind_rows(
b$ref %>%
filter(FRENTE == "REPRODUÇÃO") %>%
mutate(safra_ref = ifelse(FRENTE == "REPRODUÇÃO", stringr::str_extract(INDICADOR,"\\d\\d.\\d\\d"), NA)) %>%
mutate(ano1 = ifelse(!is.na(safra_ref), as.numeric(stringr::str_extract(safra_ref, "\\d\\d")), NA)) %>%
mutate(INDICADOR = ifelse(FRENTE == "REPRODUÇÃO", gsub("\\(SF([^\\)]+)\\)","", INDICADOR), INDICADOR)) %>%
mutate(INDICADOR = ifelse(FRENTE == "REPRODUÇÃO", gsub("SF(.*)","", INDICADOR), INDICADOR)) %>%
mutate(INDICADOR = trimws(INDICADOR, which  = "both")) %>%
arrange(INDICADOR, ano1) %>%
group_by(INDICADOR) %>%
filter(ano1 == last(ano1) | is.na(ano1))
) %>%
left_join(
ref, by = c("TÉCNICO", "franqueado", "safra", "nome")
) %>%
right_join(
b$ind_top, by = c("FRENTE", "INDICADOR")
) %>%
group_by(INDICADOR, safra, FRENTE, TOP) %>%
nest %>%
mutate(regua = purrr::pmap(.l=list(data, top = TOP), .f = calc_regua)) %>%
unnest(regua) %>%
arrange(safra, FRENTE, INDICADOR) -> regua
regua %>% ungroup
}
# a <- carregarDados( bd = 'bd/banco.rData', indtop = '../../ind_top_old.xlsx',
#                     franqueados = '../../franqueados_old.xlsx')
#
#
# #todos
# rTodasOld <- processarReguas(a)
# #CRIA
# rCOld <- processarReguas(a, sis = "CRIA")
# #CRIA-RECRIA-ENGORDA
# rCREOld <- processarReguas(a, sis = "CRIA-RECRIA-ENGORDA")
# #RECRIA-ENGORDA
# rREOld <- processarReguas(a, sis = "RECRIA-ENGORDA")
#
# load('bd/reguaTodas.rData')
#
b <- carregarDados( bd = 'bd/banco.rData', indtop = 'bd/ind_top.xlsx',
franqueados = 'bd/franqueados.xlsx')
#todos
reguaTodas <- processarReguas(b)
View(reguaTodas)
source('~/apps/farmview/R/preProcessamento.R')
View(reguaCompleto)
shiny::runApp()
library(farmView)
shiny::runApp()
shiny::runApp()
installed.packages()
uninstall.packages('farmview')
remove.packages('farmview')
remove.packages('farmView')
shiny::runApp()
shiny::runApp()
shiny::runApp()
installed.packages()
shiny::runApp()
shiny::runApp()
golem::add_shinyserver_file()
golem::add_shinyserver_file()
library(farmView)
library(farmView)
installed.packages()
shiny::runApp()
build()
pkgbuild::build()
runApp()
sessionInfo()
sessionInfo() %>% sort()
library(dplyr)
sessionInfo() %>% sort()
sessionInfo() %>% arrange()
sessionInfo()
runApp()
runApp()
runApp()
runApp()
runApp()
library(farmView)
library(farmView)
shiny::runApp()
shiny::runApp()
library(farmView)
shiny::runApp()
sessionInfo()
devtools::install_github('rstudio/shiny@v1.5.0')
shiny::runApp()
