pointwise(result.predict, coverage = 0.99)
attach(gas)
gas.m <- loess(NOx ~ E, data = gas, span = 2/3)
gas.se <- predict(gas.m, newdata = seq(min(E), max(E),length = 7),
se.fit = TRUE)
pointwise(gas.se)
# Gives the following result:
$upper:
[1] 1.985621 4.109807 5.480230 5.566510 3.527610 1.710617
[7] 1.472049
$fit:
[1] 1.1964144 3.6794968 5.0557086 5.1352603 3.1436568 1.1969317
[7] 0.5236823
$lower:
[1] 0.4072080 3.2491865 4.6311876 4.7040105 2.7597037
[6] 0.6832464 -0.4246841