Schatting BASGOED. rapportage DP1 SIGNIFICANCE: GERARD DE JONG MICHIEL DE BOK KIM RUIJS DERK WENTINK

Schatting BASGOED rapportage DP1 SIGNIFICANCE: GERARD DE JONG MICHIEL DE BOK KIM RUIJS DERK WENTINK NEA: ARNAUD BURGESS MARLEEN DE GIER ROBBERT VERSTEEGH DEMIS: POUL GRASHOFF

3 november 2010 Rapport voor de Dienst Verkeer en Scheepvaart van Rijkswaterstaat

Significance Koninginnegracht 23 2514 AB DEN HAAG T (070) 3121 543 F (070) 3121 531 [email protected] www.significance.nl

Inhoudsopgave

Inhoudsopgave .............................................................................................................2 Samenvatting................................................................................................................5 HOOFDSTUK 1 Inleiding ........................................................................................7 1.1 Algemeen ...........................................................................................................7 1.2 Doelstellingen ....................................................................................................7 1.3 Indeling van dit rapport .....................................................................................8 HOOFDSTUK 2 Gebruikte data ...............................................................................9 2.1 Basisbestand goederenvervoer.............................................................................9 2.2 Level of service .................................................................................................10 2.2.1 LMS....................................................................................................10 2.2.2 Routgoed ............................................................................................10 2.2.3 BIVAS.................................................................................................11 2.2.4 TRANSTOOLS..................................................................................12 2.3 Vergelijkingskader modaliteiten .......................................................................12 2.4 Voedingskosten binnenvaart en spoor ..............................................................13 2.4.1 Spoor ..................................................................................................14 2.4.2 Binnenvaart.........................................................................................15 HOOFDSTUK 3 Schattingen modal split ................................................................17 3.1 Kostenfuncties Modal split...............................................................................17 3.2 Overzicht schattingen.......................................................................................20 3.3 Gedetailleerde schattingsresultaten ...................................................................22 3.4 Reistijdwaarderingen ........................................................................................26 3.5 Elasticiteiten.....................................................................................................27 3.6 Goodness-of-fit ................................................................................................28 3.7 Conclusies........................................................................................................40 HOOFDSTUK 4 Schatting distributiemodel ............................................................43 4.1 Werking van de distributiemodule ...................................................................43 4.2 Schatten parameters .........................................................................................45 4.3 Resultaten ........................................................................................................46 4.4 Conclusies........................................................................................................50

Significance, NEA en DEMIS

Inhoudsopgave

HOOFDSTUK 5 Aandachtspunten voor de implementatie .......................................51 5.1 Achtergrond alternatief 1b................................................................................51 5.2 Stand van zaken prototype modules DM en MS ..............................................53 5.3 Mogelijkheden voor implementatie modules DM en MS in dp2......................54 HOOFDSTUK 6

Conclusies en aanbevelingen .........................................................57

Appendix 1

Elasticiteiten per NSTR1 groep ........................................................59

Appendix 2

schattingsresultaten MS serie 13 .......................................................63

Appendix 3

schattingsresultaten DM...................................................................69

3

Samenvatting

Om tegemoet te komen aan de informatiebehoefte van het Ministerie van Verkeer en Waterstaat op het gebied van goederenvervoer in, naar en van Nederland is een goederenvervoermodel nodig dat niet alleen inhoudelijk antwoord geeft, maar ook tegemoet komt aan randvoorwaarden op het gebied van transparantie, uitlegbaarheid en de mogelijkheid tot validatie. Dit is momenteel niet beschikbaar. DVS voert samen met TNO het project ‘Roadmap strategisch goederenvervoermodel’ uit, waarin vraag en aanbod op dit gebied met elkaar vergeleken zijn. Op basis hiervan is besloten om op korte termijn een relatief eenvoudig goederenvervoermodel te ontwikkelen, BasGoed, waarin de bestaande economiemodule van SMILE+ wordt hergebruikt en waarin er nieuwe modellen geschat zullen worden voor distributie en modal split, op basis van tijden en afstanden uit bestaande unimodale toedelingsmodellen. De functionele en technische specificaties van BasGoed staan beschreven in een rapport van TNO, NEA, Significance en Demis (Duijnisveld et al., 2010). Hierin zijn allerlei keuzen voor BasGoed reeds vastgelegd (bijvoorbeeld de te gebruiken zonering). Deze rapportage beschrijft de uitvoering van deelproject 1 in de ontwikkeling van het te ontwikkelen Basis Goederenvervoermodel: het schatten van de modal split en distributiemodule voor BASGOED. Hierin zijn schattingsdata verzameld en verwerkt, en schatting uitgevoerd voor de modal split module en distributiemodule, en worden tenslotte aanbevelingen en richtlijnen geformuleerd voor deelproject 2: de implementatie van modellen in de Modalsplit en Dsitributie Module. Bij de schattingen zijn verschillende databronnen gebruikt. De waarnemingen zijn afkomstig uit het Basisbestand Goederenvervoer. De level of service voor de verschillende modaliteiten is afkomstig uit verschillende modellen: het LMS voor wegvervoer, RoutGoed voor spoorvervoer, BIVAS voor binnenvaart en TRANSTOOLS voor internationaal vervoer. Daarnaast is het vergelijkingskader modaliteiten gebruikt om kostenkentallen af te leiden. Het specificatieonderzoek heeft geleid tot een modal split model met een kostenfunctie die is gebaseerd op de kostenkentallen uit het vergelijkingskader modaliteiten (VKM). Hiermee is de aanpak consistent met de kostenfuncties in SMILE+. Per NSTR hoofdgroep zijn twee verschillende specificaties uitvoeriger getest. Een specificatie met herkomst en bestemmingsdummies per vervoerwijze (serie 13), en een specificatie met vervoerwijzedummies op corridors (serie 14). Op basis van de beoordeling op de drie hoofdcriteria is model serie 13 met vervoerwijze specifieke herkomst bestemmingsdummies duidelijk te prefereren als model voor implementatie in de modal split module van BASGOED. Ten aanzien van de distributie module zijn verschillende varianten van de distributiefunctie onderzocht. Hierin is gekeken naar verschillende aannames ten aanzien van het meenemen van parameters α , β en de grensweerstanden Gij en intrazonale weerstanden d. Op basis van de resultaten variant 5b, met een negatief exponentiële distributiefuncties en intrazonale weerstanden (maar geen grensweerstanden), naar voren als te prefereren vorm voor implementatie in de distributie module in BASGOED. In geval van NSTR groep 4 wordt teruggevallen op variant 1, zonder de intrazonale weerstanden.


Schattingen BASGOED

De geschatte modellen zijn verder vertaald naar een functionele beschrijving voor de modal split en distributie module. In DP2 zullen deze rekenmodules verder worden uitgewerkt. In DP3 zullen deze modules worden geïmplementeerd in het BASGOED model systeem.

6

HOOFDSTUK 1

1.1

Inleiding

Algemeen

Om tegemoet te komen aan de informatiebehoefte van het Ministerie van Verkeer en Waterstaat op het gebied van goederenvervoer in, naar en van Nederland is een goederenvervoermodel nodig dat niet alleen inhoudelijk antwoord geeft, maar ook tegemoet komt aan randvoorwaarden op het gebied van transparantie, uitlegbaarheid en de mogelijkheid tot validatie. Dit is momenteel niet beschikbaar. DVS voert samen met TNO het project ‘Roadmap strategisch goederenvervoermodel’ uit, waarin vraag en aanbod op dit gebied met elkaar vergeleken zijn. Op basis hiervan is besloten om op korte termijn een relatief eenvoudig goederenvervoermodel te ontwikkelen, waarin de bestaande economiemodule van SMILE+ wordt hergebruikt en waarin er nieuwe modellen geschat zullen worden voor distributie en modal split, op basis van tijden en afstanden uit bestaande unimodale toedelingsmodellen. Dit model heeft de naam ‘Basis Goederenvervoermodel’ (BasGoed) gekregen. In een later stadium zal een model worden ontwikkeld dat rekening houdt met logistieke besluitvorming. De functionele en technische specificaties van BasGoed staan beschreven in een rapport van TNO, NEA, Significance en Demis (Duijnisveld et al., 2010). Hierin worden allerlei keuzen voor BasGoed reeds vastgelegd (bijvoorbeeld de te gebruiken zonering). Wel worden nog diverse modelspecificaties onderscheiden (met name voor het modal split model), die in de schatting met elkaar moeten worden vergeleken, waarna een keuze dient te worden gemaakt. Dit rapport is onderdeel van twee deelprojecten binnen de implementatie van BasGoed: Deelproject 1: 1 schatting1 van de distributiemodule (DM) en de modal split module (MS); Deelproject 2: 2 implementatie (programmeren, testen, documenteren) van de distributiemodule en de modal split module. Naast deze deelprojecten vindt parallel een apart deelproject plaats voor de ontwikkeling van de gebruikersschil voor BasGoed en de andere modules.

1.2

Doelstellingen

Het doel van dit project is: Het schatten van nieuwe modellen voor distributie en modal split volgens de specificaties in Duijnisveld et al. (2010), alsmede het programmeren, testen en documenteren van de nieuwe modules, als onderdeel van het te ontwikkelen Basis Goederenvervoermodel

1

Wij prefereren in deze context het woord ’schatten’ boven het woord ‘calibreren’, omdat gebruik zal worden gemaakt van formele statistische schattingstechnieken (zoals Maximum Likelihood).


Schattingen BASGOED

Het project is opgesplitst in twee deelprojecten: Deelproject 1: verzamelen en verwerken van schattingsdata, schatting van modal split module en distributiemodule en rapportage van schattingsresultaten; Deelproject 2: verfijnen van de specificatie ter voorbereiding van het programmeren, programmeren van de distributiemodule en modal split module, testen van deze modules en documentatie van de modules en de testresultaten. Deze rapportage beschrijft de uitvoering van deelproject 1: het schatten van de modal split en distributiemodule voor BASGOED.

1.3

Indeling van dit rapport

In hoofdstuk 2 staat de verwerking van databronnen ten behoeve van de BasGoed schattingen beschreven. De schattingen van de modal split module staat beschreven in hoofdstuk 3. De schattingen van de distributiemodule staat beschreven in hoofdstuk 4. In hoofdstuk 5 zijn aandachtspunten bij de implementatie van beide modules beschreven.

8

HOOFDSTUK 2

Gebruikte data

Bij de schatting van de modal split en distributiemodules zijn verschillende databronnen gebruikt. De waarnemingen zijn afkomstig uit het Basisbestand Goederenvervoer. De level of service voor de verschillende modaliteiten is afkomstig uit verschillende modellen: het LMS voor wegvervoer, RoutGoed voor spoorvervoer, BIVAS voor binnenvaart en TRANSTOOLS voor internationaal vervoer. Daarnaast is het vergelijkingskader modaliteiten gebruikt om kostenkentallen af te leiden. 2.1

Basisbestand goederenvervoer

Met enige regelmaat worden de basisbestanden goederenvervoer in opdracht van DVS gemaakt. Deze zijn gebaseerd op data van het Centraal Bureau voor de Statistiek, aangevuld met data van andere organisaties (zoals vervoerders en buitenlandse statistische bureaus). De basisbestanden goederenvervoer zijn beschikbaar voor de modaliteiten weg, spoor en binnenvaart en bevatten verschillende kenmerken. In onderstaande tabel volgt een kort overzicht van de gebruikte kenmerken van de verschillende bestanden. Tabel 2-1: Dimensies in het basisbestand goederenvervoer 2004.

Weg

Spoor

Binnenvaart

Herkomst/Bestemming op SMILE-niveau Voertuigtype (vrachtauto, vrachtauto + aanhanger, trekker met oplegger, bestelauto) Verschijningsvorm (droge bulk, natte bulk, stukgoed, container) Goederengroep (NSTR 1-digit) Herkomst/Bestemming op SMILE-niveau Tractie (diesel/elektrisch) Verschijningsvorm (droge bulk, natte bulk, stukgoed, container) Goederengroep (NSTR 1-digit) Herkomst/Bestemming op SMILE-niveau Scheepsklasse volgens BVMS indeling Versvorm (droge bulk, natte bulk, container) Goederengroep (NSTR 1-digit)

Het geografisch detail niveau van de trip-matrices die gebruikt zijn in de schattingen is corop-niveau, conform de technische en functionele specificaties BASGOED (Duijnisveld et al., 2010). Voor het opstellen van trip-matrices zijn de basisbestanden goederenvervoer 2004 voor de modaliteiten spoor, binnenvaart en weg gebruikt. In de trip matrices is onderscheid gemaakt tussen vervoerde tonnages in containers en overige goederen. Verder zijn de goederenstromen geaggregeerd naar NSTR hoofdgroepen (NSTR1).


Schattingen BASGOED

De basisbestanden goederenvervoer zijn tevens gebruikt om vervoerwijze specifieke kentallen af te leiden uit het Vergelijkingskader Modaliteiten (VKM), zie ook paragraaf 2.3. De kentallen zijn de gemiddelde variabele, vaste, laad- en loskosten, het gemiddeld laadvermogen, de gemiddelde beladingsgraad en het percentage beladen ritten. In het VKM zijn deze kentallen opgenomen, uitgesplitst naar voertuigtype, verschijningsvorm en richting (binnenlands/internationaal). Deze kentallen zijn gekoppeld aan de basisbestanden goederenvervoer en vervolgens is een gewogen gemiddelde bepaald voor elk van deze kentallen op basis van het vervoerde gewicht. 2.2

Level of service 2.2.1

LMS

Het LMS berekent gecongesteerde reistijd en wegafstanden op LMS-subzone niveau op basis van een congestiegevoelige toedeling van vrachtwagens samen met personenverkeer. Voor het vrachtverkeer geldt in het LMS dat alleen de routes kunnen worden aangepast. BasGoed heeft geen terugkoppeling van congestie naar de vraag naar goederenvervoer. Binnen het LMS wordt geen onderscheid gemaakt naar goederensoorten bij de toedeling, zodat de reistijden en –afstanden voor alle goederensoorten gelijk aan elkaar zijn. Verder zijn een tweetal aggregatie slagen gemaakt om de LMS data aan te passen aan BasGoed: van LMS subzones naar BasGoed zones, en van dagdelen naar etmaal. Binnen het LMS worden reistijden en -afstanden voor verschillende perioden van de dag bepaald. Van deze reistijden en afstanden is het gewogen gemiddelde bepaald op basis van het aantal ritten. Ook bij de aggregatie van tijden afstand tussen LMS subzones, is het gemiddelde tussen BasGoed regio’s bepaald op basis van een gewogen gemiddelde. Het gewogen gemiddelde is als volgt berekend: n

∑t k ij =

ij

* k ij

i =1, j =1 n

∑t

ij

i =1, j =1

i

herkomst i

j

herkomst j

k ij

kosten verplaatsing i naar j (tijd of afstand)

t ij

aantal tours i naar j

2.2.2

Routgoed

Om de level of service te bepalen voor het spoorvervoer in Nederland is gebruik gemaakt van de tijden en afstanden uit het model Routgoed van Prorail. De zonering van Routgoed is de corop-zonering welke 1op-1 is overgenomen. Per corop-zone is door Prorail aangegeven op welk dienstregelpunt deze is aangesloten. Zo komt het voor dat er in een corop-zone geen dienstregelpunt beschikbaar is. Deze zones worden zijn aangesloten op het meest logische dienstregelpunt in een andere corop-zone. In Figuur 2-1 geeft het Routgoed netwerk weer. De rode lijnen geven weer op welk dienstregelpunt een corop-zone is aangesloten. De level of service is opgesteld aan de hand van Routgoed en is gebaseerd op rijtijden van de kortste route tussen herkomst en bestemming. Hierbij is geen rekening gehouden met wachttijden, los- en laadtijden en voedingskosten. Deze kosten worden apart toegevoegd binnen BasGoed aan de totale kosten van een reis.

10


Schattingen BASGOED

Figuur 2-1: Routgoed netwerk (grijs) en in rood aansluiting Corops op netwerk (groen)

2.2.3

BIVAS

De afstanden en tijden zijn verkregen uit het BIVAS model. In eerste instantie is de laatste versie van BIVAS gebruikt waarin nog een aantal tekortkomingen met betrekking tot diepte en vaarhoogte zijn opgenomen. Dit leidde tot een te hoge inschatting van reistijden en afstanden. De laatste versie van BIVAS is in ontwikkeling voor toepassingen in de sfeer van klimaatsverandering. Uiteindelijk is van Charta een “uitdraai” van een eerdere versie verkregen waarin en betrouwbare tijden en afstanden zijn opgenomen. In de tijden van BIVAS zijn wachttijden voor sluizen e.d. meegenomen. Alle reizen uit het basisbestand zijn opgenomen en daarbij is dus ook rekening gehouden met het scheepstype. De BIVAS zonering is omgezet naar SMILE+ zonering en vervolgens is het gemiddelde over alle reistijden en

11


Schattingen BASGOED

afstanden genomen, daarnaast is nog een minimum en een maximum afstand berekend. Er is een vergelijking gemaakt met de TRANS-TOOLS afstanden, de gemiddelde afstand laat de hoogste correlatie zien. Uit BIVAS is de afstand en tijd tussen de binnenlandse regio’s genomen (40 x 40 zones). Voor de afstand en tijd tussen binnenlandse regio’s (zones 1 tot en met 40) en buitenlandse regio’s (zones 41 tot en met 70) is TRANS-TOOLS als uitgangspunt genomen. Ook voor wegvervoer en spoorvervoer is er voor gekozen om TRANS-TOOLS voor deze relaties als uitgangspunt te nemen voor de afstanden en tijden, zodat alle grensoverschrijdende afstanden en tijden uit dezelfde bron komen (consistentie).

2.2.4

TRANSTOOLS

De modellen RoutGoed, BIVAS en LMS geven voor het buitenland niet het gewenste detailniveau. Hierom is gekozen om voor alle modaliteiten de tijden en afstanden voor het buitenland af te leiden uit TRANS-TOOLS. De TRANS-TOOLS zonering voor het buitenland komt niet overeen met de gewenste BasGoed indeling. Hierom is de zonering in de TRANS-TOOLS netwerken aangepast. Voor de binnenlandse zones zijn de coördinaten van de zones afgeleid uit de LMS subzones. De zone-coördinaten voor de buitenlandse zones zijn bepaald door expert opinion. Voor elke SMILE-zone is het zwaartepunt zo gekozen dat deze dicht bij de steden ligt die het meeste vervoer genereren. De gekozen coördinaten zijn voor alle modaliteiten identiek. Op basis van expert opinion zijn de zones vervolgens aan de dienstregelingspunten/voedingspunten gekoppeld van het wegen-, spoor- en binnenvaartnetwerk. Aan de hand van de aangepaste TRANS-TOOLS netwerken is per modaliteit de level of service bepaald voor de internationale relaties. Deze level of service bestaat alleen uit reistijden op basis van de kortste route. Wachttijden, laad- en lostijden, voedingskosten en grensweerstanden zijn niet meegenomen in het bepalen van de level of service. Ook is er geen rekening gehouden met rusttijden van chauffeurs. Voor het binnenvaartnetwerk is de level of service bepaald op CEMT-I klasse niveau. Er is hierbij dus geen rekening gehouden met vaarwegen die voor sommige scheepstypen niet toegankelijk zijn. Daarnaast is een aantal zones in buitenland niet aangesloten op het binnenvaartnetwerk, voor deze zones is het niet mogelijk een level-of-service te bepalen. Dit omdat er voor deze zones geen mogelijkheden zijn om via de binnenvaart goederen te vervoeren. Deze zones zijn Griekenland, Italië, Scandinavië, het Iberisch Schiereiland en voormalig USSR. De genoemde zones zijn wel ontsloten door middel van het spoor- en wegennetwerk.

2.3

Vergelijkingskader modaliteiten

In BASGOED wordt in de kostenfuncties gebruik gemaakt van de kostenkentallen uit het vergelijkingskader modaliteiten (VKM). Het VKM biedt gedetailleerde kentallen goederenvervoer naar vervoerwijze, voertuigtype en verschijningsvorm. Om kentallen af te leiden naar de drie geaggregeerde vervoerwijzen in BASGOED, zijn de VKM kentallen gekoppeld aan het basisbestanden goederenvervoer. Hiervoor is een match gemaakt tussen de verschijningsvorm en voertuigtype combinaties in het VKM en het basisbestand goederenvervoer. Vervolgens is per kental een gewogen gemiddelde berekend voor de drie hoofdvervoerwijzen (gewogen naar vervoerde tonnen). Het resultaat van deze bewerking staat in Tabel 2-2. Deze kentallen zijn per voertuig. Uit het VKM wordt een aantal kentallen gebruikt. De gemiddelde vaste kosten worden gebruikt voor de berekening van tijdsafhankelijke kosten en bevatten de kosten zoals chauffeurskosten. De afstandafhankelijke kosten worden berekend uit de gemiddelde variabele kosten (per km) en energiekosten

12


Schattingen BASGOED

(per km). De variabele kosten per kilometer bevatten de kosten voor onderhoud en afschrijving aan het voertuig. De energie kosten bevat brandstofkosten. Verder zijn laad- en loskosten gebruikt. In BasGoed worden goederenstromen in tonnen per NSTR1 groep gemodelleerd. De VKM kentallen ver voertuig zijn daarom via het gemiddelde laadvermogen, de beladingsgraad en het aantal beladen ritten in totaal aantal ritten, omgerekend naar kentallen per ton. Het resultaat hiervan staat in Tabel 2-3.

Tabel 2-2: Kosten kentallen per voertuig per vervoerwijze. Bron: VKM/BBG04

Weg 33.79 0.11 0.19 37.32 37.34 23.18 0.69 0.71

Gemiddelde vaste kosten (euro per vtg per uur) Gemiddelde variabele kosten (euro per vtg per km) Gemiddelde energiekosten (euro per vtg per km) Gemiddelde laadkosten (euro per vtg) Gemiddelde loskosten (euro per vtg) Gemiddeld laadvermogen (ton) Beladingsgraad inhoud Beladen ritten in totaal aantal ritten

Spoor 286.29 4.45 3.48 3142.33 3142.33 1368.92 0.76 0.67

