thesis

library(reader)

## Warning: パッケージ 'reader' はバージョン 4.4.2 の R の下で造られました

## 要求されたパッケージ NCmisc をロード中です

## 
## 次のパッケージを付け加えます: 'reader'

## 以下のオブジェクトは 'package:NCmisc' からマスクされています:
## 
##     cat.path, get.ext, rmv.ext

library(ggplot2)
library(car)

## Warning: パッケージ 'car' はバージョン 4.4.2 の R の下で造られました

## 要求されたパッケージ carData をロード中です

## Warning: パッケージ 'carData' はバージョン 4.4.2 の R の下で造られました

library(stargazer)

## 
## Please cite as:

##  Hlavac, Marek (2022). stargazer: Well-Formatted Regression and Summary Statistics Tables.

##  R package version 5.2.3. https://CRAN.R-project.org/package=stargazer

thesis<-read.csv("thesis.csv")

head(thesis)

##   year Fiin  Vote    Pc     Ci   Cin  Cout Conp  Selp  Slin    Ep  Seno
## 1 1887 1.56 58.29 55179 166756 23.86 23.86 84.4 110.9 110.4 42853 10.13
## 2 1991 1.52 45.61 64810 188456 34.89 33.37 92.6 105.7 111.8 59057  9.94
## 3 1995 1.41 41.11 67558 199765 41.67 45.64 95.9 100.5 115.1 67177  9.81
## 4 1999 1.46 43.46 69910 236433 51.89 50.43 98.0  96.6 112.7 70766  9.66
## 5 2003 1.51 44.62 70494 193247 49.58 48.75 95.5  91.6 110.4 79754 10.75
## 6 2007 1.61 43.71 72460 224168 53.31 52.07 95.5  98.1 109.9 79992 12.23
##       Disp Atch     Mon  Tax Soci   Cep    Pop  Sepo Lilo Wolf Shir  Indt   Pin
## 1 24566.76 71.1  365448 26.4 10.6 30.89  98643 10477 70.2  8.1 3.36 48766 11111
## 2 27789.55 73.7  398192 24.8 11.5 33.37 111411 11079 75.6  7.3 7.86 59057 11869
## 3 27995.24 74.6  441561 22.3 13.1 45.64 121307 11890 62.5  4.4 2.75 67177 12509
## 4 28584.27 71.6  599381 20.5 13.6 50.43 128301 12390 59.1  5.1 1.38 70766 12288
## 5 27980.47 72.8 1014294 23.5 14.4 48.75 143008 15368 63.5  3.4 1.38 76788 16856
## 6 29486.89 78.1  886305 22.2 16.6 52.07 157896 19316 69.1  3.7 1.88 79754 13280
##    Pout   Lpri Livl
## 1  8935 481.96 19.5
## 2 10765 963.41 24.9
## 3 12137 442.41 23.4
## 4 12075 315.86 18.2
## 5 11880 284.61 31.4
## 6 12436 348.95 37.8

