[R] scrime Package simulatedSNP function

Claire_6700 chriseverclear at hotmail.com
Fri May 9 01:41:18 CEST 2008


Hello,

I need some help with the simulatedSNPs function from scrime package.

I am trying to simulate some genotype of a case/control disease locus. The
allele frequence are cases/controls

Sample     cases      controls
2000         .5            .10
1500         .6            .40

In each of the row, i need to simulate 100 snp and calculate the pvalue

##############Download Scrime Package###########################
library(scrime)
n.obs<-1000
n.snp<-100
vec.ia<-1
simulateSNPs(n.obs, n.snp, vec.ia, prop.explain = 1,
list.ia.val = NULL, vec.ia.num = NULL, maf = c(0.1, 0.12),
prob.val = rep(1/3, 3), list.equal = NULL, prob.equal = 0.8,
rm.redundancy = TRUE, shuffle = FALSE, shuffle.obs = FALSE, rand = NA)

What is the right parameter. I am pretty new with R.

wrong result.
  Interaction Cases Controls
1   SNP1 == 2   500        0

any help will be appreciated.

thanks
Claire
-- 
View this message in context: http://www.nabble.com/scrime-Package-simulatedSNP-function-tp17138820p17138820.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list