利用proast70.3.zip进行计算时,
f.proast(erfenlei)
What type of response data do you want to consider?
1: continuous, individual data
2: binary
3: ordinal
4: quantal
5: clustered continuous, individual data
6: clustered quantal
7: continuous, summary data
8: clustered continuous, summary data
9: quantal, CxT
10: other
Selection: 4
Do you want to fit a set of models, or choose a single model?
1: single model
2: change settings first
3: set of models
Selection: 3
Q1: Which variable do you want to consider as independent variable?
(e.g. dose, age)
1 : dose
2 : n
3 : RBC
4 : HGB
5 : HCT
6 : MCV
7 : MCH
8 : MCHC
9 : LGALT
-------- > 1
1 dose
2 n
3 RBC
4 HGB
5 HCT
6 MCV
7 MCH
8 MCHC
9 LGALT
Which response(s) you want to analyse
by set of models > 4
1 dose
2 n
3 RBC
4 HGB
5 HCT
6 MCV
7 MCH
8 MCHC
9 LGALT
Enter column number(s) with the associated sample sizes > 2
Give number of factor serving as potential covariate (e.g.sex)
-- type 0 if none ---
1: dose
2: n
3: RBC
4: HGB
5: HCT
6: MCV
7: MCH
8: MCHC
9: LGALT
Selection: 0
What type of Benchmark response do you want to consider?
type 0 if you do not need CIs
1: ED50
2: Additional risk, i.e. P[BMD] - P[0]
3: Extra risk, i.e. (P[BMD]-P[0])/(1-P[0])
4: CED for latent variable
Selection: 3
Give value for the BMR,
in terms of extra risk > 0.05
Do you want to calculate the BMD confidence interval by model averaging?
1: no
2: yes
Selection: 1
response: HGB
ANALYSIS WITH EXPONENTIAL MODELS
model converged npar loglik aic
full model 1 4 -1232.34 2472.68
得到以下错误
Error in eval(substitute(expr), data, enclos = parent.frame()) :
object 'CES' not found