Я действительно новичок в R.
Мои данные выглядят так: Каждый город принадлежит государству и имеет значения обработки и контроля для переменной 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")