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_regression_trees1.wasp

Title produced by software

Recursive Partitioning (Regression Trees)

Date of computation

Wed, 14 Dec 2011 15:14:42 -0500

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/2011/Dec/14/t1323893696bjqbcuvk0mpsdjp.htm/, Retrieved Wed, 01 May 2024 20:24:37 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=155227, Retrieved Wed, 01 May 2024 20:24:37 +0000

QR Codes:

Paste this QR Code to cite your computation.

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords

Estimated Impact

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)

-     [Kendall tau Correlation Matrix] [] [2010-12-05 18:04:16] [b98453cac15ba1066b407e146608df68]
- RMPD    [Recursive Partitioning (Regression Trees)] [regression tree] [2011-12-14 20:14:42] [8432dc408001a08517818ba7ac24bdb0] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

1418	210907	56	396	81	3	79	30	115	94	112285	24188	146283	144	145	11
869	120982	56	297	55	4	58	28	109	103	84786	18273	98364	103	101	15
1530	176508	54	559	50	12	60	38	146	93	83123	14130	86146	98	98	19
3201	385534	92	1562	63	0	121	25	96	91	119182	33251	195663	150	144	23
1583	149061	44	656	66	5	43	26	100	93	116174	27101	95757	84	84	16
1439	165446	33	511	57	0	69	25	93	60	57635	16373	85584	80	79	21
1764	237213	84	655	74	0	78	38	140	123	66198	19716	143983	130	127	24
1373	133131	55	525	52	7	44	30	99	90	57793	9028	59238	60	60	15
4041	324799	154	1436	108	0	158	47	181	168	97668	29498	151511	140	133	17
1706	230964	53	612	43	4	102	30	116	115	133824	27563	136368	151	150	19
2152	236785	119	865	75	3	77	31	116	71	101481	18293	112642	91	91	19
1036	135473	41	385	32	0	82	23	88	66	99645	22530	94728	138	132	25
1929	215147	58	639	85	0	101	36	135	117	99052	35082	121527	124	124	19
2242	344297	75	963	86	1	80	30	108	108	67654	16116	127766	119	118	28
1220	153935	33	398	56	5	50	25	89	84	65553	15849	98958	73	70	24
2515	174724	92	966	135	0	123	34	129	120	69112	26569	85646	123	119	26
2147	174415	100	801	63	0	73	31	118	114	82753	24785	98579	90	89	15
2352	225548	112	892	81	5	81	31	118	94	85323	17569	130767	116	112	21
1638	223632	73	513	52	0	105	33	125	120	72654	23825	131741	113	108	26
1222	124817	40	469	44	0	47	25	95	81	30727	7869	53907	56	52	16
1677	210767	60	643	39	3	94	35	135	133	117478	37791	146761	119	116	16
1579	170266	62	535	73	4	44	42	154	122	74007	9605	82036	129	123	20
2452	294424	77	992	59	2	107	33	127	124	101494	34461	171975	175	162	24
2662	325107	99	937	64	0	84	36	136	126	79215	24919	159676	96	92	10
186	7176	17	70	1	0	0	0	0	0	1423	603	1929	0	0	19
865	106408	30	260	32	1	33	14	46	37	31081	12558	58391	41	41	25
1793	96560	76	503	129	0	42	17	54	38	22996	7784	31580	47	47	22
2527	265769	146	927	37	2	96	32	124	120	83122	28522	136815	126	120	15
1324	149112	56	537	65	6	56	35	128	95	60578	14459	69107	80	79	21
2702	175824	107	910	107	0	57	20	80	77	39992	14526	50495	70	65	22
1383	152871	58	532	74	5	59	28	97	90	79892	22240	108016	73	70	27
1179	111665	34	345	54	4	39	28	104	80	49810	11802	46341	57	55	26
4308	362301	119	1635	715	2	76	34	125	110	100708	11912	79336	68	67	26
1831	183167	66	557	66	0	91	39	149	138	82875	18220	93176	127	127	22
1438	168809	66	452	32	0	76	28	118	100	72260	21884	127969	102	99	20
496	24188	24	218	20	0	8	4	12	7	5950	2694	15049	7	7	22
2253	329267	259	764	71	8	79	39	144	140	115762	15808	155135	148	141	21
2352	218946	41	866	112	1	76	29	108	96	80670	25239	102996	112	109	22
2144	244052	68	574	66	5	101	44	166	164	143558	29801	160604	137	133	20
4691	341570	168	1276	190	1	94	21	80	78	117105	18450	158051	135	123	21
1112	103597	43	379	66	1	27	16	60	49	23789	7132	44547	26	26	20
1973	256462	105	798	56	0	123	35	127	124	105195	35940	174141	181	166	25
2474	235800	94	921	127	8	105	23	84	62	149193	46230	184301	190	179	18
1226	196553	57	503	50	2	41	29	111	99	95260	30546	129847	107	108	22
1389	174184	53	382	52	0	72	25	98	70	55183	19746	117286	94	90	25
1496	143246	103	464	42	5	67	27	105	104	106671	15977	71180	116	114	21
2269	187559	121	717	76	8	75	36	135	116	73511	22583	109377	106	103	20
1833	187681	62	690	67	2	114	28	107	91	92945	17274	85298	143	142	20
893	73566	32	385	39	6	22	23	88	67	22618	3007	23824	26	25	18
1403	167488	45	619	77	2	69	28	104	72	83737	21113	82981	113	113	8
1425	143756	46	479	57	0	105	34	132	120	69094	17401	73815	120	118	22
1840	243199	75	752	34	3	88	28	108	105	95536	23567	132190	134	129	26
1502	182999	88	430	39	6	73	34	129	104	225920	13065	128754	54	51	18
1420	152299	53	537	63	0	62	33	122	98	61370	14587	67808	78	76	20
2970	346485	90	1000	106	0	118	38	147	111	106117	24021	131722	121	118	24
1644	193339	78	465	47	2	100	35	87	71	84651	20537	106175	145	141	17
1654	122774	45	711	162	0	24	24	90	69	15986	4527	25157	27	27	20
1054	130585	46	299	57	5	67	29	109	107	95364	30495	76669	91	91	23
937	112611	41	248	36	0	46	20	78	73	26706	7117	57283	48	48	20
3004	286468	144	1162	263	1	57	29	111	107	89691	17719	105805	68	63	22
2547	148446	91	905	63	1	135	37	141	129	126846	33473	72413	150	144	20
1626	182079	63	512	63	2	124	33	124	118	102860	21115	96971	181	168	19
1468	140344	53	472	77	6	33	25	93	73	51715	7236	71299	65	64	15
2445	220516	62	905	79	1	98	32	124	119	55801	13790	77494	97	97	20
1964	243060	63	786	110	4	58	29	112	104	111813	32902	120336	121	117	22
1381	162765	32	489	56	2	68	28	108	107	120293	25131	93913	99	100	13
1659	232138	62	617	43	0	131	31	117	90	161647	35947	181248	188	187	20
2888	265318	117	925	111	10	110	52	199	197	115929	29848	146123	138	127	17
1290	85574	34	351	71	0	37	21	78	36	24266	6943	32036	40	37	14
2845	310839	92	1144	62	9	130	24	91	85	162901	42705	186646	254	245	22
1982	225060	93	669	56	7	93	41	158	139	109825	31808	102255	87	87	24
1904	232317	54	707	74	0	118	33	126	106	129838	26675	168237	178	177	22
1391	144966	144	458	60	0	39	32	122	50	37510	8435	64219	51	49	23
1559	164709	109	572	53	0	81	31	115	63	87771	36867	115338	176	177	17
2146	220801	75	720	105	1	51	18	72	63	44418	12607	84845	66	60	23
874	99466	50	273	32	0	28	23	91	69	192565	22609	153197	56	55	25
1590	92661	61	508	133	1	40	17	45	41	35232	5892	29877	39	39	16
1590	133328	55	506	79	0	56	20	78	56	40909	17014	63506	66	64	18
1210	61361	77	451	51	0	27	12	39	25	13294	5394	22445	27	26	20
1281	100750	72	407	67	0	83	30	119	93	140867	6440	68370	58	58	18
1105	102010	53	370	66	3	28	13	50	44	44332	5818	42071	26	26	24
1272	101523	42	316	76	0	59	22	88	87	61056	18647	50517	77	76	23
1944	243511	71	603	65	0	133	42	155	110	101338	20556	103950	130	129	24
391	22938	10	154	9	0	12	1	0	0	1168	238	5841	11	11	23
1605	152474	65	577	45	0	106	32	123	83	65567	22392	84396	101	101	23
1988	99923	66	617	115	0	44	25	99	80	32334	12237	35753	36	36	13
1386	132487	41	411	97	0	71	36	136	98	40735	8388	55515	120	89	20
2395	317394	86	975	53	1	116	31	117	82	91413	22120	209056	195	193	18
387	21054	16	146	2	0	4	0	0	0	855	338	6622	4	4	21
1742	209641	42	705	52	5	62	24	88	60	97068	11727	115814	89	84	17
620	22648	19	184	44	0	12	13	39	28	44339	3704	11609	24	23	20
449	31414	19	200	22	0	18	8	25	9	14116	3988	13155	39	39	19
800	46698	45	274	35	0	14	13	52	33	10288	3030	18274	14	14	18
1684	131698	65	502	74	0	60	19	75	59	65622	13520	72875	78	78	19
2699	244749	95	964	144	2	98	33	124	115	76643	20923	142775	106	101	22
1204	128423	64	369	89	8	32	38	145	120	92696	3769	20112	37	36	22
1138	97839	38	417	42	2	25	24	87	66	94785	12252	61023	77	75	15
2158	272458	65	822	99	0	100	43	162	152	93815	28864	132432	132	131	17
1111	172494	52	389	52	0	46	43	165	139	86687	21721	112494	144	131	19
1421	108043	62	466	29	1	45	14	54	38	34553	4821	45109	40	39	20
2833	328107	65	1255	125	3	129	41	159	144	105547	33644	170875	153	144	22
2922	351067	95	1024	95	3	136	45	170	160	213688	42935	214921	220	211	21
1002	158015	29	400	40	0	59	31	119	114	71220	18864	100226	79	78	19
2186	229242	247	719	128	4	63	31	120	119	91721	17939	78876	95	90	21
1035	84207	29	356	73	11	14	30	112	101	111194	325	6940	12	12	18
1417	120445	118	457	72	0	36	16	59	56	51009	13539	49025	63	57	16
3261	324598	110	1402	128	0	113	37	136	133	135777	34538	122037	134	133	20
1587	131069	67	600	61	4	47	30	107	83	51513	12198	53782	69	69	21
1424	204271	42	480	73	0	92	35	130	116	74163	26924	127748	119	119	15
946	116048	64	230	45	0	50	20	75	50	33416	10855	77395	63	61	20
1926	250047	81	651	58	0	41	18	71	61	83305	11932	89324	55	49	23
3352	299775	95	1367	97	9	91	31	120	97	98952	14300	103300	103	101	15
1641	195838	67	564	50	1	111	31	116	98	102372	25515	112283	197	196	18
2035	173260	63	716	37	3	41	21	79	78	37238	2805	10901	16	15	22
2312	254488	83	747	50	10	120	39	150	117	103772	29402	120691	140	136	16
961	92499	32	319	57	0	25	18	71	55	21399	5250	25899	21	21	17
1900	224330	83	612	52	1	131	39	144	132	130115	28608	139296	167	163	24
1254	135781	31	433	98	2	45	14	47	44	24874	8092	52678	32	29	13
1335	74408	67	434	61	4	29	7	28	21	34988	4473	23853	36	35	23
1597	81240	66	503	89	0	58	17	68	50	45549	1572	17306	13	13	5
1645	181633	70	564	48	2	47	30	110	73	64466	14817	89455	96	96	19
2429	271856	103	824	91	1	109	37	147	86	54990	16714	147866	151	151	24
872	95227	34	239	70	0	37	32	111	48	34777	1669	14336	23	23	19
1018	98146	40	459	37	0	15	17	68	48	27114	7768	30059	21	14	20
1314	59194	31	288	247	6	7	24	80	68	37636	6387	22097	20	20	22
1335	139942	42	498	46	0	54	22	88	87	65461	18715	96841	82	72	15
1403	118612	46	454	72	2	54	12	48	43	30080	7936	41907	90	87	19
910	72880	33	376	41	0	14	19	76	67	24094	8643	27080	25	21	25
616	65475	18	225	24	2	16	13	51	46	69008	7294	35885	60	56	21
771	71965	35	252	33	1	32	15	59	56	46090	7185	28313	85	82	19
1376	135131	66	481	87	0	38	15	60	60	34029	8509	36134	41	38	17
1232	108446	60	389	90	1	22	17	68	65	46300	13275	55764	26	25	15
1544	181528	54	609	69	0	32	16	61	60	40662	10737	66956	49	47	21
1230	134019	53	422	51	0	32	18	67	54	28987	8033	47487	46	45	24
1255	121848	39	339	45	0	37	17	64	52	30594	5401	35619	41	41	22
721	81872	45	245	25	0	32	16	64	61	27913	10856	45608	23	23	19
1109	58981	36	384	38	7	0	23	91	61	42744	2154	7721	14	14	20
740	53515	28	212	52	2	5	22	88	81	12934	6117	20634	16	16	21
728	56375	30	229	74	7	10	13	49	40	41385	4820	31931	21	21	19
689	65490	22	224	38	3	27	16	62	40	18653	5615	37754	32	27	22
995	76302	31	333	26	0	29	20	76	68	30976	8702	40557	35	33	14
1613	104011	55	384	67	6	25	22	88	79	63339	15340	94238	42	42	25
2048	98104	54	636	132	2	55	17	66	47	25568	8030	44197	68	68	11
301	30989	14	93	35	0	5	17	68	41	4154	1278	4103	6	6	16
1803	135458	81	581	118	3	43	12	48	29	19474	4236	44144	68	67	19
861	63123	43	304	43	1	34	17	68	60	39067	7196	27640	84	77	17
1451	74914	30	407	64	0	35	23	90	79	65892	6371	28990	30	30	20
628	31774	23	170	48	1	0	17	66	47	4143	1574	4694	0	0	22
1161	81437	38	312	64	0	37	14	54	40	28579	9620	42648	36	36	20
979	65745	53	340	75	0	26	21	77	42	38084	8645	25836	50	48	22
675	56653	45	168	39	0	38	18	68	49	27717	8987	22779	30	29	15
1241	158399	39	443	42	0	23	18	72	57	32928	5544	40820	30	28	23
1049	73624	24	367	93	0	30	17	64	40	19499	6909	32378	33	33	20
1081	91899	35	335	60	0	18	15	59	33	36874	6745	39613	37	33	17
1688	139526	151	364	71	0	28	21	84	77	48259	16724	60865	83	80	20
617	51567	30	206	27	2	21	14	56	45	28207	7025	20107	30	30	25
1656	102538	57	490	79	1	50	15	58	45	45833	9078	48231	51	51	22
705	86678	40	238	44	0	12	15	59	50	29156	4605	39725	19	18	16
1597	150580	77	530	124	0	27	22	83	71	45588	9653	62991	41	39	25
982	99611	35	291	81	0	41	21	81	67	45097	8914	49363	54	54	18
1212	99373	63	397	92	1	12	18	72	62	28394	6700	24552	25	24	19
1143	86230	44	467	42	0	21	17	61	54	18632	5788	31493	25	24	25
435	30837	19	178	10	0	8	4	15	4	2325	593	3439	8	8	23
532	31706	13	175	24	0	26	10	32	25	25139	4506	19555	26	26	24
882	89806	42	299	64	0	27	16	62	40	27975	6382	21228	20	19	21
830	64175	42	260	48	0	37	18	72	59	21792	6928	28893	46	47	21
652	59382	49	227	49	0	29	12	41	24	26263	1514	21425	47	47	22
707	119308	30	239	48	0	32	16	61	58	23686	9238	50276	37	37	21
954	76702	49	333	62	0	35	21	67	42	49303	8204	37643	51	51	18
285	19764	12	75	19	1	10	2	8	4	5752	2416	9927	10	10	13
733	84105	20	261	45	0	17	17	66	63	20055	5432	27184	34	34	22
642	64187	27	238	36	0	10	16	61	54	20154	5576	18475	12	11	23
894	72535	14	329	44	0	17	16	64	39	19540	6095	35873	27	21	15

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	6 seconds
R Server	'Gwilym Jenkins' @ jenkins.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 & 6 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ jenkins.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=155227&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]6 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Gwilym Jenkins' @ jenkins.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=155227&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=155227&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	6 seconds
R Server	'Gwilym Jenkins' @ jenkins.wessa.net

