Displaying 1 result from an estimated 1 matches for "rdelta".
Did you mean:
delta
2002 Nov 16
0
Wishlist: allow NA values in medpolish (package 'eda') (PR#2298)
...=0.01, maxiter=10, trace.iter = TRUE)
+medpolish <- function (x, eps=0.01, maxiter=10, trace.iter = TRUE,
+ na.rm = FALSE)
{
z <- as.matrix(x)
nr <- nrow(z)
@@ -8,19 +9,19 @@
c <- numeric(nc)
oldsum <- 0
for(iter in 1:maxiter) {
- rdelta <- apply(z, 1, median)
+ rdelta <- apply(z, 1, median, na.rm = na.rm)
z <- z - matrix(rdelta, nr=nr, nc=nc)
r <- r + rdelta
- delta <- median(c)
+ delta <- median(c, na.rm = na.rm)
c <- c - delta
t <- t + delta
- cdelta...