setwd("~/Desktop")
Sweave("turboEM.Rnw")
Sweave("/Users/jmfeder/Documents/School/EM Accelaration/turboEM/inst/doc/turboEM.Rnw")
Sweave("/Users/jmfeder/Documents/School/EM Accelaration/turboEM/inst/doc/turboEM.Rnw")
q("no")
Sweave("turboEM")
Sweave("turboEM")
q("no")
Sweave("turboEM")
Sweave("turboEM")
Sweave("turboEM")
res <- turboem(par=c(0.5, 1, 3), fixptfn=fixptfn, objfn=objfn, method=c("EM", "squarem", "pem"), y=y, control.run=list(keep.objfval=TRUE))
plot(res, xlim=c(0.001, 0.02))
plot(res, xlim=c(0.001, 0.02))
Sweave("turboEM")
res
res2 <- turboem(par = c(0.9, 1, 3), fixptfn=fixptfn, objfn=objfn, boundary=boundary, pconstr=pconstr, method=c("em", "squarem", "squarem", "decme", "qn", "qn"), control.method=list(list(), list(K=2), list(K=3), list(version="v2"), list(qn=1), list(qn=2)), control.run=list(tol=1e-12, stoptype="maxtime", maxtime=1), y=y)
res2
res2 <- turboem(par = c(0.9, 1, 3), fixptfn=fixptfn, objfn=objfn, boundary=boundary, pconstr=pconstr, method=c("em", "squarem", "squarem", "decme", "qn", "qn"), control.method=list(list(), list(K=2), list(K=3), list(version="v1"), list(qn=1), list(qn=2)), control.run=list(tol=1e-12, stoptype="maxtime", maxtime=1), y=y)
res2
res2 <- turboem(par = c(0.5, 1, 3), fixptfn=fixptfn, objfn=objfn, boundary=boundary, pconstr=pconstr, method=c("em", "squarem", "squarem", "decme", "qn", "qn"), control.method=list(list(), list(K=2), list(K=3), list(version="v1"), list(qn=1), list(qn=2)), control.run=list(tol=1e-12, stoptype="maxtime", maxtime=1), y=y)
res2
q("no")
Sweave("turboEM.Rnw")
Sweave("turboEM.Rnw")
Sweave("turboEM.Rnw")
Sweave("turboEM.Rnw")
Sweave("turboEM.Rnw")
Sweave("turboEM.Rnw")
Sweave("turboEM.Rnw")
Sweave("turboEM.Rnw")
Sweave("turboEM.Rnw")
Sweave("turboEM.Rnw")
Sweave("turboEM.Rnw")
Sweave("turboEM.Rnw")
Sweave("turboEM.Rnw")
Sweave("turboEM.Rnw")
Sweave("turboEM.Rnw")
Sweave("turboEM.Rnw")
Sweave("turboEM.Rnw")
Sweave("turboEM.Rnw")
Sweave("turboEM.Rnw")
q("no")
Sweave("turboEM.Rnw")
q("no")
