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

17.98
17.83
19.45
21.04
20.03
20.01
19.64
18.52
19.59
20.09
19.82
21.09
22.64
22.11
20.42
18.58
18.24
16.87
18.64
27.17
33.69
35.92
32.30
27.34
24.96
20.52
19.86
20.82
21.24
20.20
21.42
21.69
21.86
23.23
22.47
19.52
18.82
19.00
18.92
20.24
20.94
22.38
21.76
21.35
21.90
21.69
20.34
19.41
19.08
20.05
20.35
20.27
19.94
19.07
17.87
18.01
17.51
18.15
16.70
14.51
15.00
14.78
14.66
16.38
17.88
19.07
19.65
18.38
17.46
17.71
18.10
17.16
17.99
18.53
18.55
19.87
19.74
18.42
17.30
18.03
18.23
17.44
17.99
19.04
18.88
19.07
21.36
23.57
21.25
20.45
21.32
21.96
23.99
24.90
23.71
25.39
25.17
22.21
20.99
19.72
20.83
19.17
19.63
19.93
19.79
21.26
20.17
18.32
16.71
16.06
15.02
15.44
14.86
13.66
14.08
13.36
14.95
14.39
12.85
11.28
12.47
12.01
14.66
17.34
17.75
17.89
20.07
21.26
23.88
22.64
24.97
26.08
27.18
29.35
29.89
25.74
28.78
31.83
29.77
31.22
33.88
33.08
34.40
28.46
29.58
29.61
27.24
27.41
28.64
27.60
26.45
27.47
25.88
22.21
19.67
19.33
19.67
20.74
24.42
26.27
27.02
25.52
26.94
28.38
29.67
28.85
26.27
29.42
32.94
35.87
33.55
28.25
28.14
30.72
30.76
31.59
28.29
30.33
31.09
32.15
34.27
34.74
36.76
36.69
40.28
38.02
40.69
44.94
45.95
53.13
48.46
43.33
46.84
47.97
54.31
53.04
49.83
56.26
58.70
64.97
65.57
62.37
58.30
59.43
65.51
61.63
62.90
69.69
70.94
70.96
74.41
73.05
63.87
58.88
59.37
62.03
54.57
59.26
60.56
63.97
63.46
67.48
74.18
72.39
79.93
86.20
94.62
91.73
92.95
95.35
105.56
112.57
125.39
133.93
133.44
116.61
103.90
76.65
57.44
41.02
41.74
39.16
47.98
49.79
59.16
69.68
64.09
71.06
69.46
75.82
78.08
74.30

par1 <- as.numeric(par1)
par2 <- as.numeric(par2)
x <- as.ts(x)
library(lattice)
bitmap(file='pic1.png')
plot(x,type='l',main='Run Sequence Plot',xlab='time or index',ylab='value')
grid()
dev.off()
bitmap(file='pic2.png')
hist(x)
grid()
dev.off()
bitmap(file='pic3.png')
if (par1 > 0)
{
densityplot(~x,col='black',main=paste('Density Plot   bw = ',par1),bw=par1)
} else {
densityplot(~x,col='black',main='Density Plot')
}
dev.off()
bitmap(file='pic4.png')
qqnorm(x)
qqline(x)
grid()
dev.off()
if (par2 > 0)
{
bitmap(file='lagplot1.png')
dum <- cbind(lag(x,k=1),x)
dum
dum1 <- dum[2:length(x),]
dum1
z <- as.data.frame(dum1)
z
plot(z,main='Lag plot (k=1), lowess, and regression line')
lines(lowess(z))
abline(lm(z))
dev.off()
if (par2 > 1) {
bitmap(file='lagplotpar2.png')
dum <- cbind(lag(x,k=par2),x)
dum
dum1 <- dum[(par2+1):length(x),]
dum1
z <- as.data.frame(dum1)
z
mylagtitle <- 'Lag plot (k='
mylagtitle <- paste(mylagtitle,par2,sep='')
mylagtitle <- paste(mylagtitle,'), and lowess',sep='')
plot(z,main=mylagtitle)
lines(lowess(z))
dev.off()
}
bitmap(file='pic5.png')
acf(x,lag.max=par2,main='Autocorrelation Function')
grid()
dev.off()
}
summary(x)
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Descriptive Statistics',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'# observations',header=TRUE)
a<-table.element(a,length(x))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'minimum',header=TRUE)
a<-table.element(a,min(x))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Q1',header=TRUE)
a<-table.element(a,quantile(x,0.25))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'median',header=TRUE)
a<-table.element(a,median(x))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'mean',header=TRUE)
a<-table.element(a,mean(x))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Q3',header=TRUE)
a<-table.element(a,quantile(x,0.75))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'maximum',header=TRUE)
a<-table.element(a,max(x))
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')