如何在治疗和控制之间自动进行数据集中多个组的显着性检验

Question

我真的很陌生。

我的数据如下：每个城市都属于一个州，并具有变量Num1的处理和控制值。我需要执行两种方式的t检验或Anova（无论如何）以检查治疗组城市的Num1值和城市所属的控制状态的Num1值之间的显着性（应针对每个城市进行此操作）。每个城市治疗状态对照组的结果应为p值。

City Measure_group   State  Day        Num1
A    Treatment        X    01.01.1900   279
A    Control          X    02.01.1900   256
B    Treatment        Y    03.01.1900   71
B    Control          Y    04.01.1900   372
C    Treatment        Z    05.01.1900   127
C    Control          Z    06.01.1900   303
D    Treatment        W    07.01.1900   248

我已经尝试了基本的ANOVA来完成此操作，但是它无法将所有城市相互比较而没有任何关系。

    model <- aov( Num1 ~City+State, data= data1)
    comparison <- pairs(lsmeans(model, ~City|State))
    test(comparison, by = NULL, adjust = "fdr")