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

0,026485833	0,621741667	0,220257333
0,01789525	0,619904167	0,216473333
0,01058325	0,618554167	0,212748
0,00272225	0,6168875	0,209097667
-0,006853	0,6149125	0,205557333
-0,016679833	0,6138625	0,202085
-0,025014583	0,613270833	0,198680333
-0,0328035	0,612245833	0,195361667
-0,037829083	0,610483333	0,19211
-0,04581475	0,608066667	0,188913333
-0,05082375	0,606129167	0,185792
-0,055329	0,605254167	0,182724333
-0,061978417	0,60455	0,17971
-0,06869825	0,603958333	0,176750333
-0,075414583	0,603779167	0,173851
-0,08296875	0,6040875	0,171014333
-0,08974525	0,604204167	0,168218333
-0,095605417	0,602595833	0,165493333
-0,100392167	0,599183333	0,162822667
-0,104531333	0,595816667	0,160189
-0,108372	0,5908625	0,157577667
-0,111641583	0,586108333	0,155009667
-0,1142075	0,5825125	0,152527
-0,118222833	0,5810125	0,150069
-0,123135833	0,580516667	0,14766
-0,128149	0,580458333	0,145308
-0,13175475	0,581395833	0,143004333
-0,1347415	0,582754167	0,140716333
-0,138674083	0,5836125	0,138449667
-0,141716667	0,585216667	0,136218
-0,143839917	0,585195833	0,134024667
-0,145080917	0,584070833	0,131875
-0,147142167	0,583520833	0,129801333
-0,148233	0,583779167	0,127781
-0,14834325	0,584908333	0,125799333
-0,150195833	0,586758333	0,123864333
-0,153742833	0,588091667	0,121979333
-0,158674	0,5879125	0,120146667
-0,161297833	0,588220833	0,118356
-0,16382275	0,58785	0,116604
-0,166688	0,587608333	0,114912667
-0,169795667	0,586445833	0,113262333
-0,1726025	0,5862875	0,111657667
-0,173340583	0,585870833	0,11011
-0,172330083	0,585591667	0,108586
-0,17149225	0,585675	0,107089667
-0,171626667	0,586170833	0,105603667
-0,171026917	0,5871125	0,104188333
-0,169731583	0,588041667	0,102807667
-0,171073083	0,589145833	0,101441667
-0,17323175	0,590345833	0,100107667
-0,174959583	0,591954167	0,098819333
-0,177394083	0,593770833	0,097553333
-0,18091875	0,5952	0,096311
-0,183774167	0,597179167	0,095106333
-0,1870595	0,598733333	0,093909
-0,191677083	0,5993375	0,092729333
-0,19574675	0,599279167	0,091564333
-0,19974575	0,598745833	0,090423
-0,202696	0,597454167	0,089306
-0,20577025	0,595870833	0,08822
-0,210256	0,594158333	0,08716
-0,215462333	0,59245	0,08613
-0,21980175	0,592083333	0,08513
-0,224308	0,592845833	0,084156667
-0,228106917	0,593416667	0,08321
-0,232391833	0,5936875	0,082283333
-0,236569917	0,593633333	0,081376667
-0,24034525	0,593608333	0,080847059
-0,242414	0,593733333	0,079963025
-0,244237917	0,591941667	0,079437288
-0,2456505	0,5896	0,078909402
-0,246628667	0,586170833	0,078375862
-0,249514417	0,5830875	0,077836522
-0,255432667	0,580791667	0,07693913
-0,262324417	0,578520833	0,076382456
-0,268100333	0,576891667	0,075812389
-0,273406417	0,574975	0,075228571
-0,27716975	0,572558333	0,074630631
-0,28102475	0,569704167	0,073686486
-0,284668167	0,566825	0,073072727
-0,288084167	0,563179167	0,072451376
-0,292629667	0,559779167	0,071825926
-0,296878	0,5565875	0,071196262
-0,30100175	0,554629167	0,070246729
-0,3058905	0,553766667	0,069626415
-0,3117765	0,553175	0,069009524
-0,316749583	0,552458333	0,068396154
-0,3201655	0,5515125	0,067790291
-0,321889917	0,5500125	0,066881553
-0,32190225	0,547529167	0,066290196
-0,322839583	0,545666667	0,065706931
-0,324684167	0,543225	0,065132
-0,326995083	0,540983333	0,064561616
-0,3276635	0,540329167	0,063688889
-0,327665583	0,5415125	0,063134694
-0,326595917	0,5428375	0,062585567
-0,32588075	0,5440875	0,0620375
-0,326658	0,5448	0,061494737
-0,326295667	0,544054167	0,060648421
-0,325280083	0,541679167	0,060114894
-0,32382675	0,537858333	0,059582796
-0,322155083	0,532616667	0,059052174
-0,321644333	0,528575	0,058527473
-0,324230167	0,525825	0,057687912
-0,3247365	0,525066667	0,057168889
-0,32253075	0,526145833	0,056651685
-0,320819833	0,528341667	0,056140909
-0,318928	0,531320833	0,055632184
-0,3168905	0,534745833	0,0548
-0,31547925	0,5384625	0,054306977
-0,314361167	0,541970833	0,053821176
-0,313898	0,545429167	0,053342857
-0,31417875	0,547820833	0,052872289
-0,31359175	0,548133333	0,052062651
-0,312060667	0,548591667	0,051614634
-0,312955333	0,548791667	0,051175309
-0,315607	0,550825	0,05074
-0,31803625	0,554104167	0,050313924
-0,319968	0,558625	0,049518987
-0,319647667	0,5640125	0,049102564
-0,31946425	0,569866667	0,048690909
-0,321666167	0,574883333	0,048278947
-0,323271917	0,5767875	0,047861333
-0,322705167	0,576183333	0,047034667
-0,321068417	0,574383333	0,046605405
-0,31994075	0,572333333	0,046169863
-0,319204833	0,571075	0,045727778
-0,320235333	0,5718125	0,045284507
-0,319472	0,5744875	0,044405634
-0,321291667	0,578604167	0,043954286
-0,322324083	0,581820833	0,043501449
-0,320845417	0,5840375	0,043047059
-0,319934167	0,584916667	0,042597015
-0,3198575	0,584445833	0,041701493
-0,319467083	0,583375	0,041266667
-0,319822083	0,5824125	0,040849231
-0,319455167	0,582579167	0,04044375
-0,31908525	0,583041667	0,040057143
-0,319016583	0,583891667	0,039225397
-0,319974583	0,585741667	0,038877419
-0,323622083	0,58775	0,038544262
-0,327519583	0,5894	0,038226667
-0,32819775	0,589491667	0,037925424
-0,329421083	0,589016667	0,037145763
-0,3313865	0,587558333	0,036862069
-0,335011083	0,585125	0,036589474
-0,33841175	0,582470833	0,036321429
-0,340155583	0,579804167	0,036058182
-0,342473833	0,577058333	0,035294545
-0,344401917	0,576166667	0,03502963
-0,34682125	0,5752	0,034769811
-0,347276583	0,574979167	0,034507692
-0,347961167	0,574420833	0,034243137
-0,35049325	0,572591667	0,033458824
-0,352054417	0,569829167	0,033184
-0,353268583	0,567058333	0,032914286
-0,353229833	0,564508333	0,032641667
-0,35584525	0,56205	0,032365957
-0,359304833	0,560220833	0,031574468
-0,363069333	0,559429167	0,031295652
-0,365657583	0,559766667	0,031022222
-0,368090333	0,560904167	0,030754545
-0,371147833	0,562445833	0,030483721
-0,376708	0,5634625	0,029693023
-0,38153025	0,563029167	0,029419048
-0,385772	0,56145	0,02915122
-0,389845	0,558920833	0,02888
-0,39304525	0,556433333	0,028605128
-0,39575	0,554370833	0,027815385
-0,398166667	0,5523375	0,027536842
-0,398333333	0,5504875	0,027264865
-0,396416667	0,5495125	0,027
-0,394	0,549666667	0,026731429
-0,39175	0,549045833	0,025942857
-0,390583333	0,548420833	0,025682353
-0,390833333	0,547933333	0,025418182
-0,390583333	0,547316667	0,02515
-0,389166667	0,546495833	0,024877419
-0,388083333	0,544704167	0,024051613
-0,38825	0,541883333	0,023733333
-0,38875	0,539133333	0,023393103
-0,388416667	0,5354625	0,023014286
-0,387333333	0,531325	0,022577778
-0,386416667	0,528041667	0,02157037
-0,3865	0,526	0,021015385
-0,38725	0,524391667	0,0204
-0,388	0,522725	0,0197
-0,389166667	0,520591667	0,0197
-0,390916667	0,518066667	0,01845
-0,393833333	0,516245833	0,01845
-0,398666667	0,515254167	0,01845
-0,403083333	0,515158333	0,01845
-0,4045	0,515854167	0,01845
-0,403916667	0,51785	0,0172
-0,403333333	0,520495833	0,0172
-0,403166667	0,5227375	0,0172
-0,403166667	0,524470833	0,0172
-0,403583333	0,525433333	0,0172
-0,40425	0,525904167	0,01595
-0,405666667	0,525258333	0,01595
-0,407833333	0,5241875	0,01595
-0,409666667	0,522491667	0,01595
-0,41	0,520220833	0,01595
-0,409	0,517704167	0,014716667
-0,406833333	0,515866667	0,014716667
-0,403833333	0,5141125	0,014716667
-0,402166667	0,512525	0,014716667
-0,40375	0,5110875	0,014716667
-0,40775	0,50975	0,0135
-0,411333333	0,508445833	0,0135
-0,413416667	0,507079167	0,0135
-0,415666667	0,50555	0,0135
-0,419166667	0,503491667	0,0135
-0,423	0,500958333	0,012333333
-0,426166667	0,498	0,012333333
-0,428416667	0,4953	0,012333333
-0,430333333	0,4927125	0,012333333
-0,4325	0,490375	0,012333333
-0,435083333	0,488745833	0,01125
-0,4375	0,487654167	0,01125
-0,439833333	0,4865875	0,01125
-0,442666667	0,4861625	0,01125
-0,445	0,4858375	0,01125
-0,445916667	0,4839625	0,010266667
-0,446833333	0,481295833	0,010266667
-0,448	0,4780875	0,010266667
-0,449	0,474108333	0,010266667
-0,45025	0,469916667	0,010266667
-0,452583333	0,464716667	0,0094
-0,455833333	0,458316667	0,0094
-0,458583333	0,451695833	0,0094
-0,459583333	0,444575	0,0094
-0,460083333	0,438633333	0,0094
-0,4615	0,434279167	0,008683333
-0,462416667	0,430791667	0,008683333
-0,461416667	0,427508333	0,008683333
-0,459166667	0,423604167	0,008683333
-0,457916667	0,418875	0,008683333
-0,45875	0,413345833	0,008083333
-0,460833333	0,406141667	0,008083333
-0,462333333	0,398033333	0,008083333
-0,46325	0,390445833	0,008083333
-0,465916667	0,385775	0,008083333
-0,47	0,382408333	0,007616667
-0,4735	0,379804167	0,007616667
-0,476083333	0,377204167	0,007616667
-0,47775	0,374570833	0,007616667
-0,478166667	0,370175	0,007616667
-0,477916667	0,3649625	0,00725
-0,477833333	0,3587875	0,00725
-0,47925	0,353466667	0,00725
-0,48275	0,349004167	0,00725
-0,486583333	0,345683333	0,00725
-0,48975	0,342895833	0,006966667
-0,4925	0,340233333	0,006966667
-0,496	0,3373	0,006966667
-0,501	0,334620833	0,006966667
-0,50525	0,33115	0,006966667
-0,50675	0,327116667	0,006716667
-0,50725	0,3247	0,006716667
-0,508583333	0,323495833	0,006716667
-0,50975	0,321295833	0,006716667
-0,510166667	0,319233333	0,006716667
-0,510166667	0,317654167	0,0065
-0,50925	0,316170833	0,0065
-0,508333333	0,3144125	0,0065
-0,50975	0,311591667	0,0065
-0,51275	0,308229167	0,0065
-0,515	0,303966667	0,006266667
-0,516	0,299308333	0,006266667
-0,5175	0,295945833	0,006266667
-0,52075	0,293083333	0,006266667
-0,523583333	0,291133333	0,006266667
-0,524166667	0,290295833	0,006016667
-0,524666667	0,289583333	0,006016667
-0,525416667	0,288408333	0,006016667
-0,525583333	0,286920833	0,006016667
-0,524833333	0,284720833	0,006016667
-0,522916667	0,2813625	0,00575
-0,52075	0,276920833	0,00575
-0,51875	0,271791667	0,00575
-0,516416667	0,2669	0,00575
-0,513416667	0,262979167	0,00575
-0,510416667	0,261116667	0,005483333
-0,509333333	0,261733333	0,005483333
-0,509583333	0,263270833	0,005483333
-0,509833333	0,264404167	0,005483333
-0,5095	0,263816667	0,005483333
-0,5085	0,261570833	0,005233333
-0,508416667	0,257295833	0,005233333

