Scaling of Cmin [Regulatives / Guidelines]

posted by VStus – Poland, 2017-02-06 09:55 (2634 d 16:25 ago) – Posting: # 17017
Views: 10,549

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
22,993 posts in 4,828 threads, 1,656 registered users;
118 visitors (0 registered, 118 guests [including 3 identified bots]).
Forum time: 03:21 CEST (Europe/Vienna)

Never never never never use Excel.
Not even for calculation of arithmetic means.    Martin Wolfsegger

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