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

35.18
35.24
35.00
35.50
36.32
36.04
36.71
36.70
36.17
36.50
36.94
36.42
36.40
36.68
36.54
36.36
36.23
36.39
36.44
36.43
36.47
36.42
36.44
36.44
36.51
36.51
36.26
36.41
36.19
36.45
36.38
36.70
36.69
37.00
37.20
37.28
37.18
37.13
36.87
36.88
36.85
36.87
37.40
37.21
37.40
36.94
36.66
36.59
37.15
37.00
36.47
36.51
36.17
36.62
35.97
36.48
36.34
37.05
37.11
36.92
36.89
36.94
37.19
36.78
36.25
36.67
36.84
36.54
37.09
37.02
37.04
37.47
37.36
37.38
37.18
37.19
37.35
37.33
37.98
37.72
37.75
37.94
37.82
38.07
38.00
38.00
38.00
37.94
38.24
38.52
38.82
38.73
38.58
37.67
37.79
37.65
38.23
38.10
38.46
38.18
38.38
38.78
38.99
38.90
38.88
38.82
38.84
39.16
39.34
39.92
39.54
38.64
38.47
38.01
37.81
38.20
38.31
38.41
38.46
38.59
38.71
38.90
38.65
38.95
38.84
38.97
39.20
39.01
38.78
38.80
39.60
39.42
39.32
39.28
39.49
39.12
39.08
39.45
39.73
39.44
39.34
39.40
39.52
39.60
39.42
39.30
39.43
39.75
39.23
39.39
39.63
39.56
39.31
39.48
39.65
39.27
38.71
38.66
38.78
38.83
38.82
38.62
38.39
38.57
38.60
38.11
38.16
38.40
39.18
39.03
39.12
39.22
39.39
39.77
39.65
39.81
39.79
40.32
40.33
40.48
41.12
41.24
40.82
41.06
40.29
40.18
39.91
39.81
40.25
39.93
39.90
39.95

56,81
57,20
56,69
56,91
57,33
56,78
56,87
57,93
58,50
58,57
59,45
59,02
59,61
59,62
59,62
59,52
59,60
59,48
60,48
61,71
61,63
61,15
60,86
60,58
60,45
60,52
59,91
59,98
59,95
59,50
60,72
61,20
61,74
60,97
62,29
61,95
62,04
62,48
61,97
61,79
61,70
61,23
62,23
62,95
63,38
63,87
62,04
61,70
62,30
62,29
62,67
62,75
61,83
61,90
60,37
61,95
62,66
63,62
64,89
64,73
63,93
66,13
65,97
65,86
63,24
64,14
63,34
65,47
65,31
65,28
64,41
66,43
66,50
66,67
66,00
66,29
66,04
66,31
66,65
64,96
65,69
65,71
65,36
66,40
65,27
65,45
65,20
62,75
61,51
61,02
61,75
59,66
59,73
58,92
58,59
58,99
59,30
59,44
59,72
60,33
60,75
61,64
62,00
62,03
60,40
60,14
59,82
60,50
59,60
58,87
58,36
57,93
56,96
56,67
56,67
57,79
58,30
58,02
57,79
57,06
57,00
57,99
57,45
58,28
51,75
52,38
53,12
53,02
53,30
53,25
53,31
53,10
52,84
53,52
53,26
53,03
52,92
53,70
54,08
54,15
54,18
54,22
53,76
53,54
53,27
53,91
54,00
54,00
54,02
54,45
54,53
54,80
54,67
54,87
54,44
54,44
54,10
53,46
53,20
52,03
51,52
51,60
51,44
50,43
50,35
50,18
49,85
52,50
50,03
49,44
50,16
50,11
49,21
48,51
49,18
49,90
49,87
49,82
50,20
50,19
50,67
51,35
50,74
50,37
50,40
51,44
50,89
54,62
53,23
53,10
53,91
53,91

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