####目次####
１．分析用データセット作成
２．STATAへのパッケージインストール（1度やればよい）
３．データの下処理（データセットを開くたびに実行）
４．３章のデータ処理と分析・図作成コマンド
５．４章のデータ処理と分析・図作成コマンド
６．５章のデータ処理と分析コマンド
############

####１．分析用データセットの作成####
##本書の分析時にはESS9波は統合化されておらず，ESS8波までのデータと統合したファイルを作る必要があります。本書刊行直前にESS9エディション3(ESSe03_1.dta)アップデートがリリースされましたが，本書分析時点ではエディション2(ESS9e02.dta)用いていました。エディション３ファイルや，将来公開されるだろうESS9も含めた統合ファイルを用いても分析可能と思いますが，完全な再現性を目指す場合はエディション２ファイルを用意し以下の作業を行ってください。(また統合ファイルはサイズが大きいため，８－９波データだけを結合して作ると処理が軽いという効用もあります)

/*ESS8波ファイル(ESS8e02_1)を開きESS9波を加えて統合データ作成*/
use "C:\あなたのフォルダ\ESS8e02_1.dta"
append using "C:\あなたのフォルダ\ESS9e02.dta"


####２．パッケージのインストール####
/*分析結果を表示する際SORを表示するlistcoefコマンドはSTATAに標準装備されていません。listcoefコマンドも実行したい場合にはspost13.adoをインストールしてください*/
search spost13_ado
net install spost13_ado.pkg


####３．データの処理####
/*欠損値処理*/
recode imsmetn (7 8 9 = .)
recode imdfetn (7 8 9 =.)
recode impcntr (7 8 9 = .)
recode imbgeco (77 88 99 = .)
recode imueclt (77 88 99 = .)
recode imwbcnt (77 88 99 = .)
recode edulvla (55= .)
recode edulvlb (100/199 = 1)
recode edulvlb (200/299 = 2)
recode edulvlb (300/399 = 3)
recode edulvlb (400/499 = 4)
recode edulvlb (500/599 = 5)
recode edulvlb (600/699 = 6)
recode edulvlb (700/799 = 7)
recode edulvlb (800 = 8)
recode edulvlb (5555 = .)
recode eisced (55=.)

/*職業リコード*/
recode isco08 (   0/ 999 = 0)
recode isco08 (1000/1999 = 1)
recode isco08 (2000/2999 = 2)
recode isco08 (3000/3999 = 3)
recode isco08 (4000/4999 = 4)
recode isco08 (5000/5999 = 5)
recode isco08 (6000/6999 = 6)
recode isco08 (7000/7999 = 7)
recode isco08 (8000/8999 = 8)
recode isco08 (9000/9999 = 9)

label define isco08 0 "army", modify
label define isco08 1 "manager", modify
label define isco08 2 "specialist", modify
label define isco08 3 "technician", modify
label define isco08 4 "clerk", modify
label define isco08 5 "servicesales", modify
label define isco08 6 "skilledagri", modify
label define isco08 7 "crafter", modify
label define isco08 8 "engineer", modify
label define isco08 9 "elementary", modify

/*NUTSへ統合*/
gen nuts2 = substr(region,1,4)

/*NUTS別のデータを追加する際，略記が競合すことを回避するため設定オフ*/
set varabbrev off

