SCodebits42013-03-08 10:28:25localhostconverted to 1.6 markup32007-07-30 17:15:43IanNimmoSmith22007-07-30 17:09:31IanNimmoSmith12007-07-30 17:08:06IanNimmoSmithBits of SPlus Code for analyses in 2003print.summary.gls 1) {
corr <- x$corBeta
class(corr) <- "correlation"
print(corr, title = "\n Correlation:", ...)
}
cat("\nStandardized residuals:\n")
print(x$residuals)
cat("\n")
cat("Residual standard error:", format(x$sigma), "\n")
cat("Degrees of freedom:", dd[["N"]], "total;", dd[["N"]] - dd[["p"]], "residual\n")
}]]>dwstat DW[lag, ]))/reps
}
}
else {
for(lag in 1:max.lag) {
p[lag] <- (sum(dw[lag] < DW[lag, ]))/reps
}
}
result <- list(r = r, dw = dw, p = p, alternative =
alternative)
class(result) <- "durbin.watson"
result
}
}]]>durbin.watson.lm DW[lag, ]))/reps
}
}
else {
for(lag in 1:max.lag) {
p[lag] <- (sum(dw[lag] < DW[lag, ]))/reps
}
}
result <- list(r = r, dw = dw, p = p, alternative =
alternative)
class(result) <- "durbin.watson"
result
}
}]]>