similar to: Propensity score and three treatments

Dear List, I am having difficulty running the ps() function when variables are stored as factors and was hoping someone could provide some advice on how to proceed. I am running propensity score matching as outlined in: Greg Ridgeway, Dan McCarey, Andrew Morral, Lane Burgette and Beth Ann Grin (May 3, 2013) Toolkit for Weighting and Analysis of Nonequivalent Groups: A tutorial for the twang

"Matching" version 0.48 is now available on CRAN. Matching provides functions for estimating causal effects by multivariate and propensity score matching. The package includes a variety of univariate and multivariate tests to determine if balance has been obtained by the matching procedure. These tests can also be used to determine if an experiment or quasi-experiment is balanced on

"Matching" version 0.48 is now available on CRAN. Matching provides functions for estimating causal effects by multivariate and propensity score matching. The package includes a variety of univariate and multivariate tests to determine if balance has been obtained by the matching procedure. These tests can also be used to determine if an experiment or quasi-experiment is balanced on

Hi all, i am wondering if there?s any other method to adjust for selection related bias of estimates except propensity scoring and heckit / mills ratio approach? i also read documentation of Match and twang package so far, so i don?t speak of any ATE / ATT related methods, respectively any methods that match or stratify... Is there something else ? thx in advance

Hola chic en s, alguien con experiencia en propensión score matching? Planteo duda: Clasicamente el PSM se ha utilizado en un intento de homogeneizar cohortes de enfermos quienes han estado ?expuestos? a un tratamiento x Vs aquellos que no han estado expuestos (no expuestos). Esto aplica para medicamentos o procedimientos quirúrgicos o no. Bien, En algún articulo he leído que el PSM se puede

I'm using the "Match" package to do propensity score matching. Here's some example code that shows the problem that I'm having (much of this code is taken from the Match package documentation): *data(lalonde) glm1 <- glm(treat~age + I(age^2) + educ + I(educ^2) + black + hisp + married + nodegr + re74 + I(re74^2) + re75 + I(re75^2) + u74 + u75,

Dear R-list, i have 6 different sets of samples. Each sample has about 5000 observations, with each observation comprised of 150 baseline covariates (X), 125 of which are dichotomous. Roughly 20% of the observations in each sample are "treatment" and the rest are "control" units. i am doing propensity score matching, i have already estimated propensity scores(predicted

Dear All, My question is more a statistical question than a R question. The reason I am posting here is that there are lots of excellent statistician on this list, who can always give me good advices. Per my understanding, the purpose of propensity score is to reduce the bias while estimating the treatment effect and its implementation is a 2-stage model. 1) First of all, if we assume that T =

There may be benefits to having a machine learning method that explicitly targets covariate balance. We have experimented with optimizing the weights directly to obtain the best covariate balance, but got some strange solutions for simple cases that made us wary of such methods. Machine learning methods that yield calibrated probability estimates should do well (e.g. those that optimize the

Hello all. I am doing one part of an evaluation of a mandatory welfare-to-work programme in the UK. As with all evaluations, the problem is to determine what would have happened if the initiative had not taken place. In our case, we have a number of pilot areas and no possibility of random assignment. Therefore we have been given control areas. My problem is to select for survey individuals in

Hi all, i am looking to built a simple example of a very basic propensity score adjustment, just using the estimated propensity scores as inverse probability weights (respectively 1-estimated weights for the non-treated). As far as i understood, MLE predictions of a logit model can directly be used as to estimates of the propensity score. I already considered the twang package and the

I have a longitudinal competing risk data of the form: ID COMPL SEX HEREDITY 1 0 1 2 1 0 1 2 1 3 1 2 2 0 0 1 2 1 0 1 2 2 0 1 2 2 0 1 3 0 0 1 3 0 0 1 3 0 0 1 3 0 0 1 3 2 0 1 4 0 1 2 4 0 1

Hola chic en s, alguien con experiencia en propensión score matching? Planteo duda: Clasicamente el PSM se ha utilizado en un intento de homogeneizar cohortes de enfermos quienes han estado ?expuestos? a un tratamiento x Vs aquellos que no han estado expuestos (no expuestos). Esto aplica para medicamentos o procedimientos quirúrgicos o no. Bien, En algún articulo he leído que el PSM se puede

Dear all, I am comparing logistic regression models to evaluate if one predictor explains additional variance that is not yet explained by another predictor. As far as I understand Baron and Li describe how to do this, but my question is now: how do I report this in an article? Can anyone recommend a particular article that shows a concrete example of how the results from te following simple

Hi I am getting the following errors when i try to use yum to install the net-snmp paclages. [root at sc1 yum.repos.d]# yum install net-snmp Loaded plugins: fastestmirror Determining fastest mirrors Traceback (most recent call last): File "/usr/bin/yum", line 29, in ? yummain.user_main(sys.argv[1:], exit_code=True) File "/usr/share/yum-cli/yummain.py", line 229,

Hello R listing, I did two-way anova on lm. Further question the investigator interested in is what two treatments are different? I am looking for a command which could do pairwise comparison for every treatment. Could anyone help me out? Thanks a bunch Kevin

Hi R-users I'd like to know if there is some packages or functions to do comparison of treatments after that manova pointed differences between them. Any suggestions are much appreciated. __________________________________________________________ Eng. Agr., M.Sc. Eduardo Dutra de Armas __________________________________________________________ Centro de Energia Nuclear na Agricultura

Hi, I am comparing three treatments. The data come from a longitudinal study, and for each treatment I have only 1 case, on which the observations across 20 years are made. My main aim is to compare the three treatments, but the effect of time is of interest as well. Which of the many R functions should I use for that purpose? I applied the aov function and received this:

Hi Folks, I'm a bit puzzled by the following (example): N<-factor(sample(c(1,2,3),1000,replace=TRUE)) unique(N) # [1] 3 2 1 # Levels: 1 2 3 So far so good. Now: contrasts(N)<-contr.treatment(3, base=1, contrasts=FALSE) contrasts(N) # 1 2 # 1 1 0 # 2 0 1 # 3 0 0 whereas: contr.treatment(3, base=1, contrasts=FALSE) # 1 2 3 # 1 1 0 0 # 2 0 1 0 # 3 0 0 1 contr.treatment(3, base=1,

Hi there, I need some help to figure out what is the proper model in survival analysis for my data. Subjects were randomized to 3 treatments in trial 1, some of them experience the event during the trial; After period of time those subjects were randomized to 3 treatments again in trial 2, but different from what they got in 1st trial, some of them experience the event during the 2nd trial (I