From 14a5af60a36071f6a9b4471fdf183fd91e89e1cd Mon Sep 17 00:00:00 2001
From: ranke <ranke@5fad18fb-23f0-0310-ab10-e59a3bee62b4>
Date: Mon, 1 Oct 2007 19:44:04 +0000
Subject: Moved everything into the trunk directory, in order to enable
 branching

git-svn-id: http://kriemhild.uft.uni-bremen.de/svn/chemCal@22 5fad18fb-23f0-0310-ab10-e59a3bee62b4
---
 tests/massart97.Rout.save | 123 ----------------------------------------------
 1 file changed, 123 deletions(-)
 delete mode 100644 tests/massart97.Rout.save

(limited to 'tests/massart97.Rout.save')

diff --git a/tests/massart97.Rout.save b/tests/massart97.Rout.save
deleted file mode 100644
index cb113d0..0000000
--- a/tests/massart97.Rout.save
+++ /dev/null
@@ -1,123 +0,0 @@
-
-R : Copyright 2006, The R Foundation for Statistical Computing
-Version 2.3.1 (2006-06-01)
-ISBN 3-900051-07-0
-
-R is free software and comes with ABSOLUTELY NO WARRANTY.
-You are welcome to redistribute it under certain conditions.
-Type 'license()' or 'licence()' for distribution details.
-
-R is a collaborative project with many contributors.
-Type 'contributors()' for more information and
-'citation()' on how to cite R or R packages in publications.
-
-Type 'demo()' for some demos, 'help()' for on-line help, or
-'help.start()' for an HTML browser interface to help.
-Type 'q()' to quit R.
-
-> require(chemCal)
-Loading required package: chemCal
-[1] TRUE
-> data(massart97ex1)
-> m <- lm(y ~ x, data = massart97ex1)
-> inverse.predict(m, 15)        #  6.1 +- 4.9
-$Prediction
-[1] 6.09381
-
-$`Standard Error`
-[1] 1.767278
-
-$Confidence
-[1] 4.906751
-
-$`Confidence Limits`
-[1]  1.187059 11.000561
-
-> inverse.predict(m, 90)        # 43.9 +- 4.9
-$Prediction
-[1] 43.93983
-
-$`Standard Error`
-[1] 1.767747
-
-$Confidence
-[1] 4.908053
-
-$`Confidence Limits`
-[1] 39.03178 48.84788
-
-> inverse.predict(m, rep(90,5)) # 43.9 +- 3.2
-$Prediction
-[1] 43.93983
-
-$`Standard Error`
-[1] 1.141204
-
-$Confidence
-[1] 3.168489
-
-$`Confidence Limits`
-[1] 40.77134 47.10832
-
-> 
-> data(massart97ex3)
-> attach(massart97ex3)
-> yx <- split(y, x)
-> ybar <- sapply(yx, mean)
-> s <- round(sapply(yx, sd), digits = 2)
-> w <- round(1 / (s^2), digits = 3)
-> weights <- w[factor(x)]
-> m <- lm(y ~ x, w = weights)
-> #calplot(m)
-> 
-> inverse.predict(m, 15, ws = 1.67)  # 5.9 +- 2.5
-$Prediction
-[1] 5.865367
-
-$`Standard Error`
-[1] 0.892611
-
-$Confidence
-[1] 2.478285
-
-$`Confidence Limits`
-[1] 3.387082 8.343652
-
-> inverse.predict(m, 90, ws = 0.145) # 44.1 +- 7.9
-$Prediction
-[1] 44.06025
-
-$`Standard Error`
-[1] 2.829162
-
-$Confidence
-[1] 7.855012
-
-$`Confidence Limits`
-[1] 36.20523 51.91526
-
-> 
-> m0 <- lm(y ~ x) 
-> lod(m0)
-$x
-[1] 5.406637
-
-$y
-[1] 13.63822
-
-> 
-> loq(m0)
-$x
-[1] 13.97767
-
-$y
-[1] 30.62355
-
-> loq(m, w.loq = 1.67)
-$x
-[1] 7.346231
-
-$y
-[1] 17.90784
-
-> 
-- 
cgit v1.2.1