#=================parta
##each pair of observations occurs in the same block one time, lambda=1
(maybe 2)


#=================partb
##Efficiency factor
##E = [(lambda)(t)]/[rk]
lambda=1 ##how many times each pair is repeated  
t = 7	##number of treatments, total of columns
r = 4 ##number of item w/in columns
k = 4	##number of items w/in each row
E = ((lambda)*(t))/(r*k)
E = 2*7/4*4 =14/16
#=================partC
#y_ij = \mu _ \beta_j + \tau_i + \epsilon_ij  i=1, 2, 3 , 4, 5, 6, 7  (number of colms) j = 1-7 (number of rows)
#
#\beta_j = jth block as block effect 
#\tau_i = ith treatment 
#\mu is general mean
#\epsilon = experimental error, mean=0 variance \sigma^2
#

#==========part c
BIBProb2=read.table(file="E:/FSU-FALL07-School/Fall07-school/STA5166/hw5/BIBprob2.dat", header=TRUE)
BIBProb2

fit1= aov(Data~factor(Run)+factor(Voltage), data=BIBProb2)
fit1





#==========part d
plot(fit1$fitted, fit1$resid)


#==========part e
lndata = log(BIBProb2$Data)
lndata 

BIBProb2_2= data.frame(lndata, data=BIBProb2)
BIBProb2_2
fit2= aov(lndata~factor(data.Voltage)+ factor(data.Run), data=BIBProb2_2)
fit2
summary(fit2)
par(mfrow=c(1,1))
plot(fit2$fitted, fit2$resid)





fit_BIBprob2  = aov(Data~factor(Run)+factor(Voltage), data=BIBProb2 )
fit_BIBprob2

summary(fit_BIBprob2)



#=================partD
par(mfrow=c(2,4))
##may need my fuction see other script
source("C:\\Documents and Settings\\Jaime\\Desktop\\FALL07\\STA5166\\hw1eda1-2.txt")

eda1(fit_BIBprob2$resid)



fit_BIBprob2_2  = aov(lndata~factor(Run)+factor(Voltage), data=BIBProb2 )
fit_BIBprob2_2

summary(fit_BIBprob2_2)

#par(mfrow=c(4,4))
##may need my fuction see other script
source("C:\\Documents and Settings\\Jaime\\Desktop\\FALL07\\STA5166\\hw1eda1-2.txt")

eda1(fit_BIBprob2_2$resid)