Hello everyone, I have a question about the commands “psmatch2” and “pstest”.
For propensity score matching, the result that I get with Stata version 17 is different from what is shown on a video using Stata version 14.
I summarize the whole process here ( I copied the log file into this rtf file [Array ] ):
-I have a data set of medical records;
-There are 4 variables in the dataset:
age
male
var1
outcomevar
-I generated a matching variable based on var1.
gen matchvar=0
replace matchvar =1 if var1 >14.0
-Then I used psmatch2 and paired the data.
psmatch2 matchvar age male, outcome(outcomevar) caliper(.001) noreplace neighbor(1)logit
gen pair = _id if _treated==0
replace pair = _n1 if _treated==1.
bysort pair: egen paircount = count(pair)
drop if paircount !=2
-Up to this point everything is as it is supposed to be.
-Then I used “pstest“
pstest age male var1 outcomevar
-pstest command without var1 included, works well. But if I include var1 it gives a r(2000) error.
-These are the important pieces of information:
1. var1 (a continuous variable) is the variable that we used for generating the matching variable.
2. In stata version 17, the command < pstest.. var1…> has no output just headings of the table and r(2000) error. [the result is on log file]
3. In stata version 14, the command <pstest …var1…> generate “the variable, mean and bias table” and “the summary of the distribution of the bias table” and a r(2000) error [ I can send a screenshot from the video]
4. command <pstest2 …var1…> works perfectly in stata version 17.
I checked the ado-file of pstest and pstest2, however, I couldn’t find the reason.
I would appreciate it if you could shed some light on it.
Related Posts with “psmatch2” and “pstest”: different results in version 14 and 17
About pre-treatment balance and difference-in-differences methodHello, I am a PhD student and I am planning to apply natural experiment methods (I have a proper sh…
Calculating rate and confidence interval for a series of data.Hello all I am doing some time series stuff (using R to do interrupted time series analysis - long …
Using local variable in forvalues loopHi everyone, I am trying to run the following code: count if productcode == 21 scalar nobs = r(N) …
Simple table export question: loop or by function?Hi everyone I am currently working on a monitoring and evaluation program where I need to create a …
Add string to empty cellHi all, I need add string (N/A) to empty cells for a string Variable, I used "subinstr", but lament…
Subscribe to:
Post Comments (Atom)
0 Response to “psmatch2” and “pstest”: different results in version 14 and 17
Post a Comment