From: Carnell, Rob C <CarnellR_at_BATTELLE.ORG>

Date: Sun 30 Jul 2006 - 23:42:29 GMT

R-devel@r-project.org mailing list

https://stat.ethz.ch/mailman/listinfo/r-devel Received on Mon Jul 31 09:46:05 2006

Date: Sun 30 Jul 2006 - 23:42:29 GMT

NIST maintains a repository of Statistical Reference Datasets at
http://www.itl.nist.gov/div898/strd/. I have been working through the
datasets to compare R's results to their references with the hope that
if all works well, this could become a validation package.

I am using R 2.3.1 for Windows XP.

rm(list=ls())

require(gsubfn)

data.base <- "http://www.itl.nist.gov/div898/strd/lls/data/LINKS/DATA"

reg.data <- paste(data.base, "/Filip.dat", sep="")

model <-

"V1~V2+I(V2^2)+I(V2^3)+I(V2^4)+I(V2^5)+I(V2^6)+I(V2^7)+I(V2^8)+I(V2^9)+I
(V2^10)"

filePath <- paste(defaultPath, "//NISTtest.dat", sep="") download.file(reg.data, filePath, quiet=TRUE)

A <- read.table(filePath, skip=60, strip.white=TRUE) lm.data <- lm(formula(model), A)

lm.data

Rob Carnell

R-devel@r-project.org mailing list

https://stat.ethz.ch/mailman/listinfo/r-devel Received on Mon Jul 31 09:46:05 2006

Archive maintained by Robert King, hosted by
the discipline of
statistics at the
University of Newcastle,
Australia.

Archive generated by hypermail 2.1.8, at Mon 31 Jul 2006 - 08:32:12 GMT.

*
Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help.
Please read the posting
guide before posting to the list.
*