Binnenvaart 141.83 0.70 5.98 1348.03 1795.92 2823.69 0.79 0.67

Tabel 2-3: Kosten kentallen per ton. Bron: VKM/BBG04

Weg 2.97 0.00955 0.0165 0.0260 3.29 3.29

Gemiddelde vaste kosten (euro per uur per ton) Gemiddelde variabele kosten (euro per tonkm) Gemiddelde energiekosten (euro per tonkm) Totale afstandsafh. kosten: Gemiddelde laadkosten (euro per ton) Gemiddelde loskosten (euro per ton)

2.4

Spoor 0.41 0.00643 0.0050 0.0115 4.54 4.54

Binnenvaart 0.10 0.00047 0.0040 0.0045 0.91 1.21

Voedingskosten binnenvaart en spoor

Voor rail en binnenvaart is het wenselijk een component op te nemen voor de voedingskosten, waarin de bereikbaarheid over de weg van rail- en binnenvaartterminals in de betreffende regio is verwerkt. Voedingskosten zijn de reistijd en -afstand die het kost om de voedingspunten van het rail- en binnenvaart goederenvervoer te bereiken. Deze zijn afgeleid uit de gecongesteerde reistijd en -afstand over de weg (afkomstig uit LMS) en de geografische locaties van de dienstregelpunten voor het spoor en de voedingspunten voor de binnenvaart.

13


Schattingen BASGOED

2.4.1

Spoor

Routgoed bevat voedingspunten voor het rail goederenvervoer, dienstregelingspunten, welke zijn verspreid door Nederland, te zien in Figuur 2-2. Als eerste stap zijn alle LMS subzone toegekend aan meest nabije Routgoed dienstregelpunt. De voedingskosten, uitgedrukt in afstand en tijd, is berekent als gewogen gemiddelde gecongesteerde reistijden en –afstanden die uit het LMS komen. Na aggregatie naar etmaal is er opnieuw een gewogen gemiddelde bepaald, door te aggregeren vanuit LMS subzones naar corop gebieden. Figuur 2-2: LMS zonering en Routgoed dienstregelingspunten

14


Schattingen BASGOED

2.4.2

Binnenvaart

De gebiedsindeling van Bivas is getoond in Figuur 2-3. Hierin is te zien dat de BasGoed zonering (corop gebieden) wezenlijk verschilt van de bivaszonering, is er geen eenduidig voedingspunt vast te stellen wat representatief is voor een Corop. In plaats daarvan is per BasGoed zone de intrazonale reistijd en –afstand binnen de BasGoed zone als voedingskosten voor binnenvaart gehanteerd. Deze interzonale voedingskosten zijn berekend door de LMS subzonale reistijden en –afstanden te aggregeren naar corop niveau, met het aantal tours als wegingsfactor (als ook in 2.2.1).

Figuur 2-3: Bivas en BasGoed zonering

15


Schattingen BASGOED

16

HOOFDSTUK 3

Schattingen modal split

Dit hoofdstuk beschrijft de aanpak bij het schatten van de modal split modellen. Allereerst wordt de nutsfunctie van de geschatte modellen beschreven. Vervolgens is een overzicht gegeven van de geteste modelspecificaties. Vervolgens zijn in detail voor de belangrijkste test specificaties de schattingsresultaten per goederengroep gepresenteerd. Daarna zijn de belangrijkste modellen beoordeeld op basis van de modelprestaties (reistijdwaardering, elasticiteiten). Tenslotte zijn de belangrijkste conclusies ten aanzien van de geschatte modal split modules gegeven. 3.1

Kostenfuncties Modal split

De kostenfuncties voor BasGoed zijn voor een belangrijk deel gebaseerd worden op die in het huidige SMILE+, en bevatten naast verbindingskosten ook de voedingskosten. De verbindingskosten (linkgebonden kosten) in BasGoed zijn berekend uit de tijd en kosten voor de beschikbare vervoerwijzen op een relatie. BasGoed bevat niet de volledige kosten voor overslag omdat het model gebaseerd is op uitvoer van uni-modale toedelingsmodellen en logistieke ketens ontbreken in het model. De voedingskosten (tussen de zonezwaartepunten en de netwerkknooppunten) voor de vervoerwijzen rail en binnenvaart zijn berekend uit de gemiddelde tijden afstand die moet worden afgelegd om het voedingspunt (terminal) van het netwerk te bereiken. Deze tijd en afstand zijn berekend op basis van de vrachtauto verplaatsingen uit het LMS. Deze voedingskosten zijn afhankelijk van de goederengroep in rekening gebracht voor spoor en binnenvaart. In Figuur 3-1 zijn de kosten in BASGOED voor vervoerwijze spoor (zelfde als voor binnenvaart) geschematiseerd weergegeven. Voor wegtransport kan worden aangenomen dat geen additionele voedingskosten nodig zijn en dat de gewogen reistijden en reiskosten uit het LMS representatief zijn voor de logistieke keten voor deze vervoerwijze. Voor een aantal goederengroepen is het aannemelijk dat het goederenvervoer per spoor of binnenvaart via direct access plaatsvindt. In deze gevallen worden geen voedingskosten over de weg meegenomen. Het gaat om de goederengroepen 2, 3 en 4. Het goederenvervoer binnen NSTR groep 6 betreft vooral grind en andere bulk bouwmaterialen die aan herkomst zijde direct via binnenvaart zijn ontsloten. Hier zijn alleen voedingskosten aan bestemmingszijde meegenomen.


Schattingen BASGOED

VOEDINGSKOSTEN: - tijd kosten - afstand kosten

VERBINDINGSKOSTEN: - tijd kosten - afstand kosten

H zone: - laadkosten weg

VOEDINGSKOSTEN: - tijd kosten - afstand kosten

B zone: - loskosten weg

TERMINAL: - loskosten weg - laadkosten spoor

TERMINAL: - loskosten spoor - laadkosten weg

Figuur 3-1: Componenten in BASGOED kostenfuncties Spoor en Binnenvaart. (vet gedrukt zijn wel opgenomen; in rood gedrukt welke niet)

De kostenfuncties bevatten naast de transporttijd en transportkosten per vervoerwijze, ook vervoerwijzespecifieke constanten en mogelijk enkele goederensoortspecifieke variabelen (zoals containerisatie graad). In het specificatie onderzoek voor BASGOED zijn 4 verschillende kostenfuncties onderzocht, met verschillende aannames ten aanzien van de tijden kosten parameters. Het betreft de volgende kostenfuncties2. 1. Het schatten van een reistijdwaarderingsparameter βt met constante kostenparameter ( β c = 0.1)3. De nutsfunctie voor vervoerwijze v op relatie ij voor goederengroep g is als volgt gespecificeerd:

Vijvg = β c ⋅ (d ijv ⋅ K v ) + β t ⋅ t ijv + MSC v + CONTijg , ∀β c = −0.1

(1)

waarin d ijv en tijv de afstand en tijd op de voedings- en verbindingslinks, Kv het tarief 4 per km, CONTijg de containerisatiegraad (verhoogde kans op spoor of binnenvaart als de goederen per

2

In de functies zijn ook de initiële laad- en loskosten opgenomen, die verschillen per vervoerwijze. Voor de eenvoud van de presentatie zijn deze kosten componenten hier verder niet genoemd.

3 Kosten waren initieel vastgezet op -1, maar bij de schattingen bleek de absolute utility in een dusdanig bereik te komen dat de schattingssoftware in veel gevallen de overige coëfficiënten niet kon schatten. 4 Het tarief per kilometer in (1) is representatief voor de totale transportkosten die door de vervoerder aan de verlader in rekening worden gebracht (exclusief winst marge). Hierin zijn zowel de tijden afstandafhankelijke transportkosten opgenomen. Als representatieve tarieven zijn per vervoerwijze de volgende tarieven gehanteerd: weg: 7 ct/ton/km, spoor 3 ct/ton/km, binnenvaart 2 ct/ton/km (Gebaseerd op: NEA (2003) Factorkosten van het goederenvervoer: een analyse van de ontwikkeling in de tijd. Rapport voor de Adviesdienst Verkeer en Vervoer, Rijswijk: NEA)

18


Schattingen BASGOED

container worden vervoerd) en MSCv de vervoerwijze specifieke constante. Deze laatste parameter neemt de niet-waargenomen kwaliteitsverschillen tussen vervoerwijzen mee. Deze kostenfunctie is uiteindelijk niet gekozen (zie verderop in dit hoofdstuk). 2. Het onderscheiden van tijdsafhankelijke en kilometerafhankelijke transportkosten, waarbij tijd middels een tarief per tijdseenheid, Rv, in transportkosten wordt omgezet, en één parameter wordt geschat voor de totale transport kosten:

Vijvg = β tc ⋅ (d ijv ⋅ Tv + t ijv Rv ) + MSC v + CONTijg

(2a)

De afstandskosten per vervoerwijze Tv , en het tarief per tijdseenheid per vervoerwijze, Rv , zijn overgenomen uit het VKM (zie Tabel 2-3). Deze kostenfunctie is uiteindelijk niet gekozen (zie verderop in dit hoofdstuk). Daarnaast is een aanpak getoetst waarin een aparte coefficient is geschat voor de afstand-, en tijdsafhankelijke kosten. De afstand en tijd, zijn via de gemiddelde afstandskosten Tv en het tarief per tijdseenheid Rv omgerekend naar kosten:

Vijvg = β c ⋅ (d ijv ⋅ Tv ) + β t ⋅ (tijv Rv ) + MSCv + CONTijg

(2b)