model1<-lm(Fiin~log(Cout)+log(Cin)+log(Pc)+log(Ci)+log(Indt)+log(Pop)+log(Ep)+year, data=thesis)
model2<-lm(log(Cout)~log(Vote)+Seno+log(Pin)+log(Pout)+log(Pop)+Lilo+year, data=thesis)
model3<-lm(log(Cin)~log(Pc)+log(Ci)+log(Indt)+year, data=thesis)
model4<-lm(log(Pc)~log(Ep)+log(Conp)+log(Disp)+year, data=thesis)
model5<-lm(log(Ci)~log(Cout)+Fiin+log(Ep*Slin)+Shir+log(Indt)+year, data=thesis)
model6<-lm(Vote~Fiin+Lilo+Atch+log(Cout)+log(Pin)+Seno+year, data=thesis)
model7<-lm(log(Disp)~log((Slin)*((Tax+Soci)/100))+year, data=thesis)
model8<-lm(Lilo~Fiin+log(Cout)+Atch+log(Pin)+log(Pout)+year,data=thesis)
model9<-lm(log(Pop)~Fiin+log(Cout)+log(Pin)+log(Pout)+Seno+year,data=thesis)
#model10<-lm(Atch~Livl+Seno+Fiin+log(Lpri)+log(Ep)+log(Pin)+log(Pout)+year,data=thesis)
#model10<-lm(Atch~Livl+Seno+log(Lpri)+log(Pin)+log(Pout)+year,data=thesis)
model10<-lm(Atch~Seno+log(Lpri)+log(Pop)+Slin+year,data=thesis)
model11<-lm(Atch~Livl+year,data=thesis)
model12<-lm(Livl~Seno+year,data=thesis)
model13<-lm(log(Pin)~Fiin+log(Ep)+log(Indt)+log(Disp)+year,data=thesis)
model14<-lm(log(Pout)~Fiin+log(Cout)+Wolf+log(Ep)+log(Indt)+year,data=thesis)


#model13<-lm(log(Pout)~(log(Lpri)*Bubl)+log((Slin)*((Tax+Soci)/100))+year,data=thesis)

#model1<-lm(log(Cin)~log(Pc)+log(Ci)+year, data=thesis)
#model2<-lm(log(Cout)~log(Pc)+log(Ci)+log(Vote)+year, data=thesis)
#model3<-lm(log(Pc)~log(Ep)+log(Conp)+log(Slin)+log(Disp)+year, data=thesis)
#model4<-lm(log(Ci)~log(Cout)+Fiin+log(Ep)+Slin+Shir+log(Indt)+year, data=thesis)
#model5<-lm(Vote~Lilo+Atch+log(Cout)+log(Pin)+Seno+year, data=thesis)
#model6<-lm(log(Disp)~log(Selp)+log(Slin)+log(Conp)+log(Soci)+Tax+year, data=thesis)
#model7<-lm(Lilo~Atch+year,data=thesis)
#model8<-lm(Atch~Seno+year,data=thesis)
#model9<-lm(log(Pop)~Fiin+log(Cout)+year,data=thesis)

stargazer(model1, model2,model3,model4,model5,type="html")

	Dependent variable:
	Fiin	log(Cout)	log(Cin)	log(Pc)	log(Ci)
	(1)	(2)	(3)	(4)	(5)
log(Cout)	-0.391**				2.203***
	(0.010)				(0.352)
log(Cin)	0.099*
	(0.012)
log(Pc)	-1.521*		14.046**
	(0.128)		(4.474)
log(Ci)	0.423**		-1.580**
	(0.017)		(0.528)
Fiin					5.170***
					(0.734)
log(Ep * Slin)					5.698**
					(1.010)
Shir					-0.117**
					(0.033)
log(Indt)	-1.624**		-2.023		-12.345***
	(0.086)		(1.479)		(1.503)
log(Vote)		-1.720**
		(0.247)
Seno		0.061*
		(0.017)
log(Pin)		-0.193**
		(0.031)
log(Pout)		-1.527**
		(0.285)
log(Pop)	1.542**	0.896*
	(0.026)	(0.271)
log(Ep)	0.833*			0.216**
	(0.075)			(0.070)
Lilo		-0.022**
		(0.004)
log(Conp)				0.747**
				(0.270)
log(Disp)				0.217
				(0.272)
year	0.001*	0.002	-0.012*	-0.0001	0.006
	(0.0001)	(0.001)	(0.005)	(0.0005)	(0.004)
Constant	3.070	12.407**	-86.449**	3.275	32.526***
	(0.843)	(2.819)	(22.306)	(1.945)	(4.543)
