Displaying 1 result from an estimated 1 matches for "d1xtil".
2013 Jan 18
0
problem that arises after using the new version of "BRugs"
...ic(names(table(x)))
cnts <- as.numeric(table(x))
nTil <- length(xTil)
deltaTil <- rep(NA, nTil)
for (iTil in 1:nTil) deltaTil[iTil] <- sum(delta[(1:n)[x ==
xTil[iTil]]])
sd.xTil <- sd(xTil)
xTil <- xTil/sd.xTil
range.x <- range.x/sd.xTil
d1xTil <- xTil[2:nTil] - xTil[1:(nTil - 1)]
d2xTil <- xTil[3:nTil] - xTil[1:(nTil - 2)]
rcsrCnts <- rev(cumsum(rev(cnts)))
trapLef <- 2 * cnts[1] * xTil[1] + (sum(cnts) - cnts[1]) *
(xTil[2] + xTil[1])
trapMid <- cnts[2:(nTil - 1)] * d1xTil[-length(d1xTil)] +...