使用plm对智能电表进行面板分析

Question

我有大约1000户伦敦家庭的智能电表消费面板数据。提前告知住户电价是低，中还是高。我的数据如下：

> head(tof_clean)
             DateTime        id       acorn consumption_Kwh temp_max Tariff tariff_numeric
1 2013-01-01 00:00:00 MAC000146   Adversity           0.155     7.49 Normal              0
2 2013-01-01 00:00:00 MAC000147 Comfortable           0.276     7.49 Normal              0
3 2013-01-01 00:00:00 MAC000158 Comfortable           0.152     7.49 Normal              0
4 2013-01-01 00:00:00 MAC000165    Affluent           0.401     7.49 Normal              0
5 2013-01-01 00:00:00 MAC000170 Comfortable           0.640     7.49 Normal              0
6 2013-01-01 00:00:00 MAC000173    Affluent           0.072     7.49 Normal              0

我想测试一下关税低，中或高时消费与时间之间是否有意义（我创建了一个数字变量0,1,2或对应于关税变量）。由于消费与季节性高度相关，因此我无法简单地比较消费值的均值，因此我希望使用plm软件包进行面板回归。

我已经运行了这段代码：

plm_random <- plm(formula = consumption_Kwh ~ tariff_numeric, data = tof_clean, model="random",index = c("DateTime"))

但是会收到此错误：

Error in plm.fit(data, model = models[1], effect = effect) : empty model

您能提供关于我做错了什么以及如何有效实现所需结果的建议吗？我是面板回归的新手。

谢谢！