Dear Peter,
Thanks for your kind help. I got the solution now.
Best Regards,
Shyam
>________________________________
> From: Peter Ehlers <ehl...@ucalgary.ca>
>To: shyam basnet <shyamabc2...@yahoo.com>
>Cc: "r-help@R-project.org" <r-help@r-project.org>
>Sent: Monday, July 16, 2012 1:16 AM
>Subject: Re: [R] how to extract p-value in GenMatch function
>
>On 2012-07-15 14:37, shyam basnet wrote:
>> Dear R-Users,
>>
>> I have a problem on extracting T-Stat and P-Value. I have written R-code
>> below
>>
>> library("Matching")
>>
>> data("lalonde")
>> attach(lalonde)
>> names(lalonde)
>> Y<- lalonde$re78
>>
>> Tr<- lalonde$treat
>> glm1<-
>> glm(Tr~age+educ+black+hisp+married+nodegr+re74+re75,family=binomial,data=lalonde)
>>
>> pscore.predicted<- predict(glm1)
>>
>> rr1<- Match(Y=Y,Tr=Tr,X=glm1$fitted,estimand="ATT",
>> M=1,ties=TRUE,replace=TRUE)
>>
>> summary(rr1)
>>
>>> summary(rr1)
>>
>> Estimate... 2624.3
>> AI SE...... 802.19
>> T-stat..... 3.2714
>> p.val...... 0.0010702
>>
>> Original number of observations.............. 445
>> Original number of treated obs............... 185
>> Matched number of observations............... 185
>> Matched number of observations (unweighted). 344
>>
>> In above output, I can extract "Estimate" and "AI SE" with below code:
>> rr1$est
>>
>> rr1$se
>>
>> But the problem is I could not extract T-statistic and P-value from the
>> above output.
>>
>>
>> Could you please someone help me to resolve this problem?
>
>You could look at the code for summary.Match to see that
>T-stat (not surprisingly) is calculated as est/se and
>p.val is calculated as (1 - pnorm(abs(est/se))) * 2.
>summary.Match() doesn't return these values, it just
>prints them.
>
>Peter Ehlers
>
>> Thanking you,
>>
>> Best Regards,
>>
>> Shyam Basnet
>> SLU, Uppsala, Sweden
>>
>> [[alternative HTML version deleted]]
>
>
>
>
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
