Dear All,
Can anyone provide me with reference notes(or steps) towards analysis ofÂ
(un)balanced panel data in R.
Thank you!
Kind regards,
Tabitha Mundia , Project Management Office, Equity Bank
Limited,P.O. Box 75104-00200
Head Office, Upper hill, NHIF BLDG, 14th Floor, Nairobi, Kenya
Direct Extension : +254732112721Â Â Mobile: +254722309538
Email: tabitha.mun...@equitybank.co.keskype ID: twamaeYahoo ID:
tabygath...@yahoo.com
An idea not coupled with action will never get any bigger than the brain cell
it occupied.
Arnold Glasgow
......
"Attempt something large enough that failure is guaranteedâ¦unless God steps
in!"
--- On Wed, 12/29/10, Manoj Aravind <aravin...@gmail.com> wrote:
From: Manoj Aravind <aravin...@gmail.com>
Subject: Re: [R] Problem applying McNemar's - Different values in SPSS and R
To: "Marc Schwartz" <marc_schwa...@me.com>
Cc: r-help@r-project.org
Date: Wednesday, December 29, 2010, 2:48 PM
Thank you Marc :)
It Certainly helped me to get the exact value of P.
How to understand when to apply mcnemar.exact or just mcnemar.test?
I'm a beginner to biostatistics.
Manoj Aravind
On Tue, Dec 28, 2010 at 11:00 PM, Marc Schwartz <marc_schwa...@me.com>wrote:
>
> On Dec 28, 2010, at 11:05 AM, Manoj Aravind wrote:
>
> > Hi friends,
> > I get different values for McNemar's test in R and SPSS. Which one should
> i
> > rely on when the p values differ.
> > I came across this problem when i started learning R and seriously give
> up
> > on SPSS or any other proprietary software.
> > Thank u in advance
> >
> > Output in SPSS follows
> >
> > *Crosstab*
> >
> >
> >Â Â Â Â Â Â Â Â Â hsc
> >
> > Total
> >
> >Â Â Â Â ABN
> >
> > NE
> >
> > ABN
> >
> > tvs
> >
> > ABN
> >
> > Count
> >
> > 40
> >
> > 3
> >
> > 43
> >
> >Â Â Â Â Row %
> >
> > 93.0%
> >
> > 7.0%
> >
> > 100.0%
> >
> >Â Â Â Â COL%
> >
> > 78.4%
> >
> > 30.0%
> >
> > 70.5%
> >
> >Â NE
> >
> > Count
> >
> > 11
> >
> > 7
> >
> > 18
> >
> >Â Â Â Â Row %
> >
> > 61.1%
> >
> > 38.9%
> >
> > 100.0%
> >
> >Â Â Â Â COL%
> >
> > 21.6%
> >
> > 70.0%
> >
> > 29.5%
> >
> > Total
> >
> > Count
> >
> > 51
> >
> > 10
> >
> > 61
> >
> >Â Row %
> >
> > 83.6%
> >
> > 16.4%
> >
> > 100.0%
> >
> >Â COL%
> >
> > 100.0%
> >
> > 100.0%
> >
> > 100.0%
> >
> >
> >
> > * Chi-Square Tests*
> >
> >
> >Â Â Â Value
> >
> > Exact Sig. (2-sided)
> >
> > McNemar Test
> >
> >Â .057(a)
> >
> > N of Valid Cases
> >
> > 61
> >
> >Â Â Â a Binomial distribution used.
> >
> > Output from R is as follows....
> >
> >> tvshsc<-
> >
> > + matrix(c(40,11,3,7),
> >
> > + nrow=2,
> >
> > + dimnames=list("TVS"=c("ABN","NE"),
> >
> > + "HSC"=c("ABN","NE")))
> >
> >> tvshsc
> >
> >Â Â Â Â HSC
> >
> > TVSÂ Â Â ABN NE
> >
> >Â ABNÂ 40Â 3
> >
> >Â NEÂ Â Â 11Â 7
> >
> >> mcnemar.test(tvshsc)
> >
> >
> > McNemar's Chi-squared test with continuity correction
> >
> >
> > data:Â tvshsc
> >
> > McNemar's chi-squared = 3.5, df = 1, p-value = 0.06137
> >
> > Regards
> >
> > Dr. B Manoj Aravind
>
>
> The SPSS test appears to be an exact test, whereas the default R function
> does not perform an exact test, so you are not comparing Apples to Apples...
>
> Try this using the 'exact2x2' CRAN package:
>
> > require(exact2x2)
> Loading required package: exact2x2
> Loading required package: exactci
>
> > mcnemar.exact(matrix(c(40, 11, 3, 7), 2, 2))
>
>Â Â Â Â Exact McNemar test (with central confidence intervals)
>
> data:Â matrix(c(40, 11, 3, 7), 2, 2)
> b = 3, c = 11, p-value = 0.05737
> alternative hypothesis: true odds ratio is not equal to 1
> 95 percent confidence interval:
>Â 0.04885492 1.03241985
> sample estimates:
> odds ratio
>Â 0.2727273
>
>
> HTH,
>
> Marc Schwartz
>
>
   [[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.
[[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.
