```r
# only statement # (idx) as predictor: coded as 0, 1, 2, 3
> summ(feols(pfc ~ idx, d0_pfc[idx < 4], cluster = ~ResponseId))
term result
<char> <char>
1: (Intercept) b = 70.82 [69.68, 71.97], p < .001
2: idx b = 5.53 [5.05, 6.02], p < .001
```
```r
# 3-way interation model
term result
<char> <char>
1: (Intercept) b = 70.44 [69.55, 71.33], p < .001
2: idx b = 5.68 [5.23, 6.13], p < .001
3: topicC b = 8.31 [6.53, 10.09], p < .001
4: conditionC b = -34.67 [-36.44, -32.89], p < .001
5: idx × topicC b = -1.60 [-2.50, -0.71], p < .001
6: idx × conditionC b = 8.78 [7.88, 9.67], p < .001
7: topicC × conditionC b = 6.96 [3.40, 10.51], p < .001
8: idx × topicC × conditionC b = 0.31 [-1.48, 2.09], p = .737
```
![[1732478931.png]]