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

50
52.4
57.5
52.5
57.5
57.6
48.3
52
62.1
59.1
62.6
57.9
59.3
61.5
66
61.1
63.8
69.6
57
59.9
63.8
69.8
64.6
60.8
64.7
63.6
68.8
66.4
64.4
65.3
63
61.1
67.7
72.3
65.4
63.2
69.4
62.3
71
68.6
62
68.2
66.8
65.5
76.9
78.1
67.6
80.1
64.7
70.4
84.6
75.1
69.6
81.8
74.2
72.9
84.9
80.5
79.6
90.8
76.5
70.9
82.3
77.8
75.6
81.3
71
75.1
89.2
84.1
82.7
82.4
78.2
78.5
91.5
76.6
80.6
85.9
74.5
79.4
89.7
92.7
89.6
87
80.9
76.2
89.7
79.1
82.4
90.3
85.8
83.5
85.1
90.6
87.7
86
89.7
86.2
91.1
91.3
85.5
92
91.5
80
100.9
97.3
89.1
104
80.2
83.3
97.5
86.8
84.3
93.4
90.2
82.5
93.7
93.9
91.1
96.9
88.2
100.9
109.5
91
89.5
109.6
97.9
94.9
103.5
100
107.1
108
95
102.2
131.4
104.5
105.6
106.1
98
113
113.2
105.4
100.1
100.7
96.1
98.2
123.5
93.9
94.8
103.5
105.3
105.8
112
114.5
108.3
103.8
103
97.7
118.7
115.1
110
117.3
119.1
105.9
114.1
124.6
117.3
115
103.6
113.4
122
122.5
119.6
132.6
113
107.5
139.3
134.6
125.6
124
111.9
101.5
130.2
121.9
111.3
122
116.4
119.1
133
128.9
126.1
122.3
110.2
113.6
131
123.2
120.7
142.8
131.7
131.6
139
128.5
122.7
148.4
118.6
126.3
141
120.9
127
138.5
131.9
136.3

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')