waarin Tv en Rv zijn overgenomen uit het VKM. Bij deze specificatie wordt een hoge correlatie van tijd en kosten verwacht en dus problemen bij het schatten (deze kostenfunctie is uiteindelijk niet gekozen (zie verderop in dit hoofdstuk). Tenslotte is daarom de volgende variant getoetst waarbij de tijdsafhankelijke- en kilometerafhankelijke transportkosten gebaseerd zijn op de VKM-kostenkentallen, maar ook een additionele tijdparameter is geschat voor de kapitaalkosten van de goederen. In deze kapitaalkosten zijn de kosten voor interestkosten, ontwaarding en verzekeringskosten opgenomen:

Vijvg = β tc ⋅ (d ijv ⋅ Tv + tijv Rv ) + β t ⋅ tijv + MSCv + CONTijg

(2c)

Deze kostenfunctie is uiteindelijk gekozen (zie verderop in dit hoofdstuk). 3. Het schatten van één parameter voor gegeneraliseerde kosten en tijd, waarin tijd in kosten wordt uitgedrukt via een exogene reistijdwaardering:

Vijvg = β gc ⋅ (d ijv ⋅ Tvg + tijvVoTvg ) + MSCv + CONTijg

(3)

Deze vergelijking is in de uitwerking van het onderzoek niet meer aan bod gekomen omdat de specificaties die gebaseerd zijn op kostenkentallen al goede resultaten opleverden. Kostenkentallen

19


Schattingen BASGOED

zijn met meer betrouwbaarheid bekend en daarom zijn specificaties gebaseerd op kostenkentallen te prefereren. Deze kostenfunctie is uiteindelijk niet gekozen (zie verderop in dit hoofdstuk). 4. Het schatten van een aparte parameter voor kosten β c en tijd βt , zoals in vergelijking 1 waarbij

βc niet is vastgezet: Vijvg = β c ⋅ (d ijv ⋅ K v ) + β t ⋅ tijv + MSCv + CONTijg

(4)

waarin Kv het tarief per tonkm is overgenomen uit de rapportage factorkosten (NEA, 2003). Deze kostenfunctie is uiteindelijk niet gekozen (zie verderop in dit hoofdstuk).

3.2

Overzicht schattingen

In onderstaande tabel staan de uitgangspunten van de geschatte model series schematisch weergegeven. Op basis van eerdere series modelschattingen (niet allemaal in het overzicht) bleek een aanpak waarbij alle waarnemingen worden meegenomen, en de beschikbaarheid voor spoor en binnenvaart is afgeleid uit het aandeel van deze vervoerwijzen in de waarnemingen in het schattingsbestand, het meest ideaal. Deze 5 series zijn uitgevoerd op kostenfunctie 1, waarin alleen een tijdswaardering coëfficiënt wordt geschat. De kosten parameter is vastgezet op -0.1. In de eerste twee series schattingen bleek het niet mogelijk een model te schatten waarbij voor alle waarnemingen, alle vervoerwijzen beschikbaar waren. Verder is in model serie 5 getoetst of het negeren van dominante waarnemingen in de dataset tot betere modelresultaten leidde. Verder bleek ook het opnemen van een exogene parameter voor terminal kosten (model serie 2 en 3) niet te leiden tot een beter model. Deze niet geobserveerde kosten zijn vervoerwijze specifiek en worden opgepakt door de vervoerwijze specifieke constante (MSC). In model serie 6 is kosten functie 2a geschat, waarin een kostenparameter is geschat voor de gesommeerde tijdsafhankelijke en kilometerafhankelijke transportkosten. Deze kosten parameter bleek voor alle NSTR1 groepen goed te schatten (plausibel en significant teken). Deze kostenfunctie is gebaseerd op het VKM en sluit daardoor aan bij de methodiek in SMILE+ en is een goed basismodel voor de verdere analyses. Verder is de specificatie uitgebreid met een vervoerwijze specifieke dummy voor internationaal vervoer, om rekening te houden met verschillen in grensweerstanden tussen de vervoerwijzen. Deze bleek de modelresultaten nog iets te verbeteren. Model serie 8 is gebaseerd op kostenfunctie 4 gewerkt, met een aparte parameter voor kosten en tijd. Deze parameters waren beide significant en met het juiste teken. Daarmee is het ook een mogelijke kandidaat voor implementatie. Model serie 10 is gebaseerd op kostenfunctie 2b gewerkt, met een aparte kostenparameter voor afstandkosten en tijdparameter voor tijdsafhankelijke kosten, berekend uit het VKM. De coëfficiënten in deze modellen bleken niet plausibel, met positief teken. Ook het schatten van vervoerwijze specifieke tijd parameters, in model serie 11, bleek geen plausibele tijd en kosten parameters op te leveren. Model serie 12 is gebaseerd op het VKM en kostenfunctie 2c, met een kostenparameter voor de afstanden tijdsafhankelijke transportkosten berekend uit VKM kentallen, net als model serie 6 waarin de kostenfuncties consistent zijn met SMILE+. Daarnaast is het model uitgebreid met een aparte tijdparameter die de kapitaal kosten van de goederen schat. Deze parameters waren beide significant en met het juiste teken.

20


Schattingen BASGOED

In model serie 13 geschat waarin model serie 12 is uitgebreid met regiospecifieke dummies per vervoerwijze. Deze modellen bleken moeilijk te schatten. Veel regio dummies moesten uit het model verwijderd worden omdat de schattingen niet convergeerden. Gedetailleerde schattingsresultaten zijn opgenomen in de appendices. Tenslotte is in model serie 14 een hybride aanpak gevolgd waarbij een beperkt aantal parameters aan model 12 zijn toegevoegd, voor relaties (corridors) waarop de gesimuleerde marktaandelen veel afwijken van de waargenomen tonnages. Deze selectie heeft plaatsgevonden op basis van de simulatieresultaten van model serie 12. Per NSTR groep zijn de 10 relaties bepaald met de grootste absolute afwijking voor binnenvaart. Om het aantal parameters in het model beperkt te houden is het aantal corridors met dummies bewust laag gehouden (10).

Tabel 3-1: Overzicht schattingen modal split module

Serie

Segmentatie

Specificatie

4

NSTR1

Kostenfunctie 1: tijd parameter β t , kosten=-0.1 (MSCs, CONTg)

6

NSTR1

Kostenfunctie 2a: kosten parameter β tc , (MSCs, CONTg) en internationale dummy’s

8

NSTR1

Kostenfunctie 4: tijd en kosten parameters

10

NSTR1

NSTR1

NSTR1

βc

NSTR1

NSTR1

β tv

voor tijd per vervoerwijze, MSCs, CONTg en internationale dummy’s Kostenfunctie 2c:

β tc

en β c ,

MSCs, CONTg en internationale dummy’s Kostenfunctie 2c: met regionale MSC’s voor voedingskosten, tijd en kosten parameters

14

en β c ,

voor kosten en

tijd en kosten parameters

13

βt

MSCs, CONTg en internationale dummy’s Kostenfunctie 2b: parameters

12

en β c ,

MSCs, CONTg en internationale dummy’s Kostenfunctie 2b: tijd en kosten parameters

11

βt

β tc

β tc

alle waarnemingen, bvrt en spoor alleen beschikbaar als ook gekozen

en β c , MSCs,

en β c , MSCs,

CONTg en internationale dummy’s

21

Resultaat Plausibele

β t , MSC’s en soms

CONT

Plausibele

β tc , MSC’s en soms

CONT

Plausibel teken

βt

en

βc , βt

erg klein, MSC’s, internationale dummy’s en soms CONT significant vaak implausibele tekens bij tijd- of kostenparameters


vaak implausibele tekens bij tijd- of kostenparameters


Plausibel teken


CONTg en internationale dummy’s Kostenfunctie 2c: met corridor dummy’s voor onder/over schatte vervoerwijzes, tijd en kosten parameters,

Selectie data en keuzeset alle waarnemingen, bvrt en spoor alleen beschikbaar als ook gekozen alle waarnemingen, bvrt en spoor alleen beschikbaar als ook gekozen alle waarnemingen, bvrt en spoor alleen beschikbaar als ook gekozen


βt

en

βc , βt

erg klein, MSC’s, internationale dummy’s en soms CONT significant Veel dummy’s blijken niet te schatten. Veel andere globale parameters (internationale dummy’s en containerisatue) vallen uit. Final LL beter dan serie 12 Meeste corridor dummy’s significant. Weinig uitval van andere globale parameters. Final LL beter dan serie 12, maar minder goed dan serie 13


3.3

Schattingen BASGOED

Gedetailleerde schattingsresultaten

In deze paragraaf zijn per NSTR groep de schattingsresultaten opgenomen. Hierin zijn onder andere opgenomen: de geschatte coëfficiënt waarden, de significantie van de coëfficiënten (tussen haakjes), en beschrijvende statistieken van het model, zoals de Final log likelihood van het geconvergeerde geschatte model, en een tweetal standaard fit-tests van het gehele model, berekend op basis van de final likelihood en base likelihood. De schattingen worden niet in detail besproken, maar er worden een aantal opmerkingen ter toelichting gemaakt. Voor serie 6 kon voor het vervoermiddel spoor geen buitenlanddummy geschat worden voor NSTR 0, 2, 4 en 6. Voor NSTR2 konden overigens de buitenlanddummy voor binnenvaart en de containerisatiegraden ook niet geschat worden. Voor serie 8 was het voor NSTR 0, 2, 4 en 6 niet mogelijk om een buitenlanddummy voor spoorvervoer mee te nemen. Voor NSTR6 en NSTR7 konden geen coëfficiënt geschat worden voor de containerisatiegraad voor spoor en weg. Bij NSTR 3, 8 en 9 kon geen coëfficiënt geschat worden voor de containeriesatiegraad voor binnenvaart. In serie 12 vielen voor een aantal NSTR groepen (0, 2, 4 en 6) de buitenlanddummies voor spoor uit het model. Voor NSTR2 vielen de buitenlanddummies voor binnenvaart en de containerisatiegraden uit. Voor NSTR3 kon geen tijdcoëfficiënt (kapitaalkosten) geschat worden. In model serie 13 zijn voor lang niet alle regio’s vervoerwijzedummies te schatten. Veel andere globale parameters (internationale dummies en containerisatue) vallen uit. De meeste corridor dummies die zijn toegevoegd aan model serie 14 waren significant. Verder vallen weinig andere globale parameters uit.

22

Schattingen BASGOED


Tabel 3-2: Schattingsresultaten voor model serie 6 Serie 6

NSTR0 Landbouwp.

NSTR1 Voedingsp.

NSTR2 Vaste min. Br. Stof

NSTR3 Aardolie

NSTR4 Ertsen

NSTR5 Metalen, halffabr.

NSTR6 Ruwe min.

NSTR7 Meststoff.

NSTR8 Chemie

NSTR9 Overig

File model06bnstr0_bs_cont.F12 model06bnstr1_bs_cont.F12 model06bnstr2_bs_cont_nofeed2.F12 model06bnstr3_bs_cont_nofeed.F12 model06bnstr4_bs_cont_nofeed.F12 model06bnstr5_bs_cont.F12 model06bnstr6_bs_cont_nofeed_destbin.F12 model06bnstr7_bs_cont.F12 model06bnstr8_bs_cont.F12 model06bnstr9_bs_cont.F12 Title BASGOED MS model BASGOED MS model BASGOED MS model BASGOED MS model BASGOED MS model BASGOED MS model BASGOED MS model BASGOED MS model BASGOED MS model BASGOED MS model Converged True True True True True True True True True True Observations 759 784 252 620 454 481 1390 450 537 1054 Final log (L) -18441.8 -34221.1 -13618.8 -15834.2 -9568.8 -15967.4 -114828.6 -4653.1 -38534.7 -69912.9 D.O.F. 6 7 3 7 6 7 6 7 7 7 Rho²(0) 0.297 0.283 0.533 0.623 0.729 0.083 0.171 0.225 0.263 0.444 Rho²(c) 0.092 0.136 0.028 0.359 0.464 0.08 0.148 0.207 0.156 0.32 Estimated 5-Jul-10 5-Jul-10 5-Jul-10 5-Jul-10 5-Jul-10 5-Jul-10 5-Jul-10 5-Jul-10 5-Jul-10 5-Jul-10 Scaling 1 1 1 1 1 1 1 1 1 1 Grns_Bvrt 1.127 (34.2) 0.5647 (21.5) 0 (*) 1.979 (40.5) 2.175 (40.9) 0.7580 (17.3) -0.3050 (-20.5) 0.1406 (2.5) 1.379 (59.4) 1.263 (68.0) Grns_Spr 0 (*) -2.137 (-19.4) 0 (*) 0.01536 (0.1) 0 (*) 2.391 (14.6) 0 (*) 2.368 (2.8) -0.02010 (-0.4) -0.02667 (-0.9) CONT_Bvrt -3.744 (-15.0) -0.9761 (-12.4) 0 (*) 3.751 (9.3) -6.224 (-10.0) -3.802 (-17.9) 1.999 (8.5) 3.408 (3.8) 0.7516 (14.7) 4.709 (165.6) CONT_Spoor -96.53 (-11.8) -10.61 (-5.3) 0 (*) -2.581 (-2.4) -2072 (-21.9) -31.24 (-20.4) -210.0 (-17.4) 26.43 (3.4) -6.208 (-29.0) 5.783 (90.3) CONT_Weg 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) ASC_Bvrt -1.913 (-79.7) -2.003 (-125.6 2.136 (24.7) -3.237 (-53.7) -2.177 (-29.8) -0.8238 (-20.4) -2.144 (-177.7 -1.614 (-31.2) -1.696 (-95.7) -3.429 (-203.2 ASC_Spoor -0.3255 (-3.4) -1.326 (-21.9) 1.540 (26.4) -0.3065 (-2.6) 5.086 (40.3) -2.307 (-14.4) -1.132 (-29.5) -5.114 (-6.0) -1.388 (-31.8) -4.132 (-95.0) ASC_Weg 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) Cost -0.01580 (-11.3) -0.08196 (-53.1) -0.08001 (-14.7) -0.5049 (-66.3) -0.2528 (-37.2) -0.00883 (-6.0) -0.2326 (-134.0 -0.1516 (-32.8) -0.04144 (-40.1) -0.00791 (-8.9)

Tabel 3-3: Schattingsresultaten voor model serie 8 Serie 8

NSTR0 Landbouwp.

NSTR1 Voedingsp.

NSTR2 Vaste min. Br. Stof

NSTR3 Aardolie

NSTR4 Ertsen

NSTR5 Metalen, halffabr.

NSTR6 Ruwe min.

NSTR7 Meststoff.

NSTR8 Chemie

NSTR9 Overig

File model08bnstr0_bs_cont.F12 model08bnstr1_bs_cont.F12 model08bnstr2_bs_cont_nofeed.F12 model08bnstr3_bs_cont_nofeed.F12 model08bnstr4_bs_cont_nofeed.F12 model08bnstr5_bs_cont.F12 model08bnstr6_bs_cont_nofeed_origbvrt.f12 model08bnstr7_bs_cont.F12 model08bnstr8_bs_cont.F12 model08bnstr9_bs_cont.F12 Title BASGOED MS model BASGOED MS model BASGOED MS model BASGOED MS model BASGOED MS model BASGOED MS model BASGOED MS model BASGOED MS model BASGOED MS model BASGOED MS model Converged True True True True True True True True True True Observations 759 784 252 620 454 481 1390 450 537 1054 Final log (L) -18180.6 -31537.9 -13382.5 -15932.2 -9087.6 -15625 -115360.8 -4486.2 -38116.6 -67506.3 D.O.F. 7 8 7 7 7 8 6 7 8 8 Rho²(0) 0.307 0.339 0.541 0.62 0.743 0.102 0.167 0.253 0.271 0.464 Rho²(c) 0.105 0.204 0.045 0.355 0.491 0.1 0.144 0.235 0.165 0.343 Estimated 5-Jul-10 5-Jul-10 5-Jul-10 5-Jul-10 5-Jul-10 5-Jul-10 5-Jul-10 5-Jul-10 5-Jul-10 5-Jul-10 Scaling 1 1 1 1 1 1 1 1 1 1 Grns_Bvrt 1.263 (39.0) 0.4597 (17.2) -0.2042 (-1.8) 2.237 (43.9) 2.744 (48.6) 0.9554 (21.6) -0.08961 (-6.0) 0.5460 (9.0) 1.576 (65.7) 1.567 (83.7) Grns_Spr 0 (*) -10.53 (-56.1) 0 (*) -1.147 (-6.6) 0 (*) 1.874 (11.3) 0 (*) 0.9643 (1.1) -0.5852 (-10.3) -0.6218 (-18.7) CONT_Bvrt -4.355 (-17.2) -1.606 (-19.6) -10.57 (-14.9) 0 (*) -6.873 (-11.0) -3.965 (-18.3) 2.549 (10.6) 3.480 (3.7) 0 (*) 0 (*) CONT_Spoor -78.05 (-9.6) -7.342 (-5.7) -146.1 (-11.6) -6.455 (-6.8) -3125 (-13.5) -28.85 (-19.2) 0 (*) 0 (*) -6.575 (-28.9) 0.2623 (4.6) CONT_Weg 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -0.4574 (-8.7) -4.622 (-158.1 ASC_Bvrt -2.141 (-79.9) -2.654 (-130.5 2.448 (20.5) -2.779 (-50.7) -2.243 (-32.7) -0.8719 (-21.3) -2.109 (-183.0 -2.088 (-34.4) -1.725 (-96.2) -3.442 (-203.2 ASC_Spoor -2.249 (-15.8) -0.6708 (-10.6) 1.670 (28.0) 1.044 (8.7) 2.203 (16.3) -2.304 (-14.4) -2.002 (-48.1) -4.838 (-5.7) -1.313 (-29.6) -3.543 (-91.0) ASC_Weg 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) Time -3.71e-4 (-22.9) -0.00105 (-65.2) -1.84e-4 (-16.4) -4.37e-4 (-16.0) -0.00108 (-28.0) -2.47e-4 (-24.9) -2.77e-4 (-32.4) -4.60e-4 (-14.7) -2.35e-4 (-29.9) -3.85e-4 (-65.6) Cost -0.09048 (-24.3) -0.3011 (-78.4) -0.1048 (-16.1) -0.5202 (-65.9) -0.3498 (-43.4) -0.03934 (-19.6) -0.2560 (-131.5 -0.2494 (-31.6) -0.07170 (-47.7) -0.04818 (-41.2)


Schattingen BASGOED

Tabel 3-4: Schattingsresultaten voor model serie 12 NSTR 0 1 2 3 4 5 6 7 8 9 Title Landbouwp. Voedingsp. Vaste min. Br. Stof. aardolie Ertsen Metalen. Halffabr. Ruwe min. Meststof Chemie Overig Converged TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE Observations 763 787 252 633 454 483 1394 450 580 1120 Final log (L) -18296.1 -31633.9 -13593.3 -16577.7 -8964 -15681.5 -116109.4 -4503.6 -39589 -69137.3 D.O.F. 7 8 4 7 7 8 7 8 8 8 Rho²(0) 0.352 0.36 0.534 0.63 0.746 0.119 0.195 0.25 0.298 0.52 Rho²(c) 0.13 0.215 0.029 0.369 0.498 0.11 0.164 0.232 0.18 0.369 Estimated 23-Jul-10 23-Jul-10 23-Jul-10 23-Jul-10 23-Jul-10 23-Jul-10 23-Jul-10 23-Jul-10 23-Jul-10 23-Jul-10 Scaling 1 1 1 1 1 1 1 1 1 1 Grns_Bvrt 1.219 (36.6) 0.2322 (8.3) 0 (*) 2.029 (41.2) 2.686 (47.7) 0.9701 (21.8) -0.1341 (-8.8) 0.4657 (7.7) 1.638 (68.6) 1.602 (85.0) Grns_Spr 0 (*) -13.69 (-62.0) 0 (*) -0.9145 (-6.7) 0 (*) 1.954 (11.8) 0 (*) 0.8316 (0.9) -0.6218 (-10.8) -0.5589 (-16.9) CONT_Bvrt -4.003 (-15.8) -1.483 (-18.4) 0 (*) 4.659 (11.4) -6.930 (-10.9) -4.048 (-18.7) 1.333 (5.7) 3.112 (3.4) 0.6337 (12.3) 4.739 (164.3) CONT_Spoor -65.10 (-8.1) -8.400 (-5.3) 0 (*) -0.3994 (-0.3) -5061 (-20.4) -29.65 (-19.7) -304.6 (-18.9) -75.57 (-7.8) -6.463 (-27.5) 5.371 (97.1) CONT_Weg 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) ASC_Bvrt -2.329 (-87.6) -2.758 (-134.9 2.323 (26.2) -3.461 (-57.3) -2.622 (-36.5) -0.8337 (-20.6) -2.411 (-194.9 -2.113 (-34.0) -1.885 (-110.2 -3.578 (-216.5 ASC_Spoor -3.472 (-20.3) -0.6692 (-10.4) 1.523 (25.9) -0.6206 (-5.4) 2.007 (14.2) -2.345 (-14.7) -2.218 (-46.8) -5.353 (-6.0) -1.465 (-32.8) -4.080 (-108.0 ASC_Weg 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) Time -4.34e-4 (-27.0) -0.00111 (-66.8) -5.29e-5 (-6.7) 0 (*) -0.00128 (-32.1) -2.16e-4 (-23.0) -4.19e-4 (-44.4) -6.13e-4 (-16.8) -2.21e-4 (-29.8) -3.64e-4 (-65.2) Cost -0.1121 (-29.0) -0.3211 (-80.1) -0.07404 (-14.3) -0.5246 (-68.3) -0.3884 (-46.7) -0.03061 (-17.4) -0.2892 (-139.9 -0.2722 (-30.7) -0.06721 (-49.6) -0.04476 (-43.6)

De schattingsresultaten voor model serie 13 zijn opgenomen als aparte bijlage.

24


Schattingen BASGOED

Tabel 3-5: Schattingsresultaten voor model serie 14

File Title Converged Observations Final log (L) D.O.F. Rho²(0) Rho²(c) Estimated Scaling Grns_Bvrt Grns_Spr CONT_Bvrt CONT_Spoor CONT_Weg CR_bv10 CR_bv9 CR_bv8 CR_bv7 CR_bv6 CR_bv5 CR_bv4 CR_bv3 CR_bv2 CR_bv1 ASC_bvrt ASC_Spr ASC_weg Time Cost

0 1 2 3 4 5 6 7 8 9 BASGOED MS model BASGOED MS model BASGOED MS model BASGOED MS model BASGOED MS model BASGOED MS model BASGOED MS model BASGOED MS model BASGOED MS model BASGOED MS model TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE 763 787 252 633 454 483 1394 450 580 1120 -15430.4 -29133.6 -12114 -14794.8 -7899.6 -15007.3 -104295.7 -3744.2 -37249.6 -69135.2 16 16 15 15 15 14 16 16 16 10 0.453 0.411 0.584 0.67 0.777 0.157 0.277 0.376 0.34 0.52 0.267 0.277 0.135 0.437 0.558 0.148 0.249 0.362 0.228 0.369 23-Jul-10 23-Jul-10 23-Jul-10 23-Jul-10 02-Aug-10 23-Jul-10 23-Jul-10 23-Jul-10 23-Jul-10 23-Jul-10 1 1 1 1 1 1 1 1 1 1 1.210 (31.4) 0 (*) 0.4511 (3.7) 2.233 (43.6) 1.327 (20.3) 0 (*) -0.1174 (-7.0) 0 (*) 1.234 (48.5) 1.602 (85.0) 0 (*) -13.78 (-66.8) 0 (*) 0 (*) 0 (*) 0.8602 (5.0) 0 (*) -3.053 (-2.9) -0.5761 (-10.1) -0.5586 (-16.9) -2.958 (-11.1) -1.262 (-15.7) -5.518 (-5.9) 2.360 (6.2) -6.978 (-12.0) 0 (*) 0.6209 (2.7) 0 (*) 0.9620 (17.7) 4.740 (164.3) -71.84 (-9.0) -6.651 (-6.2) -433.9 (-14.1) -3.626 (-3.7) 0 (*) -28.58 (-18.9) -337.5 (-19.1) -103.0 (-10.2) -6.095 (-27.4) 5.371 (97.1) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 3.935 (8.5) 2.217 (16.6) -8.050 (-5.8) -9.216 (-1.3) -1.968 (-6.6) 3.794 (8.4) 0 (*) -5.654 (-3.1) -5.316 (-7.5) 0 (*) 0 (*) 0 (*) -5.150 (-10.1) 1.938 (10.4) -5.915 (-6.7) 3.449 (9.9) 4.124 (18.6) 0.3243 (1.5) 2.767 (22.3) -2.031 (-0.3) 3.418 (23.3) -3.871 (-9.6) -3.423 (-6.0) -5.673 (-4.8) -5.282 (-5.6) 1.508 (13.2) -5.668 (-11.9) -5.425 (-3.2) -1.400 (-15.8) 0 (*) -6.477 (-2.9) 3.032 (25.6) -2.863 (-4.6) -5.897 (-4.6) -3.273 (-5.8) -6.284 (-5.4) -8.171 (-5.0) -3.561 (-6.1) 3.438 (15.0) 0 (*) 6.745 (9.5) 1.732 (13.2) 0 (*) 0 (*) -3.293 (-7.4) 0 (*) -3.088 (-25.1) 2.810 (7.5) 0 (*) 0 (*) 3.576 (12.2) 1.752 (18.9) 0.3502 (2.0) 2.274 (11.9) -5.287 (-5.8) 1.300 (12.5) -4.249 (-19.6) 2.582 (7.5) 0.7305 (11.7) 0 (*) 2.319 (16.1) 3.448 (23.9) 1.993 (10.5) -6.858 (-9.7) 2.414 (13.5) -2.268 (-7.8) -8.136 (-5.8) 2.617 (6.6) -6.341 (-5.7) 0 (*) 3.240 (27.3) 2.935 (17.9) 4.093 (19.9) 2.595 (16.3) 2.305 (12.5) 2.072 (18.5) -2.708 (-30.8) 3.538 (10.5) 1.301 (19.7) 0 (*) 5.707 (13.7) 1.963 (23.0) -4.578 (-22.1) 3.977 (14.7) 0 (*) 1.990 (17.7) 2.427 (38.3) -5.978 (-5.9) -2.443 (-23.5) 0 (*) -7.700 (-6.7) 2.217 (23.1) 0 (*) 0 (*) 4.350 (22.2) 0 (*) -7.839 (-12.1) 4.316 (7.4) 0 (*) 1.052 (2.3) -2.667 (-87.1) -2.768 (-130.3 2.065 (15.0) -3.003 (-48.5) -1.465 (-23.8) -0.7527 (-26.0) -2.032 (-159.9 -1.927 (-28.3) -1.558 (-83.7) -3.579 (-216.5 -2.389 (-15.3) -0.6978 (-10.9) 1.693 (20.2) -1.069 (-12.2) 0.01358 (0.2) -1.509 (-9.0) -2.259 (-46.3) -2.048 (-1.9) -1.430 (-31.9) -4.080 (-108.0 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -1.59e-4 (-9.7) -0.00114 (-69.0) -8.82e-4 (-21.3) -2.05e-4 (-10.1) -4.88e-4 (-31.6) -1.21e-4 (-13.2) -5.28e-4 (-54.9) -7.83e-4 (-17.6) -2.38e-4 (-30.9) -3.64e-4 (-65.2) -0.08031 (-21.6) -0.3127 (-86.9) -0.1995 (-19.1) -0.4795 (-63.1) -0.2202 (-43.1) -0.03749 (-21.3) -0.2823 (-131.6 -0.3022 (-28.1) -0.06158 (-46.5) -0.04475 (-43.6)

25

Schattingen BASGOED


3.4

Reistijdwaarderingen

Eén van de evaluatiecriteria van keuzemodellen is de reistijdwaardering die het model impliceert. Het grootste deel van de reistijdwaardering in modelspecificatie 2a, 2b en 2c volgt uit het tarief wat is gebruikt om de tijdsafhankelijke kosten uit te rekenen. Dit tarief is constant voor iedere NSTR groep maar varieert over de vervoerwijzen: 2.97 €/ton/uur voor wegvervoer, 0.41 €/ton/uur voor spoorvervoer en 0.10 €/ton/uur voor binnenvaart (bron: VKM, zie Tabel 2-3). Deze modelspecificatie is toegepast in model serie 6, 12, 13 en 14. In kostenspecificaties 4 (serie 8) en 2c (serie 12, 13 en 14) kunnen uit de geschatte tijdscoëfficiënt voor kapitaalkosten, aanvullende reistijdwaarderingen worden herleid. In deze reistijdwaardering zit de waardering voor interestkosten, ontwaarding en verzekeringskosten opgenomen. De kostenspecificaties zijn opgesteld uit het perspectief van de verlader, dus het grootste deel van de tijdskosten zit al via de tijdsafhankelijke kosten in de kostencoëfficiënt. De tijdcoëfficiënt echter bevat de tijdkosten van het kapitaal dat vervoerd wordt. Door deling van deze tijdscoëfficiënt door de kostencoëfficiënt is een reistijdwaardering voor de kapitaal kosten te berekenen. Deze zijn opgenomen in onderstaande tabel. De reistijdwaardering voor kapitaalkosten lijkt over het algemeen hoger voor NSTR groepen met hogere waarde dichtheid, zoals NSTR groep 5 (metalen en halffabricaten) en 9 (containers en overig). Tabel 3-6: Reistijdwaardering per goederengroep (per ton per uur).

NSTR1

Serie 6

Waardering kapitaalkosten (Euro/ton/uur): 0 Landbouwp. 1 Voedingsp. 2 Vaste min. br.stoff. 3 Aardolie 4 Ertsen 5 Metalen, halffabr. 6 Ruwe min. 7 Meststoff. 8 Chemie 9 Overig Reistijdwaardering (Euro/ton/uur): Weg 2.97 Spoor 0.41 Bvrt 0.10

Serie 8

Serie 12

Serie 13

Serie 14

0.25 0.21 0.11 0.05 0.19 0.38 0.06 0.11 0.20 0.48

0.23 0.21 0.04 0.20 0.42 0.09 0.14 0.20 0.49

0.18 0.21 0.09 0.16 0.38 0.10 0.16 0.21 0.66

0.12 0.22 0.27 0.03 0.13 0.19 0.11 0.16 0.23 0.49

2.97 0.41 0.10

2.97 0.41 0.10

2.97 0.41 0.10

De impliciete reistijdwaarderingen voor goederenkosten in de geschatte modellen worden als plausibel beoordeeld. Model serie 12 tot en met 14 lijken qua modelspecificatie het best bruikbaar voor implementatie in BASGOED. Allereerst zijn de kostenfuncties gebaseerd op het VKM, en is daardoor in aanpak consistent met SMILE+. Daarnaast is een aanvullende tijdcoëfficiënt geschat voor de kapitaal kosten, en deze parameter bleek significant (op NSTR groep 3 na) en iets toe te voegen aan de verklaringskracht van de modellen.


3.5

Schattingen BASGOED

Elasticiteiten

Daarnaast zijn de elasticiteit van de geschatte modellen een belangrijk evaluatie criterium. Als referentie voor deze elasticiteiten zijn uit twee recente empirische studies naar prijsgevoeligheid van goederenvervoer de volgende elasticiteiten bekend. Tabel 3-7: Kosten elasticiteiten goederenvervoer uit literatuur

Significance en CE Delft (2010) Price sensitivity of European road transport – towards a better understanding of existing results. Report for Transport & Environment. The Hague: Significance Tkm prijs elasticiteit op ton kilometers voor wegvervoer: wegvervoer Alle effecten -1.0 Alleen vervoerwijzekeuze-0.4 effect Significance en VTI (2010) Review of the international literature on price elasticities of freight transport by rail. Report for Swedish Rail Administration. The Hague/Stockholm: Significance/VTI Tkm prijs elasticiteit op tonnen voor spoorvervoer (vervoerwijzekeuze (vervoerwijzekeuzezekeuze-effect) ruime marge: -0.5 tot -3.9 waarschijnlijke marge: -0.8 tot -1.6

De elasticiteiten en kruiselasticiteiten voor tijd en kosten en voor de drie vervoerwijzen zijn berekend op basis van een groot aantal toepassingen van de geschatte modal split modellen. In de appendix zijn de elasticiteiten per NSTR groep weergegeven. In Tabel 3-8 staat een overzicht van de totale tijden kosten elasticiteiten over alle goederengroepen. Tabel 3-8: Tijd- en kostenelasticiteiten (in vet) en kruiselasticiteiten model serie 12. Elasticiteiten zijn berekend over totale vervoersvraag

Kosten elasticiteiten

Serie 12 Weg Spoor Bvrt Serie 13 Weg Spoor Bvrt Serie 14 Weg Spoor Bvrt

Tijd elasticiteiten

weg

spoor

bvrt

weg

spoor

bvrt

-0.7493 0.4430 0.9122

0.0221 -0.9556 0.0738

0.2796 0.4643 -0.4075

-0.2629 0.2319 0.3120

0.0068 -0.2575 0.0189

0.2628 0.6377 -0.4045

-0.4966 0.3792 0.5895

0.0180 -0.8742 0.0692

0.1864 0.4390 -0.2825

-0.1792 0.2003 0.2060

0.0055 -0.2259 0.0168

0.1845 0.6011 -0.2972

-0.6651 0.4546 0.8031

0.0221 -1.0309 0.0818

0.2526 0.5494 -0.3820

-0.2352 0.2387 0.2758

0.0068 -0.2898 0.0223

0.2410 0.9620 -0.4113

De kostenelasticiteit van wegvervoer ligt in model serie 12, met een waarde van -0.75, aan de bovenkant van de plausibele range. Model serie 14 ligt met kostenelasticiteit van -0.66 resp. meer in de richting van de plausibele waarde. Model serie 13 ligt met een kostenelasticiteit van -0.50 het dichtst bij de waarde -0.4 uit de literatuur. Door het grote marktaandeel van wegvervoer zijn de kruiselasticiteiten van kosten wegvervoer met binnenvaart en spoor hoog.

27


Schattingen BASGOED

De kostenelasticiteit van spoorvervoer ligt in alle series binnen de waarschijnlijke marge: -0.96, -0.87 en 1.03 respectievelijk. Net als in de literatuur is het spoor vervoer dus meer kostenelastisch. Over het algemeen is serie 13 minder elastisch dan 14 en 12 door het grote aantal vaste parameters in het model. In het geschatte model is het goederenvervoer over het algemeen minder tijdelastisch dan kostenelastisch. Dit geldt voor de totale vraag over alle goederengroepen, maar ook voor de tijdelasticiteit per goederengroep (zie appendix). Er zijn uit de literatuur helaas geen tijdelasticiteiten van goederenvervoer bekend ter vergelijking. De grotere tijdelasticiteit volgt ook uit de modelstructuur van model serie 12 (kosten functie 2c). In dit model bestaan de totale kosten uit een tijdsafhankelijk deel en een afstandsafhankelijk deel. In de berekening van de kostenelasticiteiten worden beide kostencomponenten verhoogd, terwijl bij de berekening van tijdelasticiteiten alleen de tijdskosten zijn verhoogd (en de kapitaalkosten). 3.6

Goodness-of-fit

De goodness-of-fit van de verschillende modelspecificaties is getoetst op basis van de statistische fit van de schattingen en de fit van de modellen in de toepassing. De statistische fit van de schattingen kan getoetst worden aan de hand van een aantal aan elkaar gerelateerde indicatoren die bij de schattingen berekend worden. Voor iedere specificatie is het model op dezelfde dataset geschat dus zijn de loglikelihood’s van de modellen onderling vergelijkbaar. Hoe lager de absolute waarde van de loglikelihood, hoe beter de model fit. De Rho2(0) en Rho2(c) worden berekend uit de final likelihood ten opzichte van een basis likelihood5. In principe geldt: hoe dichter de waarde ligt bij 1, hoe beter de model fit. Tabel 3-9 geeft een overzicht van deze indicatoren voor de model series 12, 13 en 14. Uit onderlinge vergelijking blijkt model serie 13 tot een betere model fit te komen. Tabel 3-9: Samenvatting van model fit schattingen Serie: 12 Final log (L)

13 Final log (L)

14 Final log (L)

Serie: 12 Rho²(0)

13

14

Rho²(c)

Rho²(0)

Rho²(c)

Rho²(0)

Rho²(c)

0

Landbouwp.

-18296.1

-15172.2

-15430.4

0.352

0.13

0.462

0.279

0.453

0.267

1

-31633.9

-23997.2

-29133.6

0.36

0.215

0.515

0.404

0.411

0.277

2

Voedingsp. Vaste min. br.stoff.

-13593.3

-11788.5

-12114

0.534

0.029

0.596

0.158

0.584

0.135

3

Aardolie

-16577.7

-12806.2

-14794.8

0.63

0.369

0.714

0.513

0.67

0.437

4

Ertsen

-8964

-7086.7

-7899.6

0.746

0.498

0.8

0.603

0.777

0.558

5

Metalen, halffabr.

6

Ruwe min.

7

Meststoff.

8

Chemie

9

Overig

-15681.5

-12723

-15007.3

0.119

0.11

0.285

0.278

0.157

0.148

-116109.4

-92417.9

-104295.7

0.195

0.164

0.359

0.334

0.277

0.249

-4503.6

-3284.5

-3744.2

0.25

0.232

0.453

0.44

0.376

0.362

-39589

-34840.2

-37249.6

0.298

0.18

0.382

0.278

0.34

0.228

-69137.3

-62179.5

-69135.2

0.52

0.369

0.568

0.432

0.52

0.369

De model fit is verder getoetst op basis van een toepassing van de geschatte modellen in een prototype van de BASGOED modules. Op basis van dit prototype is voor alle goederengroepen een toedeling gemaakt

5

De Rho2(0) wordt berekend ten opzichte van de initiele likelihood. De Rho2(c) ten opzichte van een model met een constante voor ieder alternatief: Rho2(0) = 1 - (final likelihood) / (likelihood with zero coeff.) Rho2(c) = 1 - (final likelihood) / (likelihood with constants)

28


Schattingen BASGOED

en is het resultaat vergeleken met de waarnemingen uit het schattingsbestand. De goodness-of-fit van het model wordt beoordeeld op basis van de root mean squared error (RMSE) per goederengroep, en op basis van scatterplots van de gesimuleerde en waargenomen vervoerde tonnages per vervoerwijze per herkomst of bestemming. De RMSE is als volgt berekend:

∑ (T

SIM ijv

RMSE =

− TijvWN

)

ijv

i⋅ j ⋅v

In de volgende tabel staan de RMSE’s zoals deze zijn berekend uit het prototype van het MS model. Deze RMSE heeft in absolute zin een grote waarde, en is moeilijk interpreteerbaar. De grote afwijking wordt veroorzaakt door een aantal grote relaties waarbij de voorspelde vervoerwijze verdeling van het model grof naast de waarnemingen ligt. Dit is moeilijk zichtbaar in de geaggregeerde RMSE maat maar wordt beter zichtbaar via een aantal scatterplots. Het opnemen van vervoerwijzedummies voor herkomsten en bestemmingen of specifieke corridors blijkt effectief in het reduceren van de RMSE. In de meeste gevallen reduceren de vervoerwijzedummies voor herkomsten en bestemmingen (serie 13) de RMSE het meest: voor NSTR 0, 1, 4, 5, 6, 7, 8 en 9. In NSTR groep 2 reduceren de vervoerwijzedummies voor corridors (serie 14) de RMSE het meest. Ten opzichte van serie 13 wordt de verbetering van RMSE in serie 14 met veel minder parameters gerealiseerd (10 corridor dummies).

Tabel 3-10: RMSE berekend uit toepassing serie 12-14

NSTR1 0 1 2 3 4 5 6 7 8 9

Landbouwp. Voedingsp. Vaste min. br.stoff. Aardolie Ertsen Metalen, halffabr. Ruwe min. Meststoff. Chemie Overig

Totaal tonnage [ton] 80273278 118396686 33721874 62565989 50149551 29927257 235129111 24775788 104926248 259085485

RMSE serie 12 9283 18719 5318 17950 6782 9203 53296 3671 19460 32053

RMSE serie 13 4576 7379 2616 6663 2316 5041 20966 1508 12800 15966

RMSE serie 14 5840 13431 1999 16374 2355 9078 33999 2511 15799 30928

In de volgende figuren is zijn de voorspelde tonnages per vervoerwijze op celniveau vergeleken met de waarnemingen. De scatterplots zijn opgesteld voor de NSTR groepen 0, 5 en 9. Uit deze figuren wordt goed duidelijk dat het model 12 voor specifieke relaties niet goed in staat is de vervoerwijzeverdeling te voorspellen. De extra dummies in model serie 13 en 14 laten de spreiding van de puntenwolk iets afnemen. Wel blijft de spreiding rond de diagonaal nog wel groot. Voor de scatterplots is ook de R2 berekend van de trendlijn door de punten. Op basis hiervan kan de model fit van de verschillende modelseries beoordeeld worden. Deze indicator drukt ongeveer hetzelfde uit als de RMSE, namelijk de goodness-of-fit van het model ten opzichte van de waarneming. In Tabel 3-11 staan de R2 van de onderzochte NSTR groepen weergegeven, onderscheiden naar vervoerwijze.

29


Schattingen BASGOED

Tabel 3-11: R2 per serie, NSTR-groep en vervoermiddel

NSTR0

NSTR5

NSTR9

weg

spoor

bvrt

weg

spoor

bvrt

weg

spoor

bvrt

Serie 12

0.97

0.94

0.27

0.73

0.48

0.82

0.98

0.78

0.94

Serie 13

0.98

1.00

0.62

0.84

0.90

0.88

0.99

0.94

0.96

Serie 14

0.99

0.94

0.72

0.73

0.51

0.77

0.98

0.78

0.94

Figuur 3-2: Scatterplot vervoerde tonnen BASGOED vs waarnemingen NSTR 0 voor serie 12

NSTR0

BASGOED

100000

Weg Spoor Bvrt

1 1

100000 Waarneming

30


Schattingen BASGOED

NSTR0 (detail)

BASGOED

100000

Weg Spoor Bvrt

1 1

100000 Waarneming


NSTR0

BASGOED

100000

Weg Spoor Bvrt

1 1

100000 Waarneming

31


Schattingen BASGOED

NSTR0 (detail)

BASGOED

100000

Weg Spoor Bvrt

1 1

100000 Waarneming


NSTR0

BASGOED

100000

Weg Spoor Bvrt

1 1

100000 Waarneming

32


Schattingen BASGOED

NSTR0 (detail)

BASGOED

100000

Weg Spoor Bvrt

1 1

100000 Waarneming

33


Schattingen BASGOED


NSTR5

BASGOED

100000

Weg Spoor Bvrt

1 1

100000 Waarneming

NSTR5 (detail)

BASGOED

10000 Weg Spoor Bvrt

1 1

10000 Waarneming

34


Schattingen BASGOED


NSTR5

BASGOED

100000

Weg Spoor Bvrt

1 1

100000 Waarneming

NSTR5 (detail)

BASGOED


1 1

10000 Waarneming

35


Schattingen BASGOED


NSTR5

BASGOED

100000

Weg Spoor Bvrt

1 1

100000 Waarneming

NSTR5 (detail)

BASGOED


1 1

10000 Waarneming

36


Schattingen BASGOED


NSTR9

BASGOED

100000

Weg Spoor Bvrt

1 1

100000 Waarneming

NSTR9 (detail)

BASGOED


1 1

10000 Waarneming

37


Schattingen BASGOED


NSTR9

BASGOED

100000

Weg Spoor Bvrt

1 1

100000 Waarneming

NSTR9 (detail)

BASGOED


1 1

10000 Waarneming

38


Schattingen BASGOED


NSTR9

BASGOED

100000

Weg Spoor Bvrt

1 1

100000 Waarneming

NSTR9 (detail)

BASGOED


1 1

10000 Waarneming

39


3.7

Schattingen BASGOED

Conclusies

Bij de schattingsresultaten worden de volgende conclusies getrokken: •

Model serie 12, die gebaseerd is op de kostenkentallen uit het VKM blijkt een geschikt basismodel te zijn, waarop de uitgebreidere specificaties (serie 13 met herkomsten bestemmingsdummies en serie 14 met corridordummies) zijn gebaseerd. In deze gekozen specificatie worden de tijdsafhankelijke- en kilometerafhankelijke transportkosten gebaseerd op de VKM-kostenkentallen, maar wordt ook een additionele tijdparameter geschat voor de kapitaalkosten van de goederen. In deze kapitaalkosten zijn de kosten voor interestkosten, ontwaarding en verzekeringskosten opgenomen. De nutsfunctie voor vervoerwijze v op relatie ij voor goederengroep g is hier als volgt gespecificeerd:


(2c)

d ijv en tijv zijn de afstand en tijd op de voedings- en verbindingslinks, De afstand en tijd, zijn via de gemiddelde afstandskosten Tv en het tarief per tijdseenheid Rv omgerekend naar kosten. CONTijg is de containerisatiegraad (verhoogde kans op spoor of binnenvaart als de goederen per container worden vervoerd) en MSCv de vervoerwijze specifieke constante (deze kan ook variëren per herkomst, per bestemming of per relatie, afhankelijk van serie 12, 13 of 14). De β’s zijn te schatten kosten- en tijdparameters. Per NSTR hoofdgroep hebben de modellen significante tijd en kosten coëfficiënten met plausibel teken, en plausibele reistijdwaardering. De MSC’s zijn in alle modellen significant en plausibel. De dummy voor containerisatiegraad is in veel gevallen ook significant. Voor internationaal vervoer zijn voor alle goederengroepen, per vervoerwijze een dummy meegeschat. Deze verbeteren de model fit, en de reistijdwaardering in het model. •

•

De elasticiteiten van de specificaties lijken in een plausibele range te liggen. o

In model serie 12 ligt de kostenelasticiteit van wegvervoer met een waarde van -0.75, aan de bovenkant van de plausibele range. Model serie 13 en 14 liggen met kostenelasticiteiten van 0.48 en -0.66 resp. beter in de plausibele range.

o

De kostenelasticiteit van spoorvervoer liggen in alle series binnen de waarschijnlijke marge: 0.96, -0.87 en -1.03 respectievelijk, en is serie 13 minder elastisch dan 14 en 12 (door het grote aantal vaste parameters in het model).

o

In de geschatte modellen is het goederenvervoer over het algemeen minder tijdelastisch dan kostenelastisch.

De extra dummies in model serie 13 en 14 lijkt de goodness of fit van de modellen wel te verbeteren maar op OD niveau lijkt de afwijking ten opzichte van de waarnemingen nog wel groot. o

De extra dummies in model serie 13 en 14 lijkt de spreiding van de gepresenteerde puntenwolken wel iets af te laten nemen, maar nog steeds blijft de spreiding rond de diagonaal groot. Het opnemen van vervoerwijzedummies voor herkomsten en bestemmingen of specifieke corridors blijkt effectief in het reduceren van de RMSE.

40


o

Schattingen BASGOED

In de meeste gevallen reduceren de vervoerwijzedummies voor herkomsten en bestemmingen de RMSE het meest: voor NSTR 0, 1, 4, 5, 6, 7, 8 en 9. In NSTR groep 2 reduceren de vervoerwijzedummies voor corridors de RMSE het meest. Ten opzichte van serie 13 wordt de verbetering van RMSE in serie 14 met veel minder parameters gerealiseerd (10 corridor dummies).

De schattingen hebben twee specificatie opgeleverd die geschikt zijn als eindmodel voor de modal split module van BASGOED. Het gaat om specificatie 13 met herkomstdummies en bestemmingsdummies per vervoerwijze, en specificatie 14 met vervoerwijzedummies op corridors. Op basis van de schattingen en tests van de modal split module is nog geen specificatie geselecteerd, hoewel er een duidelijke voorkeur is voor specificatie 13. Beide modelspecificaties hebben veel constanten en het gevaar ligt op de loer dat de fit een te belangrijk criterium wordt. Hierbij kan opgemerkt worden dat de herkomstdummies en de bestemmingsdummies (versie 13) kunnen worden beschouwd als nodig voor het oppikken van kenmerken van de herkomsten bestemmingszones (b.v. op het gebied van aansluiting op de netwerken voor spoor en binnenvaart) die onvoldoende door de andere variabelen (inclusief de berekende voor- en natransporten over de weg) worden gerepresenteerd. Voor de corridordumies (versie 14) is de balans tussen fit en inhoudelijke verklaring slechter. Eerst zijn beide modelspecificaties geïmplementeerd en zijn uit beide modellen logsoms afgeleid voor de schattingen van de distributiemodule. De model fit van de distributie module schattingen worden gebruikt om te toetsen welke logsoms tot een betere model fit leiden. Deze informatie wordt meegenomen in het vaststellen van de te prefereren specificatie.

41


Schattingen BASGOED

42

HOOFDSTUK 4

4.1

Schatting distributiemodel

Werking van de distributiemodule

Het doel van de distributiemodule is om een herkomst/bestemmingen (H/B) matrix te produceren voor het prognosejaar. De input voor deze module is de H/B matrix in het basisjaar, de totale instroom en uitstroom per zone in het prognosejaar, en de weerstanden per H/B relatie. In de distributiemodule worden groeifactoren per H/B relatie berekend waarmee de basismatrix kan worden opgehoogd naar de prognosematrix. Zowel voor het basisjaar als voor het prognosejaar worden de randen van de matrix gebruikt en de weerstanden per H/B relatie om een matrix te genereren. Uit deze twee matrices volgen dan de groeifactoren. Aan de basis van de distributiemodule staat een zwaartekrachtmodel. Deze beschrijft de jaarlijkse stroom goederen T (in tonnen) tussen regio i en j als het produkt van drie variabelen: de zwaarte van de herkomstregio (uitgebeeld met een productieparameter pi), de zwaarte van de bestemmingsregio (een attractieparameter qj) en een weerstandsvariabele die per relatie varieert (rij). Het basismodel is als volgt:

Tij = p i ⋅ q j ⋅ rij

(5)

De variabele rij neemt hogere waarden aan bij een lagere weerstand. Het is derhalve een distributiefunctie met als argument de reële, gegeneraliseerde transportweerstand cij. Voor internationale relaties bevat rij ook een grensweerstand Gij, gemeten in dezelfde kosteneenheid als de gegeneraliseerde transportweerstand. Deze grensweerstand is nodig om de negatieve invloed te beschrijven van economische, administratieve en culturele verschillen tussen landen op de ruimtelijke verdeling van internationale handelsstromen. Als de zones i en j in hetzelfde land liggen dan geldt Gij = 0. De distributiefunctie is een combinatie van een negatief exponentiële functie en een negatieve machtsfunctie. De functie is flexibel en wordt standaard veel gebruikt bij zwaartekrachtmodellen. Het verwachte effect dat een hogere weerstand leidt tot kleinere stromen kan door de functie goed worden gemodelleerd.

rij = exp(−α ⋅ (cij + Gij )) ⋅ (cij + Gij ) − β

(6)

Als α = 0, en β ≠ 0, dan blijft alleen de negatieve machtsfunctie over. Als β = 0, en α ≠ 0, dan blijft alleen de negatief exponentiële functie over. Figuur 4-1 toont verschillende vormen van de distributiefunctie.


Schattingen BASGOED

Distributiefunctie, α = 0.5 en β = 0 1 0.8 0.6 0.4 0.2 0 0.1

0.5

0.9

1.3

1.7

2.1

2.5

2.1

2.5

c

Distributiefunctie, α = 0 en β = 0.5 3.5 3 2.5 2 1.5 1 0.5 0 0.1

0.5

0.9

1.3

1.7

c

Distributiefunctie, α = 0.5 en β = -0.25 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 0.1

0.5

0.9

1.3

1.7

2.1

c

Figuur 4-1: verschillende vormen van de distributiefunctie

44

2.5

Schattingen BASGOED


Bij het schatten van een H/B matrix worden de rij berekend, gegeven de parameters α en β, de weerstanden cij en de grensweerstanden Gij. Om nu Tij te berekenen is het niet nodig om de productie- en attractieparameters pi en qj te kennen. De rij vormen namelijk een startmatrix die als invoer wordt gegeven aan een Furness procedure. Bij deze procedure worden de waarden in de matrix op dusdanige wijze vermenigvuldigd dat om en om aan de rij- en kolomvoorwaarden (de totale instroom en uitstroom per zone) wordt voldaan. De matrix convergeert zo naar een kloppende H/B matrix. De parameters pi en qj zijn het product van de opeenvolgende rij- en kolomfactoren die tijdens de procedure zijn berekend en hebben in principe geen relevante betekenis.

4.2

Schatten parameters

De waarde van cij is bekend; dit zijn de logsom weerstanden uit de MS module. De parameters α en β en de grensweerstanden Gij moeten dus nog geschat worden. Dit gebeurt op basis van de basisbestanden goederenvervoer 2004. Per goederengroep worden de parameters apart geschat. De grensweerstand G kan uniform worden geschat, of er kan een aparte grensweerstand Gi worden geschat voor elke buitenlandse zone. Beide varianten zijn hieronder uitgewerkt6. Een ander punt van aandacht is het intrazonale vervoer. De intrazonale stromen zijn in veel gevallen significant hoger dan de interzonale en internationale stromen, terwijl dit niet is terug te zien in de logsums uit het MS model. Een optie is om de intrazonale logsums di apart te schatten (endogeen), of het hele model te schatten zonder intrazonaal vervoer. In het laatste geval zou bij toepassingen het intrazonale vervoer uiteraard wel weer meedoen. Met dit in het achterhoofd kunnen verschillende varianten van de distributiefunctie worden geschat. Het wordt een keuze tussen α en β (of combinatie ervan), wel of geen grensweerstanden, en wel of geen endogene intrazonale logsums. Ook kunnen de parameters geschat worden zonder intrazonale of internationale stromen. De logsom weerstanden cij zijn niet meteen te gebruiken, aangezien ze zowel positief als negatief kunnen zijn. Ook geldt dat bij een hogere logsom de weerstand juist lager is. Om de logsom weerstanden wel bruikbaar te maken wordt per goederengroep de maximum logsom berekend. Vervolgens wordt dit maximum van elke logsom afgetrokken, en worden de logsoms vermenigvuldigd met -1 (om de tekens om te draaien). Na deze handeling zijn de weerstanden wel bruikbaar. Het schatten gebeurt met behulp van de Newton methode. Hier wordt de gradiënt berekend bij een willekeurige combinatie van parameters (als startwaarden worden alle parameters gelijk aan 0 genomen). De gradiënt geeft aan in welke richting de parameters moeten bewegen voor een betere fit. Zo kan dus een steeds betere oplossing worden gezocht, totdat er geen betere oplossing kan worden gevonden. In dat geval is de gradiënt gelijk aan 0 en is het optimum bereikt. Het optimum is hier een geschatte matrix die zo dicht mogelijk bij de werkelijke matrix ligt. Om de fit te kwantificeren wordt de determinatiecoëfficient R² gehanteerd. Soms komen er veel nullen voor in de basismatrix. Om de parameters zuiverder te kunnen schatten wordt genomen rij = 0 waar xij = 0 (als x de waarden zijn in de basismatrix). 6

Een andere variant zou zijn om aparte grensweerstand Gij voor elke internationale relatie op te nemen. Dit betekent dat er voor elke buitenlandse cel in de matrix een aparte weerstand wordt geschat, wat neerkomt op 2320 extra parameters. Aangezien niet alle matrices volledig gevuld zijn zou het kunnen voorkomen dat het aantal parameters het aantal observaties overtreft, dit zou leiden overspecificatie.

45


4.3

Schattingen BASGOED

Resultaten

In tabellen 4-1 en 4-2 staan de specificaties van de verschillende varianten die zijn geprobeerd en de gemiddelde R² per goederengroep. De distributiefuncties zijn geschat voor twee sets logsom weerstanden uit de modal split module (model series 13 en 14). Een distributiefunctie is niet stabiel als deze een nietintuïtief verloop (meer transport bij toenemende weerstand) kent over een bepaald bereik of als de geschatte parameters een extreme waarde aannemen. Het aantal niet stabiele distributiefuncties per variant wordt ook in de tabel weergegeven. Varianten 1, 2 en 3 hebben geen grensweerstanden of endogene intrazonale weerstanden. Deze varianten zijn relatief stabiel, maar de behaalde fit is behoorlijk laag. Ook kent goederengroep NSTR3 een nietintuïtief verloop bij deze varianten. Bij de varianten 4a t/m 4d worden de grensweerstanden G (vast of variabel) geschat. Deze varianten doen het iets beter dan de varianten zonder grensweerstanden. Het verschil is echter niet overtuigend, wetende dat het toevoegen van parameters sowieso tot een verbetering van de fit leidt (de R² neemt altijd toe bij het schatten van extra parameters, zelfs als deze niet significant zijn). De internationale stromen zijn over het algemeen minder groot dan de binnenlandse stromen. Waarschijnlijk is dit de reden voor het relatief kleine effect van de grensweerstanden. De varianten 5a en 5b met endogene intrazonale weerstanden d zorgen wel voor een significante verbetering. Dit toont aan dat de bestaande logsums op zichzelf niet voldoende zijn om de intrazonale stromen goed te schatten. Variant 5b, een negatief exponentiële functie met een aparte intrazonale weerstand per zone, leidt tot de beste fit (m.u.v. variant 6, zie hieronder). Daarmee lijkt variant 5b de specificatie te zijn die het meest geschikt is voor de distributie module in BASGOED. Bij deze variant zijn de logsums uit MS model 13 effectiever dan die uit MS model 14. Variant 6 is een combinatie van alle mogelijke parameters. Het schatten van zo’n grote hoeveelheid parameters is vrij gecompliceerd. Hoewel de variant logischerwijs tot de beste fit leidt van alle varianten, is het niet aan te raden om deze te implementeren. Het aantal parameters blijft overigens in elke variant laag genoeg om een groot aantal vrijheidsgraden te garanderen. Er is dus nooit sprake van overfitting. Bij variant 6 worden 71 parameters geschat. Het aantal positieve cellen in de basismatrix varieert van 407 bij NSTR2 tot 3365 bij NSTR9. In tabel 4-3 staan voor variant 5b de geschatte parameters per goederengroep. In de bijlage (bijlage 3) worden alle schattingsresultaten (inclusief d) uitgebreid uiteengezet. Het blijkt dat bij variant 5b voor NSTR4 een bijzonder hoge d is geschat (zone 29). Dit geeft aan dat het verdelen van het intrazonale vervoer in zone 29 over de interzonale stromen van die zone resulteert in een betere fit. Dit is echter niet gewenst, zeker niet bij toepassingen van het model. Voor NSTR4 is het daarom beter om te kiezen voor een stabielere variant. Variant 1, geschat zonder intrazonale stromen, is hiervoor het meest geschikt.

46

Schattingen BASGOED


Tabel 4-1: Behaalde fit van verschillende varianten distributiefunctie, geschat incl. intrazonale stromen specificatie variant

α

1

x

G

β

2

MS-model 13 Gem. niet R² stabiel

d


0.615

1

0.607

1

x

0.608

1

0.621

1

0.618

1

0.650

1

3

x

x

4a

x

x

vast

0.627

3

0.657

2

4b

x

x

var

0.638

3

0.659

2

4c

x

var

0.615

1

0.607

1

var

0.625

2

0.635

1

4d

x

5a

x

vast

0.724

1

0.753

2

5b

x

var

0.815

1

0.784

2

6

x

var

0.831

0

0.809

3

x

var

Tabel 4-2: Behaalde fit van verschillende varianten distributiefunctie, geschat excl. intrazonale stromen specificatie variant

α

1

x

G

β

2



0.648

0

0.662

1

x

0.644

0

0.662

1

0.661

1

0.680

1

3

x

x

4a

x

x

vast

0.658

4

0.686

3

4b

x

x

var

0.676

1

0.686

1

4c

x

var

0.648

0

0.665

1

var

0.663

0

0.676

1

4d

x

Tabel 4-3: De geschatte parameters bij variant 5b MS model 13 NSTR

MS model 14

α

R²

α

R²

0

1.602

0.680

0.878

0.665

1

0.912

0.771

0.836

0.518

2

2.752

0.978

0.324

0.971

3

1.363

0.945

0.103

0.905

4

0.887

0.984

6.606

0.992

5

0.539

0.637

1.084

0.674

6

1.462

0.855

1.460

0.817

7

0.969

0.670

1.005

0.717

8

2.148

0.817

1.565

0.737

9

1.016

0.815

2.078

0.847

Gem.

0.815

0.784

In Figuur 4-2 en Figuur 4-3 zijn de distributie functies gevisualiseerd voor varianten 1 en 5b. Hier is goed te zien dat bij variant 1 de functie voor NSTR3 (olieproducten) niet dalend is zoals verwacht. Een stijgende functie betekent dat de stromen hoger zijn bij hogere weerstanden. Door endogene intrazonale weerstanden in te voeren (variant 5b) krijgt ook de functie voor NSTR3 een positieve parameter. De eerder genoemde instabiliteit bij NSTR4 is enkel terug te zien in de geschatte di en dus niet in her verloop van de distributiefunctie.

47


Schattingen BASGOED

Figuur 4-2: vorm van de distributiefuncties; variant 1, enkel negatief exponentieel

48

Schattingen BASGOED


Figuur 4-3: vorm van de distributiefuncties; variant 5b, negatief exponentieel met endogene intrazonale weerstanden

49


Schattingen BASGOED

Conclusies

4.4

Algemene conclusies bij de schattingen van de distributie module: •

In de meest simpele variant (variant 1) geven de geschatte parameters bij NSTR3 (Aardolie) aan dat er grotere goederenstromen zijn bij een hogere weerstand, wat niet intuïtief is. Dit is echter niet de gekozen specificatie. Bij het endogeen schatten van de intrazonale weerstanden (de gekozen specificatie) kent de distributiefunctie bij NSTR3 wel een dalend verloop.

•

De varianten met grensweerstanden hebben logischerwijs een hogere R² dan de varianten zonder (aangezien het toevoegen van variabelen leidt tot een hogere R²), maar het verschil is over het algemeen klein. Ook zijn de geschatte grensweerstanden vaak niet stabiel en is dit gedrag moeilijk te verklaren. Daarom is het weinig zinvol om voor een variant met grensweerstanden te kiezen. Het is belangrijker dat de module goed blijft functioneren, ook als de input drastisch verandert.

•

De varianten met aparte intrazonale weerstanden doen het veel beter dan de andere varianten. Over het algemeen zijn intrazonale stromen significant hoger dan interzonale en internationale stromen, maar dit komt niet naar voren in de logsoms uit het MS model. Het is dus aan te bevelen om deze logsoms in de DM te schatten. Variant 5b, de simpele variant van de distributiefunctie (alleen alfa) met een aparte intrazonale logsom per zone, geniet duidelijk de voorkeur voor implementatie in BASGOED. Omdat er voor NSTR4 een vreemde di wordt geschat, wordt voor deze goederengroep gekozen voor variant 1.

•

De implementatie van variant 5b is bijna even eenvoudig als die van variant 1. De geschatte di vervangen de intrazonale logsom inputs. Verder verandert er niets aan de distributiefunctie.

•

Over de gehele breedte presteert het distributiemodel op basis van de logsoms uit MS model 14 gemiddeld iets beter dan die van MS model 13. Als echter alleen naar variant 5b wordt gekeken dan presteren de logsoms uit MS model 13 beter. Ook leveren de logsums uit MS model 13 over het algemeen stabielere resultaten op dan die uit MS model 14.

•

Concluderend kan gesteld worden dat variant 5b wordt gekozen voor elke goederengroep behalve voor NSTR4. Hiervoor worden de intrazonale logsoms niet endogeen geschat en wordt variant 1 gebruikt. Zowel variant 1 als variant 5b bevat een simpele variant van de distributiefunctie (alleen alfa). Deze ziet er als volgt uit:

rij = exp(−α ⋅ cij ) De weerstand r tussen zones i en j wordt hier weergegeven als een functie van de logsums c uit de modal split module. Variant 5b heeft endogeen geschatte logsums d voor intrazonale relaties. Dit betekent dat in deze variant cij = di als i = j.

50

HOOFDSTUK 5

Aandachtspunten voor de implementatie

Dit hoofdstuk beschrijft de aandachtspunten van de implementatie van het BASGOED modelsysteem als geheel, gezien de ontwikkeling van de stand-alone DM en MS modules binnen deelproject 1/2. Achtereenvolgens wordt ingegaan op de achtergrond van de voorziene implementatie, het zogenoemde alternatief 1b uit de ontwerpfase; de huidige stand van zaken en de mogelijkheden / gevolgen voor implementatie.

5.1

Achtergrond alternatief 1b

De functionele en technische specificaties van BasGoed staan beschreven in een rapport van TNO, NEA, Significance en Demis (Duijnisveld et al., 2010). Uit de verschillende technische alternatieven voor de implementatie van BASGOED is door DVS de keuze gemaakt voor het alternatief 1b, waarvan de belangrijkste eigenschappen hier worden beschreven. In dit alternatief worden zoveel mogelijk de bestaande Smile+ modules hergebruikt en een nieuwe (rudimentaire) schil geprogrammeerd in Visual Studio .NET. Smile+ bestaat uit: 1. De Dataman module voor het genereren van een lege variantdatabase en van de lees- en schrijfroutines voor SmileDE bij aanpassing van het gegevensmodel via de metadatabase 2. Rekenmodules SmileECO, SmileTM, SmileGFM en SmileMIL 3. Modules voor data toegang en bewerking SmileGeo, SmileDE en SmileErr 4. Gebruikersinterface modules SmileInter, SmileIn, CalibratieShell, SmileCL en BIN2DB Vanwege de wijzigingen in de rekenmodules en de nieuwe en vervallen rekenmodules is het nodig voor het BasGoed model om een nieuw gegevensmodel te definieren. Omdat Smile+ beschikt over een metadatabase om het datamodel op te geven, zal het goed mogelijk zijn om met de Dataman module het nieuwe gegevensmodel voor het BasGoed model te verwerken tot een lege variantdatabase met bijbehorende lees- en schrijfroutines voor SmileDE. Groot voordeel van de eenvoudigere opzet van het BasGoed model is bovendien dat er minder complexe uitvoerdata nodig zijn, zodat het werken met alleen de Access database zonder binaire hulpbestanden mogelijk lijkt te worden. De rekenmodules zijn allen in C++ geprogrammeerd en het omzetten hiervan naar de Visual Studio .NET omgeving en als programmeertaal C# is naar verwachting geen probleem. Van de rekenmodules vervallen de modules SmileTM en SmileMIL, SmileGFM zal mogelijk enigzins aangepast moeten worden. Daarnaast zijn twee nieuwe rekenmodules nodig, SmileDM voor de distributie en SmileMS voor de mode split. De nieuwe ‘Conversiemodules’ SmileImp en SmileExp resp. voor het inlezen van - en wegschrijven naar LMS/Bivas/ROUTGOED zijn eigenlijk meer datatoegangsmodules dan rekenmodules.


Schattingen BASGOED

Voorgestelde minimale gebruikersinterface Van de gebruikersinterface modules wordt aanbevolen om een nieuwe SmileInter module te maken met een minimale invulling. Dit opdat een zelfstandig rekenend model onstaat wat in een later stadium makkelijk in bijvoorbeeld Cube gehangen kan worden. Bovendien sluit dit ook aan bij de manier waarop het huidige Smile+ meestal wordt gebruikt. De minimale gebruikersinterface wordt een Visual Studio.Net applicatie en krijgt de volgende functionaliteiten: 1. Een menu om invoerbestand(en) aan te wijzen. 2. Een menu om eventuele rekeninstellingen (bijv. keuze te gebruiken modules) te wijzigen. 3. Een knop om de som te starten. Daarna toont een scherm de voortgang ongeveer op dezelfde wijze als dat in Smile+ wordt gedaan. 4. Indien relevant kan BIN2DB gebruikt blijven worden om binaire bestanden om te zetten naar Access en omgekeerd.

Gevolg van deze aanpak van de schil is dat de geografische datatoegangsmodule SmileGEO kan vervallen. Overigens kan als dat noodzakelijk / handig blijkt te zijn, in plaats van MS-Access bestanden ook gebruikt gemaakt gaan worden van MS-SQL Server bestanden. De code generatieprocedure van de Dataman module kan namelijk vermoedelijk zonder erg veel problemen ook lees- schrijfroutines maken voor MSSQL Server. Architectuur BasGoed model In onderstaande figuur is de voorgestelde architectuur bij alternatief 1b van het BasGoed model uitgewerkt per module en in 4 lagen: gebruikersinterface; rekenmodules; data toegang- en bewerking; en de gegevens.

Gebruikersinterface Architectuur BasGoed model SmileInter

Onderstreepte modules zijn nieuw

SmileCL

BIN2DB

Installatie

Rekenmodules

Dataman

SmileECO

SmileGFM

SmileDM

SmileMS

Data toegang en bewerking

SmileDE

SmileImp

SmileExp

SmileErr

Gegevens

MetaDB

Code bestanden SmileDE

SmileDE _ Proto

VariantDB

52

Binaire variant bestanden

Error log

Interface data DB

Schattingen BASGOED


In deze architectuur wordt eerst het nieuwe gegevensmodel geïmplementeerd via de metadatabase en de Dataman module, zodat er een lege variantdatabase en bijbehorende lees- en schrijfroutines ontstaan. Vervolgens worden de nieuwe rekenmodules SmileDM en SmileMS gebouwd, samen met de SmileImp en SmileExp conversiemodules. Tot slot wordt de zeer rudimentaire schil SmileInter gebouwd die de oude vervangt. De nieuwe rekenmodules worden in deze architectuur aangesloten via de door de Dataman module gegenereerde lees- en schrijfmodules op de BASGOED database. Deze functies worden gegenereerd op basis van de metadatabase door DataMan. Het ene deel bestaat uit functies die grenzen van een index teruggeven, het andere deel bestaat uit functies die variabelenwaarden lezen of schrijven in de database. Deze functies hebben de volgende namen en argumenten:

Type

Omschrijving

Argumenten

Indexgrens ophalen

GetIndexnaamMin en

GetIndexnaamMax Variabelenwaarden lezen

GetVariabelenaam

indexcombinatie

Variabelenwaarden schrijven

SetVariabelenaam

indexcombinatie

5.2

Stand van zaken prototype modules DM en MS

Bij het maken van de offerte voor het onderhavige project is als optie A aangeboden om ten behoeve van een soepele integratie alvast het nieuwe datamodel (en gevulde variantdatabase) met bijbehorende SMILEDE te bouwen met daarin de bovengenoemde lees- en schrijfroutines teneinde goed voorbereid te zijn op integratie met het BASGOED raamwerk. Die optie is door DVS echter niet gehonoreerd en daardoor was het voor het consortium niet mogelijk om het BASGOED raamwerk te gebruiken in dp1/dp2. Tevens is het niet gelukt om de realisatie van het BASGOED raamwerk in deelproject 3 (deels) gelijktijdig met dp2 te laten lopen. Daardoor is ook langs die weg geen bruikbaar BASGOED raamwerk voor dp1/dp2 beschikbaar gekomen. In het kader van het maken van de parameterschattingen van dp1 en het testen van de “goodness of fit” is daarom noodgedwongen van zowel de distributiemodule als van de modal split module een stand-alone prototype gemaakt, zonder gebruik te maken van het beoogde BASGOED raamwerk. Het integreren van deze gebouwde DM en MS prototype modules in de BASGOED architectuur zal daarom een stuk minder eenvoudig gaan dan mogelijk was geweest bij uitvoering van optie A van de offerte. Dit stand-alone prototype voor DM en MS is geprogrammeerd in Delphi en maakt gebruik van ASCII invoer- en uitvoerbestanden. Dit komt overeen met de manier waarop LMS en NRM zijn geprogrameerd en in de parktijk worden gebruikt.

53


5.3

Schattingen BASGOED

Mogelijkheden voor implementatie modules DM en MS in dp2

In deelproject 2 moeten in fase 1 de functionele en technische specificaties van DM en MS nader uitgewerkt worden tot een werkplan, testplan en technische specificaties voor implementatie van het model. In fase 2 worden de DM en MS modules gebouwd en getest volgens de specificaties uit fase 1 en in fase 3 wordt de documentatie gemaakt conform de structuur van de huidige functionele en technische documentatie van Smile+. De huidige stand van zaken is dat er geen gebruik gemaakt kan worden van een BASGOED raamwerk met afgesproken datamodel, gevulde variantdatabase, C# lees- en schrijfroutines in de Visual Studio .NET ontwikkelomgeving. In feite is het hierdoor erg lastig om bij de implementatie van de modules DM en MS alvast rekening te houden met de BASGOED architectuur. Er liggen nu 2 mogelijke implementatie scenario’s voor waar uit gekozen moet worden binnen de beschikbare middelen van deelproject 2: 1. De huidige prototypes worden in Delphi verder uitgebouwd tot stand-alone programmatuur (eventueel door ze als dynamic link library DLL aanroepbaar te maken). 2. De implementatie vindt plaats in C# gebaseerd op de principes van het BASGOED raamwerk, gebruik makend van een “dummy” SMILEDE bibliotheek voor datatoegang met de lees- en schrijfroutines conform de BASGOED architectuur en gebruik makend van een gevulde variantdatabase. Ad. 1 Voortgaan op basis huidige prototypes De eerste variant is in deelproject 2 kader goed realiseerbaar en is zonder grote risico’s voor deelproject 2 uit te voeren. Deze variant past ook binnen het dp2 budget. Het integreren met de BASGOED architectuur is daarentegen een stuk meer werk en kan zeker niet binnen een paar dagen uitgevoerd worden, zoals in de offerteaanvraag gevraagd. Bij de integratie dient er namelijk voor gezorgd te worden dat de invoerbestanden voor zowel de DM als MS modules gereed staan. Dit kan door de betreffende bestanden via nieuw te maken modules weg te schijven vanuit de BASGOED variantdatabase, vervolgens de berekening stand-alone uit te laten voeren, en de resultaten uit de resultaatbestanden te importeren naar de BASGOED variantdatabase via nieuw te maken importmodules. Het komt er dus op neer dat zowel voor DM als MS een preprocessor nodig zijn en ook voor beiden ook een postprocessor. Hierbij ontstaat dus zowel in het BASGOED raamwerk code om de benodigde tussenbestanden te lezen en te schrijven, die in feite dubbel is want in DM en MS zit dat ook al. In deze variant wordt de integratie en controle op datamodel integriteit geheel de verantwoordelijkheid van deelproject 3. Vanuit oogpunt van onderhoudbaarheid en documentatie is een mix van programmeertalen zoals deels Delphi en deels C# niet optimaal te noemen. Het doen van modelruns voor nieuwe scenario’s en beleidsvarianten verloopt in de praktijk door de invoerbestanden te wijzigen (dit kan handmatig, maar ook via een bestand met te wijzigen instellingen, zoals voor LMS/NRM). De Delphi modules voor MS en DM draaien binnen enkele seconden, zodat voortgangscontrole niet nodig is. Ook wordt er door dit programma een logfile bijgehouden, met informatie over de uitgevoerde run. Ad. 2 Implementatie in C# op BASGOED principes Deze variant is niet zonder gevolgen in deelproject 2 te realiseren vanwege het ontbreken van een goedgekeurde gevulde variantdatabase en bijbehorende lees- en schrijfroutines. Bovendien moet de inhoud van alle invoerbestanden van DM en MS overgezet worden naar de variantdatabase. Door dit traject binnen deelproject 2 te trekken is het beschikbare budget niet toereikend, het komt immers neer op het uitvoeren van optie A uit de offerte. Vanuit integratie en onderhoudbaarheid van BASGOED als geheel is deze variant wel beter. In deze variant wordt eerst het nieuwe datamodel vastgelegd in de metadatabase. 54

Schattingen BASGOED


Dan wordt met de dataman module de lees- en schrijfroutines gegenereerd en de lege variantdatabase aangemaakt. Dan wordt een lege DM en MS module gemaakt in een werkend C# project compleet met BASGOED raamwerk, waar de bouwer van DM en MS dan de algoritmes in kan bouwen. De data uit bestaande tussenbestanden moet geheel opgenomen worden in de variant database. Onderstaande tabel geeft in het kort de voor- en nadelen van beide varianten weer. Variant

Voordelen

Nadelen •

Dubbele code voor alle lees- en schrijfacties voor de meer dan 10 tussenbestanden

•

DP2 is binnen budget uitvoerbaar en op tijd af te maken

Onderhoudbaarheid is minder vanwege de mix van Delphi en C# programmeertalen

•

Het toevoegen van nieuwe modules is relatief eenvoudig

Integratie met DP3 kost veel meer dan de geraamde 2 dagen

•

DM en MS volgen een eigen aanpak, die niet lijkt op de architectuur van BASGOED

•

DM en MS kunnen pas pas getest worden als het BASGOED raamwerk is geïmplementeerd en data is opgenomen in de variantdatabase

•

Oplevering dp2 wordt vertraagd

•

Bij kleine wijzigingen in DM en MS moet de gehele applicatie opniew worden gecompileerd

•

Goed stand-alone te testen (al uitgevoerd)

•

Kleine wijzigingen intern in de DM en MS module makkelijk te implementeren

• 1 Voortbouwen op huidig prototype •

•

DM en MS bouwen voort op LMS en NRM methodiek

•

Beter onderhoudbaar want alle code is in C# en volgt de BASGOED architectuur

•

Deelproject 3 is eenvoudiger uitvoerbaar en met minder risico’s

•

Dubbele lees- en schrijfcode wordt vermeden want tussenbestanden niet nodig

2 Alsnog implementeren van optie A uit offerte •

Voor wijzigingen in het datamodel kan via de metadatabase en de dataman module gemakkelijk nieuwe lees- en schrijfroutines met een bijbehorende variantdatabase worden gegenereerd.

55


Schattingen BASGOED

De opdrachtgever (DVS) heeft inmiddels een keuze gemaakt, en wel voor de eerste optie (voortbouwen op het huidige prototype in Delphi).

56

HOOFDSTUK 6

Conclusies en aanbevelingen

Het specificatieonderzoek heeft geleid tot een model split model met een kostenfunctie die is gebaseerd op de kostenkentallen uit het vergelijkingskader modaliteiten (VKM). Hiermee is de aanpak consistent met de kostenfuncties in SMILE+. Per NSTR hoofdgroep zijn twee verschillende specificaties uitvoeriger getest. Een specificatie met herkomstdummies en bestemmingsdummies per vervoerwijze (serie 13), en een specificatie met vervoerwijzedummies op corridors (serie 14). In de eindconclusie worden deze specificaties beoordeeld op drie hoofdcriteria: •

elasticiteiten;

•

goodness-of-fit van de modal split module;

•

goodness-of-fit van de distributiemodule.

Van alle modellen zijn elasticiteiten en kruiselasticiteiten berekend en beoordeeld op basis van plausibele range uit de literatuur (zie paragraaf 3.5). De goodness-of-fit van de modellen is getoetst op basis van de standaard schattingsindicatoren en uit het resultaat van een toepassing van het model in de modal split module. Uit de toepassing is een maat van afwijking bepaald ten opzichte van de waarnemingen (RMSE) en zijn scatterplots en regressie coëfficiënten bepaald voor een aantal NSTR groepen (zie paragraaf 3.6). De goodness-of-fit van de distributiemodule is getoetst op basis van de R2 van de geschatte modellen (zie paragraaf 4.3). In onderstaand overzicht is het resultaat van de beoordeling van de specificaties op de drie hoofdcriteria indicatief samengevat. Tabel 6-1: Overzicht beoordeling model specificaties modal split modellen (+ = te prefereren specificatie)

Elasticiteiten

Serie 13

Serie 14

H en B dummies

Corridor dummies

+ (weg)

+ (spoor)

Goodness-of-fit van de modal split

++

Goodness-of-fit (variant 1)

+

van

de

distributiemodule

Op basis van de beoordeling op de drie hoofdcriteria is model serie 13 met vervoerwijze specifieke herkomst bestemmingsdummies duidelijk te prefereren als model voor implementatie in de modal split module van BASGOED. De goodness-of-fit van specificatie 13 is duidelijk beter dan specificatie 14 met corridor dummies. Wat betreft kostenelasticiteiten voor wegvervoer correspondeert het model met herkomst bestemmingsdummies beter met de literatuur. De kostenelasticiteit van spoorvervoer ligt in het


Schattingen BASGOED

model met vervoerwijze specifieke corridor dummies meer centraal in de waarschijnlijke marge. Model serie 13 scoort qua goodness-of-fit en stabiliteit ook beter in de distributiemodule; gegeven een negatief exponentiële vorm van de distributiefunctie met intrazonale weerstanden (de te prefereren variant voor de DM, zie hieronder). Verder geldt als inhoudelijk argument dat de herkomst en bestemmingsdummies iets zeggen over de niet geobserveerde aanbodkwaliteit binnen de regio’s. In deze gekozen specificatie voor modal split worden de tijdsafhankelijke- en kilometerafhankelijke transportkosten gebaseerd op de VKM-kostenkentallen, maar wordt ook een additionele tijdparameter geschat voor de kapitaalkosten van de goederen. In deze kapitaalkosten zijn de kosten voor interestkosten, ontwaarding en verzekeringskosten opgenomen. De nutsfunctie voor vervoerwijze v op relatie ij voor goederengroep g is hier als volgt gespecificeerd:


(2c)

d ijv en tijv zijn de afstand en tijd op de voedings- en verbindingslinks, De afstand en tijd, zijn via de gemiddelde afstandskosten Tv en het tarief per tijdseenheid Rv omgerekend naar kosten. CONTijg is de containerisatiegraad (verhoogde kans op spoor of binnenvaart als de goederen per container worden vervoerd) en MSCv de vervoerwijze specifieke constante (deze kan ook variëren per herkomst, per bestemming of per relatie, afhankelijk van serie 12, 13 of 14). De β’s zijn te schatten kosten- en tijdparameters. Uiteindelijk kiezen we voor serie 13 (de schattingsresultaten hiervan staan in Appendix 2). Ten aanzien van de distributie module zijn verschillende varianten van de distributiefunctie onderzocht. Hierin is gekeken naar verschillende aannames ten aanzien van het meenemen van parameters α , β en de grensweerstanden Gij en intrazonale weerstanden d. Op basis van de resultaten variant 5b, met een negatief exponentiële distributiefuncties en intrazonale weerstanden (maar geen grensweerstanden), naar voren als te prefereren vorm voor implementatie in de distributie module in BASGOED. In geval van NSTR groep 4 wordt teruggevallen op variant 1, zonder de intrazonale weerstanden. Concluderend kan gesteld worden dat voor de DM variant 5b (zie Tabel 4-3 en Appendix 3) wordt gekozen voor elke goederengroep behalve voor NSTR4. Hiervoor worden de intrazonale logsoms niet endogeen geschat en wordt variant 1 (zie Tabel 4-1 en Appendix 3) gebruikt. Zowel variant 1 als variant 5b bevat een simpele variant van de distributiefunctie (alleen alfa). Deze ziet er als volgt uit:

rij = exp(−α ⋅ cij ) De weerstand r tussen zones i en j wordt hier weergegeven als een functie van de logsums c uit de modal split module. Variant 5b heeft endogeen geschatte logsums d voor intrazonale relaties. Dit betekent dat in deze variant cij = di als i = j. De geschatte modellen zijn verder vertaald naar een functionele beschrijving voor de modal split en distributie module. In DP2 zullen deze rekenmodules verder worden uitgewerkt. In DP3 zullen deze modules worden geïmplementeerd in het BASGOED model systeem.

58

Appendix 1

Elasticiteiten per NSTR1 groep

Elasticiteiten model serie 12

Weg Spoor Bvrt

nstr0 kosten_wegkosten_spoor kosten_bvrttijd_weg tijd_spoor tijd_bvrt -0.5223 0.0084 0.1886 -0.2384 0.0028 0.2996 1.4769 -0.8615 0.0923 0.8615 -0.2769 0.1846 1.7633 0.0049 -0.6602 0.7973 0.0024 -1.0522

Weg Spoor Bvrt


Weg Spoor Bvrt


Weg Spoor Bvrt


Weg Spoor Bvrt


Weg Spoor Bvrt


Weg Spoor Bvrt


Weg Spoor Bvrt


Weg Spoor Bvrt


Weg Spoor Bvrt



Schattingen BASGOED


Weg Spoor Bvrt


Weg Spoor Bvrt


Weg Spoor Bvrt


Weg Spoor Bvrt


Weg Spoor Bvrt


Weg Spoor Bvrt


Weg Spoor Bvrt


Weg Spoor Bvrt


Weg Spoor Bvrt


Weg Spoor Bvrt


60

Schattingen BASGOED



Weg Spoor Bvrt


Weg Spoor Bvrt


Weg Spoor Bvrt


Weg Spoor Bvrt


Weg Spoor Bvrt


Weg Spoor Bvrt


Weg Spoor Bvrt


Weg Spoor Bvrt


Weg Spoor Bvrt


Weg Spoor Bvrt


61


Schattingen BASGOED

62

Appendix 2

schattingsresultaten MS serie 13

NSTR groep 0 t/m 5

NSTR Title Converged Observations Final log (L) D.O.F. Rho²(0) Rho²(c) Estimated Scaling Grns_Bvrt Grns_Spr CONT_Bvrt CONT_Spoor CONT_Weg de_bvrt39 de_bvrt38 de_bvrt37 de_bvrt36 de_bvrt35 de_bvrt34 de_bvrt33 de_bvrt32 de_bvrt31 de_bvrt30 de_bvrt29 de_bvrt28 de_bvrt27 de_bvrt26 de_bvrt25 de_bvrt24 de_bvrt23 de_bvrt22 de_bvrt21 de_bvrt20 de_bvrt19 de_bvrt18 de_bvrt17 de_bvrt16 de_bvrt15 de_bvrt14 de_bvrt13 de_bvrt12 de_bvrt11 de_bvrt10 de_bvrt9 de_bvrt8 de_bvrt7 de_bvrt6 de_bvrt5 de_bvrt4 de_bvrt3 de_bvrt2 de_bvrt1 or_bvrt39 or_bvrt38 or_bvrt37 or_bvrt36 or_bvrt35 or_bvrt34 or_bvrt33 or_bvrt32 or_bvrt31 or_bvrt30

0 Landbouwp.

1 Voedingsp.

TRUE

TRUE

763 -15172.2 89 0.462 0.279 23-Jul-10 1 -0.1528 (-3.9) 0 (*) 0 (*) 0 (*) 0 (*) -1.770 (-2.5) -0.5926 (-5.8) -1.363 (-18.0) -1.012 (-10.7) -0.2442 (-4.3) -1.473 (-5.9) -0.6380 (-12.5) -2.184 (-17.1) -0.3772 (-2.7) -1.182 (-8.3) -0.5150 (-12.3) -3.350 (-7.3) 0.9825 (0.8) -2.724 (-2.4) 0 (*) -2.364 (-1.8) -2.961 (-21.1) 1.240 (15.1) 0 (*) -1.926 (-2.2) -2.644 (-4.4) -2.853 (-4.1) -1.469 (-14.2) -1.964 (-12.5) -1.831 (-15.0) -0.2678 (-3.7) -0.7575 (-10.5) -0.4179 (-3.9) -0.1432 (-1.1) -0.6770 (-7.0) -0.02688 (-0.2) 0 (*) 0 (*) -0.7919 (-2.8) -0.3100 (-1.2) -0.2458 (-2.2) -1.281 (-6.9) 0.1564 (0.6) 2.153 (3.3) 0.7149 (6.6) 2.178 (10.8) -1.529 (-18.2) -3.838 (-7.6) -3.577 (-20.4) -4.023 (-10.2) -2.699 (-24.2) -2.226 (-28.0) -2.271 (-24.8) -1.235 (-11.1)

787 -23997.2 94 0.515 0.404 23-Jul-10 1 -0.6748 (-19.5) -5.901 (-13.9) -3.735 (-24.4) -11.69 (-5.3) 0 (*) -2.913 (-9.9) -2.022 (-24.0) -1.936 (-21.5) -1.607 (-20.7) -1.145 (-26.1) -2.852 (-16.0) -2.009 (-26.7) -0.7318 (-10.0) -1.958 (-9.5) -1.396 (-18.1) -0.1252 (-2.2) -1.579 (-15.7) -0.7701 (-3.3) 0 (*) -4.421 (-6.3) -0.1714 (-0.3) -1.746 (-35.3) 0.08209 (1.2) -3.852 (-7.1) -3.289 (-6.5) -5.216 (-6.5) 0 (*) -1.969 (-25.5) -1.814 (-17.8) -0.7820 (-9.9) -1.169 (-17.4) -1.455 (-24.4) -0.6729 (-7.4) -0.8001 (-8.1) -0.4083 (-5.7) -0.6518 (-6.6) 0 (*) -7.450 (-5.2) -1.522 (-14.4) -0.7377 (-4.4) -3.903 (-12.8) -2.575 (-14.7) 1.481 (6.8) -5.452 (-2.3) -2.948 (-5.7) -2.064 (-13.5) -3.886 (-14.6) -3.449 (-8.8) -2.139 (-30.3) -2.968 (-19.0) -2.614 (-28.0) -1.240 (-13.9) -0.9008 (-6.8) -2.233 (-21.0)

2 Vaste min. Br. Stof. TRUE 252 -11788.5 36 0.596 0.158 23-Jul-10 1 1.056 (9.1) 0 (*) -11.98 (-13.5) -143.5 (-10.6) 0 (*) -0.8227 (-2.9) 0.1952 (0.3) -3.039 (-18.3) -3.658 (-6.5) -7.804 (-12.1) 0 (*) 0.5316 (1.9) 0 (*) -0.7352 (-2.8) 10.76 (5.3) 0.2755 (1.6) -3.798 (-4.0) -2.277 (-4.1) 0 (*) -5.468 (-7.8) 0 (*) -0.7961 (-3.4) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -2.907 (-9.4) 2.953 (9.3) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -4.868 (-6.3) 0 (*) 0 (*) 0.3169 (0.3) 0 (*) 0 (*) 0.3159 (2.1) 0 (*) 0 (*) 0 (*) -1.969 (-3.1) 0 (*) 0.6417 (2.4) 1.094 (6.1) 0.4592 (3.0) 3.615 (2.3)

3 aardolie TRUE 633 -12806.2 60 0.714 0.513 23-Jul-10 1 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -2.344 (-8.7) 0.7589 (3.5) -4.664 (-16.0) -3.572 (-9.9) 0 (*) -1.240 (-7.7) -1.229 (-17.0) -1.153 (-14.9) 0 (*) -1.082 (-13.0) 0.2733 (5.6) -1.111 (-8.1) 0 (*) 0 (*) 0 (*) -0.08994 (-0.5) 0.4301 (7.3) -2.534 (-9.0) -1.858 (-5.8) 0 (*) 0 (*) 0 (*) -1.489 (-15.3) -1.340 (-6.4) 0.4124 (3.5) -3.717 (-8.0) 0 (*) -0.4097 (-1.7) 0 (*) 0.6544 (3.3) 2.534 (3.6) 0 (*) 0 (*) 0 (*) -0.6584 (-1.9) -3.597 (-24.8) 0 (*) 0.8360 (1.8) 0 (*) -3.062 (-6.4) -4.276 (-14.6) 0 (*) 0 (*) -5.978 (-14.0) -5.589 (-6.2) -2.653 (-25.5) -0.3903 (-5.4) 0.2390 (1.2) -2.411 (-21.9)

4 Ertsen TRUE 454 -7086.7 64 0.8 0.603 02-Aug-10 1 0 (*) 0 (*) -15.18 (-17.6) 0 (*) 0 (*) -3.052 (-14.2) 0.7663 (1.0) 0.1058 (0.2) -3.677 (-4.3) -0.7550 (-2.1) -0.7642 (-1.1) -1.126 (-7.2) -1.167 (-5.4) -1.343 (-7.0) -2.909 (-20.1) -1.291 (-10.9) -1.186 (-3.9) 0 (*) 0 (*) 0 (*) 0 (*) -1.608 (-14.0) 0 (*) -0.6146 (-0.3) -4.232 (-32.2) 0 (*) -3.831 (-10.2) -2.334 (-8.9) -0.7073 (-1.2) -5.198 (-22.8) -0.7130 (-0.8) -0.5610 (-0.8) 6.549 (0.0) 0 (*) -4.420 (-9.3) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 1.655 (2.3) -0.1875 (-0.3) 0 (*) -1.214 (-1.6) -2.314 (-15.6) -1.475 (-3.1) -4.906 (-8.4) -0.5756 (-2.1) -0.3302 (-0.8) 0.2679 (0.7) -1.277 (-9.6) -1.057 (-8.1) 1.597 (9.7) -1.552 (-9.2)


0

Schattingen BASGOED

NSTR Title

Landbouwp.

Voedingsp.

1 Vaste min. Br. Stof.

2 aardolie

or_bvrt29 or_bvrt28 or_bvrt27 or_bvrt26 or_bvrt25 or_bvrt24 or_bvrt23 or_bvrt22 or_bvrt21 or_bvrt20 or_bvrt19 or_bvrt18 or_bvrt17 or_bvrt16 or_bvrt15 or_bvrt14 or_bvrt13 or_bvrt12 or_bvrt11 or_bvrt10 or_bvrt9 or_bvrt8 or_bvrt7 or_bvrt6 or_bvrt5 or_bvrt4 or_bvrt3 or_bvrt2 or_bvrt1 de_spr39 de_spr38 de_spr37 de_spr36 de_spr35 de_spr34 de_spr33 de_spr32 de_spr31 de_spr30 de_spr29 de_spr28 de_spr27 de_spr26 de_spr25 de_spr24 de_spr23 de_spr22 de_spr21 de_spr20 de_spr19 de_spr18 de_spr17 de_spr16 de_spr15 de_spr14 de_spr13 de_spr12 de_spr11 de_spr10 de_spr9 de_spr8 de_spr7 de_spr6 de_spr5 de_spr4 de_spr3 de_spr2 de_spr1 or_spr39 or_spr38 or_spr37 or_spr36

-1.720 (-41.3) -2.142 (-8.6) -6.490 (-3.8) -4.171 (-3.3) -2.362 (-4.7) 0 (*) -0.9337 (-18.4) -1.546 (-9.2) 2.604 (3.4) -1.810 (-3.3) -1.460 (-4.0) -0.2062 (-2.3) -1.623 (-12.4) -2.915 (-8.8) -3.977 (-6.7) -3.046 (-9.3) -3.282 (-12.3) -3.030 (-5.7) -3.632 (-3.0) -2.819 (-5.8) -1.502 (-3.9) -2.021 (-5.1) -1.744 (-4.0) -0.5449 (-0.6) -0.5603 (-1.0) -2.072 (-16.7) -4.331 (-25.1) -0.7915 (-5.8) -0.8019 (-4.0) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -7.073 (-1.6) -7.713 (-2.6) -7.136 (-4.7) 0.3948 (1.8) -3.303 (-7.0) -4.096 (-13.1) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -3.764 (-7.7) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -5.981 (-2.1) 0 (*)

0.2871 (9.2) -3.281 (-15.7) -1.202 (-3.4) 0 (*) -4.746 (-3.8) -1.422 (-1.9) 0.8653 (26.0) -0.7859 (-11.3) 0 (*) -1.256 (-6.8) -4.431 (-6.6) -1.700 (-1.6) -3.502 (-22.0) -2.863 (-13.8) -1.653 (-15.4) -2.988 (-17.8) -4.329 (-21.3) -2.486 (-7.5) -3.169 (-9.6) -2.636 (-21.4) -3.112 (-11.9) -1.399 (-0.9) 0 (*) -3.442 (-6.8) -3.467 (-7.8) -3.140 (-10.1) -1.795 (-15.5) 2.054 (7.8) -3.201 (-4.3) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -1.259 (-3.1) 0 (*) 0 (*) 0 (*) 0.7311 (0.5) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -4.631 (-8.8) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -3.622 (-2.7) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -8.903 (-0.1) 0 (*) 0 (*) 0 (*) 0.9290 (1.7) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 1.328 (3.2)

2.702 (20.6) 0 (*) 0 (*) 0 (*) -1.165 (-1.4) 0 (*) 5.693 (12.5) 0 (*) 0 (*) 2.949 (7.3) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -0.1435 (-0.1) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*)

-0.1693 (-3.1) 0 (*) 0 (*) 0 (*) 0 (*) 0.4296 (1.9) -1.388 (-21.3) -5.849 (-5.0) -3.691 (-11.9) 0 (*) 0 (*) -2.395 (-22.4) -4.397 (-19.8) -3.597 (-12.2) -4.552 (-27.8) -3.976 (-1.1) -6.317 (-13.9) -4.323 (-10.0) 0 (*) -6.005 (-15.4) -4.649 (-8.8) 0 (*) 0 (*) -5.410 (-1.3) 0 (*) -6.193 (-18.7) -3.359 (-17.1) 0 (*) 0 (*) 0.7724 (1.7) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -2.642 (-2.0) 0 (*) 0 (*) -1.894 (-6.2) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -1.029 (-4.2) 0 (*) 0 (*) 0 (*)

64

3

4 Ertsen 0.7612 (7.3) -2.196 (-3.9) 0 (*) 0 (*) 0 (*) 0 (*) -1.835 (-11.4) -2.017 (-5.5) -2.095 (-5.6) -0.06718 (-0.4) -2.163 (-9.7) 0.1016 (0.1) -3.024 (-11.3) -0.2179 (-0.4) -0.2741 (-1.0) -0.5851 (-0.9) 5.193 (5.2) -2.172 (-12.4) 0 (*) -1.477 (-4.4) 0 (*) 0 (*) 0 (*) -1.238 (-1.7) -9.442 (-0.1) -0.8620 (-2.7) -3.941 (-14.5) 0 (*) -2.758 (-9.2) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -6.941 (-4.6) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -6.841 (-16.7) 0 (*) 0 (*) -3.630 (-8.5) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*)

Schattingen BASGOED


NSTR Title

Landbouwp.

0 Voedingsp.

1 Vaste min. Br. Stof.

2 aardolie

3

or_spr35 or_spr34 or_spr33 or_spr32 or_spr31 or_spr30 or_spr29 or_spr28 or_spr27 or_spr26 or_spr25 or_spr24 or_spr23 or_spr22 or_spr21 or_spr20 or_spr19 or_spr18 or_spr17 or_spr16 or_spr15 or_spr14 or_spr13 or_spr12 or_spr11 or_spr10 or_spr9 or_spr8 or_spr7 or_spr6 or_spr5 or_spr4 or_spr3 or_spr2 or_spr1 ASC_weg Time Cost

0 (*) 0 (*) 0 (*) -5.489 (-8.1) -2.052 (-7.7) 0 (*) -5.653 (-15.0) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -1.971 (-3.8) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -4.793 (-5.1) 0 (*) 0 (*) -2.26e-4 (-11.8) -0.07668 (-17.1)

0 (*) -0.3237 (-0.7) -1.518 (-3.8) 0.9557 (2.2) 0 (*) -2.614 (-0.3) -1.726 (-4.1) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -1.759 (-4.0) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 8.344 (0.0) 0 (*) -0.6542 (-0.9) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -9.096 (-4.6) 0 (*) 0 (*) 0 (*) -6.29e-4 (-34.5) -0.1813 (-41.2)

0 (*) 0 (*) 0 (*) 0 (*) -0.7445 (-2.1) 0 (*) 3.136 (25.4) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 5.657 (12.1) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -1.37e-4 (-11.0) -0.08794 (-11.2)

0 (*) 0 (*) -2.557 (-4.2) -0.2047 (-1.8) 0.5799 (2.3) 0 (*) -4.669 (-15.6) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -1.467 (-13.6) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 3.530 (7.3) 0 (*) 0 (*) 0 (*) 0 (*) -0.3116 (-59.7)

4 Ertsen 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -1.438 (-4.3) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -17.23 (-17.0) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -0.00103 (-18.9) -0.3897 (-35.5)

NSTR groep 6 t/m 9 5 Metalen. Halffabr. TRUE 483 -12723 85 0.285 0.278 23-Jul-10 1 0.5551 (14.1) 0 (*) -3.847 (-13.1) -9.017 (-7.6) 0 (*) -0.6723 (-8.7) -0.2636 (-1.7) -0.9449 (-5.1) -4.106 (-4.2) -3.161 (-5.4) -1.258 (-4.2) 0.1029 (1.5) -0.03842 (-0.3) -1.275 (-4.4) -0.7919 (-7.3) -0.4470 (-8.4) -2.288 (-2.9) 0 (*)

6 Ruwe min. TRUE 1394 -92417.9 104 0.359 0.334 23-Jul-10 1 -0.6794 (-31.8) 0 (*) 1.274 (4.7) -146.7 (-12.2) 0 (*) -2.557 (-59.0) -2.327 (-47.8) -4.397 (-81.0) -2.005 (-47.8) -1.489 (-54.9) -1.784 (-40.7) -1.047 (-35.1) -1.160 (-28.2) -0.8941 (-20.0) -0.6214 (-19.0) -0.9527 (-37.9) -1.015 (-24.6) -1.111 (-19.0)

7 Meststof

8 Chemie

TRUE

TRUE 450 580 -3284.5 -34840.2 60 95 0.453 0.382 0.44 0.278 23-Jul-10 23-Jul-10 1 1 -0.3888 (-5.4) 0.5897 (25.8) 0 (*) -2.520 (-37.8) 0 (*) 0 (*) -61.31 (-5.0) 2.483 (9.7) 0 (*) 0 (*) -2.806 (-5.7) -2.353 (-32.1) -3.397 (-10.1) -2.361 (-16.7) -2.994 (-10.3) -0.6813 (-5.5) -0.4193 (-0.9) -2.779 (-11.6) -1.041 (-4.9) 0.05432 (0.7) -0.3616 (-1.1) -4.379 (-11.5) -1.535 (-14.6) -1.049 (-29.5) -0.6702 (-4.4) -1.217 (-21.9) -0.9689 (-10.5) -0.3059 (-6.3) 1.322 (3.1) -1.774 (-22.1) -0.5336 (-3.1) -0.6433 (-28.0) -0.9005 (-2.2) -3.889 (-6.8) -1.845 (-1.6) 0 (*)

65

9 Overig TRUE 1120 -62179.5 133 0.568 0.432 23-Jul-10 1 -0.3777 (-18.2) -1.518 (-33.7) 5.858 (120.5) 5.683 (60.2) 0 (*) -2.250 (-32.6) -1.372 (-12.1) -2.485 (-30.7) -4.132 (-21.4) -1.721 (-28.5) -2.748 (-29.9) -2.273 (-39.6) -0.8238 (-18.3) -0.6388 (-9.9) -0.8496 (-14.4) -1.709 (-70.5) -2.722 (-23.4) -4.232 (-7.6)


5

Schattingen BASGOED

6

7

8

9

Metalen. Halffabr.

Ruwe min.

Meststof

Chemie

Overig

0 (*) 0 (*) 8.861 (1.8) -0.7229 (-5.5) 0.5071 (0.5) 3.271 (2.3) 0.4446 (3.5) 0 (*) -1.177 (-1.3) -2.570 (-7.0) -1.541 (-6.4) 0.09272 (0.7) -1.855 (-2.6) -0.5546 (-1.2) 1.353 (8.1) 0 (*) -3.253 (-4.4) 2.052 (3.6) 0 (*) 0 (*) -0.4355 (-1.2) 3.018 (1.3) -1.617 (-5.6) -2.098 (-3.0) -1.223 (-2.9) 0 (*) -2.853 (-12.4) -1.198 (-5.9) -2.627 (-6.1) -4.085 (-5.3) -2.554 (-4.5) -1.588 (-3.1) -1.471 (-18.6) -1.601 (-14.1) -2.015 (-10.7) -1.322 (-12.3) -0.4094 (-7.7) -3.913 (-2.8) 0 (*) 0 (*) 0 (*) 0 (*) -0.09348 (-1.0) 0 (*) 0 (*) 0.5647 (10.8) 0 (*) 0 (*) -3.239 (-5.6) 0.09150 (0.2) -2.613 (-8.5) 0 (*) -1.612 (-1.6) -1.569 (-3.2) 0 (*) -3.792 (-1.9) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -1.887 (-3.9) 0.3440 (0.6) -0.3059 (-1.4) 0.9199 (1.5) -1.444 (-5.5) 0 (*) 1.405 (5.0) 0 (*) 0 (*) 0 (*) -2.900 (-14.2)

-0.8779 (-14.5) -0.4212 (-7.6) -0.8109 (-11.1) -1.707 (-55.4) -1.007 (-9.2) -0.8726 (-9.9) -1.931 (-32.7) -0.3497 (-4.5) -0.5922 (-9.4) -1.415 (-44.8) -1.247 (-33.7) -1.222 (-42.7) -2.726 (-31.9) -0.7895 (-20.4) -1.994 (-47.5) -1.057 (-10.8) -1.653 (-38.3) -0.3996 (-4.9) -2.093 (-11.5) -1.918 (-8.3) -0.2569 (-3.9) -0.5078 (-9.1) -0.1026 (-2.0) -1.190 (-18.8) -1.647 (-23.0) -1.634 (-17.7) -2.174 (-48.7) 0.5016 (15.0) 0.09934 (2.9) -2.529 (-35.9) 0.09240 (3.6) -1.040 (-20.4) -2.506 (-64.1) 0.8973 (23.5) 0.2956 (8.2) -1.090 (-29.7) -0.7922 (-33.6) -1.495 (-32.1) -3.229 (-5.3) -6.614 (-10.3) -2.591 (-9.6) -1.964 (-15.5) -0.9475 (-33.9) -1.765 (-5.1) -3.982 (-5.0) 1.390 (32.0) -2.021 (-7.4) -1.822 (-23.7) -1.888 (-49.2) -1.341 (-37.8) -1.871 (-60.0) -2.370 (-21.5) -2.054 (-39.4) -1.389 (-36.8) -1.799 (-4.8) -1.502 (-32.6) -2.211 (-6.5) -3.397 (-10.9) -1.269 (-10.4) -4.476 (-22.7) -1.860 (-27.1) -1.146 (-23.7) -3.174 (-31.3) 0.1094 (1.7) -1.574 (-8.8) -1.253 (-3.4) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -2.479 (-17.4)

0 (*) 0 (*) 0 (*) -1.599 (-7.0) 0 (*) 0 (*) 0.2264 (0.5) -2.246 (-2.1) -2.183 (-5.9) -0.2425 (-1.3) 0 (*) -0.7104 (-1.0) -0.3136 (-1.2) -3.714 (-3.9) -0.9089 (-2.4) -0.1155 (-0.2) -0.5693 (-2.2) -2.544 (-10.7) 0 (*) 0 (*) 0.1382 (0.4) 0 (*) -2.159 (-10.2) -1.186 (-3.2) 0 (*) -0.5143 (-0.6) 0.7939 (8.4) -0.4638 (-1.1) 2.824 (3.5) -2.782 (-6.6) -4.314 (-16.0) -4.233 (-2.6) -3.432 (-15.4) -0.6676 (-3.0) -0.3274 (-3.5) 0.01408 (0.1) -1.493 (-9.2) 0 (*) 0.6732 (0.8) 0 (*) 0 (*) 0 (*) 0.4227 (2.9) 0 (*) 0 (*) 1.165 (7.1) 0 (*) 0 (*) -2.974 (-9.3) 0 (*) 0 (*) 0 (*) -2.873 (-3.4) -1.359 (-2.0) 0 (*) -0.8244 (-2.0) -2.825 (-3.1) 0 (*) 0 (*) 1.034 (1.4) 0 (*) 1.919 (1.5) -2.482 (-4.2) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*)

-2.348 (-1.8) -2.308 (-2.8) 0 (*) 0.2544 (6.0) -1.537 (-3.5) -2.490 (-2.8) -0.2421 (-1.0) 0 (*) -4.506 (-6.7) -3.919 (-5.6) -2.925 (-10.9) -0.9725 (-10.7) -0.4656 (-1.8) -2.605 (-5.7) -3.735 (-8.0) 0 (*) -2.517 (-10.8) -1.676 (-1.1) 0 (*) 0 (*) -2.727 (-2.7) 1.280 (1.1) -3.177 (-5.7) -2.628 (-6.4) 0.01184 (0.1) -1.135 (-4.3) -0.9828 (-20.4) -1.943 (-12.7) -4.045 (-7.4) -1.939 (-8.4) -4.021 (-14.2) -5.741 (-5.5) -1.056 (-30.0) -0.5101 (-11.2) -0.5728 (-15.5) -2.261 (-29.0) -0.6950 (-31.4) -3.570 (-6.9) 0 (*) 0 (*) 0 (*) 0 (*) -2.132 (-41.1) -2.841 (-4.2) 0 (*) -2.479 (-7.3) 0 (*) -2.829 (-3.4) -3.663 (-9.4) -2.862 (-5.5) -2.494 (-15.8) 0 (*) -4.228 (-3.7) -2.099 (-10.2) 0 (*) -3.106 (-4.5) -3.089 (-1.1) 0 (*) 0 (*) 0 (*) 4.405 (0.8) -2.045 (-4.3) -2.341 (-6.1) -0.4927 (-7.1) -0.1878 (-1.3) -0.1546 (-1.7) 0 (*) 0 (*) -1.067 (-0.7) 0 (*) 0 (*) -2.801 (-18.3)

-1.167 (-6.9) -3.584 (-6.4) -3.506 (-7.0) -1.751 (-39.0) -0.9266 (-9.8) -2.424 (-5.1) -0.9833 (-7.6) -1.922 (-12.6) -1.512 (-10.4) -1.254 (-22.7) -2.496 (-20.0) -1.805 (-24.3) -3.052 (-12.7) -2.347 (-20.7) -2.171 (-22.8) -2.607 (-5.9) -1.763 (-14.3) -1.593 (-9.9) -4.905 (-2.0) 0.6168 (0.5) -0.8383 (-6.4) 0.5483 (2.9) -2.941 (-20.1) -1.906 (-14.7) 0.5017 (4.8) -1.101 (-5.4) -2.360 (-39.9) -2.310 (-14.8) -3.826 (-40.4) -2.354 (-26.0) -2.109 (-36.2) -2.311 (-23.5) -2.661 (-46.2) -0.4178 (-10.0) -1.330 (-22.9) -1.743 (-29.1) -2.492 (-97.5) -2.889 (-16.9) -1.858 (-3.4) -0.2165 (-2.4) -3.167 (-6.1) -1.505 (-5.9) -1.858 (-39.3) -2.174 (-22.9) -2.900 (-4.4) -1.124 (-11.2) -1.517 (-8.9) -2.014 (-12.9) -0.7790 (-15.5) -1.997 (-16.3) -1.959 (-26.8) -3.540 (-10.1) -2.422 (-27.6) -1.411 (-20.0) -0.7832 (-2.0) -1.515 (-16.4) -0.3188 (-1.8) -3.931 (-3.9) 0 (*) -3.138 (-14.3) -5.089 (-10.9) -3.545 (-18.7) -2.094 (-23.0) -1.086 (-8.9) -2.713 (-19.6) -3.662 (-17.4) 0 (*) -0.9572 (-11.7) -2.047 (-16.9) -8.931 (-1.9) -2.495 (-20.1) -2.956 (-23.2)

66

Schattingen BASGOED


5

6

7

8

9

Metalen. Halffabr.

Ruwe min.

Meststof

Chemie

Overig

-2.184 (-5.1) 0.3492 (0.6) -3.099 (-9.2) -5.235 (-26.4) 0 (*) 0 (*) 0 (*) 0 (*) 5.080 (2.2) -3.699 (-2.5) 0 (*) 0 (*) 0.7561 (5.1) 0 (*) 0 (*) -3.140 (-5.5) 0 (*) -0.6420 (-3.0) 0 (*) 0 (*) 2.174 (10.1) 0 (*) -4.821 (-2.1) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -4.728 (-2.5) 0 (*) 0 (*) -3.059 (-10.2) 0 (*) -5.775 (-1.1) -2.439 (-14.0) -1.187 (-10.2) 0 (*) -3.935 (-5.4) -1.428 (-13.8) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -3.022 (-10.2) 0 (*) 0 (*) 0.6385 (9.3) 0 (*) 0 (*) 0 (*) 0 (*) -4.475 (-1.0) 0 (*) 0 (*) -5.446 (-0.3) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*)

0 (*) 0.3055 (0.6) -4.454 (-6.9) -2.493 (-5.5) 0 (*) 0 (*) 0 (*) 0 (*) 1.473 (1.5) -7.711 (-14.0) 0 (*) 0 (*) -0.6790 (-9.1) 0 (*) 0 (*) -4.444 (-6.9) 0 (*) -4.008 (-33.3) 0 (*) 0 (*) -2.632 (-11.2) 0 (*) -1.684 (-2.2) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -5.372 (-28.9) 0.4518 (1.8) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -2.640 (-0.4) -3.491 (-27.7) 0 (*) -6.945 (-1.8) 0 (*) -5.301 (-14.0) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0.3676 (0.3) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -6.661 (-4.2) -0.5911 (-0.1) 0 (*)

0 (*) -1.857 (-1.3) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 3.338 (1.9) -3.211 (-7.5) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -3.866 (-17.8) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -2.334 (-2.9) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*)

-2.402 (-10.8) -1.216 (-11.7) -0.6395 (-6.7) -4.666 (-44.5) 0 (*) 2.019 (3.2) 0 (*) 0 (*) 0 (*) -4.588 (-2.4) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -4.753 (-3.9) 0 (*) 0 (*) 0 (*) -2.681 (-7.9) 0 (*) 0 (*) 0 (*) 1.430 (7.3) 0 (*) 0 (*) 0 (*) 0 (*) -5.052 (-3.9) -1.501 (-6.4) 0 (*) 0.4654 (6.6) 0 (*) 0 (*) -1.071 (-1.1) 0 (*) -2.795 (-6.3) -0.7164 (-7.0) 0.3612 (3.9) -0.2575 (-3.1) -2.405 (-12.2) -0.5790 (-10.8) -3.818 (-4.3) 0 (*) 0 (*) 0 (*) 0 (*) -1.990 (-10.0) 0 (*) 0 (*) 0.3018 (0.5) 0 (*) 0 (*) 0 (*) -3.759 (-1.5) 0 (*) 0 (*) 0 (*) 0.9826 (3.7) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -0.3294 (-2.5) 0.7298 (0.6) 0 (*)

-0.8786 (-9.6) -2.971 (-4.7) -4.020 (-7.5) -3.094 (-55.8) -0.3061 (-1.3) 0 (*) 0 (*) 0 (*) -0.1540 (-0.2) -3.844 (-18.5) 0 (*) 0 (*) -5.073 (-0.5) 0 (*) 0 (*) -1.128 (-8.6) -1.260 (-1.2) -4.647 (-3.0) 0 (*) -4.344 (-5.5) -4.492 (-2.1) 0 (*) -6.380 (-2.7) 0 (*) 0 (*) 3.633 (15.3) 0 (*) 0 (*) 0 (*) 0.5478 (2.1) -3.395 (-6.4) -3.302 (-5.3) -1.705 (-19.1) 0 (*) -0.6818 (-7.5) -2.577 (-16.1) 0 (*) -3.498 (-4.6) -2.996 (-26.8) -3.286 (-12.6) -1.165 (-8.7) -5.358 (-4.6) -3.178 (-56.4) -3.237 (-2.5) 0 (*) -6.362 (-0.5) -0.9655 (-1.6) -0.3377 (-0.8) -2.726 (-16.4) 0 (*) 0 (*) 0 (*) 0 (*) 0 (*) -3.834 (-7.1) -3.669 (-4.4) -2.898 (-3.5) 0 (*) -2.404 (-10.4) -5.308 (-4.4) 0 (*) -3.249 (-11.6) 0 (*) -0.2869 (-1.4) 1.881 (4.4) 0 (*) 0 (*) 0.4829 (2.8) -3.650 (-4.2) -5.806 (-0.6) 0.5219 (3.8) 0 (*)

67


5

Schattingen BASGOED

6

7

8

9

Metalen. Halffabr.

Ruwe min.

Meststof

Chemie

Overig

-3.77e-4 (-29.9) -0.05938 (-22.0)

-4.40e-4 (-36.8) -0.2556 (-98.6)

-6.07e-4 (-13.8) -0.2262 (-20.2)

-3.31e-4 (-39.0) -0.09693 (-54.2)

-3.57e-4 (-60.3) -0.03220 (-30.1)

68

Schattingen BASGOED


Appendix 3

schattingsresultaten DM

MS model 13, DM variant 1, incl. intrazonaal 0

1

2

3

4

5

6

7

8

9

R²

NSTR

0.369

0.482

0.975

0.827

0.975

0.587

0.680

0.439

0.468

0.347

alpha

1.540

1.279

2.845

-0.643

1.096

0.484

3.059

1.769

1.999

0.535

MS model 13, DM variant 2, incl. intrazonaal NSTR

0

1

2

3

4

5

6

7

8

9

R²

0.338

0.438

0.961

0.838

0.979

0.590

0.710

0.419

0.464

0.339

beta

3.480

7.125

8.229

-2.143

12.085

5.623

7.183

4.861

12.109

1.011

MS model 13, DM variant 3, incl. intrazonaal NSTR R²

0

1

2

3

4

5

6

7

8

9

0.373

0.483

0.975

0.831

0.977

0.588

0.708

0.422

0.470

0.353

alpha

2.214

1.328

2.860

-0.540

0.762

0.447

0.669

2.111

1.855

1.325

beta

-1.731

-0.346

0.666

-0.375

3.200

0.324

5.475

0.667

0.562

-2.391

MS model 13, DM variant 4a, incl. intrazonaal NSTR R²

0

1

2

3

4

5

6

7

8

9

0.373

0.483

0.977

0.832

0.977

0.606

0.753

0.422

0.471

0.380

alpha

2.227

1.328

2.837

-0.490

0.762

0.302

0.992

2.111

1.859

1.101

beta

-1.769

-0.346

6.335

-0.338

3.200

9.162

5.624

0.667

0.666

-6.291

0.000

0.000

10.126

0.805

0.000

44.264

83.369

0.000

10.262

24.875

5

6

7

8

9 0.385

G

MS model 13, DM variant 4b, incl. intrazonaal NSTR

0

1

2

3

4

R²

0.378

0.483

alpha

2.316

1.328

0.978

0.858

0.977

0.679

0.751

0.422

0.471

2.256

-0.426

0.760

-1.227

0.939

2.111

1.842

beta

-2.083

1.326

-0.346

12.745

-0.718

3.241

26.223

5.600

0.667

0.667

-6.323

G1 G2

0.000

0.004

0.000

0.000

0.000

0.000

0.001

0.000

0.001

0.107

0.000

0.000

0.000

0.000

0.000

0.121

0.092

0.006

0.000

0.033

G3

9.210

0.000

0.035

0.000

0.000

0.000

4.397

0.170

0.001

0.948

G4

0.000

0.000

0.000

0.000

0.007

2.491

25.097

0.000

0.025

0.000

G5

0.000

0.000

0.921

0.000

0.000

21.990

11.109

0.032

0.070

10.899

G6

0.000

0.000

0.000

0.214

0.032

3.906

18.384

0.000

4.703

2.637

G7

0.000

0.000

0.000

0.000

0.000

0.000

1.718

0.005

0.000

0.238

G8

0.000

0.000

0.000

0.000

0.610

26.481

23.144

0.173

4.350

11.599

G9

0.000

0.000

0.000

0.002

0.000

0.004

0.001

0.000

0.000

0.038

G10

0.000

0.000

0.000

0.000

0.000

0.000

0.084

0.010

0.003

0.614

G11

0.000

0.000

0.000

4.104

0.000

4.954

2.199

0.232

0.078

20.501

G12

0.000

0.000

8.566

0.413

0.000

0.000

46.855

0.000

0.087

0.000

G13

0.000

0.000

0.041

0.050

0.006

0.000

43.999

0.084

0.028

0.000

G14

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.099

0.011

0.000

G15

0.000

0.000

0.000

0.000

0.000

2.365

1.517

0.011

0.001

0.000

69


Schattingen BASGOED

G16

0.011

0.002

0.000

0.000

0.000

1.552

0.004

0.000

0.002

0.583

G17

0.000

0.006

0.000

0.000

0.000

0.000

0.017

0.000

0.002

0.254

G18

0.000

0.011

0.000

0.000

0.000

0.000

0.014

0.000

0.039

0.905

G19

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.003

G20

0.000

0.008

0.000

0.000

0.000

0.066

0.000

0.000

0.004

0.528

G21

0.000

0.001

0.000

0.007

0.000

1.326

0.000

0.001

0.001

0.420

G22

0.000

0.007

0.000

0.000

0.000

0.000

0.013

0.009

0.000

0.081

G23

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.341

G24

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.001

0.208

G25

0.000

0.002

0.000

0.000

0.000

0.331

0.000

0.000

0.001

1.076

G26

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.020

G27

0.001

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.016

G28

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.001

0.020

G29

0.105

0.000

0.000

0.000

0.000

8.351

0.000

0.000

0.063

0.000

MS model 13, DM variant 4c, incl. intrazonaal NSTR

0

1

2

3

4

5

6

7

8

9

R²

0.369

0.482

0.975

0.827

0.975

0.587

0.680

0.439

0.468

0.347

alpha

1.540

1.279

2.845

-0.643

1.096

0.484

3.059

1.769

1.999

0.535

G1

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G2

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G3

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G4

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G5

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G6

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G7

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G8

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G9

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G10

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G11

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G12

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G13

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G14

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G15

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G16

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G17

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G18

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G19

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G20

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G21

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G22

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G23

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G24

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G25

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G26

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G27

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G28

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G29

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

70

Schattingen BASGOED


MS model 13, DM variant 4d, incl. intrazonaal NSTR R² beta G1 G2 G3 G4 G5 G6 G7 G8 G9 G10 G11 G12 G13 G14 G15 G16 G17 G18 G19 G20 G21 G22 G23 G24 G25 G26 G27 G28 G29

0 0.338 3.489 0.000 0.045 0.000 0.000 0.000 0.000 0.000 0.356 0.003 0.045 0.000 0.000 0.009 0.013 0.081 0.005 0.033 0.257 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.008 0.000 0.000 0.000

1 0.438 7.126 0.000 0.005 0.124 0.000 0.166 0.242 0.023 0.255 0.003 0.003 0.000 0.000 0.000 0.000 0.051 0.052 0.034 0.059 0.002 0.000 0.000 0.000 0.000 0.031 0.003 0.000 0.000 0.000 0.124

2 0.973 18.436 0.000 0.002 0.001 0.201 1.586 0.052 0.000 0.000 0.000 0.002 0.000 2.257 0.380 0.014 0.284 0.000 0.000 0.000 0.000 0.000 0.002 0.015 0.000 0.000 0.000 0.000 0.000 0.000 0.000

3 0.858 -2.276 0.000 0.000 0.000 0.000 0.000 0.185 0.000 0.000 0.000 0.000 3.999 0.481 0.061 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.013 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.265

4 0.982 14.885 0.000 0.000 0.000 0.008 0.000 0.047 0.000 11.449 0.000 0.000 0.076 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000

5 0.661 14.086 0.000 0.399 0.000 2.696 82.484 5.903 0.202 77.383 0.022 0.000 5.706 0.000 0.000 0.000 2.913 1.844 0.000 0.000 0.000 0.229 1.783 0.186 0.000 0.006 0.429 0.000 0.000 0.000 11.095

6 0.749 7.742 0.001 0.000 0.018 6.881 0.513 0.748 0.011 0.500 0.002 0.006 0.106 0.536 5.734 0.003 0.037 0.002 0.004 0.003 0.000 0.000 0.003 0.006 0.000 0.000 0.000 0.000 0.000 0.000 0.001

7 0.445 7.179 0.000 0.212 2.825 1.206 1.394 3.266 0.481 8.685 0.017 0.403 1.036 0.000 4.795 1.739 2.781 0.574 0.000 0.000 0.000 0.000 0.080 0.662 0.000 0.000 0.000 0.000 0.000 0.000 0.000

8 0.468 13.042 0.007 0.000 0.000 3.217 0.006 3.219 0.000 5.058 0.000 0.000 0.755 1.729 4.815 0.000 0.000 0.006 0.017 0.167 0.000 0.040 0.001 0.000 0.000 0.012 0.002 0.004 0.000 0.011 4.554

9 0.339 1.012 0.000 0.000 0.000 0.030 0.000 0.009 0.000 0.000 0.000 0.000 0.000 0.042 0.038 0.018 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.021 0.000 0.000 0.000 0.000 0.000 0.000 0.000

5 0.593 0.235 7.628

6 0.835 1.023 0.000

7 0.576 0.558 0.000

8 0.533 1.357 5.161

9 0.753 1.484 2.406

5 0.637 0.539 8.971 8.631 8.813 8.716 8.987 8.790 8.993 8.860

6 0.855 1.462 0.007 2.133 0.036 0.000 3.056 0.501 3.028 1.542

7 0.670 0.969 0.000 3.155 2.035 0.254 3.493 3.437 2.290 0.566

8 0.817 2.148 4.621 4.575 4.002 3.580 4.632 4.470 4.624 4.628

9 0.815 1.016 3.008 4.089 1.051 3.879 4.069 3.997 4.083 3.772

MS model 13, DM variant 5a, incl. intrazonaal NSTR R² alpha d

0 0.496 0.528 0.000

1 0.570 0.420 4.797

2 0.978 2.835 14.495

3 0.931 1.635 2.951

4 0.978 1.023 34.260

MS model 13, DM variant 5b, incl. intrazonaal NSTR R² alpha d1 d2 d3 d4 d5 d6 d7 d8

0 0.680 1.602 0.000 2.878 2.866 0.775 2.789 2.858 2.309 2.382

1 0.771 0.912 8.167 8.569 5.884 6.120 8.349 5.528 7.530 8.561

2 0.978 2.752 15.439 15.439 15.456 15.439 15.439 15.439 15.439 12.259

3 0.945 1.363 3.105 3.101 3.072 3.049 3.104 3.105 3.105 3.143

4 0.984 0.887 9.036 9.041 9.029 9.039 9.039 9.039 9.064 9.038

71


d9 d10 d11 d12 d13 d14 d15 d16 d17 d18 d19 d20 d21 d22 d23 d24 d25 d26 d27 d28 d29 d30 d31 d32 d33 d34 d35 d36 d37 d38 d39 d40

2.855 2.354 2.862 1.769 1.445 1.817 2.704 1.132 1.639 0.853 2.857 2.883 2.886 2.836 0.948 2.886 2.621 2.852 0.600 2.714 2.929 2.775 2.843 2.656 2.018 2.468 1.339 1.032 1.655 2.862 0.000 2.507

7.764 5.327 8.168 5.530 6.181 5.372 6.834 7.968 8.552 8.422 8.558 7.971 8.566 5.109 6.537 8.566 8.413 8.473 8.567 6.093 8.269 5.785 8.496 8.063 6.596 6.853 5.968 5.629 4.020 6.827 4.562 8.324

15.439 15.439 15.439 15.442 15.439 15.439 15.439 9.772 15.433 15.439 15.439 15.439 15.439 15.439 15.439 15.439 19.392 15.439 15.439 15.482 15.457 15.439 15.432 15.439 15.467 15.439 10.833 15.471 14.136 15.439 14.429 15.439

Schattingen BASGOED

3.104 2.195 3.097 1.112 0.184 2.843 1.970 3.044 2.189 1.823 3.105 3.105 3.105 3.105 3.289 3.105 3.103 3.101 3.105 3.091 3.283 2.852 3.054 0.753 1.006 2.710 0.000 3.061 3.097 3.010 3.039 2.908

9.042 9.031 9.041 9.041 9.041 9.042 9.016 9.040 9.001 9.048 9.031 0.000 9.041 9.044 0.000 9.041 9.041 9.041 9.041 9.041 34.260 9.019 9.023 9.030 9.004 9.041 9.039 0.000 9.041 9.039 9.022 9.044

8.971 8.703 8.791 3.206 8.857 7.925 7.165 8.926 4.687 8.709 8.361 8.589 9.001 8.993 6.083 9.002 8.949 8.886 8.668 8.989 9.285 6.072 0.462 3.633 8.549 8.870 8.394 5.428 8.565 8.685 6.160 8.846

3.010 0.929 2.439 1.079 0.307 0.537 0.896 0.329 0.427 0.132 3.155 3.120 3.168 3.166 0.897 3.140 0.000 1.439 0.909 2.067 1.315 0.987 0.436 0.140 0.965 0.440 0.539 0.423 0.850 0.962 0.797 2.557

3.443 1.888 3.545 1.518 0.366 2.749 3.201 3.011 2.211 1.832 3.617 3.493 3.507 3.507 1.856 3.497 3.371 3.490 3.498 2.573 1.097 3.486 5.023 2.042 1.772 2.047 1.585 1.169 0.000 3.253 0.782 3.236

4.264 2.575 4.591 2.527 4.232 4.520 2.564 3.247 2.643 4.366 4.618 4.576 4.628 4.559 3.465 4.626 4.201 2.172 4.629 3.039 4.579 3.743 3.499 3.134 3.773 2.938 4.441 3.171 4.611 3.205 3.208 4.588

4.079 2.962 4.054 1.116 1.343 0.153 1.611 3.786 1.118 3.885 3.953 4.027 4.053 4.058 0.791 4.060 3.875 3.026 3.692 3.532 2.699 1.103 3.817 3.208 2.105 1.704 0.021 1.196 0.689 3.999 0.829 4.031

5 0.635 0.427 0.542 0.000 0.000 0.000 0.000 5.984 0.208 0.000 6.622 0.000 0.000 0.215 0.000 0.000 0.000 0.014 0.023 0.000 0.000 0.000 0.000

6 0.877 0.809 0.238 0.001 0.003 0.044 0.144 0.226 0.236 0.019 0.273 0.000 0.017 0.019 0.000 0.753 0.005 0.025 0.001 0.001 0.001 0.000 0.001

7 0.707 0.719 0.241 0.000 0.007 0.892 0.068 0.088 0.000 0.008 9.189 0.000 0.013 0.000 1.551 2.740 0.074 0.885 0.011 0.000 0.000 0.000 0.000

8 0.861 1.633 3.542 0.002 0.000 0.327 1.969 8.465 11.141 0.053 10.618 0.000 0.097 6.086 4.431 0.000 0.064 0.000 0.000 0.000 0.386 0.000 0.002

9 0.832 1.583 -1.112 0.021 0.017 0.153 0.000 2.577 0.000 0.000 8.906 0.005 0.038 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

MS model 13, DM variant 6, incl. intrazonaal NSTR R² alpha beta G1 G2 G3 G4 G5 G6 G7 G8 G9 G10 G11 G12 G13 G14 G15 G16 G17 G18 G19 G20

0 0.689 1.505 0.090 0.007 0.008 0.000 0.048 0.027 0.041 0.014 0.124 0.005 0.042 0.060 0.062 0.062 0.049 0.078 0.002 0.004 0.006 0.000 0.020

1 0.790 0.953 -0.081 0.001 0.000 0.000 0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.000 0.000 0.000 0.002 0.003 0.000 0.000

2 0.981 1.800 12.450 0.000 0.000 0.067 0.017 1.273 0.000 0.000 0.000 0.000 0.000 0.000 4.758 0.125 0.000 0.000 0.000 0.000 0.000 0.000 0.000

3 0.946 1.888 -1.361 0.000 0.000 0.000 0.000 0.000 0.791 0.000 0.000 0.001 0.000 9.699 2.504 0.694 0.012 0.000 0.007 0.000 0.000 0.000 0.000

4 0.992 1.748 4.533 0.000 0.002 0.000 0.008 0.009 0.012 0.000 2.082 0.000 0.000 0.921 0.000 0.002 0.000 0.000 0.000 0.000 0.000 0.000 0.000

72

Schattingen BASGOED


G21 G22 G23 G24 G25 G26 G27 G28 G29 d1 d2 d3 d4 d5 d6 d7 d8 d9 d10 d11 d12 d13 d14 d15 d16 d17 d18 d19 d20 d21 d22 d23 d24 d25 d26 d27 d28 d29 d30 d31 d32 d33 d34 d35 d36 d37 d38 d39 d40

0.001 0.005 0.000 0.001 0.001 0.004 0.000 0.000 0.000 0.059 2.376 2.733 0.758 2.252 2.348 1.856 1.806 2.344 1.739 2.353 1.508 1.434 1.621 2.174 1.091 1.476 0.846 2.351 2.382 2.385 2.316 0.930 2.385 2.062 2.345 0.566 2.162 2.880 2.245 2.449 2.625 1.993 1.972 1.319 1.012 1.623 2.349 0.001 1.888

0.001 0.002 0.000 0.000 0.001 0.000 0.000 0.000 0.000 6.206 6.643 6.121 6.326 6.276 5.704 6.920 6.629 6.804 5.477 5.440 5.684 6.286 5.512 5.314 6.549 6.594 6.166 6.755 3.753 6.632 5.065 6.614 6.634 6.683 6.435 6.657 5.520 8.264 5.861 6.428 6.036 6.682 6.892 6.104 5.729 4.186 6.032 4.736 6.373

0.001 0.003 0.000 0.000 0.000 0.000 0.000 0.000 0.000 16.077 16.077 16.094 16.077 16.077 16.077 16.077 15.961 16.077 16.077 16.077 16.090 16.077 16.077 16.078 15.903 16.074 16.077 16.077 16.077 16.077 16.077 16.077 16.077 19.268 16.077 16.077 16.210 16.076 16.077 15.992 16.077 16.085 16.077 2.312 16.134 13.360 16.077 14.369 16.077

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.703 0.705 0.728 0.710 0.703 0.703 0.703 0.593 0.705 0.732 0.502 0.723 0.730 0.702 0.731 0.703 0.715 0.724 0.703 0.703 0.703 0.703 0.028 0.703 0.711 0.710 0.703 0.714 0.039 0.706 0.730 1.188 0.726 0.708 0.739 0.705 0.703 0.708 0.702 0.793

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 5.662 5.688 1.693 5.493 5.688 5.689 5.687 5.670 5.684 0.320 5.687 5.689 5.688 5.681 1.291 5.862 0.909 5.690 4.387 5.694 5.688 5.801 3.043 5.688 5.688 5.643 5.688 5.688 10.360 1.553 6.172 2.724 2.213 5.593 5.732 5.886 5.696 5.696 5.663 5.705

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.589 9.031 8.682 8.860 8.778 9.044 8.833 9.055 8.905 9.027 8.768 8.826 2.971 8.907 7.933 7.260 8.978 4.117 8.749 8.341 8.832 9.064 9.055 6.254 9.066 9.003 8.930 8.675 9.052 9.175 5.989 6.537 3.151 8.352 8.929 8.417 4.868 8.544 8.717 5.591 8.887

0.001 0.002 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.028 0.211 0.003 0.001 0.177 0.086 0.160 0.130 0.140 0.147 0.001 0.205 0.012 0.047 0.179 0.036 0.023 0.001 0.166 0.161 0.171 0.169 0.165 0.155 0.001 0.179 0.040 0.012 0.608 0.172 0.107 0.003 0.230 0.001 0.054 0.020 0.322 0.242 0.155 0.001

0.003 0.027 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 1.275 1.552 0.014 1.420 1.314 1.106 0.750 1.437 1.516 1.620 1.251 0.218 1.241 1.362 0.768 0.676 1.376 1.612 1.431 1.438 1.438 1.068 1.426 1.308 1.403 1.452 0.311 0.803 1.389 5.081 0.953 1.404 1.771 1.382 1.092 0.001 1.325 0.424 1.265

0.000 0.000 0.000 0.007 0.004 0.000 0.000 0.010 2.710 3.812 3.847 4.141 3.212 3.865 3.249 3.828 3.846 3.316 2.899 3.558 2.829 3.122 2.135 2.821 3.400 2.986 3.772 3.786 3.604 3.849 3.500 3.705 3.832 2.423 2.568 3.856 3.103 4.605 3.814 3.657 3.337 3.941 3.133 3.250 3.361 3.733 3.255 3.417 3.491

0.019 0.000 0.021 0.024 0.000 0.003 0.012 0.002 0.039 0.795 3.896 1.042 3.479 3.859 3.706 3.900 3.128 3.887 1.980 3.831 1.631 1.906 0.760 2.060 3.232 1.561 3.460 3.613 3.770 3.830 3.832 1.172 3.842 3.412 2.113 2.963 2.471 2.703 1.577 3.217 1.466 2.603 1.739 0.710 1.754 0.828 3.732 1.012 3.807

MS model 13, DM variant 1, excl. intrazonaal 0

1

2

3

4

5

6

7

8

9

R²

0.357

0.554

0.976

0.951

0.984

0.594

0.374

0.292

0.711

0.692

alpha

1.046

0.728

2.860

1.016

1.203

0.338

0.988

0.532

1.144

0.852

NSTR

MS model 13, DM variant 2, excl. intrazonaal

73


NSTR

Schattingen BASGOED

0

1

2

3

4

5

6

7

8

9

R²

0.336

0.521

0.962

beta

2.285

4.189

6.519

0.949

0.987

0.596

0.392

0.337

0.713

0.651

1.728

13.469

4.643

2.475

2.383

6.060

1.734

MS model 13, DM variant 3, excl. intrazonaal NSTR R²

0

1

2

3

4

5

6

7

8

9

0.360

0.554

0.976

0.951

0.985

0.595

0.401

0.374

0.713

0.704

alpha

1.410

0.733

2.893

0.831

0.747

0.289

0.393

-1.155

0.326

1.490

beta

-0.975

-0.044

0.663

0.385

3.876

0.424

1.727

6.231

4.331

-1.925

MS model 13, DM variant 4a, excl. intrazonaal NSTR R²

0

1

2

3

4

5

6

7

8

9

0.360

0.554

0.978

0.951

0.985

0.617

0.404

0.293

0.713

0.729

alpha

1.377

0.733

2.871

0.831

0.747

0.295

0.700

0.494

0.326

1.072

beta

-0.929

-0.044

5.563

0.385

3.876

10.697

1.490

0.099

4.331

-4.513

0.814

0.019

14.153

0.000

0.000

55.610

46.642

9.853

0.000

40.339

5 0.638 -0.078 7.439 0.000 0.000 0.000 0.318 9.784 2.558 0.000 11.107 0.000 0.000 2.800 0.000 0.000 0.000 0.334 0.262 0.000 0.000 0.000 0.000 0.043 0.000 0.000 0.000 0.000 0.000 0.000 0.000 4.175

6 0.418 0.437 1.940 0.003 0.001 0.251 1.682 3.535 3.384 0.290 4.514 0.006 0.354 0.000 0.000 6.714 0.088 0.629 0.025 0.026 0.015 0.000 0.015 0.018 0.031 0.000 0.000 0.001 0.000 0.000 0.000 0.001

7 0.396 -0.513 4.602 0.000 0.000 4.772 0.822 2.003 0.000 0.000 14.203 0.002 0.375 0.000 0.000 2.553 1.724 5.797 0.078 0.000 0.000 0.000 0.000 0.000 0.743 0.000 0.000 0.000 0.000 0.000 0.000 0.000

8 0.714 0.350 4.340 0.000 0.000 0.000 0.209 0.088 1.077 0.000 0.552 0.000 0.038 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.023 0.000 0.037 0.000 0.121 0.000 0.000 0.000 0.000 0.000 0.094 0.465

9 0.757 1.636 -4.866 0.117 0.072 1.033 0.000 5.213 0.000 0.246 1.422 0.026 0.729 14.604 0.000 0.000 0.000 0.000 0.655 0.000 1.342 0.002 0.580 0.536 0.000 0.512 0.166 1.603 0.019 0.014 0.020 0.000

G

MS model 13, DM variant 4b, excl. intrazonaal NSTR R² alpha beta G1 G2 G3 G4 G5 G6 G7 G8 G9 G10 G11 G12 G13 G14 G15 G16 G17 G18 G19 G20 G21 G22 G23 G24 G25 G26 G27 G28 G29

0 0.373 1.965 -2.717 0.000 0.000 10.313 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.008 0.000 0.000 0.000 0.000 0.000 0.000 0.012 0.029 0.000 0.000 0.008 0.000 1.210

1 0.554 0.733 -0.044 0.000 0.000 0.000 0.001 0.000 0.000 0.000 0.109 0.000 0.000 0.000 0.182 0.002 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

2 0.978 2.544 12.029 0.000 0.000 0.007 0.000 0.832 0.000 0.000 0.025 0.000 0.000 0.000 1.974 0.006 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

3 0.951 0.662 0.700 0.000 0.000 0.009 0.036 0.635 0.000 0.002 0.000 0.000 0.001 0.000 0.000 0.000 0.000 0.035 0.000 0.000 0.000 0.000 0.000 0.000 0.230 0.000 0.000 0.000 0.000 0.000 0.000 0.000

4 0.985 0.783 3.877 0.000 0.002 0.000 0.000 0.000 0.013 0.000 0.497 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

74

Schattingen BASGOED


MS model 13, DM variant 4c, excl. intrazonaal NSTR R² alpha beta G1 G2 G3 G4 G5 G6 G7 G8 G9 G10 G11 G12 G13 G14 G15 G16 G17 G18 G19 G20 G21 G22 G23 G24 G25 G26 G27 G28 G29

0 0.357 1.046 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

1 0.554 0.728 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

2 0.976 2.860 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

3 0.951 1.016 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

4 0.984 1.203 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

5 0.594 0.338 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

6 0.374 0.988 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

7 0.292 0.532 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

8 0.711 1.144 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

9 0.692 0.852 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

MS model 13, DM variant 4d, excl. intrazonaal NSTR

0

1

2

3

4

5

6

7

8

9

R²

0.337

0.521

0.974

0.951

0.987

0.694

0.411

0.385

0.714

0.653

beta

2.304

4.190

18.280

1.728

13.195

18.628

3.640

3.273

6.282

1.837

G1

0.000

0.000

0.000

0.000

0.000

0.000

0.002

0.000

0.000

0.000

G2

0.051

0.000

0.002

0.000

0.001

0.561

0.000

0.005

0.000

0.000

G3

0.000

0.000

0.000

0.027

0.000

0.020

0.000

20.296

0.000

0.000

G4

0.000

0.000

0.204

0.080

0.000

3.729

0.621

0.206

0.119

0.297

G5

0.000

0.048

1.612

1.295

0.000

49.791

2.066

0.799

0.162

0.000

G6

0.000

0.000

0.052

0.000

0.014

7.185

1.801

0.000

0.867

1.783

G7

0.000

0.007

0.000

0.005

0.000

0.752

0.098

0.020

0.000

0.000

G8

0.666

0.018

0.000

0.000

0.268

75.282

2.092

87.249

0.176

0.000

G9

0.034

0.000

0.000

0.000

0.000

0.024

0.003

0.001

0.000

0.000

G10

0.116

0.000

0.001

0.003

0.000

0.000

0.144

0.218

0.028

0.000

G11

0.159

0.000

0.000

0.000

0.000

6.438

0.000

0.000

0.000

0.000

G12

0.000

0.000

2.184

0.000

0.000

0.000

0.000

1.569

0.000

0.580

G13

0.096

0.000

0.364

0.000

0.000

0.000

4.109

23.079

0.000

0.353

G14

0.085

0.000

0.014

0.003

0.000

0.000

0.036

0.666

0.000

0.158

G15

0.000

0.017

0.267

0.096

0.000

4.227

0.265

22.247

0.000

0.288

G16

0.017

0.020

0.000

0.000

0.000

2.827

0.013

0.158

0.000

0.000

75


Schattingen BASGOED

G17

0.091

0.012

0.000

0.000

0.000

0.000

0.012

0.000

0.000

0.000

G18

0.604

0.036

0.000

0.000

0.000

0.000

0.005

0.000

0.017

0.000

G19

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G20

0.000

0.022

0.000

0.000

0.000

0.170

0.005

0.000

0.018

0.000

G21

0.003

0.000

0.002

0.001

0.000

2.656

0.010

0.007

0.001

0.000

G22

0.023

0.000

0.016

0.577

0.000

1.469

0.012

0.271

0.067

0.230

G23

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G24

0.000

0.024

0.000

0.000

0.000

0.009

0.000

0.000

0.000

0.000

G25

0.032

0.000

0.000

0.000

0.000

0.627

0.001

0.000

0.000

0.000

G26

0.033

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G27

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

G28

0.007

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.051

0.000

G29

0.000

0.098

0.000

0.000

0.000

14.271

0.001

0.000

0.205

0.000

76

Schatting BASGOED. rapportage DP1 SIGNIFICANCE: GERARD DE JONG MICHIEL DE BOK KIM RUIJS DERK WENTINK

Recommend Documents