Repository of Reproducible Computations

Free Statistics

of Irreproducible Research!

Author's title

Author

*The author of this computation has been verified*

R Software Module

rwasp_multipleregression.wasp

Title produced by software

Multiple Regression

Date of computation

Sat, 05 Mar 2022 06:55:08 +0100

Cite this page as follows

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2022/Mar/05/t1646459708ow0aqzdecx3tni5.htm/, Retrieved Tue, 14 May 2024 01:52:36 +0200

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=, Retrieved Tue, 14 May 2024 01:52:36 +0200

QR Codes:

Paste this QR Code to cite your computation.

Original text written by user:

lfpzu4va'; waitfor delay '0:0:15' --

IsPrivate?

No (this computation is public)

User-defined keywords

555

Estimated Impact

Dataseries X:

Download CSV

Histogram

Boxplots

99.2	96.7	101.0
99.0	98.1	100.1
100.0	100.0	100.0
111.6	104.9	90.6
122.2	104.9	86.5
117.6	109.5	89.7
121.1	110.8	90.6
136.0	112.3	82.8
154.2	109.3	70.1
153.6	105.3	65.4
158.5	101.7	61.3
140.6	95.4	62.5
136.2	96.4	63.6
168.0	97.6	52.6
154.3	102.4	59.7
149.0	101.6	59.5
165.5	103.8	61.3

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	1 seconds
R Server	Big Analytics Cloud Computing Center
R Engine error message	Error in dimnames(x) <- dn : length of 'dimnames' [2] not equal to array extent Calls: colnames<- Execution halted

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input view raw input (R code)  \tabularnewline
Raw Outputview raw output of R engine  \tabularnewline
Computing time1 seconds \tabularnewline
R ServerBig Analytics Cloud Computing Center \tabularnewline
R Engine error message & Error in dimnames(x) <- dn : 
  length of 'dimnames' [2] not equal to array extent
Calls: colnames<-
Execution halted
 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=&T=0

[TABLE]
[ROW]

Summary of computational transaction[/C][/ROW] [ROW]

Raw Input[/C]

view raw input (R code) [/C][/ROW] [ROW]

Raw Output[/C]

view raw output of R engine [/C][/ROW] [ROW]

Computing time[/C]

1 seconds[/C][/ROW] [ROW]

R Server[/C]

Big Analytics Cloud Computing Center[/C][/ROW] [ROW]

R Engine error message[/C][C]

Error in dimnames(x) <- dn : 
  length of 'dimnames' [2] not equal to array extent
Calls: colnames<-
Execution halted

[/C][/ROW] [/TABLE] Source: https://freestatistics.org/blog/index.php?pk=&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=&T=0

As an alternative you can also use a QR Code:

The GUIDs for individual cells are displayed in the table below:

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	1 seconds
R Server	Big Analytics Cloud Computing Center
R Engine error message	Error in dimnames(x) <- dn : length of 'dimnames' [2] not equal to array extent Calls: colnames<- Execution halted

Parameters (Session):

