##########################################################################################
#
# Auswertungsdatei für die Durchführung eines Auswertungsgangs mit der Forschungsdatenbank
# Nichtwohngebäude (Projekt "Forschungsdatenbank Nichtwohngebäude (ENOB:dataNWG)“,
# gefördert durch das Bundesministerium für Wirtschaft und Energie;Förderkennzeichen:
# 03ET1315) 
#
#########################################
#
#Bitte folgende Angaben vervollständigen:
#
#Dateiname des Programmcodes: <Name>.txt 
#erstellt: <Datum> 
#von: <Name> 
#E-Mail: <E-Mail-Adresse> 
#Tel.: <Telefonnummer> 
# 	
#########################################


###########################################################
# 1. Laden des Headers
#	 Lädt Gebäudedaten, stichprobenrelevante Daten 
#	 inkl. der Schätzfunktionen in den Arbeitsspeicher 
###########################################################
source('//Dateiserver/energie/1075_ENOBdataNWG/02_work_IWU/AP1_Konzeptphase/AP1.4_Erhebungsmethodik/AP1.4.1_Datenbank/Umsetzung_Auswertung/05_Schaetzer_BRE/20210509_header_bre.r')
# Verzeichnis für Exportdateien
pfad <- 'C:/Users/Nutzer/Desktop/'



##########################################
# 2. Festlegen des Arbeitsverzeichnisses
##########################################
setwd('//Dateiserver/energie/1075_ENOBdataNWG/02_work_IWU/AP4_Auswertung/')
getwd()


###########################################################################################
# 3. Aufbereiten von dichotomen Merkmalen und
#    Aufbereiten der Ausprägungen von uk_geb
#    NICHT ÄNDERN!
###########################################################################################
gebaeude <- dichot_merkm_aufbereiten()



############### Hinweis ################################################################
# VARIABLENNAMEN SIND UNABHÄNGIG VON DER SCHREIBWEISE IN DER VARIABLENDOKUMENTATION
# KLEIN ZU SCHREIBEN!
########################################################################################


###########################################################################################
# 4. Ggf. neue Merkmale definieren
# Hinweis: Nur fehlende Werte dürfen negativ codiert sein.
###########################################################################################
# Aggregierte Hauptgebäudefunktion hk_geb_agg
gebaeude$hk_geb_agg <- ((gebaeude$hk_geb < 0) * -1 
                      + ((gebaeude$hk_geb >= 1 & gebaeude$hk_geb <= 7) | gebaeude$hk_geb == 9) * 1 
                      + (gebaeude$hk_geb == 8 | gebaeude$hk_geb == 10 | gebaeude$hk_geb == 11) * 2
)

# Klassifizierte oberirdische BGF ober_bgf_kl
gebaeude$ober_bgf_kl <- ((gebaeude$ober_bgf < 0 ) * -1 
                       + (gebaeude$ober_bgf >= 0 & gebaeude$ober_bgf < 100 ) * 1          
                       + (gebaeude$ober_bgf >=100 &  gebaeude$ober_bgf < 1000) * 2
                       + (gebaeude$ober_bgf >=1000 &  gebaeude$ober_bgf < 5000) * 3  
                       + (gebaeude$ober_bgf >=5000 &  gebaeude$ober_bgf < 10000) * 4
                       + (gebaeude$ober_bgf >=10000 &  gebaeude$ober_bgf < 25000) * 5
                       + (gebaeude$ober_bgf >=25000 &  gebaeude$ober_bgf < 50000) * 6
                       + (gebaeude$ober_bgf >=50000 ) * 7)


#########################################################################################
# 5. Definition von Teilmengen
#########################################################################################
grundlage1 <- subset(gebaeude,enev_rel_nwg_be == 1)
grundlage2 <- subset(grundlage1,energietraeger >= 0 & hk_geb_agg >= 0)
grundlage3 <- subset(gebaeude,enev_rel_nwg_be == 1 &
                    uk_geb %in% c(100,101,102,103,104,105,202,205,306,307,308))
 
#############################################################################################
# 6. Auswertungen
#############################################################################################
Zellenprozente(grundlage2, pfad, "energietraeger", "hk_geb_agg")
Spaltenprozente(grundlage2, pfad, "energietraeger", "hk_geb_agg")
Mittelwert(grundlage3, pfad, "ober_bgf", "ober_bgf_kl", "leer")