我正在尝试使用tau = 0.5模拟数据。我正在使用以下r code来模拟数据。
#rm(list=ls())
library(quantreg)
set.seed(12)
n = 100
tau = 0.5
u <- runif(n)
z <- rbinom(n, 1, 0.37)
v <- rnorm(n)
d <- z * (u > 0.5 * v)
x <- runif(n, 0, 1)
x1 <- (u - tau) + x * 2
y <- (1 - d) * x + d * x1
它正确吗？另外，对于给定的r和y值，是否有任何x函数可用于模拟tau？我注意到arima.sim()可用于模拟时间序列数据。