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

83
79
92
83
92
103
82
86
106
79
86
76
108
82
108
118
127
123
72
105
63
86
58
59
100
100
78
94
105
89
101
92
105
76
80
66
117
94
107
110
110
106
94
71
101
84
89
119
97
82
89
70
101
81
74
107
97
83
95
82
88
74
104
73
73
81
79
83
111
138
81
107
66
81
74
96
86
69
73
71
64
79
60
111
107
90
98
77
93
68
74
70
80
81
72
81
92
81
78
92
92
107
98
86
77
96
104
77
65
61
117
84
69
85
116
115
55
64
117
68
104
66
70
89
79
70
63
79
62
94
83
118
62
78
83
91
84
76
100
80
98
89
98
88
81
88
75
77
88
65
69
76
53
82
67
84
112
91
104
79
90
60
79
99
68
79
107
114
81
83
61
82
134
102
132
72
72
102
92
94
86
84
66
129
88
109
84
73
86
113
88
90
82
111
73
91
72
111
108
83
81
111
69
106
115
132
78
90
132
91
115
65
77
71
74
76
115
100
70
71
74
60
58
105
105
100
74
77
77
109
68
67
96
86
85
91
104
94
67
79
73
93
87
66
79
94
84
67
121
82
116
83
66
66
71
83
93
83
112
79
135
85
91
103
77
70
53
85
88
65
119
93
84
70
64
63
152
83
66
83
106
85
85
84
78
94
82
66
69
83
83
124
101
113
107
83
79
85
62
83
101
60
86
101
73
70
88
74
105
82
83
90
70
56
70
79
127
96
101
81
93
92
79
78
68
92
73
61
73
108
88
66
59
61
55
119
89
68
125
66
82
101
104
63
63
63
98
90
97
74
63
102
90
79
74
89
70
77
78
70
95
100
64
90
109
89
66
88
108
97
66
85
79
95
62
76
69
105
74
96
83
76
83
68
68
76
73
71
74
78
95
103
82
67
77
93
76
85
68
81
87
67
78
93
87
77
79
84
114
84
105
117
94
65
75
73
117
75
85
83
63
80
73
86
103
67
83
55
90
123
98
188
76
72
65
99
86
81
99
58
90
97
85
79
110
92
88
56
71
110
77
66
114
76
109
125
74
89
100
77
99
88
96
87
101
78
119
105
80
69
89
69
66
70
83
81
73
100
67
108
105
119
91
102
100
76
100
86
103
93
57
115
82
99
63
166
74
86
87
101
73
102
97
86
116
81
140
93
111
115
70
68
77
81
64
111
126
97
72
88
78
94
78
70
96
99
79
102
89
79
94
89
94
63
76
55
83
70
81
130
64
77
85
104
72
70
84
106
118
70
74
75
95
71
81
103
56
91
54
76
101
83
62
67
91
91
97
63
80
71
83
88
74
89
65
103
122
123
98
97
69
90
85
60
97
111
105
72
123
129
83
86
88
50
62
79
99
82
71
124
81
83
121
65
61
65
71
71
62
74
63
75
82
94
107
94
91
102
73
107
74
80
112
71
78
69
88
106
56
118
90
76
83
77
122
85
67
93
59
92
58
115
96
120
87
118
87
101
55
98
83
102
79

14
13
16
13
14
15
14
13
15
13
13
13
16
16
12
15
16
16
12
16
13
14
11
12
15
15
14
15
15
15
16
15
16
12
16
12
18
18
15
17
17
16
14
14
15
14
14
17
15
14
15
14
15
15
14
15
14
16
15
15
14
16
16
13
13
15
14
13
17
18
13
16
13
15
15
13
14
14
14
15
12
16
13
15
16
14
17
12
16
13
13
14
14
13
15
15
15
12
12
14
12
15
14
15
13
15
14
13
14
14
17
14
12
15
16
16
13
14
15
13
17
13
13
12
14
15
13
15
13
15
14
16
13
15
15
14
13
15
16
15
16
16
16
14
12
13
15
15
15
14
14
14
13
15
13
15
16
15
13
13
14
15
15
16
15
12
14
16
14
14
14
12
18
14
16
12
14
13
14
16
13
13
13
14
13
14
14
14
15
15
14
13
14
17
13
15
14
16
14
14
14
17
13
14
15
18
15
15
16
15
15
12
14
14
15
14
17
15
14
13
14
14
13
14
14
14
14
15
14
15
13
12
14
14
14
15
15
15
13
13
13
14
13
12
14
15
15
10
14
14
18
15
12
13
17
14
16
14
17
15
19
13
17
16
15
15
14
14
14
12
16
17
14
13
14
13
18
14
12
15
14
16
13
15
16
13
15
13
13
12
14
15
14
15
15
15
14
15
13
13
17
13
11
16
13
15
13
14
15
15
13
15
13
13
14
14
16
15
15
15
16
13
14
15
13
17
13
12
13
15
12
12
14
13
13
15
14
13
17
13
16
16
13
14
13
12
15
15
16
13
13
15
14
15
15
16
14
13
13
15
15
13
11
13
16
13
13
15
15
14
12
14
15
13
13
14
15
16
13
16
13
16
16
14
13
14
14
14
15
14
16
15
14
13
12
16
13
13
13
12
14
12
13
15
15
13
14
13
14
16
16
15
16
15
13
13
16
14
14
14
11
14
14
13
15
12
14
15
15
15
16
19
16
14
13
16
12
16
15
11
16
16
15
15
17
15
14
13
14
16
14
14
17
14
15
15
14
14
16
13
15
14
14
16
15
14
17
16
15
13
15
14
10
15
14
16
14
15
15
15
14
17
15
15
15
15
16
13
15
16
13
17
12
15
15
17
13
17
13
15
15
14
16
15
16
14
17
14
16
16
15
13
15
14
14
15
17
16
15
15
14
14
14
13
14
16
14
16
13
13
16
15
16
14
12
15
14
15
14
17
13
14
15
15
15
15
15
16
17
16
14
13
14
16
16
14
14
15
13
15
17
14
10
14
15
15
15
12
16
14
16
15
13
13
12
13
17
18
15
15
12
14
13
14
15
16
14
15
17
18
15
14
15
12
12
13
15
13
11
15
14
14
16
13
15
13
13
14
12
15
13
15
14
16
15
14
14
16
15
17
14
14
16
15
13
15
11
16
14
15
16
13
13
14
15
14
14
15
13
13
13
17
14
17
14
16
15
16
12
16
15
15
15

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