Displaying 1 result from an estimated 1 matches for "scov".
Did you mean:
scev
2011 Apr 20
3
Help needed!
...th.spline(time, x)
x <- x - fitted(fit)
}
gg <- NULL
for (zz in unique(subject)) {
if (sum(subject == zz) > 1) {
tt <- time[subject == zz]
xx <- x[subject == zz]
g <- expand.grid(t1 = tt, t2 = tt)
scov <- xx %*% t(xx)
if (noDiag)
scov <- scov + diag(rep(Inf, length(xx)))
g$z <- matrix(scov, ncol = 1)
gg <- rbind(gg, g[g$z < Inf, ])
}
}
nobs <- nrow(gg)
tt <- min(time) + (max(time) - min(time)) * (...