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
|
R version 2.9.2 (2009-08-24)
Copyright (C) 2009 The R Foundation for Statistical Computing
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(drfit)
Loading required package: MASS
Loading required package: RODBC
> data(pyrithione)
> rpyr <- drfit(pyrithione,linlogit=TRUE,linlogitWrong=c("MSPT","MSPHI"))
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...
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=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
>
|
