隨機生成Example of nominal coding of drug use...表中矩陣

（b）表中的0.56等數字，暫時按照0,1錄入

aaa=matrix(nrow=50,ncol=2)  #生成第一個對應cc.yes和cc.no的矩陣aaa
for (i in 1:50) { #按行循環50次
  j=sample(c(1:2),1) #在1和2中隨機抽取一個數字，賦值給j，做爲矩陣aaa的列數
  aaa[i,j]=1 #給第i行第j列（第一列或第二列）賦值爲1
  aaa[i,-j]=0} #其他列賦值爲0
bbb=matrix(nrow=50,ncol=2) #生成第一個對應cm.yes和cm.no的矩陣bbb
for (i in 1:50) {
  j=sample(c(1:2),1)
  bbb[i,j]=1
  bbb[i,-j]=0}
ccc=matrix(nrow=50,ncol=2) #生成第一個對應e.yes和e.no的矩陣ccc
for (i in 1:50) {
  j=sample(c(1:2),1)
 ccc[i,j]=1
 ccc[i,-j]=0}
yyy=matrix(nrow=50,ncol=3) #生成第一個對應COMT.AG、COMT.AG和COMT.GG的矩陣yyy
for (i in 1:50) {
j=sample(c(1:3),1)
yyy[i,j]=1
yyy[i,-j]=0}
xxx=matrix(nrow=50,ncol=3)#生成第一個對應FAAH.CA、FAAH.AA和COMT.CC的矩陣xxx
for (i in 1:50) {
  j=sample(c(1:3),1)
  xxx[i,j]=1
  xxx[i,-j]=0}
dg=cbind(aaa,bbb,ccc,yyy,xxx)
colnames(dg)[1:12]=c('cc.yes','cc.no','cm.yes','cm.no','F.yes','F.no',
                     'COMT.AG','COMT.AA','COMT.GG','FAAH.CA','FAAH.AA','FAAH.CC')
dgx=dg[,1:6]
dgy=dg[,7:12]