library(lattice)
library(lmtest)
library(car)
library(MASS)
n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
mywarning <- ''
par6 <- as.numeric(par6)
if(is.na(par6)) {
par6 <- 12
mywarning = 'Warning: you did not specify the seasonality. The seasonal period was set to s = 12.'
}
par1 <- as.numeric(par1)
if(is.na(par1)) {
par1 <- 1
mywarning = 'Warning: you did not specify the column number of the endogenous series! The first column was selected by default.'
}
if (par4=='') par4 <- 0
par4 <- as.numeric(par4)
if (!is.numeric(par4)) par4 <- 0
if (par5=='') par5 <- 0
par5 <- as.numeric(par5)
if (!is.numeric(par5)) par5 <- 0
x <- na.omit(t(y))
k <- length(x[1,])
n <- length(x[,1])
x1 <- cbind(x[,par1], x[,1:k!=par1])
mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
colnames(x1) <- mycolnames #colnames(x)[par1]
x <- x1
if (par3 == 'First Differences'){
(n <- n -1)
x2 <- array(0, dim=c(n,k), dimnames=list(1:n, paste('(1-B)',colnames(x),sep='')))
for (i in 1:n) {
for (j in 1:k) {
x2[i,j] <- x[i+1,j] - x[i,j]
}
}
x <- x2
}
if (par3 == 'Seasonal Differences (s)'){
(n <- n - par6)
x2 <- array(0, dim=c(n,k), dimnames=list(1:n, paste('(1-Bs)',colnames(x),sep='')))
for (i in 1:n) {
for (j in 1:k) {
x2[i,j] <- x[i+par6,j] - x[i,j]
}
}
x <- x2
}
if (par3 == 'First and Seasonal Differences (s)'){
(n <- n -1)
x2 <- array(0, dim=c(n,k), dimnames=list(1:n, paste('(1-B)',colnames(x),sep='')))
for (i in 1:n) {
for (j in 1:k) {
x2[i,j] <- x[i+1,j] - x[i,j]
}
}
x <- x2
(n <- n - par6)
x2 <- array(0, dim=c(n,k), dimnames=list(1:n, paste('(1-Bs)',colnames(x),sep='')))
for (i in 1:n) {
for (j in 1:k) {
x2[i,j] <- x[i+par6,j] - x[i,j]
}
}
x <- x2
}
if(par4 > 0) {
x2 <- array(0, dim=c(n-par4,par4), dimnames=list(1:(n-par4), paste(colnames(x)[par1],'(t-',1:par4,')',sep='')))
for (i in 1:(n-par4)) {
for (j in 1:par4) {
x2[i,j] <- x[i+par4-j,par1]
}
}
x <- cbind(x[(par4+1):n,], x2)
n <- n - par4
}
if(par5 > 0) {
x2 <- array(0, dim=c(n-par5*par6,par5), dimnames=list(1:(n-par5*par6), paste(colnames(x)[par1],'(t-',1:par5,'s)',sep='')))
for (i in 1:(n-par5*par6)) {
for (j in 1:par5) {
x2[i,j] <- x[i+par5*par6-j*par6,par1]
}
}
x <- cbind(x[(par5*par6+1):n,], x2)
n <- n - par5*par6
}
if (par2 == 'Include Seasonal Dummies'){
x2 <- array(0, dim=c(n,par6-1), dimnames=list(1:n, paste('M', seq(1:(par6-1)), sep ='')))
for (i in 1:(par6-1)){
x2[seq(i,n,par6),i] <- 1
}
x <- cbind(x, x2)
}
if (par2 == 'Include Monthly Dummies'){
x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
for (i in 1:11){
x2[seq(i,n,12),i] <- 1
}
x <- cbind(x, x2)
}
if (par2 == 'Include Quarterly Dummies'){
x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
for (i in 1:3){
x2[seq(i,n,4),i] <- 1
}
x <- cbind(x, x2)
}
(k <- length(x[n,]))
if (par3 == 'Linear Trend'){
x <- cbind(x, c(1:n))
colnames(x)[k+1] <- 't'
}
print(x)
(k <- length(x[n,]))
head(x)
df <- as.data.frame(x)
(mylm <- lm(df))
(mysum <- summary(mylm))
if (n > n25) {
kp3 <- k + 3
nmkm3 <- n - k - 3
gqarr <- array(NA, dim=c(nmkm3-kp3+1,3))
numgqtests <- 0
numsignificant1 <- 0
numsignificant5 <- 0
numsignificant10 <- 0
for (mypoint in kp3:nmkm3) {
j <- 0
numgqtests <- numgqtests + 1
for (myalt in c('greater', 'two.sided', 'less')) {
j <- j + 1
gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value
}
if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1
if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1
if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1
}
gqarr
}
bitmap(file='test0.png')
plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index')
points(x[,1]-mysum$resid)
grid()
dev.off()
bitmap(file='test1.png')
plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index')
grid()
dev.off()
bitmap(file='test2.png')
sresid <- studres(mylm)
hist(sresid, freq=FALSE, main='Distribution of Studentized Residuals')
xfit<-seq(min(sresid),max(sresid),length=40)
yfit<-dnorm(xfit)
lines(xfit, yfit)
grid()
dev.off()
bitmap(file='test3.png')
densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals')
dev.off()
bitmap(file='test4.png')
qqPlot(mylm, main='QQ Plot')
grid()
dev.off()
(myerror <- as.ts(mysum$resid))
bitmap(file='test5.png')
dum <- cbind(lag(myerror,k=1),myerror)
dum
dum1 <- dum[2:length(myerror),]
dum1
z <- as.data.frame(dum1)
print(z)
plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals')
lines(lowess(z))
abline(lm(z))
grid()
dev.off()
bitmap(file='test6.png')
acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function')
grid()
dev.off()
bitmap(file='test7.png')
pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function')
grid()
dev.off()
bitmap(file='test8.png')
opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
plot(mylm, las = 1, sub='Residual Diagnostics')
par(opar)
dev.off()
if (n > n25) {
bitmap(file='test9.png')
plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint')
grid()
dev.off()
}
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE)
a<-table.row.end(a)
myeq <- colnames(x)[1]
myeq <- paste(myeq, '[t] = ', sep='')
for (i in 1:k){
if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '')
myeq <- paste(myeq, signif(mysum$coefficients[i,1],6), sep=' ')
if (rownames(mysum$coefficients)[i] != '(Intercept)') {
myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='')
if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='')
}
}
myeq <- paste(myeq, ' + e[t]')
a<-table.row.start(a)
a<-table.element(a, myeq)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, mywarning)
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,'Multiple Linear Regression - Ordinary Least Squares', 6, TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Variable',header=TRUE)
a<-table.element(a,'Parameter',header=TRUE)
a<-table.element(a,'S.D.',header=TRUE)
a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE)
a<-table.element(a,'2-tail p-value',header=TRUE)
a<-table.element(a,'1-tail p-value',header=TRUE)
a<-table.row.end(a)
for (i in 1:k){
a<-table.row.start(a)
a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE)
a<-table.element(a,formatC(signif(mysum$coefficients[i,1],5),format='g',flag='+'))
a<-table.element(a,formatC(signif(mysum$coefficients[i,2],5),format='g',flag=' '))
a<-table.element(a,formatC(signif(mysum$coefficients[i,3],4),format='e',flag='+'))
a<-table.element(a,formatC(signif(mysum$coefficients[i,4],4),format='g',flag=' '))
a<-table.element(a,formatC(signif(mysum$coefficients[i,4]/2,4),format='g',flag=' '))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable2.tab')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'Multiple R',1,TRUE)
a<-table.element(a,formatC(signif(sqrt(mysum$r.squared),6),format='g',flag=' '))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'R-squared',1,TRUE)
a<-table.element(a,formatC(signif(mysum$r.squared,6),format='g',flag=' '))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'Adjusted R-squared',1,TRUE)
a<-table.element(a,formatC(signif(mysum$adj.r.squared,6),format='g',flag=' '))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'F-TEST (value)',1,TRUE)
a<-table.element(a,formatC(signif(mysum$fstatistic[1],6),format='g',flag=' '))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE)
a<-table.element(a, signif(mysum$fstatistic[2],6))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE)
a<-table.element(a, signif(mysum$fstatistic[3],6))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'p-value',1,TRUE)
a<-table.element(a,formatC(signif(1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3]),6),format='g',flag=' '))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'Residual Standard Deviation',1,TRUE)
a<-table.element(a,formatC(signif(mysum$sigma,6),format='g',flag=' '))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'Sum Squared Residuals',1,TRUE)
a<-table.element(a,formatC(signif(sum(myerror*myerror),6),format='g',flag=' '))
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable3.tab')
myr <- as.numeric(mysum$resid)
myr
a <-table.start()
a <- table.row.start(a)
a <- table.element(a,'Menu of Residual Diagnostics',2,TRUE)
a <- table.row.end(a)
a <- table.row.start(a)
a <- table.element(a,'Description',1,TRUE)
a <- table.element(a,'Link',1,TRUE)
a <- table.row.end(a)
a <- table.row.start(a)
a <-table.element(a,'Histogram',1,header=TRUE)
a <- table.element(a,hyperlink( paste('https://supernova.wessa.net/rwasp_histogram.wasp?convertgetintopost=1&data=',paste(as.character(mysum$resid),sep='',collapse=' '),sep='') ,'Compute','Click here to examine the Residuals.'),1)
a <- table.row.end(a)
a <- table.row.start(a)
a <-table.element(a,'Central Tendency',1,header=TRUE)
a <- table.element(a,hyperlink( paste('https://supernova.wessa.net/rwasp_centraltendency.wasp?convertgetintopost=1&data=',paste(as.character(mysum$resid),sep='',collapse=' '),sep='') ,'Compute','Click here to examine the Residuals.'),1)
a <- table.row.end(a)
a <- table.row.start(a)
a <-table.element(a,'QQ Plot',1,header=TRUE)
a <- table.element(a,hyperlink( paste('https://supernova.wessa.net/rwasp_fitdistrnorm.wasp?convertgetintopost=1&data=',paste(as.character(mysum$resid),sep='',collapse=' '),sep='') ,'Compute','Click here to examine the Residuals.'),1)
a <- table.row.end(a)
a <- table.row.start(a)
a <-table.element(a,'Kernel Density Plot',1,header=TRUE)
a <- table.element(a,hyperlink( paste('https://supernova.wessa.net/rwasp_density.wasp?convertgetintopost=1&data=',paste(as.character(mysum$resid),sep='',collapse=' '),sep='') ,'Compute','Click here to examine the Residuals.'),1)
a <- table.row.end(a)
a <- table.row.start(a)
a <-table.element(a,'Skewness/Kurtosis Test',1,header=TRUE)
a <- table.element(a,hyperlink( paste('https://supernova.wessa.net/rwasp_skewness_kurtosis.wasp?convertgetintopost=1&data=',paste(as.character(mysum$resid),sep='',collapse=' '),sep='') ,'Compute','Click here to examine the Residuals.'),1)
a <- table.row.end(a)
a <- table.row.start(a)
a <-table.element(a,'Skewness-Kurtosis Plot',1,header=TRUE)
a <- table.element(a,hyperlink( paste('https://supernova.wessa.net/rwasp_skewness_kurtosis_plot.wasp?convertgetintopost=1&data=',paste(as.character(mysum$resid),sep='',collapse=' '),sep='') ,'Compute','Click here to examine the Residuals.'),1)
a <- table.row.end(a)
a <- table.row.start(a)
a <-table.element(a,'Harrell-Davis Plot',1,header=TRUE)
a <- table.element(a,hyperlink( paste('https://supernova.wessa.net/rwasp_harrell_davis.wasp?convertgetintopost=1&data=',paste(as.character(mysum$resid),sep='',collapse=' '),sep='') ,'Compute','Click here to examine the Residuals.'),1)
a <- table.row.end(a)
a <- table.row.start(a)
a <-table.element(a,'Bootstrap Plot -- Central Tendency',1,header=TRUE)
a <- table.element(a,hyperlink( paste('https://supernova.wessa.net/rwasp_bootstrapplot1.wasp?convertgetintopost=1&data=',paste(as.character(mysum$resid),sep='',collapse=' '),sep='') ,'Compute','Click here to examine the Residuals.'),1)
a <- table.row.end(a)
a <- table.row.start(a)
a <-table.element(a,'Blocked Bootstrap Plot -- Central Tendency',1,header=TRUE)
a <- table.element(a,hyperlink( paste('https://supernova.wessa.net/rwasp_bootstrapplot.wasp?convertgetintopost=1&data=',paste(as.character(mysum$resid),sep='',collapse=' '),sep='') ,'Compute','Click here to examine the Residuals.'),1)
a <- table.row.end(a)
a <- table.row.start(a)
a <-table.element(a,'(Partial) Autocorrelation Plot',1,header=TRUE)
a <- table.element(a,hyperlink( paste('https://supernova.wessa.net/rwasp_autocorrelation.wasp?convertgetintopost=1&data=',paste(as.character(mysum$resid),sep='',collapse=' '),sep='') ,'Compute','Click here to examine the Residuals.'),1)
a <- table.row.end(a)
a <- table.row.start(a)
a <-table.element(a,'Spectral Analysis',1,header=TRUE)
a <- table.element(a,hyperlink( paste('https://supernova.wessa.net/rwasp_spectrum.wasp?convertgetintopost=1&data=',paste(as.character(mysum$resid),sep='',collapse=' '),sep='') ,'Compute','Click here to examine the Residuals.'),1)
a <- table.row.end(a)
a <- table.row.start(a)
a <-table.element(a,'Tukey lambda PPCC Plot',1,header=TRUE)
a <- table.element(a,hyperlink( paste('https://supernova.wessa.net/rwasp_tukeylambda.wasp?convertgetintopost=1&data=',paste(as.character(mysum$resid),sep='',collapse=' '),sep='') ,'Compute','Click here to examine the Residuals.'),1)
a <- table.row.end(a)
a <- table.row.start(a)
a <-table.element(a,'Box-Cox Normality Plot',1,header=TRUE)
a <- table.element(a,hyperlink( paste('https://supernova.wessa.net/rwasp_boxcoxnorm.wasp?convertgetintopost=1&data=',paste(as.character(mysum$resid),sep='',collapse=' '),sep='') ,'Compute','Click here to examine the Residuals.'),1)
a <- table.row.end(a)
a <- table.row.start(a)
a <- table.element(a,'Summary Statistics',1,header=TRUE)
a <- table.element(a,hyperlink( paste('https://supernova.wessa.net/rwasp_summary1.wasp?convertgetintopost=1&data=',paste(as.character(mysum$resid),sep='',collapse=' '),sep='') ,'Compute','Click here to examine the Residuals.'),1)
a <- table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable7.tab')
if(n < 200) {
a<-table.start()
a<-table.row.start(a)
a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'Time or Index', 1, TRUE)
a<-table.element(a, 'Actuals', 1, TRUE)
a<-table.element(a, 'Interpolation
Forecast', 1, TRUE)
a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE)
a<-table.row.end(a)
for (i in 1:n) {
a<-table.row.start(a)
a<-table.element(a,i, 1, TRUE)
a<-table.element(a,formatC(signif(x[i],6),format='g',flag=' '))
a<-table.element(a,formatC(signif(x[i]-mysum$resid[i],6),format='g',flag=' '))
a<-table.element(a,formatC(signif(mysum$resid[i],6),format='g',flag=' '))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable4.tab')
if (n > n25) {
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'p-values',header=TRUE)
a<-table.element(a,'Alternative Hypothesis',3,header=TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'breakpoint index',header=TRUE)
a<-table.element(a,'greater',header=TRUE)
a<-table.element(a,'2-sided',header=TRUE)
a<-table.element(a,'less',header=TRUE)
a<-table.row.end(a)
for (mypoint in kp3:nmkm3) {
a<-table.row.start(a)
a<-table.element(a,mypoint,header=TRUE)
a<-table.element(a,formatC(signif(gqarr[mypoint-kp3+1,1],6),format='g',flag=' '))
a<-table.element(a,formatC(signif(gqarr[mypoint-kp3+1,2],6),format='g',flag=' '))
a<-table.element(a,formatC(signif(gqarr[mypoint-kp3+1,3],6),format='g',flag=' '))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable5.tab')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Description',header=TRUE)
a<-table.element(a,'# significant tests',header=TRUE)
a<-table.element(a,'% significant tests',header=TRUE)
a<-table.element(a,'OK/NOK',header=TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'1% type I error level',header=TRUE)
a<-table.element(a,signif(numsignificant1,6))
a<-table.element(a,formatC(signif(numsignificant1/numgqtests,6),format='g',flag=' '))
if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK'
a<-table.element(a,dum)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'5% type I error level',header=TRUE)
a<-table.element(a,signif(numsignificant5,6))
a<-table.element(a,signif(numsignificant5/numgqtests,6))
if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK'
a<-table.element(a,dum)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'10% type I error level',header=TRUE)
a<-table.element(a,signif(numsignificant10,6))
a<-table.element(a,signif(numsignificant10/numgqtests,6))
if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK'
a<-table.element(a,dum)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable6.tab')
}
}
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Ramsey RESET F-Test for powers (2 and 3) of fitted values',1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
reset_test_fitted <- resettest(mylm,power=2:3,type='fitted')
a<-table.element(a,paste('
',RC.texteval('reset_test_fitted'),'
',sep=''))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Ramsey RESET F-Test for powers (2 and 3) of regressors',1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
reset_test_regressors <- resettest(mylm,power=2:3,type='regressor')
a<-table.element(a,paste('
',RC.texteval('reset_test_regressors'),'
',sep=''))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Ramsey RESET F-Test for powers (2 and 3) of principal components',1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
reset_test_principal_components <- resettest(mylm,power=2:3,type='princomp')
a<-table.element(a,paste('
',RC.texteval('reset_test_principal_components'),'
',sep=''))
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable8.tab')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Variance Inflation Factors (Multicollinearity)',1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
vif <- vif(mylm)
a<-table.element(a,paste('
',RC.texteval('vif'),'
',sep=''))
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable9.tab')