/*地域別失業率変数追加*/
gen reuemp17 =.
replace reuemp17 =7.5 if nuts2=="AT1"
replace reuemp17 =5.2 if nuts2=="AT11"
replace reuemp17 =4.8 if nuts2=="AT12"
replace reuemp17 =10.4 if nuts2=="AT13"
replace reuemp17 =4.6 if nuts2=="AT2"
replace reuemp17 =4.8 if nuts2=="AT21"
replace reuemp17 =4.5 if nuts2=="AT22"
replace reuemp17 =3.6 if nuts2=="AT3"
replace reuemp17 =4.0 if nuts2=="AT31"
replace reuemp17 =3.1 if nuts2=="AT32"
replace reuemp17 =3.3 if nuts2=="AT33"
replace reuemp17 =3.7 if nuts2=="AT34"
replace reuemp17 =14.9 if nuts2=="BE1"
replace reuemp17 =14.9 if nuts2=="BE10"
replace reuemp17 =4.4 if nuts2=="BE2"
replace reuemp17 =5.9 if nuts2=="BE21"
replace reuemp17 =4.1 if nuts2=="BE22"
replace reuemp17 =3.4 if nuts2=="BE23"
replace reuemp17 =4.7 if nuts2=="BE24"
replace reuemp17 =3.2 if nuts2=="BE25"
replace reuemp17 =9.7 if nuts2=="BE3"
replace reuemp17 =7.6 if nuts2=="BE31"
replace reuemp17 =11.4 if nuts2=="BE32"
replace reuemp17 =10.5 if nuts2=="BE33"
replace reuemp17 =6.7 if nuts2=="BE34"
replace reuemp17 =7.4 if nuts2=="BE35"
replace reuemp17 =8.5 if nuts2=="BG3"
replace reuemp17 =11.3 if nuts2=="BG31"
replace reuemp17 =6.9 if nuts2=="BG32"
replace reuemp17 =9.4 if nuts2=="BG33"
replace reuemp17 =7.0 if nuts2=="BG34"
replace reuemp17 =4.1 if nuts2=="BG4"
replace reuemp17 =3.3 if nuts2=="BG41"
replace reuemp17 =5.3 if nuts2=="BG42"
replace reuemp17 =11.1 if nuts2=="CY0"
replace reuemp17 =11.1 if nuts2=="CY00"
replace reuemp17 =2.9 if nuts2=="CZ0"
replace reuemp17 =1.7 if nuts2=="CZ01"
replace reuemp17 =2.1 if nuts2=="CZ02"
replace reuemp17 =2.1 if nuts2=="CZ03"
replace reuemp17 =3.4 if nuts2=="CZ04"
replace reuemp17 =2.8 if nuts2=="CZ05"
replace reuemp17 =3.1 if nuts2=="CZ06"
replace reuemp17 =3.3 if nuts2=="CZ07"
replace reuemp17 =4.7 if nuts2=="CZ08"
replace reuemp17 =2.9 if nuts2=="DE1"
replace reuemp17 =3.0 if nuts2=="DE11"
replace reuemp17 =3.3 if nuts2=="DE12"
replace reuemp17 =2.9 if nuts2=="DE13"
replace reuemp17 =2.2 if nuts2=="DE14"
replace reuemp17 =2.3 if nuts2=="DE2"
replace reuemp17 =2.2 if nuts2=="DE21"
replace reuemp17 =2.1 if nuts2=="DE22"
replace reuemp17 =2.2 if nuts2=="DE23"
replace reuemp17 =3.0 if nuts2=="DE24"
replace reuemp17 =2.1 if nuts2=="DE25"
replace reuemp17 =2.3 if nuts2=="DE26"
replace reuemp17 =2.6 if nuts2=="DE27"
replace reuemp17 =7.0 if nuts2=="DE3"
replace reuemp17 =7.0 if nuts2=="DE30"
replace reuemp17 =4.5 if nuts2=="DE4"
replace reuemp17 =4.5 if nuts2=="DE40"
replace reuemp17 =4.3 if nuts2=="DE5"
replace reuemp17 =4.3 if nuts2=="DE50"
replace reuemp17 =4.2 if nuts2=="DE6"
replace reuemp17 =4.2 if nuts2=="DE60"
replace reuemp17 =3.3 if nuts2=="DE7"
replace reuemp17 =3.5 if nuts2=="DE71"
replace reuemp17 =3.6 if nuts2=="DE72"
replace reuemp17 =2.5 if nuts2=="DE73"
replace reuemp17 =5.2 if nuts2=="DE8"
replace reuemp17 =5.2 if nuts2=="DE80"
replace reuemp17 =3.8 if nuts2=="DE9"
replace reuemp17 =4.4 if nuts2=="DE91"
replace reuemp17 =4.5 if nuts2=="DE92"
replace reuemp17 =3.4 if nuts2=="DE93"
replace reuemp17 =3.0 if nuts2=="DE94"
replace reuemp17 =4.1 if nuts2=="DEA"
replace reuemp17 =4.3 if nuts2=="DEA1"
replace reuemp17 =3.7 if nuts2=="DEA2"
replace reuemp17 =3.9 if nuts2=="DEA3"
replace reuemp17 =3.9 if nuts2=="DEA4"
replace reuemp17 =4.7 if nuts2=="DEA5"
replace reuemp17 =3.3 if nuts2=="DEB"
replace reuemp17 =3.1 if nuts2=="DEB1"
replace reuemp17 =2.0 if nuts2=="DEB2"
replace reuemp17 =3.8 if nuts2=="DEB3"
replace reuemp17 =4.5 if nuts2=="DEC"
replace reuemp17 =4.5 if nuts2=="DEC1"
replace reuemp17 =4.4 if nuts2=="DED"
replace reuemp17 =4.3 if nuts2=="DED1"
replace reuemp17 =4.1 if nuts2=="DED2"
replace reuemp17 =5.1 if nuts2=="DED3"
replace reuemp17 =6.9 if nuts2=="DEE"
replace reuemp17 =6.9 if nuts2=="DEE0"
replace reuemp17 =3.6 if nuts2=="DEF"
replace reuemp17 =3.6 if nuts2=="DEF0"
replace reuemp17 =4.4 if nuts2=="DEG"
replace reuemp17 =4.4 if nuts2=="DEG0"
replace reuemp17 =5.8 if nuts2=="DK0"
replace reuemp17 =6.2 if nuts2=="DK01"
replace reuemp17 =5.6 if nuts2=="DK02"
replace reuemp17 =6.1 if nuts2=="DK03"
replace reuemp17 =5.4 if nuts2=="DK04"
replace reuemp17 =5.6 if nuts2=="DK05"
replace reuemp17 =5.8 if nuts2=="EE0"
replace reuemp17 =5.8 if nuts2=="EE00"
replace reuemp17 =23.0 if nuts2=="EL3"
replace reuemp17 =19.5 if nuts2=="EL30"
replace reuemp17 =22.9 if nuts2=="EL4"
replace reuemp17 =29.1 if nuts2=="EL41"
replace reuemp17 =24.8 if nuts2=="EL42"
replace reuemp17 =21.2 if nuts2=="EL43"
replace reuemp17 =20.7 if nuts2=="EL5"
replace reuemp17 =19.7 if nuts2=="EL51"
replace reuemp17 =26.3 if nuts2=="EL52"
replace reuemp17 =20.9 if nuts2=="EL53"
replace reuemp17 =16.8 if nuts2=="EL54"
replace reuemp17 =21.6 if nuts2=="EL6"
replace reuemp17 =21.6 if nuts2=="EL61"
replace reuemp17 =18.0 if nuts2=="EL62"
replace reuemp17 =22.5 if nuts2=="EL63"
replace reuemp17 =16.0 if nuts2=="EL64"
replace reuemp17 =17.7 if nuts2=="EL65"
replace reuemp17 =14.9 if nuts2=="ES1"
replace reuemp17 =15.7 if nuts2=="ES11"
replace reuemp17 =13.7 if nuts2=="ES12"
replace reuemp17 =13.6 if nuts2=="ES13"
replace reuemp17 =11.3 if nuts2=="ES2"
replace reuemp17 =11.3 if nuts2=="ES21"
replace reuemp17 =10.2 if nuts2=="ES22"
replace reuemp17 =12.0 if nuts2=="ES23"
replace reuemp17 =11.7 if nuts2=="ES24"
replace reuemp17 =13.4 if nuts2=="ES3"
replace reuemp17 =13.4 if nuts2=="ES30"
replace reuemp17 =18.9 if nuts2=="ES4"
replace reuemp17 =14.1 if nuts2=="ES41"
replace reuemp17 =20.8 if nuts2=="ES42"
replace reuemp17 =26.3 if nuts2=="ES43"
replace reuemp17 =15.0 if nuts2=="ES5"
replace reuemp17 =13.4 if nuts2=="ES51"
replace reuemp17 =18.2 if nuts2=="ES52"
replace reuemp17 =12.4 if nuts2=="ES53"
replace reuemp17 =24.4 if nuts2=="ES6"
replace reuemp17 =25.5 if nuts2=="ES61"
replace reuemp17 =18.0 if nuts2=="ES62"
replace reuemp17 =22.4 if nuts2=="ES63"
replace reuemp17 =27.6 if nuts2=="ES64"
replace reuemp17 =23.5 if nuts2=="ES7"
replace reuemp17 =23.5 if nuts2=="ES70"
replace reuemp17 =8.7 if nuts2=="FI1"
replace reuemp17 =9.3 if nuts2=="FI19"
replace reuemp17 =7.7 if nuts2=="FI1B"
replace reuemp17 =8.5 if nuts2=="FI1C"
replace reuemp17 =9.6 if nuts2=="FI1D"
replace reuemp17 =: if nuts2=="FI2"
replace reuemp17 =: if nuts2=="FI20"
replace reuemp17 =8.7 if nuts2=="FR1"
replace reuemp17 =8.7 if nuts2=="FR10"
replace reuemp17 =8.6 if nuts2=="FRB"
replace reuemp17 =8.6 if nuts2=="FRB0"
replace reuemp17 = if nuts2=="FRC"
replace reuemp17 =10.6 if nuts2=="FRC1"
replace reuemp17 = if nuts2=="FRC2"
replace reuemp17 =10.0 if nuts2=="FRD"
replace reuemp17 =8.3 if nuts2=="FRD1"
replace reuemp17 =11.3 if nuts2=="FRD2"
replace reuemp17 =12.0 if nuts2=="FRE"
replace reuemp17 =13.0 if nuts2=="FRE1"
replace reuemp17 =10.1 if nuts2=="FRE2"
replace reuemp17 =10.0 if nuts2=="FRF"
replace reuemp17 =9.0 if nuts2=="FRF1"
replace reuemp17 =9.5 if nuts2=="FRF2"
replace reuemp17 =11.0 if nuts2=="FRF3"
replace reuemp17 =7.1 if nuts2=="FRG"
replace reuemp17 =7.1 if nuts2=="FRG0"
replace reuemp17 =7.3 if nuts2=="FRH"
replace reuemp17 =7.3 if nuts2=="FRH0"
replace reuemp17 =9.3 if nuts2=="FRI"
replace reuemp17 =10.2 if nuts2=="FRI1"
replace reuemp17 =6.2 if nuts2=="FRI2"
replace reuemp17 =9.0 if nuts2=="FRI3"
replace reuemp17 =9.5 if nuts2=="FRJ"
replace reuemp17 =12.2 if nuts2=="FRJ1"
replace reuemp17 =7.2 if nuts2=="FRJ2"
replace reuemp17 =7.4 if nuts2=="FRK"
replace reuemp17 =8.3 if nuts2=="FRK1"
replace reuemp17 =7.2 if nuts2=="FRK2"
replace reuemp17 =10.3 if nuts2=="FRL"
replace reuemp17 =10.3 if nuts2=="FRL0"
replace reuemp17 =7.8 if nuts2=="FRM"
replace reuemp17 =7.8 if nuts2=="FRM0"
replace reuemp17 =22.3 if nuts2=="FRY"
replace reuemp17 =22.4 if nuts2=="FRY1"
replace reuemp17 =17.8 if nuts2=="FRY2"
replace reuemp17 =22.4 if nuts2=="FRY3"
replace reuemp17 =22.8 if nuts2=="FRY4"
replace reuemp17 =29.6 if nuts2=="FRY5"
replace reuemp17 =11.2 if nuts2=="HR0"
replace reuemp17 =10.9 if nuts2=="HR03"
replace reuemp17 =11.4 if nuts2=="HR04"
replace reuemp17 =2.7 if nuts2=="HU1"
replace reuemp17 =2.8 if nuts2=="HU10"
replace reuemp17 =2.9 if nuts2=="HU11"
replace reuemp17 =2.6 if nuts2=="HU12"
replace reuemp17 =3.5 if nuts2=="HU2"
replace reuemp17 =2.2 if nuts2=="HU21"
replace reuemp17 =2.4 if nuts2=="HU22"
replace reuemp17 =6.3 if nuts2=="HU23"
replace reuemp17 =5.9 if nuts2=="HU3"
replace reuemp17 =5.8 if nuts2=="HU31"
replace reuemp17 =7.4 if nuts2=="HU32"
replace reuemp17 =4.1 if nuts2=="HU33"
replace reuemp17 =6.7 if nuts2=="IE0"
replace reuemp17 =6.7 if nuts2=="IE01"
replace reuemp17 =6.7 if nuts2=="IE02"
replace reuemp17 =6.2 if nuts2=="IE04"
replace reuemp17 =7.0 if nuts2=="IE05"
replace reuemp17 =6.7 if nuts2=="IE06"
replace reuemp17 =7.4 if nuts2=="ITC"
replace reuemp17 =9.1 if nuts2=="ITC1"
replace reuemp17 =7.8 if nuts2=="ITC2"
replace reuemp17 =9.5 if nuts2=="ITC3"
replace reuemp17 =6.4 if nuts2=="ITC4"
replace reuemp17 =6.3 if nuts2=="ITH"
replace reuemp17 =3.1 if nuts2=="ITH1"
replace reuemp17 =5.7 if nuts2=="ITH2"
replace reuemp17 =6.3 if nuts2=="ITH3"
replace reuemp17 =6.7 if nuts2=="ITH4"
replace reuemp17 =6.6 if nuts2=="ITH5"
replace reuemp17 =10.0 if nuts2=="ITI"
replace reuemp17 =8.6 if nuts2=="ITI1"
replace reuemp17 =10.6 if nuts2=="ITI2"
replace reuemp17 =10.6 if nuts2=="ITI3"
replace reuemp17 =10.7 if nuts2=="ITI4"
replace reuemp17 =19.0 if nuts2=="ITF"
replace reuemp17 =11.7 if nuts2=="ITF1"
replace reuemp17 =14.6 if nuts2=="ITF2"
replace reuemp17 =20.9 if nuts2=="ITF3"
replace reuemp17 =18.9 if nuts2=="ITF4"
replace reuemp17 =12.8 if nuts2=="ITF5"
replace reuemp17 =21.6 if nuts2=="ITF6"
replace reuemp17 =20.3 if nuts2=="ITG"
replace reuemp17 =21.5 if nuts2=="ITG1"
replace reuemp17 =17.0 if nuts2=="ITG2"
replace reuemp17 =7.1 if nuts2=="LT0"
replace reuemp17 =7.1 if nuts2=="LT00"
replace reuemp17 =4.8 if nuts2=="LT01"
replace reuemp17 =8.1 if nuts2=="LT02"
replace reuemp17 =5.5 if nuts2=="LU0"
replace reuemp17 =5.5 if nuts2=="LU00"
replace reuemp17 =8.7 if nuts2=="LV0"
replace reuemp17 =8.7 if nuts2=="LV00"
replace reuemp17 =4.0 if nuts2=="MT0"
replace reuemp17 =4.0 if nuts2=="MT00"
replace reuemp17 =6.0 if nuts2=="NL1"
replace reuemp17 =7.2 if nuts2=="NL11"
replace reuemp17 =5.6 if nuts2=="NL12"
replace reuemp17 =5.0 if nuts2=="NL13"
replace reuemp17 =4.8 if nuts2=="NL2"
replace reuemp17 =5.1 if nuts2=="NL21"
replace reuemp17 =4.5 if nuts2=="NL22"
replace reuemp17 =5.7 if nuts2=="NL23"
replace reuemp17 =4.9 if nuts2=="NL3"
replace reuemp17 =4.2 if nuts2=="NL31"
replace reuemp17 =4.5 if nuts2=="NL32"
replace reuemp17 =5.6 if nuts2=="NL33"
replace reuemp17 =2.9 if nuts2=="NL34"
replace reuemp17 =4.3 if nuts2=="NL4"
replace reuemp17 =4.2 if nuts2=="NL41"
replace reuemp17 =4.8 if nuts2=="NL42"
replace reuemp17 =5.4 if nuts2=="PL1"
replace reuemp17 =4.6 if nuts2=="PL11"
replace reuemp17 =6.4 if nuts2=="PL12"
replace reuemp17 =4.0 if nuts2=="PL2"
replace reuemp17 =4.2 if nuts2=="PL21"
replace reuemp17 =3.9 if nuts2=="PL22"
replace reuemp17 =7.1 if nuts2=="PL3"
replace reuemp17 =7.2 if nuts2=="PL31"
replace reuemp17 =8.4 if nuts2=="PL32"
replace reuemp17 =7 if nuts2=="PL33"
replace reuemp17 =4.7 if nuts2=="PL34"
replace reuemp17 =3.6 if nuts2=="PL4"
replace reuemp17 =3.1 if nuts2=="PL41"
replace reuemp17 =4.7 if nuts2=="PL42"
replace reuemp17 =3.7 if nuts2=="PL43"
replace reuemp17 =4.6 if nuts2=="PL5"
replace reuemp17 =4.7 if nuts2=="PL51"
replace reuemp17 =4.3 if nuts2=="PL52"
replace reuemp17 =5.4 if nuts2=="PL6"
replace reuemp17 =5.5 if nuts2=="PL61"
replace reuemp17 =7.2 if nuts2=="PL62"
replace reuemp17 =4.2 if nuts2=="PL63"
replace reuemp17 =8.9 if nuts2=="PT1"
replace reuemp17 =9.8 if nuts2=="PT11"
replace reuemp17 =7.8 if nuts2=="PT15"
replace reuemp17 =7.1 if nuts2=="PT16"
replace reuemp17 =9.5 if nuts2=="PT17"
replace reuemp17 =8.5 if nuts2=="PT18"
replace reuemp17 =9.0 if nuts2=="PT2"
replace reuemp17 =9.0 if nuts2=="PT20"
replace reuemp17 =10.6 if nuts2=="PT3"
replace reuemp17 =10.6 if nuts2=="PT30"
replace reuemp17 =3.9 if nuts2=="RO1"
replace reuemp17 =3.6 if nuts2=="RO11"
replace reuemp17 =4.2 if nuts2=="RO12"
replace reuemp17 =4.3 if nuts2=="RO2"
replace reuemp17 =2.9 if nuts2=="RO21"
replace reuemp17 =6.5 if nuts2=="RO22"
replace reuemp17 =5.7 if nuts2=="RO3"
replace reuemp17 =6.8 if nuts2=="RO31"
replace reuemp17 =4.3 if nuts2=="RO32"
replace reuemp17 =6.2 if nuts2=="RO4"
replace reuemp17 =7.7 if nuts2=="RO41"
replace reuemp17 =4.4 if nuts2=="RO42"
replace reuemp17 =6.8 if nuts2=="SE1"
replace reuemp17 =6.3 if nuts2=="SE11"
replace reuemp17 =7.6 if nuts2=="SE12"
replace reuemp17 =6.7 if nuts2=="SE2"
replace reuemp17 =5.5 if nuts2=="SE21"
replace reuemp17 =8.4 if nuts2=="SE22"
replace reuemp17 =6.0 if nuts2=="SE23"
replace reuemp17 =6.5 if nuts2=="SE3"
replace reuemp17 =6.9 if nuts2=="SE31"
replace reuemp17 =6.2 if nuts2=="SE32"
replace reuemp17 =5.9 if nuts2=="SE33"
replace reuemp17 =6.6 if nuts2=="SI0"
replace reuemp17 =6.8 if nuts2=="SI01"
replace reuemp17 =6.4 if nuts2=="SI02"
replace reuemp17 =6.8 if nuts2=="SI03"
replace reuemp17 =6.4 if nuts2=="SI04"
replace reuemp17 =8.1 if nuts2=="SK0"
replace reuemp17 =4.2 if nuts2=="SK01"
replace reuemp17 =5.5 if nuts2=="SK02"
replace reuemp17 =9.4 if nuts2=="SK03"
replace reuemp17 =12.0 if nuts2=="SK04"
replace reuemp17 =4.8 if nuts2=="CH0"
replace reuemp17 =7.9 if nuts2=="CH01"
replace reuemp17 =4.2 if nuts2=="CH02"
replace reuemp17 =4.4 if nuts2=="CH03"
replace reuemp17 =4.5 if nuts2=="CH04"
replace reuemp17 =3.8 if nuts2=="CH05"
replace reuemp17 =2.9 if nuts2=="CH06"
replace reuemp17 =6.1 if nuts2=="CH07"
replace reuemp17 =2.7 if nuts2=="IS0"
replace reuemp17 =2.7 if nuts2=="IS00"
replace reuemp17 =. if nuts2=="LI0"
replace reuemp17 =. if nuts2=="LI00"
replace reuemp17 =4.2 if nuts2=="NO0"
replace reuemp17 =4.6 if nuts2=="NO01"
replace reuemp17 =2.8 if nuts2=="NO02"
replace reuemp17 =4.4 if nuts2=="NO03"
replace reuemp17 =5.2 if nuts2=="NO04"
replace reuemp17 =4.0 if nuts2=="NO05"
replace reuemp17 =2.9 if nuts2=="NO06"
replace reuemp17 =3.2 if nuts2=="NO07"
replace reuemp17 =5.6 if nuts2=="UKC"
replace reuemp17 =4.1 if nuts2=="UKD"
replace reuemp17 =4.8 if nuts2=="UKE"
replace reuemp17 =4 if nuts2=="UKF"
replace reuemp17 =5.5 if nuts2=="UKG"
replace reuemp17 =3.9 if nuts2=="UKH"
replace reuemp17 =5.3 if nuts2=="UKI"
replace reuemp17 =3.2 if nuts2=="UKJ"
replace reuemp17 =3.6 if nuts2=="UKK"
replace reuemp17 =4.5 if nuts2=="UKL"
replace reuemp17 =4.1 if nuts2=="UKM"
replace reuemp17 =4.6 if nuts2=="UKN"

