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

94
103
93
103
51
70
91
22
38
93
60
123
148
90
124
70
168
115
71
66
134
117
108
84
156
120
114
94
120
81
110
133
122
158
109
124
39
92
126
0
70
37
38
120
93
95
77
90
80
31
110
66
138
133
113
100
7
140
61
41
96
164
78
49
102
124
99
129
62
73
114
99
70
104
116
91
74
138
67
151
72
120
115
105
104
108
98
69
111
99
71
27
69
107
73
107
93
129
69
118
73
119
104
107
99
90
197
36
85
139
106
50
64
31
63
92
106
63
69
41
56
25
65
93
114
38
44
87
110
0
27
83
30
80
98
82
0
60
28
9
33
59
49
115
140
49
120
66
21
124
152
139
38
144
120
160
114
39
78
119
141
101
56
133
83
116
90
36
50
61
97
98
78
117
148
41
105
55
132
44
21
50
0
73
86
0
13
4
57
48
46
48
32
68
87
43
67
46
46
56
48
44
60
65
55
38
52
60
54
86
24
52
49
61
61
81
43
40
40
56
68
79
47
57
41
29
3
60
30
79
47
40
48
36
42
49
57
12
40
43
33
77
43
45
47
43
45
50
35
7
71
67
0
62
54
4
25
40
38
19
17
67
14
30
54
35
59
24
58
42
46
61
3
52
25
40
32
4
49
63
67
32
23
7
54
37
35
51
39

30
28
38
30
22
26
25
18
11
26
25
38
44
30
40
34
47
30
31
23
36
36
30
25
39
34
31
31
33
25
33
35
42
43
30
33
13
32
36
0
28
14
17
32
30
35
20
28
28
39
34
26
39
39
33
28
4
39
18
14
29
44
21
16
28
35
28
38
23
36
32
29
25
27
36
28
23
40
23
40
28
34
33
28
34
30
33
22
38
26
35
8
24
29
20
29
45
37
33
33
25
32
29
28
28
31
52
21
24
41
33
32
19
20
31
31
32
18
23
17
20
12
17
30
31
10
13
22
42
1
9
32
11
25
36
31
0
24
13
8
13
19
18
33
40
22
38
24
8
35
43
43
14
41
38
45
31
13
28
31
40
30
16
37
30
35
32
27
20
18
31
31
21
39
41
13
32
18
39
14
7
17
0
30
37
0
5
1
16
32
24
17
11
24
22
12
19
13
17
15
16
24
15
17
18
20
16
16
18
22
8
17
18
16
23
22
13
13
16
16
20
22
17
18
17
12
7
17
14
23
17
14
15
17
21
18
18
17
17
16
15
21
16
14
15
17
15
15
10
6
22
21
1
18
17
4
10
16
16
9
16
17
7
15
14
14
18
12
16
21
19
16
1
16
10
19
12
2
14
17
19
14
11
4
16
20
12
15
16

par1 <- as.numeric(par1)
library(lattice)
z <- as.data.frame(cbind(x,y))
m <- lm(y~x)
summary(m)
bitmap(file='test1.png')
plot(z,main='Scatterplot, lowess, and regression line')
lines(lowess(z),col='red')
abline(m)
grid()
dev.off()
bitmap(file='test2.png')
m2 <- lm(m$fitted.values ~ x)
summary(m2)
z2 <- as.data.frame(cbind(x,m$fitted.values))
names(z2) <- list('x','Fitted')
plot(z2,main='Scatterplot, lowess, and regression line')
lines(lowess(z2),col='red')
abline(m2)
grid()
dev.off()
bitmap(file='test3.png')
m3 <- lm(m$residuals ~ x)
summary(m3)
z3 <- as.data.frame(cbind(x,m$residuals))
names(z3) <- list('x','Residuals')
plot(z3,main='Scatterplot, lowess, and regression line')
lines(lowess(z3),col='red')
abline(m3)
grid()
dev.off()
bitmap(file='test4.png')
m4 <- lm(m$fitted.values ~ m$residuals)
summary(m4)
z4 <- as.data.frame(cbind(m$residuals,m$fitted.values))
names(z4) <- list('Residuals','Fitted')
plot(z4,main='Scatterplot, lowess, and regression line')
lines(lowess(z4),col='red')
abline(m4)
grid()
dev.off()
bitmap(file='test5.png')
myr <- as.ts(m$residuals)
z5 <- as.data.frame(cbind(lag(myr,1),myr))
names(z5) <- list('Lagged Residuals','Residuals')
plot(z5,main='Lag plot')
m5 <- lm(z5)
summary(m5)
abline(m5)
grid()
dev.off()
bitmap(file='test6.png')
hist(m$residuals,main='Residual Histogram',xlab='Residuals')
dev.off()
bitmap(file='test7.png')
if (par1 > 0)
{
densityplot(~m$residuals,col='black',main=paste('Density Plot   bw = ',par1),bw=par1)
} else {
densityplot(~m$residuals,col='black',main='Density Plot')
}
dev.off()
bitmap(file='test8.png')
acf(m$residuals,main='Residual Autocorrelation Function')
dev.off()
bitmap(file='test9.png')
qqnorm(x)
qqline(x)
grid()
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Simple Linear Regression',5,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Statistics',1,TRUE)
a<-table.element(a,'Estimate',1,TRUE)
a<-table.element(a,'S.D.',1,TRUE)
a<-table.element(a,'T-STAT (H0: coeff=0)',1,TRUE)
a<-table.element(a,'P-value (two-sided)',1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'constant term',header=TRUE)
a<-table.element(a,m$coefficients[[1]])
sd <- sqrt(vcov(m)[1,1])
a<-table.element(a,sd)
tstat <- m$coefficients[[1]]/sd
a<-table.element(a,tstat)
pval <- 2*(1-pt(abs(tstat),length(x)-2))
a<-table.element(a,pval)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'slope',header=TRUE)
a<-table.element(a,m$coefficients[[2]])
sd <- sqrt(vcov(m)[2,2])
a<-table.element(a,sd)
tstat <- m$coefficients[[2]]/sd
a<-table.element(a,tstat)
pval <- 2*(1-pt(abs(tstat),length(x)-2))
a<-table.element(a,pval)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')