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 - morri...@ibms.sinica.edu.tw
______________________________________________
R-help@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.