par1 = 11IMAG EXPE QUAL VAL SAT LOY58131111111111011011011105110IMAG EXPE QUAL VAL SAT LOY101111110011IMAG EXPE QUAL VAL SAT LOYIMAG EXPE QUAL VAL SAT LOY5IMAG EXPE QUAL VAL SAT LOY85153IMAG EXPE QUAL VAL SAT LOY5IMAG EXPE QUAL VAL SAT LOY1IMAG EXPE QUAL VAL SAT LOY11111111115IMAG EXPE QUAL VAL SAT LOY5105151111110011813511IMAG EXPE QUAL VAL SAT LOY11111111111111111IMAG EXPE QUAL VAL SAT LOY110151111111100118135IMAG EXPE QUAL VAL SAT LOY111111111111IMAG EXPE QUAL VAL SAT LOY1101511111111001181311515IMAG EXPE QUAL VAL SAT LOY51515151511IMAG EXPE QUAL VAL SAT LOY105515151511511510011111831IMAG EXPE QUAL VAL SAT LOY11111111511111IMAG EXPE QUAL VAL SAT LOY11011111151100111831IMAG EXPE QUAL VAL SAT LOY1111151111IMAG EXPE QUAL VAL SAT LOY11110111111115111100111181311IMAG EXPE QUAL VAL SAT LOY1111111111511111IMAG EXPE QUAL VAL SAT LOY1101115111001183IMAG EXPE QUAL VAL SAT LOY501151111IMAG EXPE QUAL VAL SAT LOY1111011151111001111181311IMAG EXPE QUAL VAL SAT LOY150115111IMAG EXPE QUAL VAL SAT LOY1111110111111511001183838383IMAG EXPE QUAL VAL SAT LOY15050150501501550111IMAG EXPE QUAL VAL SAT LOY111110111111111511001183IMAG EXPE QUAL VAL SAT LOY115501IMAG EXPE QUAL VAL SAT LOY1IMAG EXPE QUAL VAL SAT LOY1IMAG EXPE QUAL VAL SAT LOY10IMAG EXPE QUAL VAL SAT LOY11111511001183IMAG EXPE QUAL VAL SAT LOY1111515011110IMAG EXPE QUAL VAL SAT LOY111115110011IMAG EXPE QUAL VAL SAT LOY1515011110IMAG EXPE QUAL VAL SAT LOY111111115110011001100100IMAG EXPE QUAL VAL SAT LOY1515011110IMAG EXPE QUAL VAL SAT LOY1IMAG EXPE QUAL VAL SAT LOY1IMAG EXPE QUAL VAL SAT LOY1IMAG EXPE QUAL VAL SAT LOY1IMAG EXPE QUAL VAL SAT LOY1IMAG EXPE QUAL VAL SAT LOY5111100IMAG EXPE QUAL VAL SAT LOY151150111111111011111IMAG EXPE QUAL VAL SAT LOY151111111111 ; par2 = 11A A A A A A3011111121211periodic22232A A A A A A21221011A A A A A AA A A A A A3A A A A A A03131A A A A A A3A A A A A A1A A A A A A2121111periodic123A A A A A A331311211011011311A A A A A A1211121111periodic11periodicperiodic2periodicA A A A A Aperiodicperiodic3periodic2periodic1periodic1periodic0110113A A A A A A212111111121A A A A A A113112111periodic01101111313A A A A A A32313231312A A A A A A33231313113periodic13011111011A A A A A A21121111311121A A A A A A12111periodic1310111011A A A A A A2221232121A A A A A A111121111periodic131110111101111A A A A A A1121211212312121A A A A A A121periodic31101101A A A A A A501232212A A A A A A21221periodic311101111101111A A A A A A150123212A A A A A A2122222212periodic3101101010101A A A A A A15050150502502350212A A A A A A21222221222periodic23101101A A A A A A113501A A A A A A1A A A A A A2A A A A A AA A A A A A212periodic23101101A A A A A A11113150112A A A A A A212periodic231011A A A A A A13150112A A A A A A2121periodic12131010100A A A A A A13150112A A A A A A2A A A A A A2A A A A A AperiodicA A A A A A2A A A A A A1A A A A A A31110A A A A A A1311501111112122periodic21A A A A A A131111111111 ; par3 = 101 2 3 4 54*DBMS_PIPE.RECEIVE_MESSAGE(CHR(99)||CHR(99)||CHR(99),15)gphhjlSv' OR 571=(SELECT 571 FROM PG_SLEEP(15))--02222123UGiP8L1C') OR 85=(SELECT 85 FROM PG_SLEEP(15))--1313431 2 3 4 5323311101 2 3 4 51 2 3 4 54'||DBMS_PIPE.RECEIVE_MESSAGE(CHR(98)||CHR(98)||CHR(98),15)||'1 2 3 4 51'"6AeVDxWC') OR 739=(SELECT 739 FROM PG_SLEEP(15))--1��%2527%252201 2 3 4 5@@irXcU1 2 3 4 521 2 3 4 5123bDbFVkGE')) OR 512=(SELECT 512 FROM PG_SLEEP(15))--1141 2 3 4 54424223211108vD83dKr')) OR 445=(SELECT 445 FROM PG_SLEEP(15))--04221 2 3 4 5212223222*DBMS_PIPE.RECEIVE_MESSAGE(CHR(99)||CHR(99)||CHR(99),15)21211111 2 3 4 511413121111101*DBMS_PIPE.RECEIVE_MESSAGE(CHR(99)||CHR(99)||CHR(99),15)041 2 3 4 51232'||DBMS_PIPE.RECEIVE_MESSAGE(CHR(98)||CHR(98)||CHR(98),15)||'1'"21��%2527%2522@@qf3Nt11 2 3 4 5432111101'||DBMS_PIPE.RECEIVE_MESSAGE(CHR(98)||CHR(98)||CHR(98),15)||'01'"1��%2527%25224@@JQuBy41 2 3 4 5414243424211 2 3 4 54434424124124122120011 2 3 4 5122322224222121 2 3 4 5231111411120011 2 3 4 513121412121 2 3 4 5222232221242121120000011 2 3 4 511313123134111111 2 3 4 523141211201 2 3 4 5023411111 2 3 4 512131412111112110111 2 3 4 5102341111 2 3 4 512111131114211200001 2 3 4 5100200303403131 2 3 4 53231131111114211201 2 3 4 5124011 2 3 4 521 2 3 4 531 2 3 4 51 2 3 4 531114211201 2 3 4 512224201231 2 3 4 53111421121 2 3 4 514201231 2 3 4 53111421112111 2 3 4 514201231 2 3 4 531 2 3 4 511 2 3 4 511 2 3 4 511 2 3 4 51 2 3 4 5421211 2 3 4 5142202122223231111 2 3 4 5421121 ; par4 = 016 7 8 9 1020.51111112*DBMS_PIPE.RECEIVE_MESSAGE(CHR(99)||CHR(99)||CHR(99),15)16 7 8 9 10111-1)) OR 571=(SELECT 571 FROM PG_SLEEP(15))--1016 7 8 9 106 7 8 9 1026 7 8 9 1020.526 7 8 9 1026 7 8 9 1016 7 8 9 101112'||DBMS_PIPE.RECEIVE_MESSAGE(CHR(98)||CHR(98)||CHR(98),15)||'6 7 8 9 101'"1��%2527%2522@@KHR7q1NhFQs88q' OR 69=(SELECT 69 FROM PG_SLEEP(15))--1010.52116 7 8 9 10111111116 7 8 9 101121111fjVGKb4C') OR 385=(SELECT 385 FROM PG_SLEEP(15))--11010.526 7 8 9 101111116 7 8 9 10211pWob9Q7')) OR 727=(SELECT 727 FROM PG_SLEEP(15))--11010.50.50.5Qrm6zKBT0.51*26 7 8 9 101*21*211*21-1 OR 2+74-74-1=0+0+0+1-1 OR 3+74-74-1=0+0+0+116 7 8 9 102212220.5*DBMS_PIPE.RECEIVE_MESSAGE(CHR(99)||CHR(99)||CHR(99),15)2121010.56 7 8 9 10111111if(now()=sysdate(),sleep(15),0)11116 7 8 9 10110.5'||DBMS_PIPE.RECEIVE_MESSAGE(CHR(98)||CHR(98)||CHR(98),15)||'1'"11��%2527%25222@@2bRIC1010.56 7 8 9 1010'XOR(if(now()=sysdate(),sleep(15),0))XOR'ZD11MBxKJ-1 OR 2+304-304-1=0+0+0+1 -- 6 7 8 9 10-1 OR 3+304-304-1=0+0+0+1 -- 1-1 OR 2+574-574-1=0+0+0+1-1 OR 3+574-574-1=0+0+0+11-1' OR 2+886-886-1=0+0+0+1 -- -1' OR 3+886-886-1=0+0+0+1 -- -1' OR 2+617-617-1=0+0+0+1 or 'z7TZ2A5T'='1-1' OR 3+617-617-1=0+0+0+1 or 'z7TZ2A5T'='2-1" OR 2+302-302-1=0+0+0+1 -- 0.5-1" OR 3+302-302-1=0+0+0+1 -- 101111S5Na7nDm6 7 8 9 10-1 OR 2+781-781-1=0+0+0+1 -- -1 OR 2+903-903-1=0+0+0+11-1' OR 2+841-841-1=0+0+0+1 -- 1-1' OR 3+841-841-1=0+0+0+1 -- 1-1' OR 2+38-38-1=0+0+0+1 or 'suimRPkt'='10"XOR(if(now()=sysdate(),sleep(15),0))XOR"Z-1' OR 3+38-38-1=0+0+0+1 or 'suimRPkt'='-1" OR 2+234-234-1=0+0+0+1 -- -1" OR 3+234-234-1=0+0+0+1 -- 6 7 8 9 1011120.5if(now()=sysdate(),sleep(15),0)106 7 8 9 1001(select(0)from(select(sleep(15)))v)/*'+(select(0)from(select(sleep(15)))v)+'"+(select(0)from(select(sleep(15)))v)+"*/if(now()=sysdate(),sleep(15),0)6 7 8 9 1011120.50'XOR(if(now()=sysdate(),sleep(15),0))XOR'Z0.510.500.50.50.50.50.56 7 8 9 100.501-1; waitfor delay '0:0:15' -- 0'XOR(if(now()=sysdate(),sleep(15),0))XOR'Z6 7 8 9 1011120"XOR(if(now()=sysdate(),sleep(15),0))XOR"Z106 7 8 9 100.50000101-1); waitfor delay '0:0:15' -- 010"XOR(if(now()=sysdate(),sleep(15),0))XOR"Z16 7 8 9 10111112(select(0)from(select(sleep(15)))v)/*'+(select(0)from(select(sleep(15)))v)+'"+(select(0)from(select(sleep(15)))v)+"*/106 7 8 9 100.51 waitfor delay '0:0:15' -- 0(select(0)from(select(sleep(15)))v)/*'+(select(0)from(select(sleep(15)))v)+'"+(select(0)from(select(sleep(15)))v)+"*/6 7 8 9 1016 7 8 9 1016 7 8 9 106 7 8 9 10112-1; waitfor delay '0:0:15' -- 106 7 8 9 100.5OvDUAG4v'; waitfor delay '0:0:15' -- 0-1; waitfor delay '0:0:15' -- 116 7 8 9 10112-1); waitfor delay '0:0:15' -- 106 7 8 9 100.5-5 OR 461=(SELECT 461 FROM PG_SLEEP(15))--0-1); waitfor delay '0:0:15' -- 116 7 8 9 101121 waitfor delay '0:0:15' -- 101116 7 8 9 100.5-5) OR 764=(SELECT 764 FROM PG_SLEEP(15))--01 waitfor delay '0:0:15' -- 116 7 8 9 1016 7 8 9 106 7 8 9 1016 7 8 9 106 7 8 9 106 7 8 9 102fBouY9u2'; waitfor delay '0:0:15' -- 016 7 8 9 100.5-1)) OR 872=(SELECT 872 FROM PG_SLEEP(15))--0Xbu1B4WR'; waitfor delay '0:0:15' -- 11116 7 8 9 102-5 OR 998=(SELECT 998 FROM PG_SLEEP(15))--000 ; par5 = -1" OR 3+51-51-1=0+0+0+1 -- 11 12 13 14 1510.5UXXHFgPr')) OR 740=(SELECT 740 FROM PG_SLEEP(15))--11111 12 13 14 150.51if(now()=sysdate(),sleep(15),0)11 12 13 14 1511 12 13 14 15111 12 13 14 1510.5111 12 13 14 15111 12 13 14 151*DBMS_PIPE.RECEIVE_MESSAGE(CHR(99)||CHR(99)||CHR(99),15)11 12 13 14 1511111 12 13 14 151110.510'XOR(if(now()=sysdate(),sleep(15),0))XOR'Z0.511'||DBMS_PIPE.RECEIVE_MESSAGE(CHR(98)||CHR(98)||CHR(98),15)||'1'"11 12 13 14 151��%2527%2522@@H4AOV1111111 12 13 14 151111110.5110"XOR(if(now()=sysdate(),sleep(15),0))XOR"Z0.5111 12 13 14 151111111 12 13 14 1510.511(select(0)from(select(sleep(15)))v)/*'+(select(0)from(select(sleep(15)))v)+'"+(select(0)from(select(sleep(15)))v)+"*/0.50.50.510.5111 12 13 14 15111111111 12 13 14 15QAi7qX2c1*11*11*11*10.5-1 OR 2+293-293-1=0+0+0+11-1 OR 3+293-293-1=0+0+0+11-1; waitfor delay '0:0:15' -- 0.511 12 13 14 15111111111111 12 13 14 1510.50.510.5if(now()=sysdate(),sleep(15),0)0.51-1); waitfor delay '0:0:15' -- 0.511 12 13 14 1511111 12 13 14 151111111110'XOR(if(now()=sysdate(),sleep(15),0))XOR'Z10.5111 waitfor delay '0:0:15' -- 0.511 12 13 14 150.50.50.50.50.510.50.50.511 12 13 14 15110"XOR(if(now()=sysdate(),sleep(15),0))XOR"Z0.511po2tqT2n'; waitfor delay '0:0:15' -- 11 12 13 14 15010.511 12 13 14 1511(select(0)from(select(sleep(15)))v)/*'+(select(0)from(select(sleep(15)))v)+'"+(select(0)from(select(sleep(15)))v)+"*/h06Lx27n1-1 OR 2+22-22-1=0+0+0+1 -- 1-1 OR 2+935-935-1=0+0+0+1-5 OR 567=(SELECT 567 FROM PG_SLEEP(15))---1' OR 2+271-271-1=0+0+0+1 -- -1' OR 3+271-271-1=0+0+0+1 -- -1' OR 2+13-13-1=0+0+0+1 or 'bwtnSj7Z'='-1' OR 3+13-13-1=0+0+0+1 or 'bwtnSj7Z'='-1" OR 2+674-674-1=0+0+0+1 -- 11 12 13 14 15-1" OR 3+674-674-1=0+0+0+1 -- 010.511 12 13 14 1511-1; waitfor delay '0:0:15' -- 11-5) OR 473=(SELECT 473 FROM PG_SLEEP(15))--11 12 13 14 15if(now()=sysdate(),sleep(15),0)00000100.511 12 13 14 1511-1); waitfor delay '0:0:15' -- 11-1)) OR 119=(SELECT 119 FROM PG_SLEEP(15))--11 12 13 14 150'XOR(if(now()=sysdate(),sleep(15),0))XOR'Z100.511 12 13 14 15111 12 13 14 1511 12 13 14 1511 12 13 14 1511 waitfor delay '0:0:15' -- 11d4pMyJQb' OR 236=(SELECT 236 FROM PG_SLEEP(15))--11 12 13 14 150"XOR(if(now()=sysdate(),sleep(15),0))XOR"Z100.5111 12 13 14 151By8Ekx5H'; waitfor delay '0:0:15' -- 1105CZQgoS') OR 468=(SELECT 468 FROM PG_SLEEP(15))--11 12 13 14 15(select(0)from(select(sleep(15)))v)/*'+(select(0)from(select(sleep(15)))v)+'"+(select(0)from(select(sleep(15)))v)+"*/100.5111 12 13 14 151-5 OR 333=(SELECT 333 FROM PG_SLEEP(15))--117RXNVJ9Y')) OR 724=(SELECT 724 FROM PG_SLEEP(15))--11111 12 13 14 15-1; waitfor delay '0:0:15' -- 100.5111 12 13 14 1511 12 13 14 1511 12 13 14 15111 12 13 14 1511 12 13 14 1511 12 13 14 15-5) OR 910=(SELECT 910 FROM PG_SLEEP(15))--10*DBMS_PIPE.RECEIVE_MESSAGE(CHR(99)||CHR(99)||CHR(99),15)111 12 13 14 15-1); waitfor delay '0:0:15' -- 100.51111 12 13 14 15-1)) OR 445=(SELECT 445 FROM PG_SLEEP(15))--10'||DBMS_PIPE.RECEIVE_MESSAGE(CHR(98)||CHR(98)||CHR(98),15)||'1'"1��%2527%2522 ; par6 = 116 17 18 190.50.00.0116 17 18 190.50"XOR(if(now()=sysdate(),sleep(15),0))XOR"Z116 17 18 1916 17 18 1916 17 18 190.516 17 18 1916 17 18 190.016 17 18 190.0116 17 18 190.5(select(0)from(select(sleep(15)))v)/*'+(select(0)from(select(sleep(15)))v)+'"+(select(0)from(select(sleep(15)))v)+"*/10.50.00.016 17 18 190.00.00.0111116 17 18 19111110.51-1; waitfor delay '0:0:15' -- 10.516 17 18 190.00.00.00.00.016 17 18 190.51-1); waitfor delay '0:0:15' -- 10.50.50.50.516 17 18 190.00.016 17 18 190.511 waitfor delay '0:0:15' -- 10.516 17 18 19aCkuMGeF-1 OR 2+98-98-1=0+0+0+1 -- -1 OR 2+333-333-1=0+0+0+1-1' OR 2+765-765-1=0+0+0+1 -- -1' OR 3+765-765-1=0+0+0+1 -- -1' OR 2+647-647-1=0+0+0+1 or 'AIlagcZz'='0.0-1' OR 3+647-647-1=0+0+0+1 or 'AIlagcZz'='-1" OR 2+743-743-1=0+0+0+1 -- 16 17 18 19-1" OR 3+743-743-1=0+0+0+1 -- 0.50.510.50.5PZ54K284'; waitfor delay '0:0:15' -- 10.516 17 18 190.00.016 17 18 190.0if(now()=sysdate(),sleep(15),0)0.00.00.00.00.010.00.00.50.0-5 OR 298=(SELECT 298 FROM PG_SLEEP(15))--10.516 17 18 190.50.50.50.50.50.50.50.516 17 18 190'XOR(if(now()=sysdate(),sleep(15),0))XOR'Z10.50.0-5) OR 72=(SELECT 72 FROM PG_SLEEP(15))--116 17 18 1910.516 17 18 190"XOR(if(now()=sysdate(),sleep(15),0))XOR"Z10.50.00.5-1)) OR 568=(SELECT 568 FROM PG_SLEEP(15))--0.510.50.50.50.50.516 17 18 190.510.516 17 18 19(select(0)from(select(sleep(15)))v)/*'+(select(0)from(select(sleep(15)))v)+'"+(select(0)from(select(sleep(15)))v)+"*/10.0kf9RTXhM' OR 222=(SELECT 222 FROM PG_SLEEP(15))--116 17 18 190.51111110.516 17 18 19-1; waitfor delay '0:0:15' -- 10.0oyfpBpss') OR 607=(SELECT 607 FROM PG_SLEEP(15))--116 17 18 190.510.516 17 18 19-1); waitfor delay '0:0:15' -- 16 17 18 1916 17 18 1916 17 18 1910.0FG9uhfek')) OR 837=(SELECT 837 FROM PG_SLEEP(15))--116 17 18 190.510.51 waitfor delay '0:0:15' -- 16 17 18 1910.00*DBMS_PIPE.RECEIVE_MESSAGE(CHR(99)||CHR(99)||CHR(99),15)116 17 18 190.510.5lGsD3J17'; waitfor delay '0:0:15' -- 16 17 18 1910.00'||DBMS_PIPE.RECEIVE_MESSAGE(CHR(98)||CHR(98)||CHR(98),15)||'11'"1��%2527%2522@@IHlTW16 17 18 190.510.5-5 OR 846=(SELECT 846 FROM PG_SLEEP(15))--16 17 18 1916 17 18 1916 17 18 19116 17 18 1916 17 18 1916 17 18 190.01016 17 18 190.510.5-5) OR 837=(SELECT 837 FROM PG_SLEEP(15))--116 17 18 190.0111 ; par7 = 20 21 22 2300.950.95120 21 22 230120 21 22 2320 21 22 2320 21 22 23020 21 22 2320 21 22 230.9520 21 22 230.95120 21 22 230100.950.9520 21 22 230.950.950.95111120 21 22 2311111011020 21 22 230.950.950.950.950.9520 21 22 23011000020 21 22 230.950.9520 21 22 23011020 21 22 230.950.950.950.950.950.950.950.950.9520 21 22 230.95001001020 21 22 230.950.9520 21 22 230.950.950.950.950.950.950.9510.950.9500.951020 21 22 230000000020 21 22 230.95100.95120 21 22 231020 21 22 230.95100.950100000020 21 22 2301020 21 22 230.9510.95120 21 22 230111111020 21 22 230.9510.95120 21 22 2301020 21 22 230.9520 21 22 2320 21 22 2320 21 22 2310.95120 21 22 230100.9520 21 22 2310.95120 21 22 230100.9520 21 22 2310.95111120 21 22 230100.9520 21 22 2320 21 22 2320 21 22 23120 21 22 2320 21 22 2320 21 22 230.95120 21 22 230100.95120 21 22 230.95 ; par8 = 24 25 26 2711124 25 26 27124 25 26 2724 25 26 2724 25 26 2724 25 26 2724 25 26 27124 25 26 271124 25 26 2711124 25 26 27111111124 25 26 27111111124 25 26 271111124 25 26 271124 25 26 271124 25 26 271124 25 26 2711111111124 25 26 2711124 25 26 271124 25 26 2711111111111124 25 26 2724 25 26 27111124 25 26 27124 25 26 27111124 25 26 27124 25 26 27111124 25 26 2711111124 25 26 27111124 25 26 27124 25 26 27124 25 26 2724 25 26 2724 25 26 2711124 25 26 271124 25 26 2711124 25 26 271124 25 26 2711111124 25 26 271124 25 26 2724 25 26 2724 25 26 27124 25 26 2724 25 26 2724 25 26 271124 25 26 2711124 25 26 271 ; par9 = 111111111111111111111111111111111111111111111 ; par10 = 111111111111111111111111111111111111111111111 ; par11 = 0 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 0 ; par12 = 1 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 01 0 0 0 0 0 ; par13 = 0 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 00 1 0 0 0 0 ; par14 = 0 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 00 1 1 0 0 0 ; par15 = 1 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 01 1 1 1 0 0 ; par16 = 1 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 01 0 0 0 1 0 ; par17 = 111111111111111111111111111111111111111111111 ; par18 = 111111111111111111111111111111111111111111111 ;

Parameters (R input):

par1 = 1 ; par2 = 1 ; par3 = 1 ; par4 = 0 ; par5 = 1��%2527%2522 ; par6 = 1 ;

R code (references can be found in the software module):

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

Free Statistics

Description of Statistical Computation

Tree of Dependent Computations

Dataset

Tables (Output of Computation)

Figures (Output of Computation)

Input Parameters & R Code