install.packages("haven")
install.packages("tidyverse")
install.packages("carData")
install.packages("ellipsis")
install.packages("vctrs")
install.packages("Hmisc")
install.packages("writexl")
install.packages("broom")
install.packages("sp")
install.packages("maptools")
install.packages("RANN")
install.packages("rgdal")
install.packages("car")
# install McSpatial:
# install.packages("~/Library/CloudStorage/Dropbox/Clases/TESIS/OLD/16_PedroHerrera/THESIS/Ch_4_Infrastructure/PAPER_Data&Text/McSpatial_2.0.tar.gz", repos = NULL, type = "source")
# https://rdrr.io/cran/McSpatial/

# setwd("C:/Pedro/logit peru")
library(McSpatial)
library(readxl)
library(foreign)
library(ellipsis)
library(carData)
library(rgdal)
library(vctrs)
library(Hmisc)
library(writexl)
library(broom)
library(sp)
library(maptools)
library(RANN)
library(sf)
library(car)
# memory.limit(size = 33000)

#cargando el shape para la matriz de pesos espaciales
cmap<-readOGR (getwd(), layer="Distrito_Total")
#queen
bookw <- makew(cmap)$wmat
n = nrow(bookw)

empresas<-read_excel("DataGMMSplogit v.30.06.2021.xlsx")

#modelos logit
#################

