Error in post-hoc robust two way ANOVA

[Airen19]

I am trying to do a robust two way ANOVA due to a lack of normality, and t2way runs perfectly but mcp2atm gives an error. The data is a NO3 measurement of 12 plots, 8 samplings and in 4 different cropping systems (Rotation). I have a RCBD experimental design with 3 replicates.

I do not know how to proceed, or if my data can be analysed with a robust ANOVA. Any info will be very appreciated. Thanks!

t2way(NO3 ~ Rotation * Date, data = VBSL1)

mcp2atm(NO3 ~ Rotation * Date, data = VBSL1)

Error in nfac[unique(mf[, 2]), unique(mf[, 3])] : subscript out of bounds

[FJCC]

What is the result of running

mf <-  model.frame(NO3 ~ Rotation * Date, data = VBSL1)
summary(mf)
unique(mf[, 2])
unique(mf[, 3])

[Airen19]

mf <- model.frame(NH4 ~ Rotation * Date, data = VBSL1)
summary(mf)
NH4 Rotation Date
Min. : 0.2067 Length:82 Min. :2020-02-28 00:00:00.00
1st Qu.: 1.9096 Class :character 1st Qu.:2020-05-06 00:00:00.00
Median : 3.3181 Mode :character Median :2020-10-27 00:00:00.00
Mean : 4.6105 Mean :2020-10-30 04:05:51.22
3rd Qu.: 4.7587 3rd Qu.:2021-04-20 00:00:00.00
Max. :37.5380 Max. :2021-07-01 00:00:00.00
unique(mf[, 2])
[1] "B-B" "B-W-P"
unique(mf[, 3])
[1] "2020-02-28 UTC" "2020-05-06 UTC" "2020-08-06 UTC" "2020-10-27 UTC" "2021-01-19 UTC" "2021-04-20 UTC" "2021-07-01 UTC"

[Airen19]

I am trying with a similar dataset but now the error has changed:

M2<- t2way(NH4 ~ Phase2 * Moment2, data = season2, est="median")
M2
Call:
t2way(formula = NH4 ~ Phase2 * Moment2, data = season2, est = "median")

post hoc tests

mcp2atm(NH4 ~ Phase * Moment, data = season2, est = "median")
Error in x[[j]] :
attempt to select less than one element in integerOneIndex

mf <- model.frame(NH4 ~ Rotation * Moment, data = season2)
summary(mf)
NH4 Rotation Moment
Min. : 0.2067 Length:48 Length:48
1st Qu.: 1.3557 Class :character Class :character
Median : 2.9491 Mode :character Mode :character
Mean : 5.2154
3rd Qu.: 4.8030
Max. :37.5380