On Feb 8, 2010, at 10:28 PM, bluesky...@gmail.com wrote:

I have the R code at the end. The last command gives me "1 observation
deleted due to missingness". I don't understand what this error
message. Could somebody help me understand it and how to fix the
problem?

summary(afit)
            Df  Sum Sq Mean Sq F value  Pr(>F)
A             2   0.328 0.16382  0.1899 0.82727
B             3   2.882 0.96057  1.1136 0.34644
C             4   1.676 0.41897  0.4857 0.74617
A:B           6   5.133 0.85554  0.9919 0.43404
A:C           8  14.193 1.77406  2.0568 0.04543 *
B:C          12   7.244 0.60366  0.6999 0.74896
A:B:C        24  19.849 0.82705  0.9589 0.52442
Residuals   119 102.642 0.86254
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
1 observation deleted due to missingness


#######################
a=3
b=4
c=5

A=1:a
B=1:b
C=1:c

n=3

X=matrix(nr=a*b*c*n,nc=n)
colnames(X)=LETTERS[1:n]

for(i_a in 1:a-1) {
 for(i_b in 1:b-1) {
   for(i_c in 1:c-1) {
     for(i_n in 1:n-1) {
X[((i_a * b + i_b) * c + i_c) * n + i_n + 1, ] = c(i_a+1, i_b +1, i_c+1)
     }
   }
 }
}

which(is.na(aframe$Y)  )
[1] 180



set.seed(0)
Y=matrix(nr=a*b*c*n,nc=1)
for(i in 1:(a*b*c)) {
 for(i_n in 1:n-1) {
   fa=X[i,'A']
   fb=X[i,'B']
   fc=X[i,'C']

#Y[(i-1)*n+i_n,1]= fa +fb +fc +fa*fb +fa*fc +fb*fc +fa*fb*fc + rnorm(1)
   Y[(i-1)*n+i_n,1]= rnorm(1)
 }
}

aframe = data.frame(
 A=as.factor(X[,'A'])
 , B=as.factor(X[,'B'])
 , C=as.factor(X[,'C'])
 , Y)

afit=aov(Y ~ A * B * C, aframe)

summary(afit)

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

David Winsemius, MD
Heritage Laboratories
West Hartford, CT

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to