[R] Simulation Case/Control R Beginner

From: Claire_6700 <chriseverclear_at_hotmail.com>
Date: Tue, 06 May 2008 14:37:29 -0700 (PDT)

I am looking for a way to simulate genotypes of cases and control at a disease locus in R. I am supposed to set the allele frequency as control/cases. for each of the column below simulate 200 snp dataset.

I am looking at treesim function from popgen to stimulate the genotypes in R.

Here is what i have

size         case   control 
2000        0.6        0.15
2500        0.1        0.17

what am i missing in the treesim parameter and how do i calculate the p-value of trend test?

I am a bit confused with the parameter.

N<-2000
Sn<-200
theta<-0

a<-treesim(N, Sn, theta, rho.vec, mutations=FALSE, sample=FALSE) mat<-LDmat(.....)

Here is the documentation.
http://cran.r-project.org/web/packages/popgen/popgen.pdf

any help will really be helpful.

thank you,
Claire

-- 
View this message in context: http://www.nabble.com/Simulation-Case-Control-R-Beginner-tp17092885p17092885.html
Sent from the R help mailing list archive at Nabble.com.

______________________________________________
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 Tue 06 May 2008 - 21:49:24 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 Tue 06 May 2008 - 22:30:35 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