* Encoding: UTF-8.
* Encoding: .


*Recodes and labeling.
USE ALL.
COMPUTE filter_$=(GMPSVIND > 0).
VARIABLE LABELS filter_$ 'GMPSVIND > 0 (FILTER)'.
VALUE LABELS filter_$ 0 'Not Selected' 1 'Selected'.
FORMATS filter_$ (f1.0).
FILTER BY filter_$.
EXECUTE.

VALUE LABELS GMPSVIND 1 '	Ja, exact op locatie' 2 'Ja, maar GSV maakte aanvankelijk een andere suggestie' 
3 'ja, maar zeer ver van exacte locatie ' 4 '	nee => INELADRE'.


VALUE LABELS GSVVIND 
1	'ja, exact op locatie'
2	'ja, maar verder in de straat dan aangeduid door GSV'
3	'ja, maar obv directe buren want exacte huisnummer niet zichtbaar'
4	'ja, maar obv huisnummers paar huizen verder want exacte huisnummer niet zichtbaar'
5	'Ja, maar obv googlen publiek of commercieel gebouw nabij'
6	'ja, maar op basis van sterke gelijkenissen potentiële huizen want exacte huisnummer niet zichtbaar'
7	'Ja, maar gecensureerd (huis wazig)'
8	'Ja, maar physical boundaries'
9	'nee, andere=> INELADRE'
10	'nee, geen street view'.

VALUE LABELS GSV_INELADRE
1	     'Verwaarloosd of vernield huis/adres'
2	     'Nog niet gebouwd nog niet gereed voor bewoning'
3	      ' Adres is niet bewoond leeg, tweede huis, vakantiewoning…'
4	       'Adres is geen woonadres: alleen bedrijfs-/zakelijk doel'
5	Ž     'Adres is geen woonadres: instelling bejaardenhuis, ziekenhuis, militaire eenheid, klooster…'
6	       'Adres is niet traceerbaar, onvoldoende informatie over adres'
7	     'Andere geef details'.


VALUE LABELS GSV_HUISTYPE
1	'Boerderij'
2	'Vrijstaand huis'
3	'Halfvrijstaand huis, twee-onder-één-kap of hoekwoning'
4	'Tussenwoning, rijtjeshuis'
5	'De enige wooneenheid in een gebouw met een ander doel (commercieel eigendom), bedrijfswoning'
6	'Flat, appartement, boven- of benedenwoning'
7	'Kamerbewoning, studentenkamer of studentenappartement' 
8	'Seniorenwoning, aanleunwoning of serviceflat'
9	'Woonwagen, caravan of boot'
10	'Andere'
88	'Weet niet'.

recode GSV_HUISTYPE
(1 thru 88 = copy) (sysmis =88) into GSV_HUISTYPECOR.
EXECUTE.


VALUE LABELS GSV_APPARTYP
1	'klein (tem 3 verdiepingen)'
2	'mid-rise (tem 5 verdiepingen)'
3	'high-rise (vanaf 6+ verdiepingen'
4	'misschien een appartement'.


VALUE LABELS GSV_IMPEDIM
1	'Ja – parlofoon'
2	'Ja – gesloten poort/deur'
3	'Ja – parlofoon EN gesloten poort/deur'
4	'Nee – geen van beiden'
5	'Niet zichtbaar'.


VALUE LABELS GSV_STATHUIS
1	'In een zeer goede staat'
2	'In een goede staat'
3	'In voldoende staat'
4	'Slechte staat'
5	'Zeer slechte staat' 
6	'niet zichtbaar'.


VALUE LABELS GSV_STATBUHS
1	'veel beter'
2	'beter'
3	'vergelijkbaar'
4	'slechter'
5	'veel slechter'
6	'niet zichtbaar'.


VALUE LABELS GSV_OMGAVFAL
1	'Zeer veel'
2	'Redelijk veel'
3	'Niet erg veel'
4	'Helemaal geen'.


VALUE LABELS GSV_OMGVANDA
1	'Zeer veel'
2	'Redelijk veel'
3	'Niet erg veel'
4	'Helemaal geen'.


