From f381f9a6a8a47b89ec25cd627833a7248da7932b Mon Sep 17 00:00:00 2001
From: ranke <ranke@5fad18fb-23f0-0310-ab10-e59a3bee62b4>
Date: Tue, 23 May 2006 07:33:22 +0000
Subject: Don't do calplot and lod for linear models from weighted regression
 any more, since this is not supported (PR#8877).

git-svn-id: http://kriemhild.uft.uni-bremen.de/svn/chemCal@13 5fad18fb-23f0-0310-ab10-e59a3bee62b4
---
 tests/massart97.Rout.save | 54 +++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 54 insertions(+)
 create mode 100644 tests/massart97.Rout.save

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

diff --git a/tests/massart97.Rout.save b/tests/massart97.Rout.save
new file mode 100644
index 0000000..ae50275
--- /dev/null
+++ b/tests/massart97.Rout.save
@@ -0,0 +1,54 @@
+
+R : Copyright 2006, The R Foundation for Statistical Computing
+Version 2.3.0 (2006-04-24)
+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.
+
+> library(chemCal)
+> 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)
+> # The following concords with the book
+> inverse.predict(m, 15, ws = 1.67)
+$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)
+$Prediction
+[1] 44.06025
+
+$`Standard Error`
+[1] 2.829162
+
+$Confidence
+[1] 7.855012
+
+$`Confidence Limits`
+[1] 36.20523 51.91526
+
+> 
-- 
cgit v1.2.1