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 >