R/dlogPrgcm_HW1.theta.R

In SPmlficmcm: Semiparametric Maximum Likelihood Method for Interactions Gene-Environment in Case-Mother Control-Mother Designs

dlogPrgcm_HW1.theta <-
function(vmc,p){    
                           v1<-ifelse(vmc[,1]==0 & vmc[,2]==0,-S0(-p),0)+
                           ifelse(vmc[,1]==0 & vmc[,2]==1,S0(p),0)+
                           ifelse(vmc[,1]==1 & vmc[,2]==0,-S0(-p),0)+
                           ifelse(vmc[,1]==1 & vmc[,2]==2,S0(p),0)+
                           ifelse(vmc[,1]==2 & vmc[,2]==1,-S0(-p),0)+
                           ifelse(vmc[,1]==2 & vmc[,2]==2,S0(p),0) 
                           return(v1)
                           }