困惑於理解一個rbinom輸出

Question

p1 <- c(.25,.025,.025,.1,.2,.4) 

    T <- sample(1:6,size=N,replace=TRUE, prob=someprobabilityvector) 
    Y <- rbinom(N,1,p1[c(T)]) 

嗨，夥計們，我是新來的R和編程一般，需要一些幫助理解基本。有人能向我解釋一個向量Y上面發生了什麼。我計算出p1[c(T)]在上面做了什麼。但不知道Y在做什麼向量。所有的幫助提前感謝。

Answer 1

代碼的第一行創建六個概率的矢量：

p1 <- c(.25,.025,.025,.1,.2,.4) 

在第二行，你隨機選擇從數字一N值六（帶替換）。每個值的概率在someprobabilityvector中指定。因此，該函數將返回長N包括1和6之間

T <- sample(1:6,size=N,replace=TRUE, prob=someprobabilityvector) 

在第三行中的值的矢量，生成從與一個試驗和在p1[c(T)]指定概率問題二項分佈N隨機數。 c(T)與T相同：該矢量包括從1到6的值。該矢量用於索引矢量p1。因此，p1[c(T)]將返回一個矢量，其中包括來自矢量p1的N值。

Y <- rbinom(N,1,p1[c(T)]) 

由於指定二項式分佈具有一個試驗只，載體Y將包含零和一。