10-Fold Cross Validation
	Prediction (training)			Prediction (testing)
Actual	C1	C2	CV	C1	C2	CV
C1	694	85	0.8909	93	8	0.9208
C2	74	697	0.904	12	67	0.8481
Overall	-	-	0.8974	-	-	0.8889

\begin{tabular}{lllllllll}
\hline
10-Fold Cross Validation \tabularnewline
 & Prediction (training) & Prediction (testing) \tabularnewline
Actual & C1 & C2 & CV & C1 & C2 & CV \tabularnewline
C1 & 694 & 85 & 0.8909 & 93 & 8 & 0.9208 \tabularnewline
C2 & 74 & 697 & 0.904 & 12 & 67 & 0.8481 \tabularnewline
Overall & - & - & 0.8974 & - & - & 0.8889 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=155227&T=1

[TABLE]
[ROW][C]10-Fold Cross Validation[/C][/ROW]
[ROW][C][/C][C]Prediction (training)[/C][C]Prediction (testing)[/C][/ROW]
[ROW][C]Actual[/C][C]C1[/C][C]C2[/C][C]CV[/C][C]C1[/C][C]C2[/C][C]CV[/C][/ROW]
[ROW][C]C1[/C][C]694[/C][C]85[/C][C]0.8909[/C][C]93[/C][C]8[/C][C]0.9208[/C][/ROW]
[ROW][C]C2[/C][C]74[/C][C]697[/C][C]0.904[/C][C]12[/C][C]67[/C][C]0.8481[/C][/ROW]
[ROW][C]Overall[/C][C]-[/C][C]-[/C][C]0.8974[/C][C]-[/C][C]-[/C][C]0.8889[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=155227&T=1

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

As an alternative you can also use a QR Code:

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

10-Fold Cross Validation
	Prediction (training)			Prediction (testing)
Actual	C1	C2	CV	C1	C2	CV
C1	694	85	0.8909	93	8	0.9208
C2	74	697	0.904	12	67	0.8481
Overall	-	-	0.8974	-	-	0.8889

Confusion Matrix (predicted in columns / actuals in rows)
	C1	C2
C1	78	10
C2	7	78

\begin{tabular}{lllllllll}
\hline
Confusion Matrix (predicted in columns / actuals in rows) \tabularnewline
 & C1 & C2 \tabularnewline
C1 & 78 & 10 \tabularnewline
C2 & 7 & 78 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=155227&T=2

[TABLE]
[ROW][C]Confusion Matrix (predicted in columns / actuals in rows)[/C][/ROW]
[ROW][C][/C][C]C1[/C][C]C2[/C][/ROW]
[ROW][C]C1[/C][C]78[/C][C]10[/C][/ROW]
[ROW][C]C2[/C][C]7[/C][C]78[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=155227&T=2

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

As an alternative you can also use a QR Code:

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

Confusion Matrix (predicted in columns / actuals in rows)
	C1	C2
C1	78	10
C2	7	78

Figure 1

PNG link

Postscript link

PDF link

Figure 2

PNG link

Postscript link

PDF link

Figure 3

PNG link

Postscript link

PDF link

Parameters (Session):

par1 = kendall ;

Parameters (R input):

par1 = 10 ; par2 = quantiles ; par3 = 2 ; par4 = yes ;

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

Free Statistics

Description of Statistical Computation

Tree of Dependent Computations

Dataset

Tables (Output of Computation)

Figures (Output of Computation)

Input Parameters & R Code