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
|
R : Copyright 2006, The R Foundation for Statistical Computing
Version 2.3.1 beta (2006-05-18 r38119)
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...
Error in prof$getProfile() : number of iterations exceeded maximum of 50
Pyridin: Fitting data...
PyNO: Fitting data...
(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...
Error in nls(response ~ linlogitf(dose, 1, f, logED50, b), data = tmp, :
step factor 0.000488281 reduced below 'minFactor' of 0.000976562
Waiting for profiling to be done...
Error in numericDeriv(form[[3]], names(ind), env) :
Missing value or an infinity produced when evaluating the model
In addition: Warning message:
NaNs produced in: pnorm(q, mean, sd, lower.tail, log.p)
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...
(PyS)2: Fitting data...
Error in nls(response ~ linlogitf(dose, 1, f, logED50, b), data = tmp, :
singular gradient
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=5)
Substance ndl n lld lhd mtype logED50 2.5% 97.5%
1 Na Pyrithion 20 108 -2.10721 2.0000 linlogit -0.34613 NA NA
2 Pyridin 19 161 0.59176 3.6990 inactive NA NA NA
3 PyNO 19 81 0.59176 3.6990 no fit NA NA NA
4 (PT)2 19 81 -2.40824 2.0000 linlogit -0.41970 -0.53261 NA
5 MSO2P 28 108 -1.70927 3.6990 inactive NA NA NA
6 MSPHI 19 81 0.59176 3.6990 probit 3.66250 3.55483 3.880532
7 PyS 18 80 -0.40824 3.0000 probit 2.68179 NA NA
8 Zn Pyrithion 27 81 -2.10721 2.0000 linlogit -0.41329 -0.52933 -0.299798
9 Cu Pyrithion 19 79 -2.40824 2.0000 linlogit -0.30747 -0.47535 -0.154428
10 Fe Pyrithion 19 81 -2.40824 2.0000 linlogit -0.35346 -0.50038 -0.214596
11 (PyS)2 20 81 -0.40824 3.0000 probit 1.78405 1.67268 1.888245
12 MSPT 18 108 -0.40824 3.0000 probit 2.13848 2.07093 2.205115
13 TBT 38 135 -2.70927 2.3979 linlogit -0.15810 -0.26021 -0.050455
14 NaJ 10 108 0.59176 3.3010 inactive NA NA NA
unit sigma a b c
1 然 0.208578 -0.34613 2.11883 -0.266989
2 然 NA NA NA NA
3 然 0.091889 NA NA NA
4 然 0.205091 -0.41970 1.84673 -0.515056
5 然 NA NA NA NA
6 然 0.129433 3.66250 0.43875 NA
7 然 0.134019 2.68179 0.12505 NA
8 然 0.228779 -0.41329 1.74248 0.610790
9 然 0.244945 -0.30747 1.01507 -0.048974
10 然 0.233832 -0.35346 1.18236 0.025107
11 然 0.232987 1.78405 0.40980 NA
12 然 0.228813 2.13848 0.21126 NA
13 然 0.192641 -0.15810 1.04066 -0.024161
14 然 NA NA NA NA
>
|
