lcd.marg {LogConcDEAD} | R Documentation |
This function numerically integrates the log MLE of 2-d data to
produce axis-aligned marginals for use in plot.LogConcDEAD
.
lcd.marg (out, gridlen=100) lcd.marg2 (out, gridlen=100)
out |
Output from lcd.mle (of class LogConcDEAD) |
gridlen |
Size of grid |
Does a (currently rather crude) approximation to the integral of the
log MLE, for the purpose of plotting axis-aligned marginals. The
integrate
routine is used to perform numerical integration.
Currently use lcd.marg
for X_1 and lcd.marg2
for
X_2.
x |
x-values of grid (for X_1 axis) |
y |
y-values of grid (for X_2 axis) |
marg |
Integrated MLE |
Madeleine Cule mlc40@cam.ac.uk
Robert Gramacy
Richard Samworth
#simple normal data, a couple of plots set.seed(101) x <- matrix(rnorm(200),ncol=2) out <- lcd.mle(x) m1 <- lcd.marg(out) m2 <- lcd.marg2(out) par(mfrow=c(1,2)) plot(out,marg=1,g.marg=m1) plot(out,marg=2,g.marg=m2)