version 1.4-2

DESCRIPTION

@@ -1,15 +1,19 @@
 Package: polyapost
-Version: 1.1-6
-Date: 2014-04-21
-Imports: boot
-Title: Simulating from the Polya posterior
+Version: 1.4-2
+Date: 2015-07-07
+Imports: boot, stats
+Depends: R (>= 3.0.0), rcdd (>= 1.1-5)
+Title: Simulating from the Polya Posterior
 Author: Glen Meeden <glen@stat.umn.edu> and Radu Lazar
-        <lazar@stat.umn.edu>
+        <lazar@stat.umn.edu> and Charles J. Geyer <charlie@stat.umn.edu>
 Maintainer: Glen Meeden <glen@stat.umn.edu>
-Description: Generate dependent samples from a non-full dimensional
-        polytope via a Markov Chain sampler
+ByteCompile: TRUE
+Description: Simulate via Markov chain Monte Carlo (hit-and-run algorithm)
+    a Dirichlet distribution conditioned to satisfy a finite set of linear
+    equality and inequality constraints (hence to lie in a convex polytope
+    that is a subset of the unit simplex).
 License: GPL (>= 2)
-Packaged: 2014-04-21 20:45:00 UTC; gmeeden
 NeedsCompilation: yes
+Packaged: 2015-07-07 22:44:18 UTC; geyer
 Repository: CRAN
-Date/Publication: 2014-04-22 07:25:57
+Date/Publication: 2015-07-08 10:17:05

MD5

@@ -1,30 +1,38 @@
-eadf626f64cd4373588b71e1a14de944 *DESCRIPTION
-c76df6367ec5ce07212d916fcc833e16 *NAMESPACE
+f371f7ae4cf82bcbf975ccb2f47225ab *DESCRIPTION
+b7608e4bdf9470a5d33898041a55e8e3 *NAMESPACE
 af205ca09fab4fe4f460b9c4348eb85d *R/checkconstr.R
 a24b9132c8639e5deaee307bc7fe9e23 *R/constrmat.R
 8d02e6e9a42bcb8f7ab982cedd2626c1 *R/constrppmn.R
 011a1e7a2adcf3cfb491295bd5dbe7d9 *R/constrppprob.R
-8af116268c5ae7cdd51ffcf32a488be7 *R/cwpolya.R
+05b5b6f481312e763815463282aac5da *R/cwpolya.R
 600fb1c902b0c084e6ccd9aaeea3ff7a *R/feasible.R
 c6bdabbb7c3e74638991f0aa7d0f8988 *R/feasible1.R
-e90a2fe8385423ec55bb745fe57a3dfb *R/means.R
+abce94d7aeaad527b76013a71b6a2340 *R/hitrun.R
+945cc6bf2ed3409bc123f02f5bda1479 *R/means.R
 8d8ddb268afce5ec9658b67bf362220f *R/nullspace.R
 bde8037ae45667692d2a889433a73e31 *R/polyaest.R
 a1a00e355ceb163ef53365bfca339836 *R/polyap.R
 2be02f772609a4f9913e4bf2a30c876d *R/probvect.R
-d2d0660f39acf21c1149c4945c4c8b48 *R/probvect1.R
+7901d2bcd10051987b69652cb5a12e05 *R/probvect1.R
 be39dcf61a777387307071ea97d954aa *R/wtpolyap.R
-d7f807e1d748648a1764c6a866aed137 *build/vignette.rds
-8bd03911b62d364e44db6c2396d2b524 *inst/doc/pp1.R
-33d591a767d4e338d00e677dedba1b02 *inst/doc/pp1.Rnw
-60cf84f029be4a5325c1528e9e60e3ec *inst/doc/pp1.pdf
+ad77e0565ab48ba72848c6b3636f1514 *build/vignette.rds
+a174ddebf8949bf3a09ca869b1d23db8 *inst/doc/pp1.R
+c7ee19bc09131a6685f27bcbd18e4f37 *inst/doc/pp1.Rnw
+4f30d7e099051a7eaa8105f77b515526 *inst/doc/pp1.pdf
 3b716132aa2bd8685f02914ffd67b2c1 *man/constrppmn.Rd
-bfe800f8e58196f92238b7c95ea743a1 *man/constrppprob.Rd
-5b6324e86f7901e956c84d1c78909919 *man/feasible.Rd
+dd16907048de6ac08d118282acb818ee *man/constrppprob.Rd
+9a9e1dfb73af54d6522d3901e9c7d48a *man/feasible.Rd
+9476a048a7f5aa18b5130c551cff0689 *man/hitrun.Rd
 70a66177a47ae41060719beac631cc0d *man/polyap.Rd
 f3306aed9f43d4ce80ac07816b69a2ed *man/wtpolyap.Rd
-acc74c55459b9a9eeec62c095169c1cb *src/cwpolya.c
-d067003129a24718720d6524880cdf08 *src/means.c
-3210d9b6a13447a22754e5819e7cb2dc *src/probvect.c
-33d591a767d4e338d00e677dedba1b02 *vignettes/pp1.Rnw
-4388f62ae698c61d4700c06f5e555ecc *vignettes/ref.bib
+2fa4c7011c2bc0f7449ae151d5cc44ae *src/Makevars
+03ed829b5f3b2a4ca099d56f5a5216e0 *src/cwpolya.c
+9b7bfbd067a6f235f3c683b8430b76a6 *src/hitrun.c
+83c75ee506da5303afe8842c5896b970 *src/init.c
+8de2f794828c82d8e868a9dc861a0e7d *src/means.c
+d09330139cdab8b048fab14f91fe61e2 *src/polyapost.h
+949b366c817a7df2fc83d142efe9d62d *src/probvect.c
+6c8b99f22810226994b882e98443efe3 *tests/hitrun.R
+8bc7f4c14f434edd4175ead3b450ddd0 *tests/hitrun.Rout.save
+c7ee19bc09131a6685f27bcbd18e4f37 *vignettes/pp1.Rnw
+98f6b1389b95534b6e22db7c8d213d8e *vignettes/ref.bib

NAMESPACE

@@ -1,11 +1,19 @@
-useDynLib(polyapost)
+
+# Refer to all C/Fortran routines by their name prefixed by C_
+useDynLib(polyapost, .registration = TRUE, .fixes = "C_")
 
 importFrom(boot, simplex)
+import(rcdd)
+importFrom(stats, runif)
+importFrom(stats, quantile)
 
 export(polyap)
 export(wtpolyap)
 export(feasible)
 export(constrppmn)
 export(constrppprob)
 
+export(hitrun)
+S3method(hitrun, hitrun)
+S3method(hitrun, default)
 

R/cwpolya.R

@@ -8,11 +8,10 @@ cwpolya <- function(samp, wts, k) {
 	k <- as.integer(k)
 	if (k <= 0)
 		stop("k must be positive integer")
-       out<-.C("cwpolya",
+       out<-.C(C_cwpolya,
        x=as.double(c(samp, rep(0, k))),
        w=as.double(cumsum(wts)),
        n=as.integer(length(samp)),
-       k=as.integer(k),
-       PACKAGE="polyapost")
+       k=as.integer(k))
        return(out$x)
 }