Displaying 2 results from an estimated 2 matches for "cimedian".
2001 Oct 23
0
problems with postscript device
...oblems with other functions:
boot.stat <-
function (x, B = 1000, norm = FALSE)
{
xname <- deparse(substitute(x))
means <- meds <- numeric(B)
for (i in 1:B) {
s <- sample(x, replace = TRUE)
means[i] <- mean(s)
meds[i] <- median(s)
}
CImedian <- quantile(meds, c(0.025, 0.5, 0.975))
CImean <- quantile(means, c(0.025, 0.5, 0.975))
hmean <- hist(means, freq = FALSE, plot = FALSE)
hmeds <- hist(meds, freq = FALSE, plot = FALSE)
split.screen(c(1, 2))
xlims <- range(c(hmean$breaks, hmeds$breaks))
ymax &l...
2001 Oct 26
1
postscript problem (PR#1147)
...n its head!
Kjetil Halvorsen
boot.stat <-
function (x, B = 1000, norm = FALSE)
{
xname <- deparse(substitute(x))
means <- meds <- numeric(B)
for (i in 1:B) {
s <- sample(x, replace = TRUE)
means[i] <- mean(s)
meds[i] <- median(s)
}
CImedian <- quantile(meds, c(0.025, 0.5, 0.975))
CImean <- quantile(means, c(0.025, 0.5, 0.975))
hmean <- hist(means, freq = FALSE, plot = FALSE)
hmeds <- hist(meds, freq = FALSE, plot = FALSE)
split.screen(c(1, 2))
xlims <- range(c(hmean$breaks, hmeds$breaks))
ymax &l...