Я новичок в R. Я попытался сделать множественную регрессию, используя зависимую переменную "daglig" и некоторые независимые переменные с именами OPd, NOKUSD и др. c. Также я использовал несколько фиктивных переменных с именем Inn_ut, a, b, c ... ao.
Когда я набираю сценарий и запускаю каждое предложение, я получаю «NA» для дат в моем наборе данных. Даты не включены в R, они указаны на моем листе Excel.
Может кто-нибудь сказать мне, если я сделал что-то не так? Являются ли результаты, которые я получил, «понятными»? Я никогда не делал этого раньше.
Вот мой сценарий:
library(tokenizers)
library(corpus)
library(tidygraph)
library(igraph)
library(openxlsx)
library(ggplot2)
library(stringr)
library(stargazer)
data2 <- read.xlsx("Ferdigreg2.xlsx")
head(data2)
lm(daglig~OPd,data=data2)
library(stargazer)
m4 <- lm(daglig~OPd+NOKUSD+NOKEUR+APRd+OSEBX,data=data2)
library(corrplot)
stargazer(m4, title="data2", align=TRUE,
no.space=TRUE,type="text")
m5 <- lm(daglig~OPd+NOKUSD+NOKEUR+APRd+OSEBX+Inn_ut+a+b+c+d+e+f+g+h+i+j+k+l+m+n+o+p+q+r+s+t+u+v+w+x+y+z+aa+ab+ac+ad+ae+af+ag+ah+ai+aj+ak+al+am+an+ao,data=data2)
stargazer(m5, title="data2", align=TRUE,
no.space=TRUE,type="text")
А вот мои результаты в консоли:
> library(tokenizers)
> library(corpus)
> library(tidygraph)
> library(igraph)
> library(openxlsx)
> library(ggplot2)
> library(stringr)
> library(stargazer)
> data2 <- read.xlsx("Ferdigreg2.xlsx")
> head(data2)
X1 X2 Dato daglig OPd NOKUSD NOKEUR APRd OSEBX Inn_ut a b c d e f g h i j k l m n o p q r s t u v w x y z aa ab ac ad ae af ag ah ai aj
1 577.7497 23.3 42461 - 0.030042918 36.42 8.2576 9.4401 22.6 -0.013780881 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 2 569.7878 NA 42464 -0.008849558 36.05 8.2917 9.4359 22.4 -0.005152444 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 3 566.8520 NA 42465 0.017857143 35.88 8.3571 9.4995 22.8 -0.012376776 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 4 559.8362 NA 42466 -0.008771930 37.77 8.3592 9.4760 22,6 0,005455167 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 5 562.8902 NA 42467 -0.017699115 37,15 8,3198 9,4546 22,2 - 0.004379184 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 6 560.4252 NA 42468 0.004504505 40.71 8.2886 9.4183 22.3 0.010371054 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 ak al am a a 05.2016 06.2016 07.2016 08.2016 09.2016 10.2016 11.2016 12.2016 01.2017 02.2017 03.2017 04.2017 05.2017 06.2017 07.2017 08.2017 09.2017 10.2017 1 0 0 0 0 0 НА НА НА НА НА НА НА НА НА НА НА НА НА НА НА НА 2 0 0 0 0 0 Н А Н А Н А Н А Н А Н А Н А Н А Н А Н А Н А Н А Н А Н А Н А Н А Н А Н А Н А Н А Н НА Н А Н N N NA 4 0 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 6 NA 0 0 NA NA NA NA NA NA NA NA NA NA NA NA 11.2017 12.2017 01.20178 02.2018 03.2018 04.2018 05.2018 06.2018 07.2018 08.2018 09.2018 10.2018 11.2018 12.2018 01.2019 02.2019 03.2019 04.2019 05.2019 06.2019 1 NA NA NA NA NA NA NA NA NA NA NA Нет нет нет нет нет нет нет нет больше нет нет больше нет ничего нет больше нет ничего нет больше нет ничего нет больше нет больше нет ничего нет больше нет больше нет ничего нет больше нет ничего больше нет. Нет Нет Нет Нет Нет Нет Нет Нет Нет Нет Нет Нет Нет Нет Нет Нет Нет Нет Нет Нет Нет Нет Нет Нет Нет Нет больше нет больше нет ничего нет больше нет. А Н А Н А Н А Н А Н А Н А Н А Н А Н А Н А Н А Н А Н А Н А Н А Н А Н А Н НА НА 07.2019 08.2019 09.2019 1 Н А Н А Н 2 NA NA NA 3 NA NA NA 4 NA NA NA 5 NA NA NA 6 NA NA NA
lm (daglig ~ OPd, data = data2)
Call:
lm(formula = daglig ~ OPd, data = data2)
Coefficients:
(Intercept) OPd
-0.0059818 0.0001071
> library(stargazer)
> m4 <- lm(daglig~OPd+NOKUSD+NOKEUR+APRd+OSEBX,data=data2)
> library(corrplot)
> stargazer(m4, title="data2", align=TRUE,
+ no.space=TRUE,type="text")
data2
===============================================
Dependent variable:
---------------------------
daglig
-----------------------------------------------
OPd 0.0001***
(0.00005)
NOKUSD 0.001
(0.003)
NOKEUR -0.004
(0.003)
APRd 0.0004**
(0.0002)
OSEBX 0.807***
(0.088)
Constant 0.015
(0.030)
-----------------------------------------------
Observations 874
R2 0.101
Adjusted R2 0.095
Residual Std. Error 0.022 (df = 868)
F Statistic 19.425*** (df = 5; 868)
===============================================
Note: *p<0.1; **p<0.05; ***p<0.01
> m5 <- lm(daglig~OPd+NOKUSD+NOKEUR+APRd+OSEBX+Inn_ut+a+b+c+d+e+f+g+h+i+j+k+l+m+n+o+p+q+r+s+t+u+v+w+x+y+z+aa+ab+ac+ad+ae+af+ag+ah+ai+aj+ak+al+am+an+ao,data=data2)
> stargazer(m5, title="data2", align=TRUE,
+ no.space=TRUE,type="text")
data2
===============================================
Dependent variable:
---------------------------
daglig
-----------------------------------------------
OPd 0.0002***
(0.0001)
NOKUSD 0.007
(0.011)
NOKEUR 0.009
(0.014)
APRd 0.005***
(0.001)
OSEBX 0.796***
(0.086)
Inn_ut -0.039***
(0.011)
a 0.003
(0.007)
b -0.009
(0.007)
c -0.019***
(0.007)
d -0.016**
(0.007)
e -0.027***
(0.008)
f -0.039***
(0.009)
g -0.029***
(0.009)
h 0.019
(0.013)
i 0.018
(0.013)
j 0.012
(0.014)
k 0.016
(0.012)
l 0.013
(0.011)
m 0.017*
(0.009)
n -0.021***
(0.007)
o -0.025***
(0.007)
p -0.025***
(0.008)
q -0.025***
(0.008)
r -0.023***
(0.008)
s -0.026***
(0.008)
t -0.024**
(0.010)
u -0.024**
(0.010)
v -0.013
(0.011)
w -0.006
(0.010)
x -0.002
(0.010)
y 0.005
(0.009)
z -0.018**
(0.008)
aa -0.028***
(0.008)
ab -0.025***
(0.008)
ac -0.020**
(0.008)
ad -0.017**
(0.007)
ae -0.024***
(0.008)
af -0.008
(0.009)
ag -0.023**
(0.009)
ah -0.057***
(0.010)
ai -0.055***
(0.010)
aj -0.076***
(0.012)
ak -0.117***
(0.016)
al -0.061***
(0.010)
am -0.041***
(0.009)
an -0.025***
(0.007)
ao
Constant -0.275**
(0.109)
-----------------------------------------------
Observations 874
R2 0.197
Adjusted R2 0.152
Residual Std. Error 0.022 (df = 827)
F Statistic 4.405*** (df = 46; 827)
===============================================
Note: *p<0.1; **p<0.05; ***p<0.01