/*人口密度変数追加*/
gen density16 =.
replace density16 = if nuts2=="AT1"
replace density16 =79.4 if nuts2=="AT11"
replace density16 =87.7 if nuts2=="AT12"
replace density16 =4693.4 if nuts2=="AT13"
replace density16 = if nuts2=="AT2"
replace density16 =59.9 if nuts2=="AT21"
replace density16 =76 if nuts2=="AT22"
replace density16 = if nuts2=="AT3"
replace density16 =124.6 if nuts2=="AT31"
replace density16 =77.7 if nuts2=="AT32"
replace density16 =59.3 if nuts2=="AT33"
replace density16 =152.5 if nuts2=="AT34"
replace density16 =7415.3 if nuts2=="BE1"
replace density16 =7454.6 if nuts2=="BE10"
replace density16 =483.5 if nuts2=="BE2"
replace density16 =655.9 if nuts2=="BE21"
replace density16 =363.6 if nuts2=="BE22"
replace density16 =506.7 if nuts2=="BE23"
replace density16 =536 if nuts2=="BE24"
replace density16 =379.5 if nuts2=="BE25"
replace density16 =214.6 if nuts2=="BE3"
replace density16 =365.5 if nuts2=="BE31"
replace density16 =355.4 if nuts2=="BE32"
replace density16 =287.4 if nuts2=="BE33"
replace density16 =64 if nuts2=="BE34"
replace density16 =134.9 if nuts2=="BE35"
replace density16 = if nuts2=="BG3"
replace density16 =41.2 if nuts2=="BG31"
replace density16 =55.3 if nuts2=="BG32"
replace density16 =64.7 if nuts2=="BG33"
replace density16 =53.6 if nuts2=="BG34"
replace density16 = if nuts2=="BG4"
replace density16 =104.7 if nuts2=="BG41"
replace density16 =64.7 if nuts2=="BG42"
replace density16 =92.4 if nuts2=="CY0"
replace density16 =92.4 if nuts2=="CY00"
replace density16 =136.7 if nuts2=="CZ0"
replace density16 =2626.8 if nuts2=="CZ01"
replace density16 =124.4 if nuts2=="CZ02"
replace density16 =70.9 if nuts2=="CZ03"
replace density16 =132.1 if nuts2=="CZ04"
replace density16 =123.1 if nuts2=="CZ05"
replace density16 =123 if nuts2=="CZ06"
replace density16 =133.6 if nuts2=="CZ07"
replace density16 =228 if nuts2=="CZ08"
replace density16 =307.7 if nuts2=="DE1"
replace density16 =387.1 if nuts2=="DE11"
replace density16 =403.9 if nuts2=="DE12"
replace density16 =241 if nuts2=="DE13"
replace density16 =210.6 if nuts2=="DE14"
replace density16 =183.7 if nuts2=="DE2"
replace density16 =268.3 if nuts2=="DE21"
replace density16 =118.5 if nuts2=="DE22"
replace density16 =113.6 if nuts2=="DE23"
replace density16 =147 if nuts2=="DE24"
replace density16 =241.8 if nuts2=="DE25"
replace density16 =153.9 if nuts2=="DE26"
replace density16 =187.3 if nuts2=="DE27"
replace density16 =4160.8 if nuts2=="DE3"
replace density16 =4193.2 if nuts2=="DE30"
replace density16 =85.6 if nuts2=="DE4"
replace density16 =85.8 if nuts2=="DE40"
replace density16 =1717.4 if nuts2=="DE5"
replace density16 =1726.7 if nuts2=="DE50"
replace density16 =2517.5 if nuts2=="DE6"
replace density16 =2533.7 if nuts2=="DE60"
replace density16 =293.6 if nuts2=="DE7"
replace density16 =532.4 if nuts2=="DE71"
replace density16 =193.8 if nuts2=="DE72"
replace density16 =147.1 if nuts2=="DE73"
replace density16 =71.42 if nuts2=="DE8"
replace density16 =71.4 if nuts2=="DE80"
replace density16 =167.4 if nuts2=="DE9"
replace density16 =197.5 if nuts2=="DE91"
replace density16 =237.6 if nuts2=="DE92"
replace density16 =110.6 if nuts2=="DE93"
replace density16 =167.9 if nuts2=="DE94"
replace density16 =527.9 if nuts2=="DEA"
replace density16 =1002.4 if nuts2=="DEA1"
replace density16 =607.6 if nuts2=="DEA2"
replace density16 =379.1 if nuts2=="DEA3"
replace density16 =316.4 if nuts2=="DEA4"
replace density16 =450.4 if nuts2=="DEA5"
replace density16 =205.6 if nuts2=="DEB"
replace density16 =185.8 if nuts2=="DEB1"
replace density16 =108.3 if nuts2=="DEB2"
replace density16 =300 if nuts2=="DEB3"
replace density16 =388.2 if nuts2=="DEC"
replace density16 =388.4 if nuts2=="DEC1"
replace density16 =224.7 if nuts2=="DED"
replace density16 =206.1 if nuts2=="DED1"
replace density16 =224.3 if nuts2=="DED2"
replace density16 =262.2 if nuts2=="DED3"
replace density16 =110.9 if nuts2=="DEE"
replace density16 =110.7 if nuts2=="DEE0"
replace density16 =185.5 if nuts2=="DEF"
replace density16 =186.3 if nuts2=="DEF0"
replace density16 =134.5 if nuts2=="DEG"
replace density16 =134.1 if nuts2=="DEG0"
replace density16 =135.9 if nuts2=="DK0"
replace density16 =738.5 if nuts2=="DK01"
replace density16 =117.3 if nuts2=="DK02"
replace density16 =101.3 if nuts2=="DK03"
replace density16 =101.7 if nuts2=="DK04"
replace density16 =76 if nuts2=="DK05"
replace density16 =30.3 if nuts2=="EE0"
replace density16 =30.3 if nuts2=="EE00"
replace density16 = if nuts2=="EL3"
replace density16 =992.5 if nuts2=="EL30"
replace density16 = if nuts2=="EL4"
replace density16 =52.4 if nuts2=="EL41"
replace density16 =63.7 if nuts2=="EL42"
replace density16 =75.9 if nuts2=="EL43"
replace density16 = if nuts2=="EL5"
replace density16 =42.9 if nuts2=="EL51"
replace density16 =99.6 if nuts2=="EL52"
replace density16 =29.5 if nuts2=="EL53"
replace density16 =37 if nuts2=="EL54"
replace density16 = if nuts2=="EL6"
replace density16 =52 if nuts2=="EL61"
replace density16 =89.5 if nuts2=="EL62"
replace density16 =60.2 if nuts2=="EL63"
replace density16 =36 if nuts2=="EL64"
replace density16 =37.5 if nuts2=="EL65"
replace density16 = if nuts2=="ES1"
replace density16 =92.4 if nuts2=="ES11"
replace density16 =98.2 if nuts2=="ES12"
replace density16 =110.5 if nuts2=="ES13"
replace density16 = if nuts2=="ES2"
replace density16 =301.2 if nuts2=="ES21"
replace density16 =61.7 if nuts2=="ES22"
replace density16 =62.2 if nuts2=="ES23"
replace density16 =27.8 if nuts2=="ES24"
replace density16 = if nuts2=="ES3"
replace density16 =809.2 if nuts2=="ES30"
replace density16 = if nuts2=="ES4"
replace density16 =26.1 if nuts2=="ES41"
replace density16 =25.8 if nuts2=="ES42"
replace density16 =26.4 if nuts2=="ES43"
replace density16 = if nuts2=="ES5"
replace density16 =232.3 if nuts2=="ES51"
replace density16 =213.3 if nuts2=="ES52"
replace density16 =229.9 if nuts2=="ES53"
replace density16 = if nuts2=="ES6"
replace density16 =96.9 if nuts2=="ES61"
replace density16 =130 if nuts2=="ES62"
replace density16 =4242.2 if nuts2=="ES63"
replace density16 =6061.1 if nuts2=="ES64"
replace density16 = if nuts2=="ES7"
replace density16 =289.1 if nuts2=="ES70"
replace density16 = if nuts2=="FI1"
replace density16 =23.7 if nuts2=="FI19"
replace density16 =179.1 if nuts2=="FI1B"
replace density16 =36.9 if nuts2=="FI1C"
replace density16 =6.4 if nuts2=="FI1D"
replace density16 = if nuts2=="FI2"
replace density16 =18.7 if nuts2=="FI20"
replace density16 =1013.8 if nuts2=="FR1"
replace density16 =1015.3 if nuts2=="FR10"
replace density16 =65.7 if nuts2=="FRB"
replace density16 =65.7 if nuts2=="FRB0"
replace density16 =59 if nuts2=="FRC"
replace density16 =51.8 if nuts2=="FRC1"
replace density16 =72.8 if nuts2=="FRC2"
replace density16 =111.24 if nuts2=="FRD"
replace density16 =83.2 if nuts2=="FRD1"
replace density16 =151.8 if nuts2=="FRD2"
replace density16 =188.9 if nuts2=="FRE"
replace density16 =328.1 if nuts2=="FRE1"
replace density16 =99.6 if nuts2=="FRE2"
replace density16 =96.9 if nuts2=="FRF"
replace density16 =228.1 if nuts2=="FRF1"
replace density16 =52.1 if nuts2=="FRF2"
replace density16 =99.3 if nuts2=="FRF3"
replace density16 =116.1 if nuts2=="FRG"
replace density16 =116.4 if nuts2=="FRG0"
replace density16 =120.7 if nuts2=="FRH"
replace density16 =120.9 if nuts2=="FRH0"
replace density16 =70.5 if nuts2=="FRI"
replace density16 =82.3 if nuts2=="FRI1"
replace density16 =43.3 if nuts2=="FRI2"
replace density16 =69.8 if nuts2=="FRI3"
replace density16 =79.8 if nuts2=="FRJ"
replace density16 =102.4 if nuts2=="FRJ1"
replace density16 =66.6 if nuts2=="FRJ2"
replace density16 =112.6 if nuts2=="FRK"
replace density16 =52.3 if nuts2=="FRK1"
replace density16 =148.6 if nuts2=="FRK2"
replace density16 =160.8 if nuts2=="FRL"
replace density16 =161 if nuts2=="FRL0"
replace density16 =38.1 if nuts2=="FRM"
replace density16 =38.3 if nuts2=="FRM0"
replace density16 = if nuts2=="FRY"
replace density16 =254.9 if nuts2=="FRY1"
replace density16 =339.9 if nuts2=="FRY2"
replace density16 =3.3 if nuts2=="FRY3"
replace density16 =345.2 if nuts2=="FRY4"
replace density16 =673.4 if nuts2=="FRY5"
replace density16 = if nuts2=="HR0"
replace density16 =56.7 if nuts2=="HR03"
replace density16 =88.6 if nuts2=="HR04"
replace density16 = if nuts2=="HU1"
replace density16 =433.4 if nuts2=="HU10"
replace density16 =3443.2 if nuts2=="HU11"
replace density16 =197.6 if nuts2=="HU12"
replace density16 = if nuts2=="HU2"
replace density16 =99 if nuts2=="HU21"
replace density16 =87.9 if nuts2=="HU22"
replace density16 =65.2 if nuts2=="HU23"
replace density16 = if nuts2=="HU3"
replace density16 =86.4 if nuts2=="HU31"
replace density16 =84.5 if nuts2=="HU32"
replace density16 =69.4 if nuts2=="HU33"
replace density16 = if nuts2=="IE0"
replace density16 =39.2 if nuts2=="IE01"
replace density16 =96 if nuts2=="IE02"
replace density16 =33.8 if nuts2=="IE04"
replace density16 =53.9 if nuts2=="IE05"
replace density16 =162.7 if nuts2=="IE06"
replace density16 =278 if nuts2=="ITC"
replace density16 =176.9 if nuts2=="ITC1"
replace density16 =39.2 if nuts2=="ITC2"
replace density16 =294.2 if nuts2=="ITC3"
replace density16 =439.2 if nuts2=="ITC4"
replace density16 =187 if nuts2=="ITH"
replace density16 =70.9 if nuts2=="ITH1"
replace density16 =87.7 if nuts2=="ITH2"
replace density16 =279.5 if nuts2=="ITH3"
replace density16 =161.4 if nuts2=="ITH4"
replace density16 =198.1 if nuts2=="ITH5"
replace density16 =208 if nuts2=="ITI"
replace density16 =165.2 if nuts2=="ITI1"
replace density16 =107.9 if nuts2=="ITI2"
replace density16 =163.9 if nuts2=="ITI3"
replace density16 =348.8 if nuts2=="ITI4"
replace density16 =191 if nuts2=="ITF"
replace density16 =124.5 if nuts2=="ITF1"
replace density16 =71.2 if nuts2=="ITF2"
replace density16 =436.5 if nuts2=="ITF3"
replace density16 =212.1 if nuts2=="ITF4"
replace density16 =58.9 if nuts2=="ITF5"
replace density16 =133.4 if nuts2=="ITF6"
replace density16 =135 if nuts2=="ITG"
replace density16 =199.4 if nuts2=="ITG1"
replace density16 =69.3 if nuts2=="ITG2"
replace density16 =45.8 if nuts2=="LT0"
replace density16 =45.8 if nuts2=="LT00"
replace density16 =85.5 if nuts2=="LT01"
replace density16 =38.8 if nuts2=="LT02"
replace density16 = if nuts2=="LU0"
replace density16 =225.1 if nuts2=="LU00"
replace density16 = if nuts2=="LV0"
replace density16 =31 if nuts2=="LV00"
replace density16 = if nuts2=="MT0"
replace density16 =1450.2 if nuts2=="MT00"
replace density16 = if nuts2=="NL1"
replace density16 =247.3 if nuts2=="NL11"
replace density16 =190.3 if nuts2=="NL12"
replace density16 =183.5 if nuts2=="NL13"
replace density16 = if nuts2=="NL2"
replace density16 =340.9 if nuts2=="NL21"
replace density16 =408.3 if nuts2=="NL22"
replace density16 =282.9 if nuts2=="NL23"
replace density16 = if nuts2=="NL3"
replace density16 =910.4 if nuts2=="NL31"
replace density16 =1016 if nuts2=="NL32"
replace density16 =1260.8 if nuts2=="NL33"
replace density16 =211.5 if nuts2=="NL34"
replace density16 = if nuts2=="NL4"
replace density16 =505.1 if nuts2=="NL41"
replace density16 =517.8 if nuts2=="NL42"
replace density16 = if nuts2=="PL1"
replace density16 =136.7 if nuts2=="PL11"
replace density16 =151.3 if nuts2=="PL12"
replace density16 = if nuts2=="PL2"
replace density16 =221.9 if nuts2=="PL21"
replace density16 =371.6 if nuts2=="PL22"
replace density16 = if nuts2=="PL3"
replace density16 =84.9 if nuts2=="PL31"
replace density16 =117.7 if nuts2=="PL32"
replace density16 =106.5 if nuts2=="PL33"
replace density16 =57.9 if nuts2=="PL34"
replace density16 = if nuts2=="PL4"
replace density16 =117.1 if nuts2=="PL41"
replace density16 =77 if nuts2=="PL42"
replace density16 =72.8 if nuts2=="PL43"
replace density16 = if nuts2=="PL5"
replace density16 =144.9 if nuts2=="PL51"
replace density16 =102.2 if nuts2=="PL52"
replace density16 = if nuts2=="PL6"
replace density16 =117.3 if nuts2=="PL61"
replace density16 =61.5 if nuts2=="PL62"
replace density16 =129 if nuts2=="PL63"
replace density16 = if nuts2=="PT1"
replace density16 =170.2 if nuts2=="PT11"
replace density16 =89.5 if nuts2=="PT15"
replace density16 =80.4 if nuts2=="PT16"
replace density16 =1002.5 if nuts2=="PT17"
replace density16 =23.2 if nuts2=="PT18"
replace density16 = if nuts2=="PT2"
replace density16 =106.5 if nuts2=="PT20"
replace density16 = if nuts2=="PT3"
replace density16 =320 if nuts2=="PT30"
replace density16 = if nuts2=="RO1"
replace density16 =76.4 if nuts2=="RO11"
replace density16 =69.2 if nuts2=="RO12"
replace density16 = if nuts2=="RO2"
replace density16 =89.9 if nuts2=="RO21"
replace density16 =78.9 if nuts2=="RO22"
replace density16 = if nuts2=="RO3"
replace density16 =90.2 if nuts2=="RO31"
replace density16 =1300.6 if nuts2=="RO32"
replace density16 = if nuts2=="RO4"
replace density16 =69.7 if nuts2=="RO41"
replace density16 =56.9 if nuts2=="RO42"
replace density16 = if nuts2=="SE1"
replace density16 =344.9 if nuts2=="SE11"
replace density16 =43 if nuts2=="SE12"
replace density16 = if nuts2=="SE2"
replace density16 =25.4 if nuts2=="SE21"
replace density16 =105.9 if nuts2=="SE22"
replace density16 =67.7 if nuts2=="SE23"
replace density16 = if nuts2=="SE3"
replace density16 =13.3 if nuts2=="SE31"
replace density16 =5.3 if nuts2=="SE32"
replace density16 =3.4 if nuts2=="SE33"
replace density16 = if nuts2=="SI0"
replace density16 =88.5 if nuts2=="SI01"
replace density16 =124.6 if nuts2=="SI02"
replace density16 =88.5 if nuts2=="SI03"
replace density16 =124.6 if nuts2=="SI04"
replace density16 = if nuts2=="SK0"
replace density16 =315.5 if nuts2=="SK01"
replace density16 =123.4 if nuts2=="SK02"
replace density16 =83 if nuts2=="SK03"
replace density16 =103.4 if nuts2=="SK04"
replace density16 = if nuts2=="CH0"
replace density16 =194.2 if nuts2=="CH01"
replace density16 =189.8 if nuts2=="CH02"
replace density16 =587.1 if nuts2=="CH03"
replace density16 =893.1 if nuts2=="CH04"
replace density16 =102.7 if nuts2=="CH05"
replace density16 =186.9 if nuts2=="CH06"
replace density16 =129.2 if nuts2=="CH07"
replace density16 = if nuts2=="IS0"
replace density16 =3.3 if nuts2=="IS00"
replace density16 = if nuts2=="LI0"
replace density16 =238.7 if nuts2=="LI00"
replace density16 = if nuts2=="NO0"
replace density16 =249 if nuts2=="NO01"
replace density16 =7.6 if nuts2=="NO02"
replace density16 =28.9 if nuts2=="NO03"
replace density16 =31.7 if nuts2=="NO04"
replace density16 =18.9 if nuts2=="NO05"
replace density16 =11.5 if nuts2=="NO06"
replace density16 =4.4 if nuts2=="NO07"
replace density16 =306.8 if nuts2=="UKC"
replace density16 =510.4 if nuts2=="UKD"
replace density16 =351 if nuts2=="UKE"
replace density16 =300.9 if nuts2=="UKF"
replace density16 =444.9 if nuts2=="UKG"
replace density16 =319.2 if nuts2=="UKH"
replace density16 =5549.5 if nuts2=="UKI"
replace density16 =471.3 if nuts2=="UKJ"
replace density16 =230.5 if nuts2=="UKK"
replace density16 =149.8 if nuts2=="UKL"
replace density16 =69.2 if nuts2=="UKM"
replace density16 =134.6 if nuts2=="UKN"

