drgegne.gm: causal models with three possible causal factors: G, G*E and...
In IFP: Identifying Functional Polymorphisms

Description Usage Arguments Value See Also Examples

provides concordance probabilities of relative pairs for a causal model with G, G*E and E components

1	drgegne(fdg,frg,fdge,frge,eg,e)

`fdg`	an array (size=number of dominant genes+recessive genes) of dominant gene frequencies including 0 values of recessive genes of G component
`frg`	an array (size=number of dominant genes+recessive genes) of recessive gene frequencies including 0 values of dominant genes of G component
`fdge`	an array (size=number of dominant genes+recessive genes) of dominant gene frequencies including 0 values of recessive genes of G*E component
`frge`	an array (size=number of dominant genes+recessive genes) of recessive gene frequencies including 0 values of dominant genes of G*E component
`eg`	a proportion of population who are exposed to environmental cause of GE interactiong the genetic cause of GE during their entire life
`e`	a proportion of population who are exposed to environmental cause during their entire life

matrix of NN, ND, and DD probabilities of 9 relative pairs: 1:mzt,2:parent-offspring,3:dzt,4:sibling,5:2-direct(grandparent-grandchild),6:3rd(uncle-niece),7:3-direct(great-grandparent-great-grandchild),8:4th (causin),9:4d(great-great-grandparent-great-great-grandchild)

drgn drgene

### PLI=0.01.
ppt<-0.01



pg<-0.002  # the proportion of G component in total populations
pge<-0.005  # the proportion of G*E component in total populations
e<-1-(1-ppt)/(1-pg)/(1-pge)   
  # the proportion of E component in total populations

fd<-0.001  # one dominant gene
tt<-2      # the number of recessive genes

temp<-sqrt(1-((1-pg)/(1-fd)^2)^(1/tt))
fr<-c(array(0,length(fd)),array(temp,tt))
fd<-c(fd,array(0,tt))

ppe<-0.5
ppg<-pge/ppe

fdge<-0.002
ttge<-2      # the number of recessive genes

temp<-sqrt(1-((1-ppg)/(1-fdge)^2)^(1/ttge))
frge<-c(array(0,length(fdge)),array(temp,ttge))
fdge<-c(fdge,array(0,ttge))


drgegne(fd,fr,fdge,frge,ppe,e)