Блог

Доброе утро, я пытаюсь создать линейную смешанную модель со случайным эффектом, используя plm пакет R. Я работаю над df , который вы можете скачать, нажав здесь, и структурирован следующим образом:

 > str(tab)
'data.frame':   498 obs. of  4 variables:
 $ X            : int  1 2 3 4 5 6 7 8 9 10 ...
 $ biking       : num  30.8 65.13 1.96 44.8 69.43 ...
 $ smoking      : num  10.9 2.22 17.59 2.8 15.97 ...
 $ heart.disease: num  11.77 2.85 17.18 6.82 4.06 ...
 

Когда я пытаюсь построить линейную модель следующим образом

 > lin_mod<-lm(heart.disease~smoking, data=tab)
 

это дает мне линейную модель без ошибок, но когда я пытаюсь ввести случайный эффект с помощью plm следующим образом

 random<-plm(heart.disease ~ smoking, data=tab, model="random")
 

или

 random<-plm(heart.disease ~ smoking X, data=tab, model="random")
 

или

 random<-plm(heart.disease ~ smoking biking, data=tab, model="random")
 

это дает мне эту ошибку:
Ошибка в plm.fit(данные, модель = модели[1L], эффект = эффект) : пустая модель

Я посмотрел, но, по-видимому, я не совершаю никаких заслуживающих внимания ошибок, но, конечно, я делаю что-то не так. Не могли бы вы дать мне несколько советов? Спасибо!