FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_regression_trees1.wasp
Title produced by softwareRecursive Partitioning (Regression Trees)
Date of computationFri, 09 Dec 2011 14:32:47 -0500
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2011/Dec/09/t1323459338kjshjzazoz7uq01.htm/, Retrieved Thu, 02 May 2024 17:26:16 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=153440, Retrieved Thu, 02 May 2024 17:26:16 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact111

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Recursive Partitioning (Regression Trees)] [] [2010-12-05 18:59:57] [b98453cac15ba1066b407e146608df68]
- R PD  [Recursive Partitioning (Regression Trees)] [Workshop 10 Recur...] [2011-12-09 18:45:00] [de8512d9b386046939a89973b76869e3]
- R  D    [Recursive Partitioning (Regression Trees)] [Workshop 10 Recur...] [2011-12-09 19:18:30] [de8512d9b386046939a89973b76869e3]
- R         [Recursive Partitioning (Regression Trees)] [Workshop 10 Recur...] [2011-12-09 19:22:17] [de8512d9b386046939a89973b76869e3]
- R  D          [Recursive Partitioning (Regression Trees)] [Workshop 10 Recur...] [2011-12-09 19:32:47] [5c44e6aad476a1bab98fc6774eca4c08] [Current]
- R               [Recursive Partitioning (Regression Trees)] [Workshop 10 Recur...] [2011-12-09 19:38:35] [de8512d9b386046939a89973b76869e3]
-  MP               [Recursive Partitioning (Regression Trees)] [Paper SHW Recursi...] [2011-12-16 14:43:06] [74be16979710d4c4e7c6647856088456]
- RMP               [Multiple Regression] [Paper SHW MLR 2] [2011-12-16 14:52:16] [de8512d9b386046939a89973b76869e3]
- R P               [Recursive Partitioning (Regression Trees)] [Paper SHW Recursi...] [2011-12-16 14:58:20] [de8512d9b386046939a89973b76869e3]
- RMP             [Recursive Partitioning (Regression Trees)] [Paper SHW Recursi...] [2011-12-16 14:42:20] [74be16979710d4c4e7c6647856088456]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
869	58	28	103	84786	98364	120982
2172	108	30	103	101193	96933	179321
901	49	22	51	38361	79234	123185
463	0	26	70	68504	42551	52746
371	1	18	22	22807	6853	33170
1495	86	44	148	71701	75851	173326
2187	104	40	124	80444	93163	258873
1491	63	34	70	53855	96037	180083
1036	82	23	66	99645	94728	135473
1882	115	36	134	114789	105499	202925
1220	50	25	84	65553	98958	153935
1289	83	39	156	97500	77900	132943
1812	105	33	110	77873	178812	221698
1731	114	43	158	90183	163253	260561
807	38	30	109	61542	27032	84853
1940	71	32	92	55813	86572	215641
1499	59	28	70	55461	85371	167542
2747	106	30	93	70106	120642	269651
2099	34	39	31	71570	78348	116408
918	20	26	66	33032	56968	78800
3373	115	39	133	139077	161632	277965
1713	85	33	113	71595	87850	150629
1438	76	28	100	72260	127969	168809
496	8	4	7	5950	15049	24188
744	21	18	61	32551	25109	65029
1161	30	14	41	31701	45824	101097
2694	92	28	102	120733	162647	233328
1769	75	28	99	73107	60622	206161
3148	128	38	129	132068	179566	311473
2474	105	23	62	149193	184301	235800
2084	55	36	73	46821	75661	177939
1954	56	32	114	87011	96144	207176
1389	72	25	70	55183	117286	174184
2269	75	36	116	73511	109377	187559
1268	118	23	74	78664	73631	119016
1943	77	40	138	70054	86767	182192
1762	66	40	151	74011	93487	194979
1857	116	33	115	93133	94552	275541
1502	73	34	104	225920	128754	182999
1441	99	30	108	62133	66363	135649
1416	53	22	69	43836	61724	120221
1317	30	26	99	38692	68580	145790
870	49	8	27	56622	55792	80953
2008	75	45	93	67267	129484	241066
1885	68	33	69	41140	87831	204713
1369	81	28	99	138599	136048	182613
2845	130	24	85	162901	186646	310839
1391	39	32	50	37510	64219	144966
602	13	19	64	43750	19630	43287
1743	74	20	31	40652	76825	155754
2014	109	31	92	85872	109427	201940
2143	151	32	106	89275	118168	235454
1401	54	31	114	120662	146304	224549
530	23	11	30	25162	24610	61857
387	4	0	0	855	6622	21054
1742	62	24	60	97068	115814	209641
449	18	8	9	14116	13155	31414
1606	64	40	140	110681	68847	184510
568	16	8	21	8773	13983	38214
1459	48	35	124	83209	65176	151101
1955	130	38	120	103487	180759	250579
1002	59	31	114	71220	100226	158015
956	32	28	78	56926	54454	85439
3604	95	40	141	115168	170745	351619
1035	14	30	101	111194	6940	84207
1701	70	32	90	51633	86839	165543
1249	19	27	36	75345	44830	141722
3352	91	31	97	98952	103300	299775
1369	135	41	148	123969	58106	104389
2201	87	32	105	135400	122422	199476
207	4	0	0	6023	7953	14688
151	7	0	0	1644	4245	7199
474	12	5	13	6179	21509	46660
141	0	1	4	3926	7670	17547
1318	46	24	46	73224	53608	152601

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time4 seconds
R Server'George Udny Yule' @ yule.wessa.net

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input & view raw input (R code)  \tabularnewline
Raw Output & view raw output of R engine  \tabularnewline
Computing time & 4 seconds \tabularnewline
R Server & 'George Udny Yule' @ yule.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=153440&T=0

