[R] Don't know which model in"drc" package is to be used to find EC values.
Syamkumar R
syamkumar4 at gmail.com
Wed Feb 23 16:20:19 CET 2011
Hi every one,
I am using the package 'drc' to model root elongation using dose
response data. I don't know which model I should use. Though I don't
know which model I should use, I tried the following codes given
below. But it produced the error messages.Can any one tell me the code
in 'drc' package to find out the EC (Effective Concentration) values
and Confidence Intervals. The data is given below (also attached a
".csv" file).
dat
conc Length
1 1e-04 28.55
2 1e-04 14.47
3 1e-04 29.93
4 1e-04 31.66
5 1e-04 10.78
6 1e-04 16.83
7 1e-04 10.52
8 1e-04 31.44
9 1e-04 16.56
10 1e-04 10.30
11 1e-04 20.08
12 1e-04 14.04
13 1e-04 19.93
14 1e-04 14.94
15 1e-04 38.68
16 1e-04 37.92
17 1e-04 3.54
18 1e-04 13.58
19 1e-04 7.50
20 1e-04 12.81
21 1e-04 6.80
22 1e-04 19.61
23 1e-04 26.75
24 1e-04 20.28
25 1e-04 29.95
26 1e-04 12.90
27 1e-04 11.88
28 1e-04 18.72
29 1e-04 16.14
30 1e-04 0.00
31 1e-03 45.50
32 1e-03 45.58
33 1e-03 34.32
34 1e-03 33.88
35 1e-03 26.59
36 1e-03 21.91
37 1e-03 14.08
38 1e-03 18.28
39 1e-03 11.79
40 1e-03 19.93
41 1e-03 19.03
42 1e-03 19.06
43 1e-03 26.44
44 1e-03 34.35
45 1e-03 24.25
46 1e-03 37.12
47 1e-03 46.14
48 1e-03 45.30
49 1e-03 11.51
50 1e-03 43.21
51 1e-03 30.71
52 1e-03 24.23
53 1e-03 33.86
54 1e-03 21.19
55 1e-03 12.24
56 1e-03 30.82
57 1e-03 21.50
58 1e-03 29.79
59 1e-03 44.50
60 1e-03 0.00
61 1e-02 14.44
62 1e-02 26.84
63 1e-02 34.56
64 1e-02 13.19
65 1e-02 28.46
66 1e-02 21.45
67 1e-02 25.62
68 1e-02 27.87
69 1e-02 11.19
70 1e-02 0.00
71 1e-02 13.33
72 1e-02 26.38
73 1e-02 30.35
74 1e-02 11.02
75 1e-02 21.25
76 1e-02 25.31
77 1e-02 27.26
78 1e-02 34.75
79 1e-02 12.97
80 1e-02 0.00
81 1e-02 26.45
82 1e-02 27.07
83 1e-02 59.99
84 1e-02 27.85
85 1e-02 34.53
86 1e-02 43.07
87 1e-02 18.24
88 1e-02 28.16
89 1e-02 16.03
90 1e-02 0.00
91 1e-01 9.45
92 1e-01 5.04
93 1e-01 9.60
94 1e-01 0.00
95 1e-01 0.00
96 1e-01 0.00
97 1e-01 0.00
98 1e-01 0.00
99 1e-01 0.00
100 1e-01 0.00
101 1e-01 0.00
102 1e-01 0.00
103 1e-01 0.00
104 1e-01 0.00
105 1e-01 0.00
106 1e-01 0.00
107 1e-01 0.00
108 1e-01 0.00
109 1e-01 0.00
110 1e-01 0.00
111 1e-01 3.63
112 1e-01 7.25
113 1e-01 0.00
114 1e-01 0.00
115 1e-01 0.00
116 1e-01 0.00
117 1e-01 0.00
118 1e-01 0.00
119 1e-01 0.00
120 1e-01 0.00
121 1e+00 8.35
122 1e+00 5.40
123 1e+00 0.00
124 1e+00 0.00
125 1e+00 0.00
126 1e+00 0.00
127 1e+00 0.00
128 1e+00 0.00
129 1e+00 0.00
130 1e+00 0.00
131 1e+00 5.90
132 1e+00 0.00
133 1e+00 0.00
134 1e+00 0.00
135 1e+00 0.00
136 1e+00 0.00
137 1e+00 0.00
138 1e+00 0.00
139 1e+00 0.00
140 1e+00 0.00
141 1e+00 11.00
142 1e+00 0.00
143 1e+00 0.00
144 1e+00 0.00
145 1e+00 0.00
146 1e+00 0.00
147 1e+00 0.00
148 1e+00 0.00
149 1e+00 0.00
150 1e+00 0.00
151 1e+01 0.00
152 1e+01 0.00
153 1e+01 0.00
154 1e+01 0.00
155 1e+01 0.00
156 1e+01 0.00
157 1e+01 0.00
158 1e+01 0.00
159 1e+01 0.00
160 1e+01 0.00
161 1e+01 0.00
162 1e+01 0.00
163 1e+01 0.00
164 1e+01 0.00
165 1e+01 0.00
166 1e+01 0.00
167 1e+01 0.00
168 1e+01 0.00
169 1e+01 0.00
170 1e+01 0.00
171 1e+01 0.00
172 1e+01 0.00
173 1e+01 0.00
174 1e+01 0.00
175 1e+01 0.00
176 1e+01 0.00
177 1e+01 0.00
178 1e+01 0.00
179 1e+01 0.00
180 1e+01 0.00
m1 <- drm(Length ~ conc, data = dat, fct = CRS.4c())
Error in optim(startVec, opfct, hessian = TRUE, method = optMethod,
control = list(maxit = maxIt, :
non-finite finite-difference value [1]
Error in drmOpt(opfct, opdfct1, startVecSc, optMethod, constrained, warnVal, :
Convergence failed
Thanks in advance
Syamkumar R
School of Environmental Studies
CUSAT
Kerala, India
More information about the R-help
mailing list