demoslogit1544<- splogit(SIC_1544 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit1544<-cbind(demoslogit1544$coef,demoslogit1544$se,demoslogit1544$u,demoslogit1544$gmat)
logit1544_2<-as.data.frame(logit1544)
write_xlsx(logit1544_2, path = "splog1544.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit1552<- splogit(SIC_1552 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit1552<-cbind(demoslogit1552$coef,demoslogit1552$se,demoslogit1552$u,demoslogit1552$gmat)
logit1552_2<-as.data.frame(logit1552)
write_xlsx(logit1552_2, path = "splog1552.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit1711<- splogit(SIC_1711 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F,  minp=.0001)
logit1711<-cbind(demoslogit1711$coef,demoslogit1711$se,demoslogit1711$u,demoslogit1711$gmat)
logit1711_2<-as.data.frame(logit1711)
write_xlsx(logit1711_2, path = "splog1711.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit1712<- splogit(SIC_1712 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit1712<-cbind(demoslogit1712$coef,demoslogit1712$se,demoslogit1712$u,demoslogit1712$gmat)
logit1712_2<-as.data.frame(logit1712)
write_xlsx(logit1712_2, path = "splog1712.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit1729<- splogit(SIC_1729 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit1729<-cbind(demoslogit1729$coef,demoslogit1729$se,demoslogit1729$u,demoslogit1729$gmat)
logit1729_2<-as.data.frame(logit1729)
write_xlsx(logit1729_2, path = "splog1729.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit1730<- splogit(SIC_1730 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit1730<-cbind(demoslogit1730$coef,demoslogit1730$se,demoslogit1730$u,demoslogit1730$gmat)
logit1730_2<-as.data.frame(logit1730)
write_xlsx(logit1730_2, path = "splog1730.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit1810<- splogit(SIC_1810 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit1810<-cbind(demoslogit1810$coef,demoslogit1810$se,demoslogit1810$u,demoslogit1810$gmat)
logit1810_2<-as.data.frame(logit1810)
write_xlsx(logit1810_2, path = "splog1810.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit1912<- splogit(SIC_1912 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit1912<-cbind(demoslogit1912$coef,demoslogit1912$se,demoslogit1912$u,demoslogit1912$gmat)
logit1912_2<-as.data.frame(logit1912)
write_xlsx(logit1912_2, path = "splog1912.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit2102<- splogit(SIC_2102 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit2102<-cbind(demoslogit2102$coef,demoslogit2102$se,demoslogit2102$u,demoslogit2102$gmat)
logit2102_2<-as.data.frame(logit2102)
write_xlsx(logit2102_2, path = "splog2102.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit2109<- splogit(SIC_2109 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit2109<-cbind(demoslogit2109$coef,demoslogit2109$se,demoslogit2109$u,demoslogit2109$gmat)
logit2109_2<-as.data.frame(logit2109)
write_xlsx(logit2109_2, path = "splog2109.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit2211<- splogit(SIC_2211 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit2211<-cbind(demoslogit2211$coef,demoslogit2211$se,demoslogit2211$u,demoslogit2211$gmat)
logit2211_2<-as.data.frame(logit2211)
write_xlsx(logit2211_2, path = "splog2211.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit2222<- splogit(SIC_2222 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit2222<-cbind(demoslogit2222$coef,demoslogit2222$se,demoslogit2222$u,demoslogit2222$gmat)
logit2222_2<-as.data.frame(logit2222)
write_xlsx(logit2222_2, path = "splog2222.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit2413<- splogit(SIC_2413 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit2413<-cbind(demoslogit2413$coef,demoslogit2413$se,demoslogit2413$u,demoslogit2413$gmat)
logit2413_2<-as.data.frame(logit2413)
write_xlsx(logit2413_2, path = "splog2413.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit2422<- splogit(SIC_2422 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit2422<-cbind(demoslogit2422$coef,demoslogit2422$se,demoslogit2422$u,demoslogit2422$gmat)
logit2422_2<-as.data.frame(logit2422)
write_xlsx(logit2422_2, path = "splog2422.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit2423<- splogit(SIC_2423 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit2423<-cbind(demoslogit2423$coef,demoslogit2423$se,demoslogit2423$u,demoslogit2423$gmat)
logit2423_2<-as.data.frame(logit2423)
write_xlsx(logit2423_2, path = "splog2423.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit2424<- splogit(SIC_2424 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit2424<-cbind(demoslogit2424$coef,demoslogit2424$se,demoslogit2424$u,demoslogit2424$gmat)
logit2424_2<-as.data.frame(logit2424)
write_xlsx(logit2424_2, path = "splog2424.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit2429<- splogit(SIC_2429 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit2429<-cbind(demoslogit2429$coef,demoslogit2429$se,demoslogit2429$u,demoslogit2429$gmat)
logit2429_2<-as.data.frame(logit2429)
write_xlsx(logit2429_2, path = "splog2429.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit2519<- splogit(SIC_2519 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit2519<-cbind(demoslogit2519$coef,demoslogit2519$se,demoslogit2519$u,demoslogit2519$gmat)
logit2519_2<-as.data.frame(logit2519)
write_xlsx(logit2519_2, path = "splog2519.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit2520<- splogit(SIC_2520 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit2520<-cbind(demoslogit2520$coef,demoslogit2520$se,demoslogit2520$u,demoslogit2520$gmat)
logit2520_2<-as.data.frame(logit2520)
write_xlsx(logit2520_2, path = "splog2520.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit2710<- splogit(SIC_2710 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit2710<-cbind(demoslogit2710$coef,demoslogit2710$se)
write.table(logit2710, file = "splog2710.csv", sep = ",", col.names = NA)

demoslogit2732<- splogit(SIC_2732 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit2732<-cbind(demoslogit2732$coef,demoslogit2732$se,demoslogit2732$u,demoslogit2732$gmat)
logit2732_2<-as.data.frame(logit2732)
write_xlsx(logit2732_2, path = "splog2732.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit2891<- splogit(SIC_2891 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit2891<-cbind(demoslogit2891$coef,demoslogit2891$se,demoslogit2891$u,demoslogit2891$gmat)
logit2891_2<-as.data.frame(logit2891)
write_xlsx(logit2891_2, path = "splog2891.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit2893<- splogit(SIC_2893 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit2893<-cbind(demoslogit2893$coef,demoslogit2893$se,demoslogit2893$u,demoslogit2893$gmat)
logit2893_2<-as.data.frame(logit2893)
write_xlsx(logit2893_2, path = "splog2893.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit2899<- splogit(SIC_2899 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit2899<-cbind(demoslogit2899$coef,demoslogit2899$se,demoslogit2899$u,demoslogit2899$gmat)
logit2899_2<-as.data.frame(logit2899)
write_xlsx(logit2899_2, path = "splog2899.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit2915<- splogit(SIC_2915 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit2915<-cbind(demoslogit2915$coef,demoslogit2915$se,demoslogit2915$u,demoslogit2915$gmat)
logit2915_2<-as.data.frame(logit2915)
write_xlsx(logit2915_2, path = "splog2915.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit2919<- splogit(SIC_2919 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit2919<-cbind(demoslogit2919$coef,demoslogit2919$se,demoslogit2919$u,demoslogit2919$gmat)
logit2919_2<-as.data.frame(logit2919)
write_xlsx(logit2919_2, path = "splog2919.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit2924<- splogit(SIC_2924 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit2924<-cbind(demoslogit2924$coef,demoslogit2924$se,demoslogit2924$u,demoslogit2924$gmat)
logit2924_2<-as.data.frame(logit2924)
write_xlsx(logit2924_2, path = "splog2924.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit2929<- splogit(SIC_2929 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit2929<-cbind(demoslogit2929$coef,demoslogit2929$se,demoslogit2929$u,demoslogit2929$gmat)
logit2929_2<-as.data.frame(logit2929)
write_xlsx(logit2929_2, path = "splog2929.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit3190<- splogit(SIC_3190 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit3190<-cbind(demoslogit3190$coef,demoslogit3190$se,demoslogit3190$u,demoslogit3190$gmat)
logit3190_2<-as.data.frame(logit3190)
write_xlsx(logit3190_2, path = "splog3190.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit3311<- splogit(SIC_3311 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit3311<-cbind(demoslogit3311$coef,demoslogit3311$se,demoslogit3311$u,demoslogit3311$gmat)
logit3311_2<-as.data.frame(logit3311)
write_xlsx(logit3311_2, path = "splog3311.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit3693<- splogit(SIC_3693 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit3693<-cbind(demoslogit3693$coef,demoslogit3693$se,demoslogit3693$u,demoslogit3693$gmat)
logit3693_2<-as.data.frame(logit3693)
write_xlsx(logit3693_2, path = "splog3693.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit3694<- splogit(SIC_3694 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit3694<-cbind(demoslogit3694$coef,demoslogit3694$se,demoslogit3694$u,demoslogit3694$gmat)
logit3694_2<-as.data.frame(logit3694)
write_xlsx(logit3694_2, path = "splog3694.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit3699<- splogit(SIC_3699 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit3699<-cbind(demoslogit3699$coef,demoslogit3699$se,demoslogit3699$u,demoslogit3699$gmat)
logit3699_2<-as.data.frame(logit3699)
write_xlsx(logit3699_2, path = "splog3699.xlsx",col_names = TRUE, format_headers = TRUE)

demoslogit3720<- splogit(SIC_3720 ~ Water_supply + Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
logit3720<-cbind(demoslogit3720$coef,demoslogit3720$se,demoslogit3720$u,demoslogit3720$gmat)
logit3720_2<-as.data.frame(logit3720)
write_xlsx(logit3720_2, path = "splog3720.xlsx",col_names = TRUE, format_headers = TRUE)

# List of the coefficients to export, excluding Water_supply
demoslogit1544<- splogit(SIC_1544 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit1552<- splogit(SIC_1552 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit1711<- splogit(SIC_1711 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F,  minp=.0001)
demoslogit1712<- splogit(SIC_1712 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit1729<- splogit(SIC_1729 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit1730<- splogit(SIC_1730 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit1810<- splogit(SIC_1810 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit1912<- splogit(SIC_1912 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit2102<- splogit(SIC_2102 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit2109<- splogit(SIC_2109 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit2211<- splogit(SIC_2211 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit2222<- splogit(SIC_2222 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit2413<- splogit(SIC_2413 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit2422<- splogit(SIC_2422 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit2423<- splogit(SIC_2423 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit2424<- splogit(SIC_2424 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit2429<- splogit(SIC_2429 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit2519<- splogit(SIC_2519 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit2520<- splogit(SIC_2520 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit2710<- splogit(SIC_2710 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit2732<- splogit(SIC_2732 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit2891<- splogit(SIC_2891 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit2893<- splogit(SIC_2893 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit2899<- splogit(SIC_2899 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit2915<- splogit(SIC_2915 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit2919<- splogit(SIC_2919 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit2924<- splogit(SIC_2924 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit2929<- splogit(SIC_2929 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit3190<- splogit(SIC_3190 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit3311<- splogit(SIC_3311 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit3693<- splogit(SIC_3693 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit3694<- splogit(SIC_3694 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit3699<- splogit(SIC_3699 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)
demoslogit3720<- splogit(SIC_3720 ~ Sewerage_service + Electricity_supply + Education + Bank_credit_percapita_dummy + Travel_times , wmat=bookw, data=empresas, silent=F, minp=.0001)