[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