R/ps.R

In kimihua1995/adjKMtest:

Defines functions ps

Documented in ps

# calculating propensity score
ps <- function(treat,x){
  glm.sol1<-glm((treat)~x, family=binomial(link="logit"))
  beta=glm.sol1$coef
  pre<-cbind(1,x)%*%beta
  ex<-(1-1/(1+exp(pre)))

  return(ex)
}