1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
|
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.
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(drfit)
Loading required package: MASS
Loading required package: RODBC
> data(pyrithione)
> rpyr <- drfit(pyrithione,linlogit=TRUE,linlogitWrong=c("MSPT","MSPHI","PyS"))
Na Pyrithion: Fitting data...
Waiting for profiling to be done...
Pyridin: Fitting data...
PyNO: Fitting data...
Error in nls(response ~ linlogitf(dose, 1, f, logED50, b), data = tmp, :
Convergence failure: false convergence (8)
(PT)2: Fitting data...
Waiting for profiling to be done...
MSO2P: Fitting data...
MSPHI: Fitting data...
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...
Waiting for profiling to be done...
Cu Pyrithion: Fitting data...
Waiting for profiling to be done...
Fe Pyrithion: Fitting data...
Waiting for profiling to be done...
MSPT: Fitting data...
Waiting for profiling to be done...
TBT: Fitting data...
Waiting for profiling to be done...
NaJ: Fitting data...
> 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
|