VALUE LABELS OORTTYPE
1	'verstedelijkt (nabij kern van stad of gemeente)'
2	'suburbs'
3	'dorpskern'
4	'landelijk'
5	'industriezone of KMO-zone'.


if intnum2 =  666666 MULTIIV = 1.
if intnum2 < 666666  &  intnum3 > 666665 MULTIIV = 2.
if intnum3 < 666666  MULTIIV = 3.

recode MULTIIV (1=1) (2 thru 3 =2) into MULTIIV_DICH.
execute.

recode GSVVIND (1 thru 8 = 1) (9 thru 10 = 2) into GSVVINDDICH.
execute. 
value labels GSVVINDDICH 1'found in Google' 2 'not found'.


recode GSVVIND (1 thru 3 = 1) (4 thru 6 = 2) (7 = 3) (8 thru 87 = 4)
into GSVVIND_4CAT.

recode GSVJAAR (2007 thru 2015 = copy) (missing= 9999) into GSVJAARMIS.



*Table 1.
USE ALL.
FREQUENCIES VARIABLES=ESS7OUTCOME
  /STATISTICS=MODE
  /ORDER=ANALYSIS.

*Table 2.

CROSSTABS
  /TABLES=  GSVVIND   BY GSVJAAR
  /FORMAT=AVALUE TABLES
  /STATISTICS=CHISQ GAMMA 
  /CELLS=COUNT  
  /COUNT ROUND CELL.


*Table 3.

CROSSTABS
  /TABLES=type_recode by  GSV_HUISTYPE_recode
  /FORMAT=AVALUE TABLES
  /STATISTICS=CHISQ GAMMA 
  /CELLS=COUNT  
  /COUNT ROUND CELL.



CROSSTABS
  /TABLES=type BY  GSV_HUISTYPE 
  /FORMAT=AVALUE TABLES
  /STATISTICS=CHISQ GAMMA 
  /CELLS=COUNT ROW 
  /COUNT ROUND CELL.
CROSSTABS
  /TABLES=physa BY  GSV_STATHUIS 
  /FORMAT=AVALUE TABLES
  /STATISTICS=CHISQ GAMMA 
  /CELLS=COUNT ROW 
  /COUNT ROUND CELL.
CROSSTABS
  /TABLES=littera BY  GSV_OMGAVFAL 
  /FORMAT=AVALUE TABLES
  /STATISTICS=CHISQ GAMMA 
  /CELLS=COUNT ROW 
  /COUNT ROUND CELL.
CROSSTABS
  /TABLES=vandaa BY  GSV_OMGVANDA 
  /FORMAT=AVALUE TABLES
  /STATISTICS=CHISQ GAMMA 
  /CELLS=COUNT ROW 
  /COUNT ROUND CELL.
CROSSTABS
  /TABLES=access BY   GSV_IMPEDIM  
  /FORMAT=AVALUE TABLES
  /STATISTICS=CHISQ GAMMA 
  /CELLS=COUNT ROW 
  /COUNT ROUND CELL.



recode littera GSV_OMGAVFAL vandaa GSV_OMGVANDA 
(1 thru 3 = 1) ( 4 = 2) (6=6) into
littera_dich GSV_OMGAVFAL_dich vandaa_dich GSV_OMGVANDA_dich.
execute.
VALUE LABELS littera_dich GSV_OMGAVFAL_dich vandaa_dich GSV_OMGVANDA_dich
1 'iets' 2' geen'.

recode physa   GSV_STATHUIS 
(1 thru 2 = 1) (3=2) (4 thru 5= 3) into physa_tri  GSV_STATHUIS_tri. 
execute.
VALUE LABELS physa_dich  GSV_STATHUIS_tri 1 'goed/zeer goed' 2'voldoende staat' 3'slecht/zeer slecht'.
recode physa   GSV_STATHUIS 
(1 thru 2 = 1) (3 thru 5= 0) into physa_dich  GSV_STATHUIS_dich. 
execute.
VALUE LABELS physa_dich  GSV_STATHUIS_dich 1 'goed/zeer goed' 0'slecht/zeer slecht/voldoende staat'.


