```
# NOT RUN {
fit <- kppm(redwood ~ x + y)
vcov(fit)
vcov(fit, what="corr")
# confidence interval
confint(fit)
# cross-check the confidence interval by hand:
sd <- sqrt(diag(vcov(fit)))
t(coef(fit) + 1.96 * outer(sd, c(lower=-1, upper=1)))
# }
```

Run the code above in your browser using DataLab