aboutsummaryrefslogtreecommitdiff
path: root/tests
diff options
context:
space:
mode:
Diffstat (limited to 'tests')
-rw-r--r--tests/IM1xIPC81.R2
-rw-r--r--tests/IM1xIPC81.Rout.save54
-rw-r--r--tests/IM1xVibrio.R2
-rw-r--r--tests/IM1xVibrio.Rout.save46
-rw-r--r--tests/pyrithione.R4
-rw-r--r--tests/pyrithione.Rout.save72
6 files changed, 98 insertions, 82 deletions
diff --git a/tests/IM1xIPC81.R b/tests/IM1xIPC81.R
index 52b5653..475fcae 100644
--- a/tests/IM1xIPC81.R
+++ b/tests/IM1xIPC81.R
@@ -1,4 +1,4 @@
library(drfit)
data(IM1xIPC81)
rIM1xIPC81 <- drfit(IM1xIPC81,linlogit=TRUE)
-print(rIM1xIPC81,digits=5)
+print(rIM1xIPC81,digits=4)
diff --git a/tests/IM1xIPC81.Rout.save b/tests/IM1xIPC81.Rout.save
index db37e67..203ac66 100644
--- a/tests/IM1xIPC81.Rout.save
+++ b/tests/IM1xIPC81.Rout.save
@@ -1,7 +1,8 @@
-R version 2.9.2 (2009-08-24)
-Copyright (C) 2009 The R Foundation for Statistical Computing
+R version 2.15.1 (2012-06-22) -- "Roasted Marshmallows"
+Copyright (C) 2012 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
+Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
@@ -31,7 +32,7 @@ Error in nls(response ~ linlogitf(dose, 1, f, logED50, b), data = tmp, :
Error in numericDeriv(form[[3L]], names(ind), env, ifelse(internalPars < :
Missing value or an infinity produced when evaluating the model
In addition: Warning message:
-In pnorm(q, mean, sd, lower.tail, log.p) : NaNs produced
+In pnorm(-log10(dose), -logED50, scale) : NaNs produced
IM15 BF4: Fitting data...
@@ -45,7 +46,7 @@ Waiting for profiling to be done...
Error in numericDeriv(form[[3L]], names(ind), env, ifelse(internalPars < :
Missing value or an infinity produced when evaluating the model
In addition: Warning message:
-In pnorm(q, mean, sd, lower.tail, log.p) : NaNs produced
+In pnorm(-log10(dose), -logED50, scale) : NaNs produced
IM18 BF4: Fitting data...
@@ -53,7 +54,7 @@ Waiting for profiling to be done...
Error in numericDeriv(form[[3L]], names(ind), env, ifelse(internalPars < :
Missing value or an infinity produced when evaluating the model
In addition: Warning message:
-In pnorm(q, mean, sd, lower.tail, log.p) : NaNs produced
+In pnorm(-log10(dose), -logED50, scale) : NaNs produced
IM19 BF4: Fitting data...
@@ -61,28 +62,31 @@ Waiting for profiling to be done...
Error in numericDeriv(form[[3L]], names(ind), env, ifelse(internalPars < :
Missing value or an infinity produced when evaluating the model
In addition: Warning message:
-In pnorm(q, mean, sd, lower.tail, log.p) : NaNs produced
+In pnorm(-log10(dose), -logED50, scale) : NaNs produced
IM1-10 BF4: Fitting data...
Waiting for profiling to be done...
-> print(rIM1xIPC81,digits=5)
- Substance ndl n lld lhd mtype logED50 2.5% 97.5% unit
-1 IM13 BF4 9 81 0.591760 3.0000 inactive NA NA NA microM
-2 IM14 BF4 20 216 -0.010300 3.1761 no fit NA NA NA microM
-3 IM15 BF4 9 135 0.591760 3.0000 inactive NA NA NA microM
-4 IM16 BF4 9 108 0.591760 3.0000 inactive NA NA NA microM
-5 IM17 BF4 9 81 0.591760 3.0000 linlogit 2.5786 2.50572 2.66174 microM
-6 IM18 BF4 9 135 0.591760 3.0000 linlogit 1.6806 1.62297 1.74189 microM
-7 IM19 BF4 9 81 0.591760 3.0000 linlogit 1.6496 1.59848 1.70315 microM
-8 IM1-10 BF4 11 162 -0.010300 3.0000 linlogit 0.7697 0.68701 0.85442 microM
- sigma a b c
-1 NA NA NA NA
-2 NA NA NA NA
-3 NA NA NA NA
-4 NA NA NA NA
-5 0.23759 2.5786 2.3002 0.014680
-6 0.23249 1.6806 2.2373 0.057193
-7 0.14527 1.6496 1.9770 0.109561
-8 0.29877 0.7697 1.9363 0.458090
+> print(rIM1xIPC81,digits=4)
+ Substance ndl n lld lhd mtype logED50 2.5% 97.5% unit sigma
+1 IM13 BF4 9 81 0.5918 3.000 inactive NA NA NA microM NA
+2 IM14 BF4 20 216 -0.0103 3.176 no fit NA NA NA microM NA
+3 IM15 BF4 9 135 0.5918 3.000 inactive NA NA NA microM NA
+4 IM16 BF4 9 108 0.5918 3.000 inactive NA NA NA microM NA
+5 IM17 BF4 9 81 0.5918 3.000 linlogit 2.5786 2.506 2.6617 microM 0.2376
+6 IM18 BF4 9 135 0.5918 3.000 linlogit 1.6806 1.623 1.7419 microM 0.2325
+7 IM19 BF4 9 81 0.5918 3.000 linlogit 1.6496 1.598 1.7031 microM 0.1453
+8 IM1-10 BF4 11 162 -0.0103 3.000 linlogit 0.7697 0.687 0.8544 microM 0.2988
+ a b c
+1 NA NA NA
+2 NA NA NA
+3 NA NA NA
+4 NA NA NA
+5 2.5786 2.300 0.01468
+6 1.6806 2.237 0.05719
+7 1.6496 1.977 0.10956
+8 0.7697 1.936 0.45809
>
+> proc.time()
+ user system elapsed
+ 0.852 0.024 0.869
diff --git a/tests/IM1xVibrio.R b/tests/IM1xVibrio.R
index 8fb23d2..83aa878 100644
--- a/tests/IM1xVibrio.R
+++ b/tests/IM1xVibrio.R
@@ -1,4 +1,4 @@
library(drfit)
data(IM1xVibrio)
rIM1xVibrio <- drfit(IM1xVibrio)
-print(rIM1xVibrio,digits=5)
+print(rIM1xVibrio,digits=4)
diff --git a/tests/IM1xVibrio.Rout.save b/tests/IM1xVibrio.Rout.save
index ccfdc48..e54ba59 100644
--- a/tests/IM1xVibrio.Rout.save
+++ b/tests/IM1xVibrio.Rout.save
@@ -1,7 +1,8 @@
-R version 2.9.2 (2009-08-24)
-Copyright (C) 2009 The R Foundation for Statistical Computing
+R version 2.15.1 (2012-06-22) -- "Roasted Marshmallows"
+Copyright (C) 2012 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
+Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
@@ -52,23 +53,26 @@ Waiting for profiling to be done...
IM1-10 BF4: Fitting data...
Waiting for profiling to be done...
-> print(rIM1xVibrio,digits=5)
- Substance ndl n lld lhd mtype logED50 2.5% 97.5% unit
-1 IM13 BF4 13 22 -5.30103 4.6990 probit 3.93991 3.86362 4.01601 microM
-2 IM14 BF4 12 57 -4.30103 4.6990 probit 3.54421 3.50302 3.58555 microM
-3 IM15 BF4 7 20 -0.30103 4.6990 probit 3.13985 3.11359 3.16624 microM
-4 IM16 BF4 13 30 -4.30103 4.6990 probit 3.18042 3.14358 3.21762 microM
-5 IM17 BF4 10 20 -1.30103 4.6990 probit 2.43173 2.36807 2.49512 microM
-6 IM18 BF4 10 17 -2.30103 2.6990 probit 1.40152 1.31467 1.48806 microM
-7 IM19 BF4 13 22 -5.30103 2.6990 probit 0.71579 0.67038 0.75925 microM
-8 IM1-10 BF4 11 20 -6.30103 2.6990 probit -0.17903 -0.25694 -0.10372 microM
- sigma a b
-1 0.063029 3.93991 0.57631
-2 0.059054 3.54421 0.65613
-3 0.026955 3.13985 0.55267
-4 0.045141 3.18042 0.48816
-5 0.047733 2.43173 0.63865
-6 0.056113 1.40152 0.80085
-7 0.029563 0.71579 0.69693
-8 0.040987 -0.17903 0.80194
+> print(rIM1xVibrio,digits=4)
+ Substance ndl n lld lhd mtype logED50 2.5% 97.5% unit sigma
+1 IM13 BF4 13 22 -5.301 4.699 probit 3.9399 3.8636 4.0160 microM 0.06303
+2 IM14 BF4 12 57 -4.301 4.699 probit 3.5442 3.5030 3.5855 microM 0.05905
+3 IM15 BF4 7 20 -0.301 4.699 probit 3.1398 3.1136 3.1662 microM 0.02695
+4 IM16 BF4 13 30 -4.301 4.699 probit 3.1804 3.1436 3.2176 microM 0.04514
+5 IM17 BF4 10 20 -1.301 4.699 probit 2.4317 2.3681 2.4951 microM 0.04773
+6 IM18 BF4 10 17 -2.301 2.699 probit 1.4015 1.3147 1.4881 microM 0.05611
+7 IM19 BF4 13 22 -5.301 2.699 probit 0.7158 0.6704 0.7592 microM 0.02956
+8 IM1-10 BF4 11 20 -6.301 2.699 probit -0.1790 -0.2569 -0.1037 microM 0.04099
+ a b
+1 3.9399 0.5763
+2 3.5442 0.6561
+3 3.1398 0.5527
+4 3.1804 0.4882
+5 2.4317 0.6387
+6 1.4015 0.8009
+7 0.7158 0.6969
+8 -0.1790 0.8019
>
+> proc.time()
+ user system elapsed
+ 0.864 0.024 0.879
diff --git a/tests/pyrithione.R b/tests/pyrithione.R
index baec676..3a3ff23 100644
--- a/tests/pyrithione.R
+++ b/tests/pyrithione.R
@@ -1,4 +1,4 @@
library(drfit)
data(pyrithione)
-rpyr <- drfit(pyrithione,linlogit=TRUE,linlogitWrong=c("MSPT","MSPHI"))
-print(rpyr,digits=3)
+rpyr <- drfit(pyrithione,linlogit=TRUE,linlogitWrong=c("MSPT","MSPHI","PyS"))
+print(rpyr,digits=4)
diff --git a/tests/pyrithione.Rout.save b/tests/pyrithione.Rout.save
index af67552..eb92b80 100644
--- a/tests/pyrithione.Rout.save
+++ b/tests/pyrithione.Rout.save
@@ -1,7 +1,8 @@
-R version 2.9.2 (2009-08-24)
-Copyright (C) 2009 The R Foundation for Statistical Computing
+R version 2.15.1 (2012-06-22) -- "Roasted Marshmallows"
+Copyright (C) 2012 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
+Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
@@ -19,7 +20,7 @@ Type 'q()' to quit R.
Loading required package: MASS
Loading required package: RODBC
> data(pyrithione)
-> rpyr <- drfit(pyrithione,linlogit=TRUE,linlogitWrong=c("MSPT","MSPHI"))
+> rpyr <- drfit(pyrithione,linlogit=TRUE,linlogitWrong=c("MSPT","MSPHI","PyS"))
Na Pyrithion: Fitting data...
@@ -47,6 +48,10 @@ Waiting for profiling to be done...
PyS: Fitting data...
Waiting for profiling to be done...
+Error in numericDeriv(form[[3L]], names(ind), env, ifelse(internalPars < :
+ Missing value or an infinity produced when evaluating the model
+In addition: Warning message:
+In pnorm(-log10(dose), -logED50, scale) : NaNs produced
Zn Pyrithion: Fitting data...
@@ -70,33 +75,36 @@ Waiting for profiling to be done...
NaJ: Fitting data...
-> print(rpyr,digits=3)
- Substance ndl n lld lhd mtype logED50 2.5% 97.5% unit sigma
-1 Na Pyrithion 20 108 -2.107 2.0 linlogit -0.346 -0.428 -0.2648 microM 0.209
-2 Pyridin 19 161 0.592 3.7 inactive NA NA NA microM NA
-3 PyNO 19 81 0.592 3.7 no fit NA NA NA microM 0.107
-4 (PT)2 19 81 -2.408 2.0 linlogit -0.420 -0.533 NA microM 0.205
-5 MSO2P 28 108 -1.709 3.7 inactive NA NA NA microM NA
-6 MSPHI 19 81 0.592 3.7 probit 3.663 3.555 3.8802 microM 0.129
-7 PyS 18 80 -0.408 3.0 linlogit 2.756 NA 2.8228 microM 0.146
-8 Zn Pyrithion 27 81 -2.107 2.0 linlogit -0.413 -0.529 -0.2998 microM 0.229
-9 Cu Pyrithion 19 79 -2.408 2.0 linlogit -0.307 -0.475 -0.1544 microM 0.245
-10 Fe Pyrithion 19 81 -2.408 2.0 linlogit -0.353 -0.500 -0.2146 microM 0.234
-11 MSPT 18 108 -0.408 3.0 probit 2.138 2.071 2.2051 microM 0.229
-12 TBT 38 135 -2.709 2.4 linlogit -0.158 -0.260 -0.0505 microM 0.193
-13 NaJ 10 108 0.592 3.3 inactive NA NA NA microM NA
- a b c
-1 -0.346 2.119 -0.266998
-2 NA NA NA
-3 NA NA NA
-4 -0.420 1.847 -0.515071
-5 NA NA NA
-6 3.663 0.439 NA
-7 2.756 0.743 -0.000993
-8 -0.413 1.742 0.610900
-9 -0.307 1.015 -0.048966
-10 -0.353 1.182 0.025082
-11 2.138 0.211 NA
-12 -0.158 1.041 -0.024165
-13 NA NA NA
+> print(rpyr,digits=4)
+ Substance ndl n lld lhd mtype logED50 2.5% 97.5% unit
+1 Na Pyrithion 20 108 -2.1072 2.000 linlogit -0.3461 -0.4276 -0.26479 microM
+2 Pyridin 19 161 0.5918 3.699 inactive NA NA NA microM
+3 PyNO 19 81 0.5918 3.699 no fit NA NA NA microM
+4 (PT)2 19 81 -2.4082 2.000 linlogit -0.4197 -0.5326 NA microM
+5 MSO2P 28 108 -1.7093 3.699 inactive NA NA NA microM
+6 MSPHI 19 81 0.5918 3.699 probit 3.6625 3.5546 3.88019 microM
+7 PyS 18 80 -0.4082 3.000 probit 2.6818 NA NA microM
+8 Zn Pyrithion 27 81 -2.1072 2.000 linlogit -0.4133 -0.5293 -0.29980 microM
+9 Cu Pyrithion 19 79 -2.4082 2.000 linlogit -0.3075 -0.4753 -0.15444 microM
+10 Fe Pyrithion 19 81 -2.4082 2.000 linlogit -0.3535 -0.5004 -0.21460 microM
+11 MSPT 18 108 -0.4082 3.000 probit 2.1385 2.0709 2.20511 microM
+12 TBT 38 135 -2.7093 2.398 linlogit -0.1581 -0.2602 -0.05047 microM
+13 NaJ 10 108 0.5918 3.301 inactive NA NA NA microM
+ sigma a b c
+1 0.2086 -0.3461 2.1188 -0.26700
+2 NA NA NA NA
+3 0.1073 NA NA NA
+4 0.2051 -0.4197 1.8467 -0.51507
+5 NA NA NA NA
+6 0.1294 3.6625 0.4387 NA
+7 0.1340 2.6818 0.1250 NA
+8 0.2288 -0.4133 1.7425 0.61081
+9 0.2449 -0.3075 1.0151 -0.04897
+10 0.2338 -0.3535 1.1823 0.02508
+11 0.2288 2.1385 0.2113 NA
+12 0.1926 -0.1581 1.0407 -0.02416
+13 NA NA NA NA
>
+> proc.time()
+ user system elapsed
+ 1.556 0.044 1.595

Contact - Imprint