library(actuar) # for rinvweibull
library(deamer)
set.seed(123)
RT <- rinvweibull(10000, shape=5.53861156, scale=488)/1000
RT <- RT[RT<1.5]
noise <- 1/rnorm(10000, mean=0.0023853421, sd=0.0004784688)/1000
noise <- noise[noise<1.5]

ST <- deamerSE(RT, errors=noise)
plot(ST)