Re: [R] Help on performing multiple chi-square tests

From: Jorge Ivan Velez <jorgeivanvelez_at_gmail.com>
Date: Sun, 09 Nov 2008 17:32:40 -0500

Dear Tania,
An option would be splitting the data by marker and the perform the chi-square test selecting both the statistics and p-value for each marker. Here is an example:

# Dummy data set
mydata=read.table(textConnection("
Marker Treatment Genotype1 Genotype2 Genotype3

1        A         23        57        32
1        B         43        59        12
2        A         13        27        12
2        B         23        29        22"),header=TRUE)


# Chi-square

res=lapply(split(mydata,mydata$Marker),function(x){
res=chisq.test(x[,-c(1,2)]) # Deleting columns 1 and 2 of every list
res=c(res$statistic,res$p.value)

names(res)=c('statistic','pvalue')
res
}
)

# The output
res2=do.call(rbind,res)
rownames(res2)=paste('marker_',1:length(res),sep="") res2

# statistic       pvalue
#marker_1 15.169487 0.0005081452
#marker_2  2.010403 0.3659709476


HTH, Jorge

On Sun, Nov 9, 2008 at 4:14 PM, QUESADA,TANIA <tquesada_at_ufl.edu> wrote:

> Hi,
>
> I am new to R and I need to perform multiple chi-square tests. I manage to
> perform one at a time, but is there a specific command to do multiple tests?
>
> For example, I have a table that looks like this:
>
> Marker Treatment Genotype1 Genotype2 Genotype3
> 1 A 23 57 32
> 1 B 43 59 12
> ...
> ...
> n A ## ## ##
> n B ## ## ##
>
> I can perform an individual chi-square test for marker 1, but if n is too
> large (i.e. 4000 markers), I would like to write the code to make it perform
> tests for each marker types and give me an output table with the p-values
> for each markers.
>
> I would appreciate any suggestions you may have.
>
> Thanks!
>
> Tania Quesada
>
>
>
> --
> QUESADA,TANIA
>
> ______________________________________________
> R-help_at_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_at_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. Received on Sun 09 Nov 2008 - 22:34:52 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Mon 10 Nov 2008 - 00:30:23 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive