随机生成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]