/*人口密度の千分率とログデータ用意*/
gen density16x = density16/1000
gen density16l = log(density16)

/*いったんオフにした変数名略記設定をオンに*/
set varabbrev on

/* 右翼支持（広義）変数の作成*/
gen frcl = .
replace frcl = 1 if cntry =="AT" & prtclcat == 3
replace frcl = 1 if cntry =="AT" & prtcldat == 3
replace frcl = 1 if cntry =="BE" & prtclcbe == 3
replace frcl = 1 if cntry =="BE" & prtclcbe == 7
replace frcl = 1 if cntry =="BE" & prtclcbe == 15
replace frcl = 1 if cntry =="BE" & prtcldbe == 3
replace frcl = 1 if cntry =="BE" & prtcldbe == 7
replace frcl = 1 if cntry =="BE" & prtcldbe == 15
replace frcl = 1 if cntry =="BG" & prtcldbg == 3
replace frcl = 1 if cntry =="BG" & prtcldbg == 8
replace frcl = 1 if cntry =="BG" & prtcldbg == 9
replace frcl = 1 if cntry =="BG" & prtcldbg == 11
replace frcl = 1 if cntry =="CH" & prtclfch == 1
replace frcl = 1 if cntry =="CH" & prtclfch == 10
replace frcl = 1 if cntry =="CH" & prtclfch == 12
replace frcl = 1 if cntry =="CH" & prtclgch == 1
replace frcl = 1 if cntry =="CH" & prtclgch == 10
replace frcl = 1 if cntry =="CH" & prtclgch == 12
replace frcl = 1 if cntry =="CZ" & prtcldcz == 5
replace frcl = 1 if cntry =="CZ" & prtcldcz == 7
replace frcl = 1 if cntry =="CZ" & prtclecz == 5
replace frcl = 1 if cntry =="CZ" & prtclecz == 8
replace frcl = 1 if cntry =="CY" & prtclbcy == 4
replace frcl = 1 if cntry =="CY" & prtclbcy == 5
replace frcl = 1 if cntry =="DE" & prtclede == 6
replace frcl = 1 if cntry =="DE" & prtclede == 6
replace frcl = 1 if cntry =="EE" & prtclfee == 6
replace frcl = 1 if cntry =="EE" & prtclgee == 6
replace frcl = 1 if cntry =="ES" & prtclfes == 16
replace frcl = 1 if cntry =="FI" & prtcldfi == 4
replace frcl = 1 if cntry =="FI" & prtclefi == 4
replace frcl = 1 if cntry =="FI" & prtclefi == 5
replace frcl = 1 if cntry =="FR" & prtclefr == 2
replace frcl = 1 if cntry =="FR" & prtclffr == 11
replace frcl = 1 if cntry =="HU" & prtclfhu == 3
replace frcl = 1 if cntry =="HU" & prtclfhu == 4
replace frcl = 1 if cntry =="HU" & prtclghu == 3
replace frcl = 1 if cntry =="HU" & prtclghu == 4
replace frcl = 1 if cntry =="HR" & prtclahr == 10
replace frcl = 1 if cntry =="HR" & prtclahr == 13
replace frcl = 1 if cntry =="IS" & prtclbis == 4
replace frcl = 1 if cntry =="IT" & prtclcit == 3
replace frcl = 1 if cntry =="IT" & prtclcit == 5	
replace frcl = 1 if cntry =="IT" & prtcldit == 9
replace frcl = 1 if cntry =="IT" & prtcldit == 10
replace frcl = 1 if cntry =="LV" & prtclalv == 4
replace frcl = 1 if cntry =="LT" & prtclblt == 5
replace frcl = 1 if cntry =="LT" & prtclblt == 6
replace frcl = 1 if cntry =="NL" & prtclenl == 3
replace frcl = 1 if cntry =="NL" & prtclfnl == 3
replace frcl = 1 if cntry =="NL" & prtclfnl == 13
replace frcl = 1 if cntry =="NO" & prtclbno == 8
replace frcl = 1 if cntry =="NO" & prtclbno == 8
replace frcl = 1 if cntry =="PL" & prtclgpl == 1
replace frcl = 1 if cntry =="PL" & prtclgpl == 2
replace frcl = 1 if cntry =="PL" & prtclgpl == 6
replace frcl = 1 if cntry =="PL" & prtclhpl == 1
replace frcl = 1 if cntry =="PL" & prtclhpl == 5
replace frcl = 1 if cntry =="PL" & prtclhpl == 8
replace frcl = 1 if cntry =="SE" & prtclbse ==9
replace frcl = 1 if cntry =="SE" & prtclcse ==9
replace frcl = 1 if cntry =="SI" & prtclesi ==6
replace frcl = 1 if cntry =="SI" & prtclfsi ==11
replace frcl = 1 if cntry =="SI" & prtclfsi ==8
replace frcl = 1 if cntry =="SK" & prtcldsk ==6
replace frcl = 1 if cntry =="SK" & prtcldsk ==8
replace frcl = 1 if cntry =="GB" & prtclbgb ==1
replace frcl = 1 if cntry =="GB" & prtclbgb ==7
replace frcl = 1 if cntry =="GB" & prtclbgb ==10
replace frcl = 1 if cntry =="GB" & prtclcgb ==1
replace frcl = 1 if cntry =="GB" & prtclcgb ==7
replace frcl = 1 if cntry =="GB" & prtclcgb ==10

