R：如何生成這個條件二項式隨機變量？

Question

我想從以下模型生成Xi。我知道如何從binomial distribution生成數據，我可以做的是如下，但我怎樣才能生成數據在同一時間滿足條件，Y = 1或Y =-1？

set.seed(100) 
n1<-50 
n2<-50 
p1<- 0.4 
p2<- 0.3 

y <- c(rep(1,50),rep(-1,50)) 
xi1<- c(rbinom(n1,1,p1), rbinom(n2,1,p2)) 

Answer 1

因爲rbinom()被矢量，我們可以做

n <- ifelse(y == 1, n1, n2) 
p <- ifelse(y == 1, p1, p2) 
rbinom(n, 1, p) 

你最好閱讀?rbinom密切。