[TABLE]
[ROW][C]Summary of computational transaction[/C][/ROW]
[ROW][C]Raw Input[/C][C]view raw input (R code) [/C][/ROW]
[ROW][C]Raw Output[/C][C]view raw output of R engine [/C][/ROW]
[ROW][C]Computing time[/C][C]4 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'George Udny Yule' @ yule.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=153440&T=0

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

As an alternative you can also use a QR Code:  
QR Code


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

Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time4 seconds
R Server'George Udny Yule' @ yule.wessa.net






Goodness of Fit
Correlation0.8044
R-squared0.6471
RMSE24679.7249

\begin{tabular}{lllllllll}
\hline
Goodness of Fit \tabularnewline
Correlation & 0.8044 \tabularnewline
R-squared & 0.6471 \tabularnewline
RMSE & 24679.7249 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=153440&T=1

[TABLE]
[ROW][C]Goodness of Fit[/C][/ROW]
[ROW][C]Correlation[/C][C]0.8044[/C][/ROW]
[ROW][C]R-squared[/C][C]0.6471[/C][/ROW]
[ROW][C]RMSE[/C][C]24679.7249[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=153440&T=1

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

As an alternative you can also use a QR Code:  
QR Code


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

Goodness of Fit
Correlation0.8044
R-squared0.6471
RMSE24679.7249






Actuals, Predictions, and Residuals
#ActualsForecastsResiduals
18478685163.35-377.350000000006
210119385163.3516029.65
3383612252415837
46850468840.8-336.800000000003
52280722524283
67170168840.82860.2
78044485163.35-4719.35000000001
85385585163.35-31308.35
99964585163.3514481.65
1011478985163.3529625.65
116555385163.35-19610.35
129750068840.828659.2
1377873126394.692307692-48521.6923076923
1490183126394.692307692-36211.6923076923
156154268840.8-7298.8
165581368840.8-13027.8
175546168840.8-13379.8
187010685163.35-15057.35
197157068840.82729.2
203303268840.8-35808.8
21139077126394.69230769212682.3076923077
227159568840.82754.2
237226085163.35-12903.35
24595022524-16574
25325512252410027
2631701225249177
27120733126394.692307692-5661.69230769231
287310768840.84266.2
29132068126394.6923076925673.30769230769
30149193126394.69230769222798.3076923077
314682168840.8-22019.8
328701185163.351847.64999999999
335518385163.35-29980.35
347351185163.35-11652.35
357866468840.89823.2
367005468840.81213.2
377401185163.35-11152.35
389313385163.357969.65
39225920126394.69230769299525.3076923077
406213368840.8-6707.8
41438362252421312
423869268840.8-30148.8
43566222252434098
4467267126394.692307692-59127.6923076923
454114068840.8-27700.8
46138599126394.69230769212204.3076923077
47162901126394.69230769236506.3076923077
483751068840.8-31330.8
49437502252421226
50406522252418128
518587285163.35708.649999999994
528927585163.354111.64999999999
53120662126394.692307692-5732.69230769231
5425162225242638
5585522524-21669
569706885163.3511904.65
571411622524-8408
5811068168840.841840.2
59877322524-13751
608320968840.814368.2
61103487126394.692307692-22907.6923076923
627122085163.35-13943.35
635692668840.8-11914.8
64115168126394.692307692-11226.6923076923
6511119468840.842353.2
665163368840.8-17207.8
677534568840.86504.2
689895285163.3513788.65
6912396968840.855128.2
7013540085163.3550236.65
71602322524-16501
72164422524-20880
73617922524-16345
74392622524-18598
757322468840.84383.2

\begin{tabular}{lllllllll}
\hline
Actuals, Predictions, and Residuals \tabularnewline
# & Actuals & Forecasts & Residuals \tabularnewline
1 & 84786 & 85163.35 & -377.350000000006 \tabularnewline
2 & 101193 & 85163.35 & 16029.65 \tabularnewline
3 & 38361 & 22524 & 15837 \tabularnewline
4 & 68504 & 68840.8 & -336.800000000003 \tabularnewline
5 & 22807 & 22524 & 283 \tabularnewline
6 & 71701 & 68840.8 & 2860.2 \tabularnewline
7 & 80444 & 85163.35 & -4719.35000000001 \tabularnewline
8 & 53855 & 85163.35 & -31308.35 \tabularnewline
9 & 99645 & 85163.35 & 14481.65 \tabularnewline
10 & 114789 & 85163.35 & 29625.65 \tabularnewline
11 & 65553 & 85163.35 & -19610.35 \tabularnewline
12 & 97500 & 68840.8 & 28659.2 \tabularnewline
13 & 77873 & 126394.692307692 & -48521.6923076923 \tabularnewline
14 & 90183 & 126394.692307692 & -36211.6923076923 \tabularnewline
15 & 61542 & 68840.8 & -7298.8 \tabularnewline
16 & 55813 & 68840.8 & -13027.8 \tabularnewline
17 & 55461 & 68840.8 & -13379.8 \tabularnewline
18 & 70106 & 85163.35 & -15057.35 \tabularnewline
19 & 71570 & 68840.8 & 2729.2 \tabularnewline
20 & 33032 & 68840.8 & -35808.8 \tabularnewline
21 & 139077 & 126394.692307692 & 12682.3076923077 \tabularnewline
22 & 71595 & 68840.8 & 2754.2 \tabularnewline
23 & 72260 & 85163.35 & -12903.35 \tabularnewline
24 & 5950 & 22524 & -16574 \tabularnewline
25 & 32551 & 22524 & 10027 \tabularnewline
26 & 31701 & 22524 & 9177 \tabularnewline
27 & 120733 & 126394.692307692 & -5661.69230769231 \tabularnewline
28 & 73107 & 68840.8 & 4266.2 \tabularnewline
29 & 132068 & 126394.692307692 & 5673.30769230769 \tabularnewline
30 & 149193 & 126394.692307692 & 22798.3076923077 \tabularnewline
31 & 46821 & 68840.8 & -22019.8 \tabularnewline
32 & 87011 & 85163.35 & 1847.64999999999 \tabularnewline
33 & 55183 & 85163.35 & -29980.35 \tabularnewline
34 & 73511 & 85163.35 & -11652.35 \tabularnewline
35 & 78664 & 68840.8 & 9823.2 \tabularnewline
36 & 70054 & 68840.8 & 1213.2 \tabularnewline
37 & 74011 & 85163.35 & -11152.35 \tabularnewline
38 & 93133 & 85163.35 & 7969.65 \tabularnewline
39 & 225920 & 126394.692307692 & 99525.3076923077 \tabularnewline
40 & 62133 & 68840.8 & -6707.8 \tabularnewline
41 & 43836 & 22524 & 21312 \tabularnewline
42 & 38692 & 68840.8 & -30148.8 \tabularnewline
43 & 56622 & 22524 & 34098 \tabularnewline
44 & 67267 & 126394.692307692 & -59127.6923076923 \tabularnewline
45 & 41140 & 68840.8 & -27700.8 \tabularnewline
46 & 138599 & 126394.692307692 & 12204.3076923077 \tabularnewline
47 & 162901 & 126394.692307692 & 36506.3076923077 \tabularnewline
48 & 37510 & 68840.8 & -31330.8 \tabularnewline
49 & 43750 & 22524 & 21226 \tabularnewline
50 & 40652 & 22524 & 18128 \tabularnewline
51 & 85872 & 85163.35 & 708.649999999994 \tabularnewline
52 & 89275 & 85163.35 & 4111.64999999999 \tabularnewline
53 & 120662 & 126394.692307692 & -5732.69230769231 \tabularnewline
54 & 25162 & 22524 & 2638 \tabularnewline
55 & 855 & 22524 & -21669 \tabularnewline
56 & 97068 & 85163.35 & 11904.65 \tabularnewline
57 & 14116 & 22524 & -8408 \tabularnewline
58 & 110681 & 68840.8 & 41840.2 \tabularnewline
59 & 8773 & 22524 & -13751 \tabularnewline
60 & 83209 & 68840.8 & 14368.2 \tabularnewline
61 & 103487 & 126394.692307692 & -22907.6923076923 \tabularnewline
62 & 71220 & 85163.35 & -13943.35 \tabularnewline
63 & 56926 & 68840.8 & -11914.8 \tabularnewline
64 & 115168 & 126394.692307692 & -11226.6923076923 \tabularnewline
65 & 111194 & 68840.8 & 42353.2 \tabularnewline
66 & 51633 & 68840.8 & -17207.8 \tabularnewline
67 & 75345 & 68840.8 & 6504.2 \tabularnewline
68 & 98952 & 85163.35 & 13788.65 \tabularnewline
69 & 123969 & 68840.8 & 55128.2 \tabularnewline
70 & 135400 & 85163.35 & 50236.65 \tabularnewline
71 & 6023 & 22524 & -16501 \tabularnewline
72 & 1644 & 22524 & -20880 \tabularnewline
73 & 6179 & 22524 & -16345 \tabularnewline
74 & 3926 & 22524 & -18598 \tabularnewline
75 & 73224 & 68840.8 & 4383.2 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=153440&T=2

[TABLE]
[ROW][C]Actuals, Predictions, and Residuals[/C][/ROW]
[ROW][C]#[/C][C]Actuals[/C][C]Forecasts[/C][C]Residuals[/C][/ROW]
[ROW][C]1[/C][C]84786[/C][C]85163.35[/C][C]-377.350000000006[/C][/ROW]
[ROW][C]2[/C][C]101193[/C][C]85163.35[/C][C]16029.65[/C][/ROW]
[ROW][C]3[/C][C]38361[/C][C]22524[/C][C]15837[/C][/ROW]
[ROW][C]4[/C][C]68504[/C][C]68840.8[/C][C]-336.800000000003[/C][/ROW]
[ROW][C]5[/C][C]22807[/C][C]22524[/C][C]283[/C][/ROW]
[ROW][C]6[/C][C]71701[/C][C]68840.8[/C][C]2860.2[/C][/ROW]
[ROW][C]7[/C][C]80444[/C][C]85163.35[/C][C]-4719.35000000001[/C][/ROW]
[ROW][C]8[/C][C]53855[/C][C]85163.35[/C][C]-31308.35[/C][/ROW]
[ROW][C]9[/C][C]99645[/C][C]85163.35[/C][C]14481.65[/C][/ROW]
[ROW][C]10[/C][C]114789[/C][C]85163.35[/C][C]29625.65[/C][/ROW]
[ROW][C]11[/C][C]65553[/C][C]85163.35[/C][C]-19610.35[/C][/ROW]
[ROW][C]12[/C][C]97500[/C][C]68840.8[/C][C]28659.2[/C][/ROW]
[ROW][C]13[/C][C]77873[/C][C]126394.692307692[/C][C]-48521.6923076923[/C][/ROW]
[ROW][C]14[/C][C]90183[/C][C]126394.692307692[/C][C]-36211.6923076923[/C][/ROW]
[ROW][C]15[/C][C]61542[/C][C]68840.8[/C][C]-7298.8[/C][/ROW]
[ROW][C]16[/C][C]55813[/C][C]68840.8[/C][C]-13027.8[/C][/ROW]
[ROW][C]17[/C][C]55461[/C][C]68840.8[/C][C]-13379.8[/C][/ROW]
[ROW][C]18[/C][C]70106[/C][C]85163.35[/C][C]-15057.35[/C][/ROW]
[ROW][C]19[/C][C]71570[/C][C]68840.8[/C][C]2729.2[/C][/ROW]
[ROW][C]20[/C][C]33032[/C][C]68840.8[/C][C]-35808.8[/C][/ROW]
[ROW][C]21[/C][C]139077[/C][C]126394.692307692[/C][C]12682.3076923077[/C][/ROW]
[ROW][C]22[/C][C]71595[/C][C]68840.8[/C][C]2754.2[/C][/ROW]
[ROW][C]23[/C][C]72260[/C][C]85163.35[/C][C]-12903.35[/C][/ROW]
[ROW][C]24[/C][C]5950[/C][C]22524[/C][C]-16574[/C][/ROW]
[ROW][C]25[/C][C]32551[/C][C]22524[/C][C]10027[/C][/ROW]
[ROW][C]26[/C][C]31701[/C][C]22524[/C][C]9177[/C][/ROW]
[ROW][C]27[/C][C]120733[/C][C]126394.692307692[/C][C]-5661.69230769231[/C][/ROW]
[ROW][C]28[/C][C]73107[/C][C]68840.8[/C][C]4266.2[/C][/ROW]
[ROW][C]29[/C][C]132068[/C][C]126394.692307692[/C][C]5673.30769230769[/C][/ROW]
[ROW][C]30[/C][C]149193[/C][C]126394.692307692[/C][C]22798.3076923077[/C][/ROW]
[ROW][C]31[/C][C]46821[/C][C]68840.8[/C][C]-22019.8[/C][/ROW]
[ROW][C]32[/C][C]87011[/C][C]85163.35[/C][C]1847.64999999999[/C][/ROW]
[ROW][C]33[/C][C]55183[/C][C]85163.35[/C][C]-29980.35[/C][/ROW]
[ROW][C]34[/C][C]73511[/C][C]85163.35[/C][C]-11652.35[/C][/ROW]
[ROW][C]35[/C][C]78664[/C][C]68840.8[/C][C]9823.2[/C][/ROW]
[ROW][C]36[/C][C]70054[/C][C]68840.8[/C][C]1213.2[/C][/ROW]
[ROW][C]37[/C][C]74011[/C][C]85163.35[/C][C]-11152.35[/C][/ROW]
[ROW][C]38[/C][C]93133[/C][C]85163.35[/C][C]7969.65[/C][/ROW]
[ROW][C]39[/C][C]225920[/C][C]126394.692307692[/C][C]99525.3076923077[/C][/ROW]
[ROW][C]40[/C][C]62133[/C][C]68840.8[/C][C]-6707.8[/C][/ROW]
[ROW][C]41[/C][C]43836[/C][C]22524[/C][C]21312[/C][/ROW]
[ROW][C]42[/C][C]38692[/C][C]68840.8[/C][C]-30148.8[/C][/ROW]
[ROW][C]43[/C][C]56622[/C][C]22524[/C][C]34098[/C][/ROW]
[ROW][C]44[/C][C]67267[/C][C]126394.692307692[/C][C]-59127.6923076923[/C][/ROW]
[ROW][C]45[/C][C]41140[/C][C]68840.8[/C][C]-27700.8[/C][/ROW]
[ROW][C]46[/C][C]138599[/C][C]126394.692307692[/C][C]12204.3076923077[/C][/ROW]
[ROW][C]47[/C][C]162901[/C][C]126394.692307692[/C][C]36506.3076923077[/C][/ROW]
[ROW][C]48[/C][C]37510[/C][C]68840.8[/C][C]-31330.8[/C][/ROW]
[ROW][C]49[/C][C]43750[/C][C]22524[/C][C]21226[/C][/ROW]
[ROW][C]50[/C][C]40652[/C][C]22524[/C][C]18128[/C][/ROW]
[ROW][C]51[/C][C]85872[/C][C]85163.35[/C][C]708.649999999994[/C][/ROW]
[ROW][C]52[/C][C]89275[/C][C]85163.35[/C][C]4111.64999999999[/C][/ROW]
[ROW][C]53[/C][C]120662[/C][C]126394.692307692[/C][C]-5732.69230769231[/C][/ROW]
[ROW][C]54[/C][C]25162[/C][C]22524[/C][C]2638[/C][/ROW]
[ROW][C]55[/C][C]855[/C][C]22524[/C][C]-21669[/C][/ROW]
[ROW][C]56[/C][C]97068[/C][C]85163.35[/C][C]11904.65[/C][/ROW]
[ROW][C]57[/C][C]14116[/C][C]22524[/C][C]-8408[/C][/ROW]
[ROW][C]58[/C][C]110681[/C][C]68840.8[/C][C]41840.2[/C][/ROW]
[ROW][C]59[/C][C]8773[/C][C]22524[/C][C]-13751[/C][/ROW]
[ROW][C]60[/C][C]83209[/C][C]68840.8[/C][C]14368.2[/C][/ROW]
[ROW][C]61[/C][C]103487[/C][C]126394.692307692[/C][C]-22907.6923076923[/C][/ROW]
[ROW][C]62[/C][C]71220[/C][C]85163.35[/C][C]-13943.35[/C][/ROW]
[ROW][C]63[/C][C]56926[/C][C]68840.8[/C][C]-11914.8[/C][/ROW]
[ROW][C]64[/C][C]115168[/C][C]126394.692307692[/C][C]-11226.6923076923[/C][/ROW]
[ROW][C]65[/C][C]111194[/C][C]68840.8[/C][C]42353.2[/C][/ROW]
[ROW][C]66[/C][C]51633[/C][C]68840.8[/C][C]-17207.8[/C][/ROW]
[ROW][C]67[/C][C]75345[/C][C]68840.8[/C][C]6504.2[/C][/ROW]
[ROW][C]68[/C][C]98952[/C][C]85163.35[/C][C]13788.65[/C][/ROW]
[ROW][C]69[/C][C]123969[/C][C]68840.8[/C][C]55128.2[/C][/ROW]
[ROW][C]70[/C][C]135400[/C][C]85163.35[/C][C]50236.65[/C][/ROW]
[ROW][C]71[/C][C]6023[/C][C]22524[/C][C]-16501[/C][/ROW]
[ROW][C]72[/C][C]1644[/C][C]22524[/C][C]-20880[/C][/ROW]
[ROW][C]73[/C][C]6179[/C][C]22524[/C][C]-16345[/C][/ROW]
[ROW][C]74[/C][C]3926[/C][C]22524[/C][C]-18598[/C][/ROW]
[ROW][C]75[/C][C]73224[/C][C]68840.8[/C][C]4383.2[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=153440&T=2

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

As an alternative you can also use a QR Code:  
QR Code


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

Actuals, Predictions, and Residuals
#ActualsForecastsResiduals
18478685163.35-377.350000000006
210119385163.3516029.65
3383612252415837
46850468840.8-336.800000000003
52280722524283
67170168840.82860.2
78044485163.35-4719.35000000001
85385585163.35-31308.35
99964585163.3514481.65
1011478985163.3529625.65
116555385163.35-19610.35
129750068840.828659.2
1377873126394.692307692-48521.6923076923
1490183126394.692307692-36211.6923076923
156154268840.8-7298.8
165581368840.8-13027.8
175546168840.8-13379.8
187010685163.35-15057.35
197157068840.82729.2
203303268840.8-35808.8
21139077126394.69230769212682.3076923077
227159568840.82754.2
237226085163.35-12903.35
24595022524-16574
25325512252410027
2631701225249177
27120733126394.692307692-5661.69230769231
287310768840.84266.2
29132068126394.6923076925673.30769230769
30149193126394.69230769222798.3076923077
314682168840.8-22019.8
328701185163.351847.64999999999
335518385163.35-29980.35
347351185163.35-11652.35
357866468840.89823.2
367005468840.81213.2
377401185163.35-11152.35
389313385163.357969.65
39225920126394.69230769299525.3076923077
406213368840.8-6707.8
41438362252421312
423869268840.8-30148.8
43566222252434098
4467267126394.692307692-59127.6923076923
454114068840.8-27700.8
46138599126394.69230769212204.3076923077
47162901126394.69230769236506.3076923077
483751068840.8-31330.8
49437502252421226
50406522252418128
518587285163.35708.649999999994
528927585163.354111.64999999999
53120662126394.692307692-5732.69230769231
5425162225242638
5585522524-21669
569706885163.3511904.65
571411622524-8408
5811068168840.841840.2
59877322524-13751
608320968840.814368.2
61103487126394.692307692-22907.6923076923
627122085163.35-13943.35
635692668840.8-11914.8
64115168126394.692307692-11226.6923076923
6511119468840.842353.2
665163368840.8-17207.8
677534568840.86504.2
689895285163.3513788.65
6912396968840.855128.2
7013540085163.3550236.65
71602322524-16501
72164422524-20880
73617922524-16345
74392622524-18598
757322468840.84383.2


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
par1 = pearson ; par2 = equal ; par3 = 2 ; par4 = no ;
Parameters (R input):
par1 = 5 ; par2 = none ; par3 = 3 ; par4 = no ;
R code (references can be found in the software module):
library(party)
library(Hmisc)
par1 <- as.numeric(par1)
par3 <- as.numeric(par3)
x <- data.frame(t(y))
is.data.frame(x)
x <- x[!is.na(x[,par1]),]
k <- length(x[1,])
n <- length(x[,1])
colnames(x)[par1]
x[,par1]
if (par2 == 'kmeans') {
cl <- kmeans(x[,par1], par3)
print(cl)
clm <- matrix(cbind(cl$centers,1:par3),ncol=2)
clm <- clm[sort.list(clm[,1]),]
for (i in 1:par3) {
cl$cluster[cl$cluster==clm[i,2]] <- paste('C',i,sep='')
}
cl$cluster <- as.factor(cl$cluster)
print(cl$cluster)
x[,par1] <- cl$cluster
}
if (par2 == 'quantiles') {
x[,par1] <- cut2(x[,par1],g=par3)
}
if (par2 == 'hclust') {
hc <- hclust(dist(x[,par1])^2, 'cen')
print(hc)
memb <- cutree(hc, k = par3)
dum <- c(mean(x[memb==1,par1]))
for (i in 2:par3) {
dum <- c(dum, mean(x[memb==i,par1]))
}
hcm <- matrix(cbind(dum,1:par3),ncol=2)
hcm <- hcm[sort.list(hcm[,1]),]
for (i in 1:par3) {
memb[memb==hcm[i,2]] <- paste('C',i,sep='')
}
memb <- as.factor(memb)
print(memb)
x[,par1] <- memb
}
if (par2=='equal') {
ed <- cut(as.numeric(x[,par1]),par3,labels=paste('C',1:par3,sep=''))
x[,par1] <- as.factor(ed)
}
table(x[,par1])
colnames(x)
colnames(x)[par1]
x[,par1]
if (par2 == 'none') {
m <- ctree(as.formula(paste(colnames(x)[par1],' ~ .',sep='')),data = x)
}
load(file='createtable')
if (par2 != 'none') {
m <- ctree(as.formula(paste('as.factor(',colnames(x)[par1],') ~ .',sep='')),data = x)
if (par4=='yes') {
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'10-Fold Cross Validation',3+2*par3,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'',1,TRUE)
a<-table.element(a,'Prediction (training)',par3+1,TRUE)
a<-table.element(a,'Prediction (testing)',par3+1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Actual',1,TRUE)
for (jjj in 1:par3) a<-table.element(a,paste('C',jjj,sep=''),1,TRUE)
a<-table.element(a,'CV',1,TRUE)
for (jjj in 1:par3) a<-table.element(a,paste('C',jjj,sep=''),1,TRUE)
a<-table.element(a,'CV',1,TRUE)
a<-table.row.end(a)
for (i in 1:10) {
ind <- sample(2, nrow(x), replace=T, prob=c(0.9,0.1))
m.ct <- ctree(as.formula(paste('as.factor(',colnames(x)[par1],') ~ .',sep='')),data =x[ind==1,])
if (i==1) {
m.ct.i.pred <- predict(m.ct, newdata=x[ind==1,])
m.ct.i.actu <- x[ind==1,par1]
m.ct.x.pred <- predict(m.ct, newdata=x[ind==2,])
m.ct.x.actu <- x[ind==2,par1]
} else {
m.ct.i.pred <- c(m.ct.i.pred,predict(m.ct, newdata=x[ind==1,]))
m.ct.i.actu <- c(m.ct.i.actu,x[ind==1,par1])
m.ct.x.pred <- c(m.ct.x.pred,predict(m.ct, newdata=x[ind==2,]))
m.ct.x.actu <- c(m.ct.x.actu,x[ind==2,par1])
}
}
print(m.ct.i.tab <- table(m.ct.i.actu,m.ct.i.pred))
numer <- 0
for (i in 1:par3) {
print(m.ct.i.tab[i,i] / sum(m.ct.i.tab[i,]))
numer <- numer + m.ct.i.tab[i,i]
}
print(m.ct.i.cp <- numer / sum(m.ct.i.tab))
print(m.ct.x.tab <- table(m.ct.x.actu,m.ct.x.pred))
numer <- 0
for (i in 1:par3) {
print(m.ct.x.tab[i,i] / sum(m.ct.x.tab[i,]))
numer <- numer + m.ct.x.tab[i,i]
}
print(m.ct.x.cp <- numer / sum(m.ct.x.tab))
for (i in 1:par3) {
a<-table.row.start(a)
a<-table.element(a,paste('C',i,sep=''),1,TRUE)
for (jjj in 1:par3) a<-table.element(a,m.ct.i.tab[i,jjj])
a<-table.element(a,round(m.ct.i.tab[i,i]/sum(m.ct.i.tab[i,]),4))
for (jjj in 1:par3) a<-table.element(a,m.ct.x.tab[i,jjj])
a<-table.element(a,round(m.ct.x.tab[i,i]/sum(m.ct.x.tab[i,]),4))
a<-table.row.end(a)
}
a<-table.row.start(a)
a<-table.element(a,'Overall',1,TRUE)
for (jjj in 1:par3) a<-table.element(a,'-')
a<-table.element(a,round(m.ct.i.cp,4))
for (jjj in 1:par3) a<-table.element(a,'-')
a<-table.element(a,round(m.ct.x.cp,4))
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable3.tab')
}
}
m
bitmap(file='test1.png')
plot(m)
dev.off()
bitmap(file='test1a.png')
plot(x[,par1] ~ as.factor(where(m)),main='Response by Terminal Node',xlab='Terminal Node',ylab='Response')
dev.off()
if (par2 == 'none') {
forec <- predict(m)
result <- as.data.frame(cbind(x[,par1],forec,x[,par1]-forec))
colnames(result) <- c('Actuals','Forecasts','Residuals')
print(result)
}
if (par2 != 'none') {
print(cbind(as.factor(x[,par1]),predict(m)))
myt <- table(as.factor(x[,par1]),predict(m))
print(myt)
}
bitmap(file='test2.png')
if(par2=='none') {
op <- par(mfrow=c(2,2))
plot(density(result$Actuals),main='Kernel Density Plot of Actuals')
plot(density(result$Residuals),main='Kernel Density Plot of Residuals')
plot(result$Forecasts,result$Actuals,main='Actuals versus Predictions',xlab='Predictions',ylab='Actuals')
plot(density(result$Forecasts),main='Kernel Density Plot of Predictions')
par(op)
}
if(par2!='none') {
plot(myt,main='Confusion Matrix',xlab='Actual',ylab='Predicted')
}
dev.off()
if (par2 == 'none') {
detcoef <- cor(result$Forecasts,result$Actuals)
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Goodness of Fit',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Correlation',1,TRUE)
a<-table.element(a,round(detcoef,4))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'R-squared',1,TRUE)
a<-table.element(a,round(detcoef*detcoef,4))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'RMSE',1,TRUE)
a<-table.element(a,round(sqrt(mean((result$Residuals)^2)),4))
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,'Actuals, Predictions, and Residuals',4,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'#',header=TRUE)
a<-table.element(a,'Actuals',header=TRUE)
a<-table.element(a,'Forecasts',header=TRUE)
a<-table.element(a,'Residuals',header=TRUE)
a<-table.row.end(a)
for (i in 1:length(result$Actuals)) {
a<-table.row.start(a)
a<-table.element(a,i,header=TRUE)
a<-table.element(a,result$Actuals[i])
a<-table.element(a,result$Forecasts[i])
a<-table.element(a,result$Residuals[i])
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable.tab')
}
if (par2 != 'none') {
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Confusion Matrix (predicted in columns / actuals in rows)',par3+1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'',1,TRUE)
for (i in 1:par3) {
a<-table.element(a,paste('C',i,sep=''),1,TRUE)
}
a<-table.row.end(a)
for (i in 1:par3) {
a<-table.row.start(a)
a<-table.element(a,paste('C',i,sep=''),1,TRUE)
for (j in 1:par3) {
a<-table.element(a,myt[i,j])
}
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable2.tab')
}