/*支持政党有りの回答者のみに限定した範囲での右翼政党支持変数を作成*/
gen frcl1 = frcl
replace frcl1 = 0 if clsprty==1 & frcl!=1

/*全回答者内での右翼政党支持変数を作成*/
gen frcl2 = frcl
replace frcl2 = 0 if frcl!=1

/* 極右政党＝右翼支持（狭義）変数の作成*/
gen rrcl = .
replace rrcl = 1 if cntry =="AT" & prtclcat == 3
replace rrcl = 1 if cntry =="AT" & prtcldat == 3
replace rrcl = 1 if cntry =="BE" & prtclcbe == 7
replace rrcl = 1 if cntry =="BE" & prtcldbe == 7
replace rrcl = 1 if cntry =="BG" & prtcldbg == 3
replace rrcl = 1 if cntry =="BG" & prtcldbg == 8
replace rrcl = 1 if cntry =="BG" & prtcldbg == 9
replace rrcl = 1 if cntry =="CZ" & prtcldcz == 7
replace rrcl = 1 if cntry =="CZ" & prtclecz == 8
replace rrcl = 1 if cntry =="CY" & prtclbcy == 4
replace rrcl = 1 if cntry =="DE" & prtclede == 6
replace rrcl = 1 if cntry =="DE" & prtclede == 6
replace rrcl = 1 if cntry =="EE" & prtclfee == 6
replace rrcl = 1 if cntry =="EE" & prtclgee == 6
replace rrcl = 1 if cntry =="ES" & prtclfes == 16
replace rrcl = 1 if cntry =="FI" & prtcldfi == 4
replace rrcl = 1 if cntry =="FI" & prtclefi == 4
replace rrcl = 1 if cntry =="FR" & prtclefr == 2
replace rrcl = 1 if cntry =="FR" & prtclffr == 11
replace rrcl = 1 if cntry =="HR" & prtclahr == 10
replace rrcl = 1 if cntry =="HU" & prtclfhu == 4
replace rrcl = 1 if cntry =="HU" & prtclghu == 4
replace rrcl = 1 if cntry =="IT" & prtclcit == 3
replace frcl = 1 if cntry =="IT" & prtclcit == 5
replace rrcl = 1 if cntry =="IT" & prtcldit == 9
replace rrcl = 1 if cntry =="IT" & prtcldit == 10
replace frcl = 1 if cntry =="LV" & prtclalv == 4
replace frcl = 1 if cntry =="NL" & prtclenl == 3
replace rrcl = 1 if cntry =="NL" & prtclfnl == 3
replace frcl = 1 if cntry =="PL" & prtclgpl == 1
replace frcl = 1 if cntry =="PL" & prtclgpl == 2
replace rrcl = 1 if cntry =="PL" & prtclhpl == 8
replace rrcl = 1 if cntry =="SI" & prtclfsi ==11
replace rrcl = 1 if cntry =="SK" & prtcldsk ==8

