• QR Code (size=2)
• QR Code (size=3)
• QR Code (size=4)
• QR Code (size=5)
• QR Code (size=6)

62
67,1
75,9
67
74,2
72,2
60,2
65,8
76,2
76,6
76,8
70,6
74,5
73,5
80,2
71,5
76,6
79,6
65,5
69,2
74,8
79,4
75
67,7
72,5
71,2
78,3
76,6
74,9
76,5
69,4
67,4
77,2
82,2
75,1
70,6
75,6
73,5
79,4
77,5
72,9
78
71,5
66,6
81,8
83,5
74,6
79,8
73,9
76,6
88,9
81,7
76,5
88,8
75,5
75,2
89
87,9
85,7
89,2
82,7
81
90,3
86,3
81,5
91,1
73,1
76,4
91
86,9
89,6
90,5
86,3
86,5
98,8
84,3
91,2
95,5
78,1
81,5
94,4
98,5
95,3
91,6
92,8
90,5
102,2
91,5
94,9
102,1
88,8
89,4
97,8
108,8
100,8
95
101
101
102,5
105,6
98,3
105,5
96,4
88
108,1
107,2
92,5
95,7
84,8
85,4
94,6
86
88,6
93,3
83,1
82,6
96,7
96,2
92,6
92,7
89,9
95,4
108,4
96,2
95
109
91,9
92,2
107,1
105,6
105,4
103,9
99,2
102,4
121,8
102,3
110,1
106
91,9
100,1
112
105
103,3
101,8
100,9
104,2
116,8
97,8
100,7
107,2
96,3
95,9
104,6
107,5
102,5
94,9
98,7
96,8
108,3
103,9
102,4
107,3
101,9
92,5
105,4
113,2
105,7
101,7
101,8
102,9
109,2
105,6
103,4
108,8
98,1
90
112,8
112,2
102,2
102,5
101,8
98,8
114,3
105,2
98,3
110,1
96,4
92,1
112,2
111,6
107,6
103,4
103,6
107,7
117,9
110,4
104,4
116,2
98,9
102,1
113,7
109,5
110,3
114,5
107
109,4
124,6
104,8
112
119,2
103
106,5

library(car)
par2 <- abs(as.numeric(par2)*100)
par3 <- as.numeric(par3)*100
if(par4=='') par4 <- 0
par4 <- as.numeric(par4)
numlam <- par2 + par3 + 1
x <- x + par4
n <- length(x)
c <- array(NA,dim=c(numlam))
l <- array(NA,dim=c(numlam))
mx <- -1
mxli <- -999
for (i in 1:numlam)
{
l[i] <- (i-par2-1)/100
if (l[i] != 0)
{
if (par1 == 'Full Box-Cox transform') x1 <- (x^l[i] - 1) / l[i]
if (par1 == 'Simple Box-Cox transform') x1 <- x^l[i]
} else {
x1 <- log(x)
}
c[i] <- cor(qnorm(ppoints(x), mean=0, sd=1),sort(x1))
if (mx < c[i])
{
mx <- c[i]
mxli <- l[i]
x1.best <- x1
}
}
print(c)
print(mx)
print(mxli)
print(x1.best)
if (mxli != 0)
{
if (par1 == 'Full Box-Cox transform') x1 <- (x^mxli - 1) / mxli
if (par1 == 'Simple Box-Cox transform') x1 <- x^mxli
} else {
x1 <- log(x)
}
mypT <- powerTransform(x)
summary(mypT)
bitmap(file='test1.png')
plot(l,c,main='Box-Cox Normality Plot', xlab='Lambda',ylab='correlation')
mtext(paste('Optimal Lambda =',mxli))
grid()
dev.off()
bitmap(file='test2.png')
hist(x,main='Histogram of Original Data',xlab='X',ylab='frequency')
grid()
dev.off()
bitmap(file='test3.png')
hist(x1,main='Histogram of Transformed Data', xlab='X',ylab='frequency')
grid()
dev.off()
bitmap(file='test4.png')
qqPlot(x)
grid()
mtext('Original Data')
dev.off()
bitmap(file='test5.png')
qqPlot(x1)
grid()
mtext('Transformed Data')
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Box-Cox Normality Plot',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'# observations x',header=TRUE)
a<-table.element(a,n)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'maximum correlation',header=TRUE)
a<-table.element(a,mx)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'optimal lambda',header=TRUE)
a<-table.element(a,mxli)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'transformation formula',header=TRUE)
if (par1 == 'Full Box-Cox transform') {
a<-table.element(a,'for all lambda <> 0 : T(Y) = (Y^lambda - 1) / lambda')
} else {
a<-table.element(a,'for all lambda <> 0 : T(Y) = Y^lambda')
}
a<-table.row.end(a)
if(mx<0) {
a<-table.row.start(a)
a<-table.element(a,'Warning: maximum correlation is negative! The Box-Cox transformation must not be used.',2)
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable.tab')
if(par5=='Yes') {
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Obs.',header=T)
a<-table.element(a,'Original',header=T)
a<-table.element(a,'Transformed',header=T)
a<-table.row.end(a)
for (i in 1:n) {
a<-table.row.start(a)
a<-table.element(a,i)
a<-table.element(a,x[i])
a<-table.element(a,x1.best[i])
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable1.tab')
}
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Maximum Likelihood Estimation of Lambda',1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,paste('
',RC.texteval('summary(mypT)'),'
',sep=''))
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable3.tab')