[R] How to stop loop inside status ? ( haplo.stats package )
Nash
morrison at ibms.sinica.edu.tw
Tue Mar 10 17:48:51 CET 2009
How to stop loop inside status ?
require(haplo.stats)
## normal status
y=rep(c(0,1),each=50)
geno=as.data.frame(matrix(sample(c("A","G","T","C"),600,replace=T),100,6))
hs<- haplo.score(y, geno, trait.type="binomial", offset = NA, x.adj = NA,
min.count=5,
locus.label=NA, miss.val=c(0,NA), haplo.effect="additive",
eps.svd=1e-5, simulate=TRUE,
sim.control=score.sim.control(min.sim=200,max.sim=500))
hs$score.global.p.sim
## find loop inside ,and I can't stop it!
geno=as.data.frame(matrix("G",100,6))
hs<- haplo.score(y, geno, trait.type="binomial", offset = NA, x.adj = NA,
min.count=5,
locus.label=NA, miss.val=c(0,NA), haplo.effect="additive",
eps.svd=1e-5, simulate=TRUE,
sim.control=score.sim.control(min.sim=200,max.sim=500))
hs$score.global.p.sim
How to stop loop inside status ?
--
Nash - morrison at ibms.sinica.edu.tw
More information about the R-help
mailing list