/*支持政党有りの回答者のみに限定した範囲での極右政党支持変数を作成*/
gen rrcl1 = rrcl
replace rrcl1 = 0 if clsprty==1 & rrcl!=1

/*全回答者内での極右政党支持変数を作成*/
gen rrcl2 = rrcl
replace rrcl2 = 0 if rrcl!=1

/* 国コードの作成追加 */
gen ccode = 0
replace ccode = 1 if cntry=="AT"
replace ccode = 2 if cntry=="BE"
replace ccode = 3 if cntry=="BG"
replace ccode = 4 if cntry=="CH"
replace ccode = 5 if cntry=="CY"
replace ccode = 6 if cntry=="CZ"
replace ccode = 7 if cntry=="DE"
replace ccode = 8 if cntry=="DK"
replace ccode = 9 if cntry=="EE"
replace ccode = 10 if cntry=="ES"
replace ccode = 11 if cntry=="FI"
replace ccode = 12 if cntry=="FR"
replace ccode = 13 if cntry=="GB"
replace ccode = 14 if cntry=="GR"
replace ccode = 15 if cntry=="HR"
replace ccode = 16 if cntry=="HU"
replace ccode = 17 if cntry=="IE"
replace ccode = 18 if cntry=="IL"
replace ccode = 19 if cntry=="IS"
replace ccode = 20 if cntry=="IT"
replace ccode = 21 if cntry=="LT"
replace ccode = 22 if cntry=="LU"
replace ccode = 23 if cntry=="NL"
replace ccode = 24 if cntry=="NO"
replace ccode = 25 if cntry=="PL"
replace ccode = 26 if cntry=="PT"
replace ccode = 27 if cntry=="RU"
replace ccode = 28 if cntry=="SE"
replace ccode = 29 if cntry=="SI"
replace ccode = 30 if cntry=="SK"
replace ccode = 31 if cntry=="TR"
replace ccode = 32 if cntry=="UA"
replace ccode = 33 if cntry=="RS"

/* 国イヤー変数の作成 */
generate cy = cntry + string(essround)

