R version 2.6.0 (2007-10-03)
Copyright (C) 2007 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(IM1xIPC81)
> rIM1xIPC81 <- drfit(IM1xIPC81,linlogit=TRUE)
IM13 BF4: Fitting data...
IM14 BF4: Fitting data...
Error in nls(response ~ linlogitf(dose, 1, f, logED50, b), data = tmp, :
Convergence failure: false convergence (8)
Error in numericDeriv(form[[3]], names(ind), env, ifelse(internalPars < :
Missing value or an infinity produced when evaluating the model
In addition: Warning message:
In pnorm(q, mean, sd, lower.tail, log.p) : NaNs produced
IM15 BF4: Fitting data...
IM16 BF4: Fitting data...
IM17 BF4: Fitting data...
Waiting for profiling to be done...
Error in numericDeriv(form[[3]], names(ind), env, ifelse(internalPars < :
Missing value or an infinity produced when evaluating the model
In addition: Warning message:
In pnorm(q, mean, sd, lower.tail, log.p) : NaNs produced
IM18 BF4: Fitting data...
Waiting for profiling to be done...
Error in numericDeriv(form[[3]], names(ind), env, ifelse(internalPars < :
Missing value or an infinity produced when evaluating the model
In addition: Warning message:
In pnorm(q, mean, sd, lower.tail, log.p) : NaNs produced
IM19 BF4: Fitting data...
Waiting for profiling to be done...
Error in numericDeriv(form[[3]], names(ind), env, ifelse(internalPars < :
Missing value or an infinity produced when evaluating the model
In addition: Warning message:
In pnorm(q, mean, sd, lower.tail, log.p) : NaNs produced
IM1-10 BF4: Fitting data...
Waiting for profiling to be done...
> print(rIM1xIPC81,digits=5)
Substance ndl n lld lhd mtype logED50 2.5% 97.5% unit
1 IM13 BF4 9 81 0.591760 3.0000 inactive NA NA NA microM
2 IM14 BF4 20 216 -0.010300 3.1761 no fit NA NA NA microM
3 IM15 BF4 9 135 0.591760 3.0000 inactive NA NA NA microM
4 IM16 BF4 9 108 0.591760 3.0000 inactive NA NA NA microM
5 IM17 BF4 9 81 0.591760 3.0000 linlogit 2.5786 2.50572 2.66174 microM
6 IM18 BF4 9 135 0.591760 3.0000 linlogit 1.6806 1.62298 1.74189 microM
7 IM19 BF4 9 81 0.591760 3.0000 linlogit 1.6496 1.59848 1.70315 microM
8 IM1-10 BF4 11 162 -0.010300 3.0000 linlogit 0.7697 0.68701 0.85442 microM
sigma a b c
1 NA NA NA NA
2 NA NA NA NA
3 NA NA NA NA
4 NA NA NA NA
5 0.23759 2.5786 2.3002 0.014680
6 0.23249 1.6806 2.2373 0.057193
7 0.14527 1.6496 1.9770 0.109561
8 0.29877 0.7697 1.9363 0.458090
>
