<?
function preprocess () { ?>
require(xtable)
<?        }

function calculate () {
?>
n=5
results <- summary.lm (lm (<? getRK ("y"); ?> ~ <? getRK ("x"); ?>) )
quarlm<-matrix(quantile(residuals(results)),nrow=1, ncol=5)
colnames(quarlm)<-c("Min","1Q","Median","2Q","Max")
coeflm<-coef(results)
siglm<-results$sigma
maxdflm<-(max(results$df))
rsq<-results$r.squared
arsq<-results$adj.r.squared
fslm<-matrix(results$fstatistic,nrow=1,ncol=3)
colnames(fslm)<-names(results$fstatistic)
<?
}

function printout () {
?>
rk.header ("Linear Regression")
rk.print("Residuals")
rk.results(xtable(quarlm))
rk.print("Coefficients")
rk.results(xtable(coeflm))
rk.print(c("Residual standard error: ", round(siglm,n), " on ", maxdflm, " degrees of fredom"))
rk.print(c("Multiple R-Squared: ", round(rsq,n), ", Adjusted R-squared: ", round(arsq,n)))
rk.print("F-statistic")
rk.results(xtable(fslm))
<?
}
?>