/* 分析除外変数の作成（非民主国とイスラエル）*/
gen nondemo = 0
replace nondemo = 1 if cntry=="UA"
replace nondemo = 1 if cntry=="RU"
replace nondemo = 1 if cntry=="IL"
replace nondemo = 1 if cntry=="RS"
replace nondemo = 1 if cntry=="ME"

/* 東欧ダミー変数の作成 */
gen east = 0
replace east = 1 if cntry=="EE"
replace east = 1 if cntry=="LV"
replace east = 1 if cntry=="LT"
replace east = 1 if cntry=="PL"
replace east = 1 if cntry=="CZ"
replace east = 1 if cntry=="SK"
replace east = 1 if cntry=="HU"
replace east = 1 if cntry=="SI"
replace east = 1 if cntry=="HR"
replace east = 1 if cntry=="BG"

/* ウェイトの作成 */
gen newweight = dweight*pweight

/* 解釈簡単化のために変数反転 */
gen imbgeco1 = 10-imbgeco
gen imueclt1 = 10-imueclt
gen imptrad1 = 6-imptrad
gen ipfrule1 = 6-ipfrule
gen impfree1 = 6-impfree
gen gincdif1 = 6-gincdif
gen trstprl1 = 10-trstprl
gen trstep1 = 10-trstep
gen euftf1 = 10-euftf
gen ipeqopt1 = 10-ipeqopt
gen stflife1 = 10-stflife
gen stfdem1 = 10-stfdem
gen stfeco1 = 10-stfeco


####第３章の分析一覧####

/*メインの右翼政党支持分析（表2.1 p72）*/
/*モデル１*/
logistic c.frcl2 c.agea c.gndr c.hinctnta c.edulvlb i.mnactic c.density16x c. reuemp17 i.ccode [pweight = newweight]if nondemo==0, vce(cluster ccode)
listcoef
/*モデル２*/
logistic c.frcl2 c.agea c.gndr c.hinctnta c.edulvlb i.mnactic c.density16x c. reuemp17 stflife1 stfeco1 stfdem1 c.euftf1 c.imptrad1 c.ipfrule1 c.atchctr c.impfree1 c.imbgeco1 c.imueclt1 i.ccode [pweight = newweight]if nondemo==0, vce(cluster ccode)
listcoef
/*モデル３*/
logistic c.frcl1 c.agea c.gndr c.hinctnta c.edulvlb i.mnactic c.density16x c. reuemp17 stflife1 stfeco1 stfdem1 c.euftf1 c.imptrad1 c.ipfrule1 c.atchctr c.impfree1 c.imbgeco1 c.imueclt1 i.ccode [pweight = newweight]if nondemo==0, vce(cluster ccode)
listcoef

/*補遺の極右政党分析（補遺表2.1 p.259）*/
/*モデル１*/
logistic c.rrcl2 c.agea c.gndr c.hinctnta c.edulvlb i.mnactic c.density16x c. reuemp17 i.ccode [pweight = newweight]if nondemo==0, vce(cluster ccode)
listcoef
/*モデル２*/
logistic c.rrcl2 c.agea c.gndr c.hinctnta c.edulvlb i.mnactic c.density16x c. reuemp17 stflife1 stfeco1 stfdem1 c.euftf1 c.imptrad1 c.ipfrule1 c.atchctr c.impfree1 c.imbgeco1 c.imueclt1 i.ccode [pweight = newweight]if nondemo==0, vce(cluster ccode)
listcoef
/*モデル３*/
logistic c.rrcl1 c.agea c.gndr c.hinctnta c.edulvlb i.mnactic c.density16x c. reuemp17 stflife1 stfeco1 stfdem1 c.euftf1 c.imptrad1 c.ipfrule1 c.atchctr c.impfree1 c.imbgeco1 c.imueclt1 i.ccode [pweight = newweight]if nondemo==0, vce(cluster ccode)
listcoef

/*補遺の右翼政党分析（職種別分析）（補遺表2.2 p.260）*/
/*モデル１*/
logistic c.frcl2 c.agea c.gndr c.hinctnta c.edulvlb ib1.isco08 c.density16x c. reuemp17 i.ccode [pweight = newweight]if nondemo==0, vce(cluster ccode)
listcoef
/*モデル２*/
logistic c.frcl2 c.agea c.gndr c.hinctnta c.edulvlb ib1.isco08 c.density16x c. reuemp17 stflife1 stfeco1 stfdem1 c.euftf1 c.imptrad1 c.ipfrule1 c.atchctr c.impfree1 c.imbgeco1 c.imueclt1 i.ccode [pweight = newweight]if nondemo==0, vce(cluster ccode)
listcoef

/*図の作成*/
/*(図2.2)*/
graph bar (mean) frcl2 [pweight = newweight] if nondemo==0, over(stflife, relabel(1"大いに不満" 2" " 3" " 4" " 5" " 6"どちらでもない" 7" " 8" " 9" " 10" " 11"大いに満足")) ytitle("") ylabel(0(0.05)0.2) title(人生への満足度と右翼政党支持確率) scheme(s1mono)

graph bar (mean) frcl2 [pweight = newweight] if nondemo==0, over(hinctnta, relabel(1"最貧困層" 2" " 3" " 4" " 5" " 6" " 7" " 8" " 9" " 10"最富裕層")) ytitle("") ylabel(0(0.05)0.2) title(家計所得と右翼政党支持確率) scheme(s1mono)

/*(図2.3)*/
graph bar (mean) frcl2 [pweight = newweight] if nondemo==0, over(imptrad1, relabel(1"伝統を重視しない" 2" " 3" " 4" " 5" " 6"伝統を重視する")) ytitle("") ylabel(0(0.05)0.2) title(伝統や習慣への追随志向と右翼政党支持確率) scheme(s1mono)

graph bar (mean) frcl2 [pweight = newweight] if nondemo==0, over(ipfrule1, relabel(1"規則を重視しない" 2" " 3" " 4" " 5" " 6"規則を重視する")) ytitle("") ylabel(0(0.05)0.2) title(命令や規則への服従志向と右翼政党支持確率) scheme(s1mono)

graph bar (mean) frcl2 [pweight = newweight] if nondemo==0, over(imueclt1, relabel(1"自国文化は豊かになる" 2" " 3" " 4" " 5" " 6"どちらでもない" 7" " 8" " 9" " 10" " 11"自国文化が破壊される")) ytitle("") ylabel(0(0.05)0.2) title(移民による自国文化への影響の認識と右翼政党支持確率) scheme(s1mono)

graph bar (mean) frcl2 [pweight = newweight] if nondemo==0, over(imbgeco1, relabel(1"自国経済に良い" 2" " 3" " 4" " 5" " 6"どちらでもない" 7" " 8" " 9" " 10" " 11"自国経済に悪い")) ytitle("") ylabel(0(0.05)0.2) title(移民による自国経済への影響の認識と右翼政党支持確率) scheme(s1mono)

/*(図2.4)*/
graph bar (mean) frcl2 [pweight = newweight] if nondemo==0, over(stfdem, relabel(1"大いに不満" 2" " 3" " 4" " 5" " 6"どちらでもない" 7" " 8" " 9" " 10" " 11"大いに満足")) ytitle("") ylabel(0(0.05)0.2) title(自国民主政治への満足度と右翼政党支持確率) scheme(s1mono)

graph bar (mean) frcl2 [pweight = newweight] if nondemo==0, over(euftf, relabel(1"行き過ぎである" 2" " 3" " 4" " 5" " 6"どちらでもない" 7" " 8" " 9" " 10" " 11"より推進すべき")) ytitle("") ylabel(0(0.05)0.2) title(欧州統合政治への満足度と右翼政党支持確率) scheme(s1mono)

graph bar (mean) frcl2 [pweight = newweight] if nondemo==0, over(trstprl, relabel(1"信頼していない" 2" " 3" " 4" " 5" " 6"どちらでもない" 7" " 8" " 9" " 10" " 11"信頼している")) ytitle("") ylabel(0(0.05)0.2) title(自国議会への信頼度と右翼政党支持確率) scheme(s1mono)

graph bar (mean) frcl2 [pweight = newweight] if nondemo==0, over(trstep, relabel(1"信頼していない" 2" " 3" " 4" " 5" " 6"どちらでもない" 7" " 8" " 9" " 10" " 11"信頼している")) ytitle("") ylabel(0(0.05)0.2) title(欧州議会への信頼度と右翼政党支持確率) scheme(s1mono)


####第４章の分析一覧####

/*文化侵蝕懸念ダミー変数の作成*/
gen imueclt2 = imueclt1
recode imueclt2 (0/5 = 0)
recode imueclt2 (6/10 = 1)

/*文化侵蝕懸念の回帰分析(表3.1 p.98)*/
/*モデル１*/
logistic imueclt2 c.agea c.gndr c.hinctnta c.edulvlb i.mnactic c.density16x c.reuemp17 i.ccode [pweight = newweight]if nondemo==0, vce(cluster ccode)
listcoef
/*モデル２*/
logistic imueclt2 c.agea c.gndr c.hinctnta c.edulvlb i.mnactic c.density16x c.reuemp17 stflife1 stfeco1 stfdem1 c.euftf1 c.imptrad1 c.ipfrule1 c.atchctr c.impfree1 i.ccode [pweight = newweight]if nondemo==0, vce(cluster ccode)
listcoef

