[ESS-bugs] ESS 12.09-2: *R* buffer font-locking
Rodney Sparapani
rsparapa at mcw.edu
Wed Feb 13 17:03:50 CET 2013
Here's something that I have been noticing recently. Check out the
attached graphic. Notice that there is no fontlocking for rows
1, 8, 9, 13 and 14. I don't see what the pattern is. You can
generate this output as follows:
require(DPpackage)
####################################
# Neal demo example
####################################
# Data
neal <- as.data.frame(x=c(-1.48,-1.40,-1.16,-1.08,-1.02, 0.14, 0.51,
0.53, 0.78))
dimnames(neal)[[2]] <- "neal"
attach(neal)
# Initial state
state <- NULL
# MCMC parameters
nburn <- 100
nsave <- 1000
nskip <- 0
ndisplay <- 1000
mcmc <- list(nburn=nburn, nsave=nsave, nskip=nskip, ndisplay=ndisplay)
prior <- list(alpha=1, m1=0, psiinv1=1, nu1=3, k0=0.01)
fit <- DPdensity(y=neal, prior=prior, mcmc=mcmc, state=state, status=TRUE)
# Posterior means
fit
# Plot the estimated density
plot(fit, ask=FALSE)
# Extracting the density estimate
# funky output comes from this statement
cbind(fit$x1, fit$dens)
--
Rodney Sparapani, PhD Center for Patient Care and Outcomes Research
Sr. Biostatistician http://www.mcw.edu/pcor
4 wheels good, 2 wheels better! Medical College of Wisconsin (MCW)
WWLD?: What Would Lombardi Do? Milwaukee, WI, USA
More information about the ESS-bugs
mailing list