recode access GSV_IMPEDIM  
(1 thru 3 =1) (4 =2) (5 thru 6=3) into access_dich GSV_IMPED_dich.
execute.
VALUE LABELS access_dich GSV_IMPED_dich 1'ja' 2'nee' 3'niet zichtbaar'.



recode type  GSV_HUISTYPE (1 thru 2 = 1) (6 thru 7=6) (else=copy) 
into type_recode  GSV_HUISTYPE_recode.
execute.

recode type  GSV_HUISTYPE (1 thru 3 = 1) (6 thru 8 =6) (else=copy) 
into type_reco2  GSV_HUISTYPE_reco2.
execute.



USE ALL.
COMPUTE filter_$=(missing (VAL_AFVAL)=0 & missing (VAL_VANDA)=0 & missing (VAL_PHYSA)=0 & missing (VAL_HTYP2)=0 ).
VARIABLE LABELS filter_$ 'missing (VAL_AFVAL)=0 & missing (VAL_VANDA)=0 & missing (VAL_PHYSA)=0 & missing (VAL_HTYP2)=0  (FILTER)'.
VALUE LABELS filter_$ 0 'Not Selected' 1 'Selected'.
FORMATS filter_$ (f1.0).
FILTER BY filter_$.
EXECUTE.

CROSSTABS
  /TABLES=type_recode by  GSV_HUISTYPE_recode
  /FORMAT=AVALUE TABLES
  /STATISTICS=CHISQ GAMMA 
  /CELLS=COUNT column 
  /COUNT ROUND CELL.
CROSSTABS
  /TABLES=type_reco2 by  GSV_HUISTYPE_reco2
  /FORMAT=AVALUE TABLES
  /STATISTICS=CHISQ GAMMA 
  /CELLS=COUNT  
  /COUNT ROUND CELL.



CROSSTABS
  /TABLES=physa_dich BY  GSV_STATHUIS_dich 
  /FORMAT=AVALUE TABLES
  /STATISTICS=CHISQ GAMMA 
  /CELLS=COUNT column 
  /COUNT ROUND CELL.
CROSSTABS
  /TABLES=littera_dich BY  GSV_OMGAVFAL_dich 
  /FORMAT=AVALUE TABLES
  /STATISTICS=CHISQ GAMMA 
  /CELLS=COUNT column 
  /COUNT ROUND CELL.
CROSSTABS
  /TABLES=vandaa_dich BY  GSV_OMGVANDA_dich 
  /FORMAT=AVALUE TABLES
  /STATISTICS=CHISQ GAMMA 
  /CELLS=COUNT column 
  /COUNT ROUND CELL.
CROSSTABS
  /TABLES=access_dich BY   GSV_IMPED_dich  
  /FORMAT=AVALUE TABLES
  /STATISTICS=CHISQ GAMMA 
  /CELLS=COUNT column 
  /COUNT ROUND CELL.

use all.



DO IF  VAL_AFVAL < 99 &  VAL_VANDA < 99 & VAL_PHYSA < 99 & VAL_HTYP2 <99.
COUNT VAL_GSV4=VAL_AFVAL VAL_VANDA VAL_PHYSA HUISMATCH4(1).
VARIABLE LABELS  VAL_GSV4 'validated 4 gsv obs'.
END IF.
EXECUTE.

DO IF VAL_AFVAL < 99 &  VAL_VANDA < 99 & VAL_PHYSA < 99 & VAL_HTYP2 < 99 & (missing (VAL_IMPED)=0).
COUNT VAL_GSV5=VAL_AFVAL VAL_VANDA VAL_PHYSA VAL_HTYP2 VAL_IMPED (1).
VARIABLE LABELS  VAL_GSV5 'validated 5 gsv obs including imped'.
END IF.
EXECUTE.



USE ALL.
COMPUTE filter_$=(VAL_AFVAL < 99 &  VAL_VANDA < 99 & VAL_PHYSA < 99 & VAL_HTYP2 < 99  ).
VARIABLE LABELS filter_$ 'missing 5GSV=0 (FILTER)'.
VALUE LABELS filter_$ 0 'Not Selected' 1 'Selected'.
FORMATS filter_$ (f1.0).
FILTER BY filter_$.
EXECUTE.

FREQUENCIES VARIABLES=  VAL_GSV4  
  /STATISTICS=MODE
  /ORDER=ANALYSIS.