/*補遺の職種別：文化侵蝕懸念の回帰分析(補遺表３ p.261)*/
/*モデル１*/
logistic c. imueclt2 c.agea c.gndr c.hinctnta c.edulvlb ib1.isco08 c.density16x c. reuemp17 i.ccode [pweight = newweight]if nondemo==0, vce(cluster ccode)
listcoef
/*モデル２*/
logistic c. imueclt2 c.agea c.gndr c.hinctnta c.edulvlb ib1.isco08 c.density16x c. reuemp17 stflife1 stfeco1 stfdem1 c.euftf1 c.imptrad1 c.ipfrule1 c.atchctr c.impfree1 i.ccode [pweight = newweight]if nondemo==0, vce(cluster ccode)
listcoef

/*図の作成*/
/*(図3.2)*/
graph bar (mean) imueclt2 [pweight = newweight] if nondemo==0, over(hinctnta, relabel(1"最貧困層" 2" " 3" " 4" " 5" " 6"　" 7" " 8" " 9" " 10"最富裕層")) ytitle("") ylabel(0(0.1)0.5) title(家計所得と文化的反移民感情保持確率) scheme(s1mono)

graph bar (mean) imueclt2 [pweight = newweight] if nondemo==0, over(stflife, relabel(1"大いに不満" 2" " 3" " 4" " 5" " 6"どちらでもない" 7" " 8" " 9" " 10" " 11"大いに満足")) ytitle("") ylabel(0(0.1)0.5) title(人生への満足度と文化的反移民感情保持確率) scheme(s1mono)

/*(図3.3)*/
graph bar (mean) imueclt2 [pweight = newweight] if nondemo==0, over(stfdem, relabel(1"大いに不満" 2" " 3" " 4" " 5" " 6"どちらでもない" 7" " 8" " 9" " 10" " 11"大いに満足")) ytitle("") ylabel(0(0.1)0.5) title(自国民主政治への満足度と文化的反移民感情保持確率) scheme(s1mono)

graph bar (mean) imueclt2 [pweight = newweight] if nondemo==0, over(euftf, relabel(1"行き過ぎである" 2" " 3" " 4" " 5" " 6"どちらでもない" 7" " 8" " 9" " 10" " 11"より推進すべき")) ytitle("") ylabel(0(0.1)0.5) title(欧州統合政治への満足度と文化的反移民感情保持確率) scheme(s1mono)

/*(本文にない割愛した図)*/
graph bar (mean) imueclt2 [pweight = newweight] if nondemo==0, over(trstprl, relabel(1"信頼していない" 2" " 3" " 4" " 5" " 6"どちらでもない" 7" " 8" " 9" " 10" " 11"信頼している")) ytitle("") ylabel(0(0.1)0.5) title(自国議会への信頼度と文化的反移民感情保持確率) scheme(s1mono)

graph bar (mean) imueclt2 [pweight = newweight] if nondemo==0, over(trstep, relabel(1"信頼していない" 2" " 3" " 4" " 5" " 6"どちらでもない" 7" " 8" " 9" " 10" " 11"信頼している")) ytitle("") ylabel(0(0.1)0.5) title(欧州議会への信頼度と文化的反移民感情保持確率) scheme(s1mono)

graph bar (mean) imueclt2 [pweight = newweight] if nondemo==0, over(imptrad1, relabel(1"伝統を重視しない" 2" " 3" " 4" " 5" " 6"伝統を重視する")) ytitle("") ylabel(0(0.1)0.5) title(伝統や習慣への追随志向と文化的反移民感情保持確率) scheme(s1mono)

graph bar (mean) imueclt2 [pweight = newweight] if nondemo==0, over(ipfrule1, relabel(1"規則を重視しない" 2" " 3" " 4" " 5" " 6"規則を重視する")) ytitle("") ylabel(0(0.1)0.5) title(命令や規則への服従志向と文化的反移民感情保持確率) scheme(s1mono)


/*移民の人種的異同に伴う受け入れタイプ変数の作成*/
gen type1 = 0
replace type1 = 1 if imsmetn<=2 & imdfetn<=2
replace type1 = 2 if imsmetn<=2 & imdfetn>=3
replace type1 = 3 if imsmetn>=3 & imdfetn>=3
replace type1 = 4 if imsmetn>=3 & imdfetn<=2
replace type1 = 0 if imsmetn==.
replace type1 = 0 if imdfetn==.

/* 
type1 1 = 全人種移民受入
type1 2 = 異人種移民拒否racist
type1 3 = 全人種移民拒否xenophobia
type1 4 = 異人種移民受入xenophile 
*/

/*図の作成*/
tabulate cntry type1 if nondemo==0, row

/*多項ロジット分析に関係しない回答者を欠損値に*/
recode type1 ( 0 = .)
recode type1 ( 4 = .)

/*多項ロジット分析 (表3.2 p111)*/
mlogit type1 c.agea c.gndr c.hinctnta c.edulvlb i.mnactic c.density16x c. reuemp17 i.ccode [pweight = newweight]if nondemo==0, vce(cluster ccode)
listcoef
mlogit type1 c.agea c.gndr c.hinctnta c.edulvlb i.mnactic c.density16x c. reuemp17 stflife1 stfeco1 stfdem1 c.euftf1 c.imptrad1 c.ipfrule1 c.atchctr c.impfree1 i.ccode [pweight = newweight]if nondemo==0, vce(cluster ccode)
listcoef


####第５章の分析一覧####
/*反欧州統合ダミー変数の作成*/
gen euftf2 = euftf1
recode euftf2 (0/5 = 0)
recode euftf2 (6/10 = 1)

/*右翼支持・反移民・反EUの三角関係の分析（3章用・4章用のデータ処理が必要です）(補遺4 p.262)*/
tetrachoric euftf2 imueclt2 frcl2 if nondemo==0 & cntry=="AT", stats(rho p)
tetrachoric euftf2 imueclt2 frcl2 if nondemo==0 & cntry=="BE", stats(rho p)
tetrachoric euftf2 imueclt2 frcl2 if nondemo==0 & cntry=="BG", stats(rho p)
tetrachoric euftf2 imueclt2 frcl2 if nondemo==0 & cntry=="CY", stats(rho p)
tetrachoric euftf2 imueclt2 frcl2 if nondemo==0 & cntry=="CH", stats(rho p)
tetrachoric euftf2 imueclt2 frcl2 if nondemo==0 & cntry=="CZ", stats(rho p)
tetrachoric euftf2 imueclt2 frcl2 if nondemo==0 & cntry=="DE", stats(rho p)
tetrachoric euftf2 imueclt2 frcl2 if nondemo==0 & cntry=="EE", stats(rho p)
tetrachoric euftf2 imueclt2 frcl2 if nondemo==0 & cntry=="ES", stats(rho p)
tetrachoric euftf2 imueclt2 frcl2 if nondemo==0 & cntry=="FI", stats(rho p)
tetrachoric euftf2 imueclt2 frcl2 if nondemo==0 & cntry=="FR", stats(rho p)
tetrachoric euftf2 imueclt2 frcl2 if nondemo==0 & cntry=="GB", stats(rho p)
tetrachoric euftf2 imueclt2 frcl2 if nondemo==0 & cntry=="HR", stats(rho p)
tetrachoric euftf2 imueclt2 frcl2 if nondemo==0 & cntry=="HU", stats(rho p)
tetrachoric euftf2 imueclt2 frcl2 if nondemo==0 & cntry=="IT", stats(rho p)
tetrachoric euftf2 imueclt2 frcl2 if nondemo==0 & cntry=="LT", stats(rho p)
tetrachoric euftf2 imueclt2 frcl2 if nondemo==0 & cntry=="LV", stats(rho p)
tetrachoric euftf2 imueclt2 frcl2 if nondemo==0 & cntry=="NL", stats(rho p)
tetrachoric euftf2 imueclt2 frcl2 if nondemo==0 & cntry=="NO", stats(rho p)
tetrachoric euftf2 imueclt2 frcl2 if nondemo==0 & cntry=="PL", stats(rho p)
tetrachoric euftf2 imueclt2 frcl2 if nondemo==0 & cntry=="SE", stats(rho p)
tetrachoric euftf2 imueclt2 frcl2 if nondemo==0 & cntry=="SI", stats(rho p)
tetrachoric euftf2 imueclt2 frcl2 if nondemo==0 & cntry=="SK", stats(rho p)

/*国別の右翼政党支持分析（以下それぞれの国について country=="AT"の所を任意の国コードに代えて分析してください）(表4.1)*/
logistic c.frcl2 c.agea c.gndr c.hinctnta c.edulvlb i.mnactic c.density16x c. reuemp17 stflife1 stfeco1 stfdem1 c.euftf1 c.imptrad1 c.ipfrule1 c.atchctr c.impfree1 c.imbgeco1 c.imueclt1 [pweight = newweight]if cntry=="AT", vce(robust)
listcoef


/*国別の移民文化棄損懸念分析（以下それぞれの国について country=="AT"の所を任意の国コードに代えて分析してください）(表4.1)*/
logistic c.imueclt1 c.agea c.gndr c.hinctnta c.edulvlb i.mnactic c.density16x c. reuemp17 stflife1 stfeco1 stfdem1 c.euftf1 c.imptrad1 c.ipfrule1 c.atchctr c.impfree1 [pweight = newweight]if cntry=="AT", vce(robust)
listcoef