模拟趋势

Question

所以我想创建一个模拟数据集，其中我的结果（y）随着时间的推移而增加，并且与两个协变量x1和x2相关联。这是一个例子：

expit<-function(z){1/(1+exp(-(z)))}
year<- factor(sample(c(2005:2014), 1e6, replace = T, 
       prob = rep(0.10, 10)))    
x1<-rbinom(1e6, 1, .15)
x2<- rbinom(1e6, 1, expit(-2+log(2)*x1))
y<- rbinom(1e6, 1, .15 + .35*x1 + .15*x2) 

我想知道如何随着时间的增加y增加x1和x2的年增长概率。换句话说，在十年的时间里，我希望x1和x2的概率每年增加y的概率.02。

有什么建议吗？

最好，

→