Observations	10	10	10	10	10
R2	1.000	0.998	0.944	0.987	0.962
Adjusted R2	0.999	0.993	0.900	0.977	0.887
Residual Std. Error	0.002 (df = 1)	0.034 (df = 2)	0.137 (df = 5)	0.013 (df = 5)	0.048 (df = 3)
F Statistic	1,479.669** (df = 8; 1)	172.325*** (df = 7; 2)	21.207*** (df = 4; 5)	96.653*** (df = 4; 5)	12.758** (df = 6; 3)
Note:	p<0.1; p<0.05; p<0.01

stargazer(model6,model7,model8,model9,model10,type="html")

	Dependent variable:
	Vote	log(Disp)	Lilo	log(Pop)	Atch
	(1)	(2)	(3)	(4)	(5)
Fiin	23.239**		0.974	0.932**
	(4.549)		(19.252)	(0.277)
Lilo	-0.283*
	(0.066)
Atch	-0.975**		1.971**
	(0.160)		(0.466)
log(Cout)	-15.359**		-8.355	0.333*
	(2.558)		(6.282)	(0.122)
log(Pin)	-3.945**		-4.459**	0.001
	(0.562)		(1.343)	(0.020)
Seno	2.471**			0.018	1.654***
	(0.355)			(0.010)	(0.237)
log((Slin) * ((Tax + Soci)/100))		0.173*
		(0.077)
log(Pout)			-72.831**	0.514*
			(13.703)	(0.218)
log(Lpri)					4.834***
					(0.705)
log(Pop)					21.612***
					(3.455)
Slin					0.860***
					(0.124)
year	-0.015	0.001***	0.149*	-0.0002	-0.039**
	(0.015)	(0.0002)	(0.060)	(0.001)	(0.009)
Constant	200.870***	6.959***	377.890**	4.517*	-244.563***
	(18.987)	(0.316)	(94.856)	(1.537)	(29.439)
Observations	10	10	10	10	10
R2	0.997	0.945	0.961	0.988	0.994
Adjusted R2	0.988	0.930	0.883	0.963	0.986
Residual Std. Error	0.612 (df = 2)	0.016 (df = 7)	2.636 (df = 3)	0.037 (df = 3)	0.433 (df = 4)
F Statistic	107.551*** (df = 7; 2)	60.335*** (df = 2; 7)	12.286** (df = 6; 3)	40.065*** (df = 6; 3)	125.792*** (df = 5; 4)
Note:	p<0.1; p<0.05; p<0.01

stargazer(model11,model12,model13,model14,type="html")

	Dependent variable:
	Atch	Livl	log(Pin)	log(Pout)
	(1)	(2)	(3)	(4)
Livl	0.255***
	(0.067)
Seno		2.890***
		(0.459)
Fiin			6.340	-1.364**
			(3.948)	(0.338)
log(Cout)				-0.836**
				(0.152)
Wolf				0.012
				(0.011)
log(Ep)			-107.534**	-3.810*
			(30.297)	(1.211)
log(Indt)			115.542**	5.816**
			(31.958)	(1.488)
log(Disp)			-107.487**
			(30.058)
year	0.013	0.056	0.249**	0.005*
	(0.020)	(0.039)	(0.070)	(0.002)
Constant	40.794	-115.641	515.662**	-17.872**
	(37.768)	(75.744)	(149.864)	(5.331)
Observations	10	10	10	10
R2	0.808	0.906	0.784	0.966
Adjusted R2	0.753	0.880	0.514	0.898
Residual Std. Error	1.805 (df = 7)	3.971 (df = 7)	0.494 (df = 4)	0.034 (df = 3)
F Statistic	14.708*** (df = 2; 7)	33.886*** (df = 2; 7)	2.906 (df = 5; 4)	14.246** (df = 6; 3)
Note:	p<0.1; p<0.05; p<0.01

#stargazer(model1,model2,model3, model4,model5,model6,model7, model8,model9,model10,type="html")



#print(model1, type="html")
#print(model2, type="html")
#print(model3, type="html")
#print(model4, type="html")