## Not run: 
# ##See svmmatch() for a full implementation
# ##Load data
# 	data("LaLonde")
# 	Data1<-LaLonde
# 	Data1<-Data1[Data1$exper==0|Data1$treat==1,]
# 	attach(Data1)	
# 
# ##Format X matrix
# 	varnames<-c("age","educ","black","married","nodegr","hisp",
# 		"re75","re74")
# 	X<-cbind(Data1[,varnames],Data1$re75==0,Data1$re74==0)
# 	X<-as.matrix(X)
# 
# ##Fit model
# 	set.seed(1)
# 	m1.param<-svmmatch(treat, X, dv=re78, burnin=100, gibbs=100, thin=5)
# 
# ##Assess ignorability:	
# 	#Assess control overlap--are there enough
# 	#controls to balance the treated?
# 		cont1<-control.overlap(m1.param)	
# 
# 	## End(Not run)
Run the code above in your browser using DataLab