我已经使用dplyr函数为两个类别生成了线性模型。现在想比较系数
首先，我不明白为什么在这种情况下为什么需要两个不同的数据集df和dfTrend，否则我在两个国家的结果中的值将相同。
calculation <- dlply (df, "categorie", function(df) lm( df$response ~ df$input , data = dfTrend))
然后执行：
ldply(calculation, coef)
我的结局是
categorie (Intercept) df$input
1 country 1 -6431798 3231.953
2 country 2 -4650245 2334.810
其次，现在我想检查一下Anova的系数是否与我正在努力解决的差异显着
anova(calculation)
结果：
UseMethod（“ anova”）中的错误： 没有适用于“方差”的适用方法应用于“列表”类的对象 或不同的方法：线性假设不能应用
linearHypothesis(calculation, country 1 = country 2)
结果：
vcov.default（model，complete = FALSE）中的错误： 没有用于类列表模型的vcov（）方法