USE ALL.
COMPUTE filter_$=(VAL_AFVAL < 99 &  VAL_VANDA < 99 & VAL_PHYSA < 99 & VAL_HTYP2 < 99 & missing (GSV_IMPEDIM)=0 ).
VARIABLE LABELS filter_$ 'missing 5GSV=0 (FILTER)'.
VALUE LABELS filter_$ 0 'Not Selected' 1 'Selected'.
FORMATS filter_$ (f1.0).
FILTER BY filter_$.
EXECUTE.
FREQUENCIES VARIABLES=    VAL_GSV5
  /STATISTICS=MODE
  /ORDER=ANALYSIS.
use all.


*Table 4 and 5.

CROSSTABS
  /TABLES= VAL_HTYP2  by  GSVVIND GSVJAAR
  /FORMAT=AVALUE TABLES
  /STATISTICS=CHISQ GAMMA 
  /CELLS=COUNT column 
  /COUNT ROUND CELL.

*Table 6.

*filter conditional on contact = IV or refusal.
USE ALL.
COMPUTE filter_$=( missing (GSV_STATHUIS_dich)=0 & missing (GSV_OMGAVFAL_dich)=0 & missing (GSV_OMGVANDA_dich)=0  & ESS7OUTCOME =1 | ESS7OUTCOME =3 ).
VARIABLE LABELS filter_$ 'NO missing 3GSV  (FILTER)'.
VALUE LABELS filter_$ 0 'Not Selected' 1 'Selected'.
FORMATS filter_$ (f1.0).
FILTER BY filter_$.
EXECUTE.



LOGISTIC REGRESSION VARIABLES  RR1AAPOR
   /METHOD=ENTER  GSV_OMGAVFAL_dich GSV_OMGVANDA_dich GSV_STATHUIS_dich 
 /CONTRAST (GSV_OMGAVFAL_dich)=Indicator
 /CONTRAST (GSV_OMGVANDA_dich)=Indicator
 /CONTRAST (GSV_STATHUIS_dich)=Indicator
 /CONTRAST (GSV_IMPED_dich)=Indicator
  /CASEWISE OUTLIER(3)
  /CRITERIA=PIN(0.05) POUT(0.10) ITERATE(20) CUT(0.5).
LOGISTIC REGRESSION VARIABLES  RR1AAPOR
    /METHOD=ENTER  littera_dich vandaa_dich physa_dich  
 /CONTRAST (access_dich)=Indicator
 /CONTRAST (littera_dich)=Indicator
 /CONTRAST (vandaa_dich)=Indicator
 /CONTRAST (physa_dich)=Indicator
  /CASEWISE OUTLIER(3)
  /CRITERIA=PIN(0.05) POUT(0.10) ITERATE(20) CUT(0.5).



USE ALL.
COMPUTE filter_$=( missing (GSV_STATHUIS_dich)=0 & missing (GSV_OMGAVFAL_dich)=0 & missing (GSV_OMGVANDA_dich)=0 & missing (GSV_IMPED_dich)=0  ).
VARIABLE LABELS filter_$ 'NO missing 3GSV  (FILTER)'.
VALUE LABELS filter_$ 0 'Not Selected' 1 'Selected'.
FORMATS filter_$ (f1.0).
FILTER BY filter_$.
EXECUTE.
LOGISTIC REGRESSION VARIABLES  RR1AAPOR
   /METHOD=ENTER  GSV_OMGAVFAL_dich GSV_OMGVANDA_dich GSV_STATHUIS_dich GSV_IMPED_dich
 /CONTRAST (GSV_OMGAVFAL_dich)=Indicator
 /CONTRAST (GSV_OMGVANDA_dich)=Indicator
 /CONTRAST (GSV_STATHUIS_dich)=Indicator
 /CONTRAST (GSV_IMPED_dich)=Indicator
  /CASEWISE OUTLIER(3)
  /CRITERIA=PIN(0.05) POUT(0.10) ITERATE(20) CUT(0.5).
