ANOVAPts1 <- scan(what = list(obs1=0.,obs2=0.,obs3=0.,group=""), multi.line=F)
 1 1 1 O  2 2 2 O  3 3 3 O  .5 .5 .5 O  4 4 4 O 11 11 11 O  15 15 15 O
 7 7 7 O  -1 -1 -1 O  -2 -2 -2 O  5 5 5 O  6 6 6 O 	-11 -11 -11 O  -15 -15 -15 O
 12 10 20 Y 
 12 15 15 LS 
 1 0 0 D1  2 0 0 D1  4 0 0 D1  7 0 0 D1  9 0 0 D1  11 0 0 D1  14 0 0 D1
 -1 0 0 D1  -2 0 0 D1  -4 0 0 D1  -7 0 0 D1  -9 0 0 D1  -11 0 0 D1  -14 0 0 D1
 0 1 1 D2  0 3 3 D2  0 5 5 D2  0 7 7 D2 0 10 10 D2  0 14 14 D2 0 17 17 D2 
 0 -1 -1 D2  0 -3 -3 D2  0 -5 -5 D2  0 -7 -7 D2 0 -10 -10 D2  0 -14 -14 D2 0 -17 -17 D2 
-20    -40    -40      C  -20    -36    -36      C
-20    -32    -32      C  -20    -28    -28      C
-20    -24    -24      C  -20    -20    -20      C
-20    -16    -16      C  -20    -12    -12      C
-20     -8     -8      C  -20     -4     -4      C
-20      0      0      C  -16    -36    -36      C
-16    -32    -32      C  -16    -28    -28      C
-16    -24    -24      C  -16    -20    -20      C
-16    -16    -16      C  -16    -12    -12      C
-16     -8     -8      C  -16     -4     -4      C
-16      0      0      C  -16      4      4      C
-12    -32    -32      C  -12    -28    -28      C
-12    -24    -24      C  -12    -20    -20      C
-12    -16    -16      C  -12    -12    -12      C
-12     -8     -8      C  -12     -4     -4      C
-12      0      0      C  -12      4      4      C
-12      8      8      C   -8    -28    -28      C
 -8    -24    -24      C   -8    -20    -20      C
 -8    -16    -16      C   -8    -12    -12      C
 -8     -8     -8      C   -8     -4     -4      C
 -8      0      0      C   -8      4      4      C
 -8      8      8      C   -8     12     12      C
 -4    -24    -24      C   -4    -20    -20      C
 -4    -16    -16      C   -4    -12    -12      C
 -4     -8     -8      C
 -4     -4     -4      C   -4      0      0      C
 -4      4      4      C   -4      8      8      C
 -4     12     12      C   -4     16     16      C
  0    -20    -20      C    0    -16    -16      C
  0    -12    -12      C    0     -8     -8      C
  0     -4     -4      C    0      0      0      C
  0      4      4      C    0      8      8      C
  0     12     12      C    0     16     16      C
  0     20     20      C    4    -16    -16      C
  4    -12    -12      C    4     -8     -8      C
  4     -4     -4      C    4      0      0      C
  4      4      4      C    4      8      8      C
  4     12     12      C    4     16     16      C
  4     20     20      C    4     24     24      C
  8    -12    -12      C    8     -8     -8      C
  8     -4     -4      C    8      0      0      C
  8      4      4      C    8      8      8      C
  8     12     12      C    8     16     16      C
  8     20     20      C    8     24     24      C
  8     28     28      C   12     -8     -8      C
 12     -4     -4      C   12      0      0      C
 12      4      4      C   12      8      8      C
 12     12     12      C   12     16     16      C
 12     20     20      C   12     24     24      C
 12     28     28      C
12     32     32      C
16     -4     -4      C
16      0      0      C
16      4      4      C
16      8      8      C
16     12     12      C
16     16     16      C
16     20     20      C
16     24     24      C
16     28     28      C
16     32     32      C
16     36     36      C
20      0      0      C
20      4      4      C
20      8      8      C
20     12     12      C
20     16     16      C
20     20     20      C
20     24     24      C
20     28     28      C
20     32     32      C
20     36     36      C
20     40     40      C

dim(ANOVAPts1)


ANOVAPts1  <- data.frame(ANOVAPts1)

library(rgl)

ANOVAPts1$colors   <- NA
ANOVAPts1$colors[ANOVAPts1$group == "C"]   <- "blue"
ANOVAPts1$colors[ANOVAPts1$group == "LS"]   <- "limegreen"
ANOVAPts1$colors[ANOVAPts1$group == "O"]   <- "red"
ANOVAPts1$colors[ANOVAPts1$group == "Y"]   <- "violetred"
ANOVAPts1$colors[ANOVAPts1$group == "D1"]   <- "magenta"
ANOVAPts1$colors[ANOVAPts1$group == "D2"]   <- "orangered"
    
plot3d(ANOVAPts1$obs1, ANOVAPts1$obs2, ANOVAPts1$obs3, col = ANOVAPts1$colors)