[R] chisq.test using amalgamation automatically (possible ?!?)

Mohammad Ehsanul Karim wildscop at yahoo.com
Sun Jun 26 12:00:15 CEST 2005


Dear List,


If any of observed and/or expected data has less than
5 frequencies, then  chisq.test (Pearson's Chi-squared
Test for Count Data from package:stats) gives warning
messages. For example,

x<-c(10, 14, 10, 11, 11, 7, 8, 4, 1, 4, 4, 2, 1, 1, 2,
1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1)
y<-c(9.13112391745095, 13.1626482033341,
12.6623267638188, 11.0130706413029, 9.16415925139016,
7.47441794889028, 6.03743388141852, 4.85350508692505,
3.89248001363859, 3.11803140037476, 2.49617540962629,
1.99774139023269, 1.5985926374167, 1.27909653584089,
1.02341602646530, 0.818828097315106,
0.655132353196336, 0.524159229418155,
0.418022824890164, 0.335528136508225,
0.268448671671046, 0.214779801990545,
0.171840507806838, 0.137485729582785,
0.109999238967747, 0.0880079144684513,
0.070413150156564)

Chi.Sq<-sum((c(x[1:7], sum(x[8:9]), sum(x[10:11]),
sum(x[12:27]))-c(y[1:7], sum(y[8:9]), sum(y[10:11]),
sum(y[12:27])))^2/c(y[1:7], sum(y[8:9]),
sum(y[10:11]), sum(y[12:27]))) # using amalgamation
pchisq(Chi.Sq, df=9, ncp=0, lower.tail = FALSE, log.p
= FALSE) # result being 0.8830207

but chisq.test(x,y) gives the following output with
incorrect df:

        Pearson's Chi-squared test

data:  x and y 
X-squared = 216, df = 208, p-value = 0.3373

Warning message:
Chi-squared approximation may be incorrect in:
chisq.test(x, y) 



Is there any way that we can use directly chisq.test
without having warning message in such cases (that is,
using amalgamation conveniently so that we don't have
to check each elements if they are less than 5 or not
- the whole process being automatic, may be by means
of programming)?



Any hint, help, support, references will be highly
appreciated.
Thank you for your time. 

----------------------------------

Mohammad Ehsanul Karim 

Web: http://snipurl.com/ehsan 
ISRT, University of Dhaka, BD 

----------------------------------




		
____________________________________________________ 

Rekindle the Rivalries. Sign up for Fantasy Football




More information about the R-help mailing list