LOGISTIC REGRESSION VARIABLES  RR1AAPOR
    /METHOD=ENTER  littera_dich vandaa_dich physa_dich  access_dich
 /CONTRAST (access_dich)=Indicator
 /CONTRAST (littera_dich)=Indicator
 /CONTRAST (vandaa_dich)=Indicator
 /CONTRAST (physa_dich)=Indicator
  /CASEWISE OUTLIER(3)
  /CRITERIA=PIN(0.05) POUT(0.10) ITERATE(20) CUT(0.5).

LOGISTIC REGRESSION VARIABLES  REFUSALAAPOR
   /METHOD=ENTER  GSV_OMGAVFAL_dich GSV_OMGVANDA_dich GSV_STATHUIS_dich /METHOD=ENTER GSV_IMPED_dich
 /CONTRAST (GSV_OMGAVFAL_dich)=Indicator
 /CONTRAST (GSV_OMGVANDA_dich)=Indicator
 /CONTRAST (GSV_STATHUIS_dich)=Indicator
 /CONTRAST (GSV_IMPED_dich)=Indicator
  /CASEWISE OUTLIER(3)
  /CRITERIA=PIN(0.05) POUT(0.10) ITERATE(20) CUT(0.5).
LOGISTIC REGRESSION VARIABLES  REFUSALAAPOR
    /METHOD=ENTER  littera_dich vandaa_dich physa_dich  /METHOD=ENTER access_dich
 /CONTRAST (access_dich)=Indicator
 /CONTRAST (littera_dich)=Indicator
 /CONTRAST (vandaa_dich)=Indicator
 /CONTRAST (physa_dich)=Indicator
  /CASEWISE OUTLIER(3)
  /CRITERIA=PIN(0.05) POUT(0.10) ITERATE(20) CUT(0.5).

USE ALL.
COMPUTE filter_$=( missing (GSV_STATHUIS_dich)=0 & missing (GSV_OMGAVFAL_dich)=0 & missing (GSV_OMGVANDA_dich)=0 & missing (GSV_IMPED_dich)=0 ).
VARIABLE LABELS filter_$ 'NO missing 3GSV  (FILTER)'.
VALUE LABELS filter_$ 0 'Not Selected' 1 'Selected'.
FORMATS filter_$ (f1.0).
FILTER BY filter_$.
EXECUTE.

LOGISTIC REGRESSION VARIABLES  CONTACTAAPOR
   /METHOD=ENTER  GSV_OMGAVFAL_dich GSV_OMGVANDA_dich GSV_STATHUIS_dich /METHOD=ENTER GSV_IMPED_dich 
 /CONTRAST (GSV_OMGAVFAL_dich)=Indicator
 /CONTRAST (GSV_OMGVANDA_dich)=Indicator
 /CONTRAST (GSV_STATHUIS_dich)=Indicator
 /CONTRAST (GSV_IMPED_dich)=Indicator
  /CASEWISE OUTLIER(3)
  /CRITERIA=PIN(0.05) POUT(0.10) ITERATE(20) CUT(0.5).

LOGISTIC REGRESSION VARIABLES  REFUSALAAPOR
   /METHOD=ENTER  GSV_OMGAVFAL_dich GSV_OMGVANDA_dich GSV_STATHUIS_dich  /METHOD=ENTER GSV_IMPED_dich
 /CONTRAST (GSV_OMGAVFAL_dich)=Indicator
 /CONTRAST (GSV_OMGVANDA_dich)=Indicator
 /CONTRAST (GSV_STATHUIS_dich)=Indicator
 /CONTRAST (GSV_IMPED_dich)=Indicator
  /CASEWISE OUTLIER(3)
  /CRITERIA=PIN(0.05) POUT(0.10) ITERATE(20) CUT(0.5).

LOGISTIC REGRESSION VARIABLES  CONTACTAAPOR
    /METHOD=ENTER  littera_dich vandaa_dich physa_dich  /METHOD=ENTER access_dich
 /CONTRAST (access_dich)=Indicator
 /CONTRAST (littera_dich)=Indicator
 /CONTRAST (vandaa_dich)=Indicator
 /CONTRAST (physa_dich)=Indicator
  /CASEWISE OUTLIER(3)
  /CRITERIA=PIN(0.05) POUT(0.10) ITERATE(20) CUT(0.5).