Scaling of Cmin [Regulatives / Guidelines]

posted by VStus – Poland, 2017-02-06 09:55 (3055 d 11:55 ago) – Posting: # 17017
Views: 12,402

Dear zizou,

No change in residual variance if I remove 'period' from the model.

anovadata <- lm(log(AUC) ~ prd + seq + rep + subj:seq, data=data, na.action=na.exclude)
anova(anovadata)

## Analysis of Variance Table
##
## Response: log(AUC)
##           Df  Sum Sq  Mean Sq F value  Pr(>F) 
## prd        1 0.20989 0.209885  6.5803 0.01894 *
## rep        1 0.00507 0.005067  0.1589 0.69466 
## seq:subj  18 1.55131 0.086184  2.7020 0.01873 *
## Residuals 19 0.60602 0.031896                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

MSE       <- anova(anovadata)["Residuals","Mean Sq"] 100*sqrt(exp(MSE)-1)

## [1] 18.00281

anovadata <- lm(log(AUC) ~ seq + rep + subj:seq, data=data, na.action=na.exclude)
anova(anovadata)
## Analysis of Variance Table
##
## Response: log(AUC)
##           Df  Sum Sq  Mean Sq F value  Pr(>F) 
## seq        1 0.20989 0.209885  6.5803 0.01894 *
## rep        1 0.00507 0.005067  0.1589 0.69466 
## seq:subj  18 1.55131 0.086184  2.7020 0.01873 *
## Residuals 19 0.60602 0.031896                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1



Maybe MIXED procedure should be used?

Regards, VStus

Complete thread:

UA Flag
Activity
 Admin contact
23,424 posts in 4,927 threads, 1,675 registered users;
33 visitors (0 registered, 33 guests [including 11 identified bots]).
Forum time: 22:51 CEST (Europe/Vienna)

Philosophy, like medicine, has plenty of drugs, few good remedies,
and hardly any specific cures.    Sebastien-Roch Nicolas de Chamfort

The Bioequivalence and Bioavailability Forum is hosted by
BEBAC Ing. Helmut Schütz
HTML5