Fasit, Nasjonalregnskap
Johannes Mauritzen
Høsten 2017
Labb 5: Nasjonalregnskap
Norge, som de fleste land, samler inn data og beregner et nasjonalregnskap. Dette kan sammenlignes med vanlig regnskap. Det skal gi innblikk i den økonomiske situasjon til landet. Hvis vi har en resesjon, vet vi det på grunn av nasjonalregnskapet.
Vi kan enkelt få tilgang til nasjonalregnskapet direkt via SSB sin webside. Som vanlig, begynner vi med å laste inn noen pakker som vi trenger:
library(ggplot2)
library(dplyr)
library(reshape2)
library(zoo)Data kommer fra SSB, per kvartal og vi laster det inn direkt. Du kan lese mer om nasjonalregnskapet her: https://www.ssb.no/nasjonalregnskap-og-konjunkturer
NR<-read.csv("http://data.ssb.no/api/v0/dataset/59022.csv",
sep=";", header=FALSE, dec=",", na.strings = c(".", ".."), skip=1)
colnames(NR) = c("variabel", "kvartal", "enhet", "verdi")
kvartal <- as.yearqtr(NR$kvartal, format = "%YK%q") #fra zoo
NR["kvartal"] <- as.Date(kvartal)
str(NR)## 'data.frame': 7140 obs. of 4 variables:
## $ variabel: Factor w/ 51 levels "bif.nr8_ Bruttoinvestering i alt",..: 41 41 41 41 41 41 41 41 41 41 ...
## $ kvartal : Date, format: "2000-01-01" "2000-01-01" ...
## $ enhet : Factor w/ 2 levels "Faste 2015-priser (mill. kr)",..: 1 2 1 2 1 2 1 2 1 2 ...
## $ verdi : int 203931 210362 210180 213896 215569 213891 221117 212606 207273 215122 ...levels(NR$enhet)## [1] "Faste 2015-priser (mill. kr)"
## [2] "Faste 2015-priser, sesongjustert (mill. kr)"Vi kan velge mellom sesongjustert og ikke sesongjustert. Hva betyr sesongjustert? Er det grunner for å velge den ene eller den andre?
Vi velger å bruke sesongjusterte tall:
NRSA <- filter(NR, enhet == "Faste 2015-priser, sesongjustert (mill. kr)")
NRSA <- NRSA[,c("variabel", "kvartal", "verdi")]Her ser vi at variabelen består av en kategori (brutto-investering, privatkonsum, osv) og deretter underkategorier. Vi vil splitte variablen inn til to deler - kategori og underkategori.
levels(NRSA$variabel)## [1] "bif.nr8_ Bruttoinvestering i alt"
## [2] "bif.nr83_6 Bruttoinvestering i fast realkapital"
## [3] "bif.nr83_6fn \xac Fastlands-Norge (bruttoinvestering)"
## [4] "bif.nr83_6fnxof \xac\xac Fastlands-Norge utenom offentlig forvaltning (bruttoinvestering)"
## [5] "bif.nr8307 \xac\xac\xac\xac Tjenester tilknyttet utvinning (bruttoinvestering)"
## [6] "bif.nr8368 \xac\xac\xac Boliger (husholdninger) (bruttoinvestering)"
## [7] "bif.nr83indberg \xac\xac\xac\xac Industri og bergverk (bruttoinvestering)"
## [8] "bif.nr83naer \xac\xac\xac N\xe6ringer (bruttoinvestering)"
## [9] "bif.nr83oljroer \xac Utvinning og r\xf8rtransport (bruttoinvestering)"
## [10] "bif.nr83rest \xac\xac\xac\xac Andre tjenester (bruttoinvestering)"
## [11] "bif.nr83vare \xac\xac\xac\xac Annen vareproduksjon (bruttoinvestering)"
## [12] "bif.nr84_5 \xac\xac\xac Offentlig forvaltning (bruttoinvestering)"
## [13] "bif.pub8X50A \xac Utenriks sj\xf8fart (bruttoinvestering)"
## [14] "bil.nr8lu_ Lagerendring og statistiske avvik"
## [15] "biv.nr0890 Anskaffelser minus avhendelse av verdigjenstander"
## [16] "bnpb.nr23_6fn Bruttoprodukt Fastlands-Norge, basisverdi"
## [17] "bnpb.nr23_6fnxof \xac Fastlands-Norge utenom offentlig forvaltning (bruttoprodukt)"
## [18] "bnpb.nr23_9 Bruttonasjonalprodukt, markedsverdi"
## [19] "bnpb.nr23_9fn Bruttonasjonalprodukt Fastlands-Norge, markedsverdi"
## [20] "bnpb.nr23indberg \xac\xac Industri og bergverk (bruttoprodukt)"
## [21] "bnpb.nr23oljsj Oljevirksomhet og utenriks sj\xf8fart"
## [22] "bnpb.nr23ptjfn \xac\xac Tjenester inkl. boligtjenester (bruttoprodukt)"
## [23] "bnpb.nr23vare \xac\xac Annen vareproduksjon (bruttoprodukt)"
## [24] "bnpb.nr24_5 \xac Offentlig forvaltning (bruttoprodukt)"
## [25] "bnpb.nr29korr Produktavgifter og -subsidier"
## [26] "eks.nrskiplfly \xac Skip, plattformer og fly (eksport)"
## [27] "eks.nrtjen \xac Tjenester (eksport)"
## [28] "eks.nrtot Eksport i alt"
## [29] "eks.nrtradvare \xac Tradisjonelle varer (eksport)"
## [30] "eks.puboljegs \xac R\xe5olje og naturgass (eksport)"
## [31] "imp.nrskiplfly \xac Skip, plattformer og fly (import)"
## [32] "imp.nrtjen \xac Tjenester (import)"
## [33] "imp.nrtot Import i alt"
## [34] "imp.nrtradvare \xac Tradisjonelle varer (import)"
## [35] "imp.puboljegs \xac R\xe5olje og naturgass (import)"
## [36] "koh.nr61_ \xac Konsum i husholdninger"
## [37] "koh.nr61L8 \xac\xac Husholdningenes kj\xf8p i utlandet"
## [38] "koh.nr61L9 \xac\xac Utlendingers kj\xf8p i Norge"
## [39] "koh.nr61tjen \xac\xac Tjenestekonsum"
## [40] "koh.nr61vare \xac\xac Varekonsum"
## [41] "koh.nrpriv Konsum i husholdninger og ideelle organisasjoner"
## [42] "koi.nr66_ \xac Konsum i ideelle organisasjoner"
## [43] "koo.nr64_ \xac Konsum i statsforvaltningen"
## [44] "koo.nr6401 \xac\xac\xac Konsum i statsforvaltningen, forsvar"
## [45] "koo.nr64sivil \xac\xac\xac Konsum i statsforvaltningen, sivilt"
## [46] "koo.nr65_ \xac Konsum i kommuneforvaltningen"
## [47] "koo.nroff Konsum i offentlig forvaltning"
## [48] "makrok.nrefn Ettersp\xf8rsel fra Fastlands-Norge (ekskl. lagerendring)"
## [49] "makrok.nreoff Ettersp\xf8rsel fra offentlig forvaltning"
## [50] "makrok.nrianv Innenlandsk sluttanvendelse"
## [51] "makrok.nrsanv Samlet sluttanvendelse"Vi bruker funksjonen “separate” fra pakken tidyr til å gjøre dette:
library(tidyr) #Problemer? Prøv å installere pakken først: install.package("tidyr")##
## Attaching package: 'tidyr'## The following object is masked from 'package:reshape2':
##
## smithsNRSA <- separate(NRSA, variabel, into = c("kateg", "variabel"), sep="\\.", extra = "merge")Man kan lese mer om separate og tidyr her
Vi definerer de nye kolomene kateg og variabel som faktor:
NRSA["kateg"]<-factor(NRSA$kateg)
NRSA["variabel"]<-factor(NRSA$variabel)Nå vil vi se litt nærmere på investering.
SSB kaller kategorien for Bruttoinvestering: Hvorfor brutto? Hva ville netto-investering vært?
Er investering et flyt eller en lager?
investering = NRSA[NRSA$kateg =="bif",]
# eller: investering=filter(NRSA, kateg=="bif")
investering["variabel"] <- factor(investering$variabel)
levels(investering$variabel)## [1] "nr8_ Bruttoinvestering i alt"
## [2] "nr83_6 Bruttoinvestering i fast realkapital"
## [3] "nr83_6fn \xac Fastlands-Norge (bruttoinvestering)"
## [4] "nr83_6fnxof \xac\xac Fastlands-Norge utenom offentlig forvaltning (bruttoinvestering)"
## [5] "nr8307 \xac\xac\xac\xac Tjenester tilknyttet utvinning (bruttoinvestering)"
## [6] "nr8368 \xac\xac\xac Boliger (husholdninger) (bruttoinvestering)"
## [7] "nr83indberg \xac\xac\xac\xac Industri og bergverk (bruttoinvestering)"
## [8] "nr83naer \xac\xac\xac N\xe6ringer (bruttoinvestering)"
## [9] "nr83oljroer \xac Utvinning og r\xf8rtransport (bruttoinvestering)"
## [10] "nr83rest \xac\xac\xac\xac Andre tjenester (bruttoinvestering)"
## [11] "nr83vare \xac\xac\xac\xac Annen vareproduksjon (bruttoinvestering)"
## [12] "nr84_5 \xac\xac\xac Offentlig forvaltning (bruttoinvestering)"
## [13] "pub8X50A \xac Utenriks sj\xf8fart (bruttoinvestering)"Vi rydder opp i navnene:
invest_var <- c("Bruttoinvestering i alt",
"fast realkapital",
"Fastlands-Norge",
"Fastlands-Norge utenom offentlig forvaltning",
"Tjenester tilknyttet utvinning" ,
"Boliger",
"Industri og bergverk",
"Naeringer",
"Utvinning og roertransport",
"Andre tjenester",
"Annen vareproduksjon",
"Offentlig forvaltning",
"Utenriks sjoefart")
#her skal vi endre kategori-variablene, som blir kalt "levels"
levels(investering$variabel) <- invest_varOg da kan vi bruke ggplot og facet_wrap til å se på alle subkategoriene:
ggplot(investering, aes(x=kvartal, y=verdi)) +
geom_line() +
facet_wrap(~variabel, ncol=2, scales="free")
Kanskje dette er litt for mye informasjon - la oss begrense det litt:
invest_l <- c("Bruttoinvestering i alt",
"fast realkapital",
"Fastlands-Norge utenom offentlig forvaltning",
"Boliger", "Industri og bergverk",
"Naeringer", "Utvinning og roertransport",
"Offentlig forvaltning")Med denne kommandoen skal vi begrense vår datasett til kategoriene i variabel som vi har samlet in invest_l
investering_l <- investering[investering$variabel %in% invest_l,]
ggplot(investering_l, aes(x=kvartal, y=verdi)) +
geom_line() +
facet_wrap(~variabel, ncol=2, scales="free")
Kan man se noen trender her?
Hvilken typer investering har gått ned de siste årene? Hvilken typer har gått opp? Kan du forklare hvorfor?
Vi har hatt veldig lave renter de siste par årene. Hvordan har det påvirket investering og i hvilken kategori?
Handelsbalansen - eksport og import
Nå ser vi litt på handelsbalansen - eksport og import
eksport <- NRSA[NRSA$kateg == "eks", ]
levels(eksport$variabel)## [1] "nr0890 Anskaffelser minus avhendelse av verdigjenstander"
## [2] "nr23_6fn Bruttoprodukt Fastlands-Norge, basisverdi"
## [3] "nr23_6fnxof \xac Fastlands-Norge utenom offentlig forvaltning (bruttoprodukt)"
## [4] "nr23_9 Bruttonasjonalprodukt, markedsverdi"
## [5] "nr23_9fn Bruttonasjonalprodukt Fastlands-Norge, markedsverdi"
## [6] "nr23indberg \xac\xac Industri og bergverk (bruttoprodukt)"
## [7] "nr23oljsj Oljevirksomhet og utenriks sj\xf8fart"
## [8] "nr23ptjfn \xac\xac Tjenester inkl. boligtjenester (bruttoprodukt)"
## [9] "nr23vare \xac\xac Annen vareproduksjon (bruttoprodukt)"
## [10] "nr24_5 \xac Offentlig forvaltning (bruttoprodukt)"
## [11] "nr29korr Produktavgifter og -subsidier"
## [12] "nr61_ \xac Konsum i husholdninger"
## [13] "nr61L8 \xac\xac Husholdningenes kj\xf8p i utlandet"
## [14] "nr61L9 \xac\xac Utlendingers kj\xf8p i Norge"
## [15] "nr61tjen \xac\xac Tjenestekonsum"
## [16] "nr61vare \xac\xac Varekonsum"
## [17] "nr64_ \xac Konsum i statsforvaltningen"
## [18] "nr6401 \xac\xac\xac Konsum i statsforvaltningen, forsvar"
## [19] "nr64sivil \xac\xac\xac Konsum i statsforvaltningen, sivilt"
## [20] "nr65_ \xac Konsum i kommuneforvaltningen"
## [21] "nr66_ \xac Konsum i ideelle organisasjoner"
## [22] "nr8_ Bruttoinvestering i alt"
## [23] "nr83_6 Bruttoinvestering i fast realkapital"
## [24] "nr83_6fn \xac Fastlands-Norge (bruttoinvestering)"
## [25] "nr83_6fnxof \xac\xac Fastlands-Norge utenom offentlig forvaltning (bruttoinvestering)"
## [26] "nr8307 \xac\xac\xac\xac Tjenester tilknyttet utvinning (bruttoinvestering)"
## [27] "nr8368 \xac\xac\xac Boliger (husholdninger) (bruttoinvestering)"
## [28] "nr83indberg \xac\xac\xac\xac Industri og bergverk (bruttoinvestering)"
## [29] "nr83naer \xac\xac\xac N\xe6ringer (bruttoinvestering)"
## [30] "nr83oljroer \xac Utvinning og r\xf8rtransport (bruttoinvestering)"
## [31] "nr83rest \xac\xac\xac\xac Andre tjenester (bruttoinvestering)"
## [32] "nr83vare \xac\xac\xac\xac Annen vareproduksjon (bruttoinvestering)"
## [33] "nr84_5 \xac\xac\xac Offentlig forvaltning (bruttoinvestering)"
## [34] "nr8lu_ Lagerendring og statistiske avvik"
## [35] "nrefn Ettersp\xf8rsel fra Fastlands-Norge (ekskl. lagerendring)"
## [36] "nreoff Ettersp\xf8rsel fra offentlig forvaltning"
## [37] "nrianv Innenlandsk sluttanvendelse"
## [38] "nroff Konsum i offentlig forvaltning"
## [39] "nrpriv Konsum i husholdninger og ideelle organisasjoner"
## [40] "nrsanv Samlet sluttanvendelse"
## [41] "nrskiplfly \xac Skip, plattformer og fly (eksport)"
## [42] "nrskiplfly \xac Skip, plattformer og fly (import)"
## [43] "nrtjen \xac Tjenester (eksport)"
## [44] "nrtjen \xac Tjenester (import)"
## [45] "nrtot Eksport i alt"
## [46] "nrtot Import i alt"
## [47] "nrtradvare \xac Tradisjonelle varer (eksport)"
## [48] "nrtradvare \xac Tradisjonelle varer (import)"
## [49] "pub8X50A \xac Utenriks sj\xf8fart (bruttoinvestering)"
## [50] "puboljegs \xac R\xe5olje og naturgass (eksport)"
## [51] "puboljegs \xac R\xe5olje og naturgass (import)"eksport["variabel"] <- factor(eksport$variabel)levels(eksport$variabel) <- c("Skip, plattformer og fly",
"Tjenester",
"Eksport i alt",
"Tradisjonelle varer",
"Raaolje og naturgass")Oppgaver:
Oppgave 1.)
Prøv å plotte eksport-kategoriene med ggplot og facet_wrap.
Kan man gjøre det samme for import?
Hvilken trender ser man i disse dataseriene?
Har norge en positiv handelsbalanse?
Om eksport av olje og gass skulle forsvinne, hva ville skjedd med handelsbalansen?
Vi kan også finne data om privatkonsum (“koh”) og offentligkonsum (“okon”). Dra ut disse dataene og bruk ggplot til å se på alle subkategoriene. Finner man noen trender?
import <- NRSA[NRSA$kateg == "imp",]
import["variabel"] <- factor(import$variabel)
levels(import$variabel) <- c("Skip, plattformer og fly",
"Tjenester",
"Import i alt",
"Tradisjonelle varer",
"Raaolje og naturgass")ggplot(import, aes(x=kvartal, y=verdi)) +
geom_line() +
facet_wrap(~variabel, ncol=2, scales="free")
Privat og offentlig konsum
pkon <- NRSA[NRSA$kateg=="koh",]
okon <- NRSA[NRSA$kateg=="koo",]
pkon$variabel <- factor(pkon$variabel)
okon$variabel <- factor(okon$variabel)
levels(pkon$variabel)<-c("Konsum i husholdninger",
"Husholdningenes kjoep i utlandet",
"Utlendingers kjoep i Norge",
"Tjenestekonsum",
"Varekonsum",
"Konsum i husholdninger, totalt")
levels(okon$variabel) <- c("Statsforvaltningen",
"Statsforvaltningen, forsvar",
"Statsforvaltningen, sivilt",
"Kommuneforvaltningen",
"Offentlig forvaltning")ggplot(pkon, aes(x=kvartal, y=verdi)) +
geom_line() +
facet_wrap(~variabel, ncol=2, scales="free")
ggplot(okon, aes(x=kvartal, y=verdi)) +
geom_line() +
facet_wrap(~variabel, ncol=2, scales="free")
Oppgave 2: Intro til økonometri
I denne oppgaven skal vi begynne å øve på litt formell økonometri. Hovedverktøyet er linear-regresjon.
Vi skal prøve å svare på et enkelt spørsmål: Hvor mye påvirker investeringer i boliger BNP?
Vi kan begynne med å finne fram til data og rydde litt:
bnp <- NRSA[NRSA$kateg == "bnpb",]
bnp$variabel <- factor(bnp$variabel)
bnp_var <- c("Bruttoprodukt Fastlands-Norge, basisverdi",
"Fastlands-Norge utenom offentlig forvaltning (bruttoprodukt)",
"Bruttonasjonalprodukt, markedsverdi",
"Bruttonasjonalprodukt Fastlands-Norge, markedsverdi",
"Industri og bergverk (bruttoprodukt)",
"Oljevirksomhet og utenriks sjoefart",
"Tjenester inkl. boligtjenester (bruttoprodukt)",
"nnen vareproduksjon (bruttoprodukt)",
"Offentlig forvaltning (bruttoprodukt)",
"Produktavgifter og -subsidier")
levels(bnp$variabel) <- bnp_varNå er det flere måter å måle BNP på. Hvilken metrikk skal vi velge? Total BNP - “Bruttonasjonalprodukt, markedsverdi” - inkluderer også olje og gass. Er det grunn for å inkludere dette, eller burde vi se på kun BNP, Fastlands-Norge.
Jeg velger å bruke fastlands-Norge:
privat_bnp_fastland = bnp[bnp$variabel == "Bruttonasjonalprodukt Fastlands-Norge, markedsverdi",]
privat_bnp_fastland = privat_bnp_fastland[,c("kvartal", "verdi")]
colnames(privat_bnp_fastland)[2]<-"privat_bnp_fastland"Nå får vi også tall for investering:
investering = NRSA[NRSA$kateg =="bif",]
investering["variabel"] <- factor(investering$variabel)
levels(investering$variabel)## [1] "nr8_ Bruttoinvestering i alt"
## [2] "nr83_6 Bruttoinvestering i fast realkapital"
## [3] "nr83_6fn \xac Fastlands-Norge (bruttoinvestering)"
## [4] "nr83_6fnxof \xac\xac Fastlands-Norge utenom offentlig forvaltning (bruttoinvestering)"
## [5] "nr8307 \xac\xac\xac\xac Tjenester tilknyttet utvinning (bruttoinvestering)"
## [6] "nr8368 \xac\xac\xac Boliger (husholdninger) (bruttoinvestering)"
## [7] "nr83indberg \xac\xac\xac\xac Industri og bergverk (bruttoinvestering)"
## [8] "nr83naer \xac\xac\xac N\xe6ringer (bruttoinvestering)"
## [9] "nr83oljroer \xac Utvinning og r\xf8rtransport (bruttoinvestering)"
## [10] "nr83rest \xac\xac\xac\xac Andre tjenester (bruttoinvestering)"
## [11] "nr83vare \xac\xac\xac\xac Annen vareproduksjon (bruttoinvestering)"
## [12] "nr84_5 \xac\xac\xac Offentlig forvaltning (bruttoinvestering)"
## [13] "pub8X50A \xac Utenriks sj\xf8fart (bruttoinvestering)"invest_var <- c("Bruttoinvestering i alt",
"fast realkapital",
"Fastlands-Norge",
"Fastlands-Norge utenom offentlig forvaltning",
"Tjenester tilknyttet utvinning" ,
"Boliger (husholdninger)",
"Industri og bergverk ",
"Naeringer ",
"Utvinning og roertransport",
"Andre tjenester",
"Annen vareproduksjon",
"Offentlig forvaltning",
"Utenriks sjoefart")
levels(investering$variabel) <- invest_var
inv_bolig <- investering[investering$variabel=="Boliger (husholdninger)",]
inv_bolig <- inv_bolig[,c("kvartal", "verdi")]
colnames(inv_bolig)[2] <- "bolig_investering"
head(inv_bolig)## kvartal bolig_investering
## 1541 2000-01-01 24011
## 1542 2000-04-01 23634
## 1543 2000-07-01 24494
## 1544 2000-10-01 25184
## 1545 2001-01-01 26151
## 1546 2001-04-01 26332Nå setter vi de to dataseriene vi er interessert i en dataframe ved å bruke merge
model_data <- merge(privat_bnp_fastland, inv_bolig, by="kvartal")
head(model_data)## kvartal privat_bnp_fastland bolig_investering
## 1 2000-01-01 448789 24011
## 2 2000-04-01 445405 23634
## 3 2000-07-01 449328 24494
## 4 2000-10-01 450443 25184
## 5 2001-01-01 459190 26151
## 6 2001-04-01 455303 26332Vi kan bruke ggplot til å visualisere dataen:
model_data_long <- melt(model_data, id.var = "kvartal")
ggplot(model_data_long, aes(x=kvartal, y=value)) +
geom_line() +
facet_wrap(~variable, nrow=2, scales="free")
Er det noe som er verdt å legge merke til her?
Svar
En ting som vi burde være obs på er at det finnes en trend i begge disse seriene. Dette er problematisk. Vi kan tenke oss at det finnes en manglende variabel som forårsaker at begge bolig investering of bnp øker over tid. Når vi kjører en regresjon uten å ta hensikt til dette vil vi da få resultater som gjenspeiler dette og ikke nødvendigvis forholdet mellom bolig investering og BNP:
Nå kan vi kjøre en regresjon. Vi kan begynne med å kjøre den enkleste regresjonen mulig:
\[bnp = \alpha + \beta bolig\_investering + \epsilon\]
Dette er veldig enkelt å gjøre i R:
model1<- lm(privat_bnp_fastland ~ bolig_investering, data = model_data)Her har vi kjørt regresjonen med funksjonen lm som står for “linear model”, og lagret resultatene i en objekt som heter model1. Ellers kan man se på de ulike delene (eller “attributes”) ved å bruke str
Vi kan oppsummere resultatene ved å bruke summary.
summary(model1)##
## Call:
## lm(formula = privat_bnp_fastland ~ bolig_investering, data = model_data)
##
## Residuals:
## Min 1Q Median 3Q Max
## -55382 -26748 -4381 23341 73262
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 2.174e+05 2.113e+04 10.29 1.66e-15 ***
## bolig_investering 1.012e+01 6.080e-01 16.64 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 32700 on 68 degrees of freedom
## Multiple R-squared: 0.8028, Adjusted R-squared: 0.7999
## F-statistic: 276.8 on 1 and 68 DF, p-value: < 2.2e-16- Hvordan kan vi tolke disse resultatene?
Svar
I følge resultatene, vil en økning av 1.000 kr i boliginvestering, føre til en økning av 10.000 kr i BNP.
- Gir de noen mening? Er det riktig?
Nei, vi kan ikke tolke dette som at 1.000 kr i bolig investering forårsaker 10.000 i høyere bnp. Vi har ikke tatt hensyn til felles trenden i datasettene. Dette er en form av det som vi skal kalle “ikke-stasjonæritet”, som gjør vår regresjon nærmest meningsløst.
Et annet probem er at BNP (Y) og investering i boliger (I) er endogene. De bestemmer hverandre. Derfor kan vi ikke si “investering i bolig forårsaker endring i BNP”. Det ville vært like riktig og like feil å skrive “En endring i BNP forårsaker en endring i BNP.”
- Hvordan kunne man endre regresjonen til å få resultater med mer mening?
To forslag - vi kunne sette in en tidstrend. Det vil si at vi satt in data for kvartal, som kontrollerer for trenden:
model1<- lm(privat_bnp_fastland ~ kvartal + bolig_investering, data = model_data)Et annet ide er å transformere datasettene til endringer eller prosent endringer. Da kan vi tolke resultatene som “en endring i bolig investering fører til en endring i BNP”
Men problemet med endogenitet er fortsatt der. Dette er et mye vanskeligere problem å løse.