blob: a3f84c0ce06ab695197a8a4dce9f4514c32367d7 (
plain) (
blame)
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
|
R version 3.1.1 (2014-07-10) -- "Sock it to Me"
Copyright (C) 2014 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
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)
> data(XY)
> rXY <- drfit(XY,logit=TRUE,weibull=TRUE,chooseone=FALSE)
Control: Fitting data...
Substance X: Fitting data...
Waiting for profiling to be done...
Waiting for profiling to be done...
Substance Y: Fitting data...
> print(rXY,digits=5)
Substance ndl n lld lhd mtype logED50 2.5% 97.5% unit sigma
1 Control 1 6 -Inf -Inf inactive NA NA NA mg/L NA
2 Substance X 4 12 1 3 probit 2.1619 2.1059 2.2144 mg/L 0.041317
3 Substance X 4 12 1 3 logit 2.1609 2.1121 2.2079 mg/L 0.037587
4 Substance X 4 12 1 3 weibull 2.1744 NA NA mg/L 0.057239
5 Substance Y 4 12 1 3 active NA NA NA mg/L NA
a b
1 NA NA
2 2.1619 0.50613
3 2.1609 0.29914
4 3.0096 2.03624
5 NA NA
>
> proc.time()
user system elapsed
0.609 0.379 0.491
|
