Не удалось сопоставить перехваты через emmeans в R

Я хотел бы проверить симетрию в реакции наблюдателя на контрастные стимулы разной полярности, положительные (белый) и отрицательные (черные). Я взял время реакции (RT) как зависимую переменную по четырем различным контрастам. Известно, что время отклика следует кривой Пьерона, асимптоты которой помещены (1) на пороге наблюдателя (Inf) и (2) на базовом RT, расположенном где-то между 250 и 450 мсек. Эти знания позволяют нам линеаризовать взаимосвязь, преобразуя независимую переменную (эффективный контраст EC) как 1 / EC ^ 2 (tEC), поэтому уравнение, связывающее RT с EC, становится:

RT = m * tEC + RT0

Чтобы проверить симметрию, я установил критерии: одинаковый наклон и одинаковый отрезок в двух полярностях подразумевают симметрию. Чтобы получить коэффициенты, я сделал линейную модель с взаимодействием (кодирование через фиктивную переменную для полярности: положительная или отрицательная). Для меня результат lm ясен, но некоторые коллеги предпочитают что-то более похожее на результат ANOVA. Поэтому я решил использовать emmeans, чтобы создать контраст. С наклоном все в порядке, но при вычислении пересечений начинаются проблемы. Перехваты, вычисленные lm, сильно отличаются от вывода emmeans, и выводы также разные. Ниже я воспроизведу этот пример. Вопрос двоякий: можно ли использовать emmeans для решения моей проблемы? Если нет, можно ли контрастировать с помощью других пакетов (какой)?

Данные

Код

# Model
model <- lm(RT1000 ~ tEC * Polarity, data = Data)

# emmeans
library(emmeans)
# Slopes
m.slopes <- lstrends(model, "Polarity", var="tEC")
# Intercepts
m.intercept <- lsmeans(model, "Polarity")

# Contrasts
pairs(m.slopes)
pairs(m.intercept)

Выходы

Вычисленные перехваты с помощью emmeans отличаются от тех, которые вычисляются с помощью lm. Думаю, проблема в том, что модель не определена для EC = 0. Но я не уверен.