setwd("c:/glmed/data/")

rm(list=ls(all=TRUE))

# forma wide
(df.a <- read.table("Table_12_17a.txt",header=TRUE))

# forma long
df.b <- read.delim("Table_12_17b.txt",header=TRUE)

fit <- aov(Opinion ~ Influence + Error(as.factor(Subject)/Influence), data=df.b)
summary(fit)

#
X <- df.a[,2:4]
X <- as.matrix(X)

mcol <- colMeans(X)
mrig <- rowMeans(X)
mgen <- mean(X)

n <- nrow(X)
y <- as.vector(X)
yr <- rep(mrig,ncol(X))
yc <- rep(mcol, each=n)
yg <- rep(mgen, length(y))

tmp <- (y - yr) - (yc - yg)
sum(tmp)
sum(tmp^2)
summary(fit)