Relationship between Brain Mass and Body Mass
/****************************************************************/
/* S A S S A M P L E L I B R A R Y */
/* */
/* NAME: ROBUST4 */
/* TITLE: Relationship between Brain Mass and Body Mass */
/* PRODUCT: IML */
/* DATA: */
/* */
/* SUPPORT: Rick Wicklin UPDATE: SEP 2013 */
/* REF: */
/* MISC: */
/* */
/****************************************************************/
%include sampsrc(robustmc); /* define graphing modules */
proc iml;
load module=_all_; /* load graphing modules */
/* Log(Body Mass) Log(Brain Mass) */
mass={ 0.1303338 0.9084851, 2.6674530 2.6263400,
1.5602650 2.0773680, 1.4418520 2.0606980,
0.0170333 0.7403627, 4.0681860 1.6989700,
3.4060290 3.6630410, 2.2720740 2.6222140,
2.7168380 2.8162410, 1.0000000 2.0606980,
0.5185139 1.4082400, 2.7234560 2.8325090,
2.3159700 2.6085260, 1.7923920 3.1205740,
3.8230830 3.7567880, 3.9731280 1.8450980,
0.8325089 2.2528530, 1.5440680 1.7481880,
-0.9208187 0.0000000, -1.6382720 -0.3979400,
0.3979400 1.0827850, 1.7442930 2.2430380,
2.0000000 2.1959000, 1.7173380 2.6434530,
4.9395190 2.1889280, -0.5528420 0.2787536,
-0.9136401 0.4771213, 2.2833010 2.2552720};
create BrainLog from mass[c={"Body" "Brain"}];
append from mass;
close BrainLog;
optn = j(5,1,.);
optn[1] = 1; /* print basic output */
optn[2] = 1; /* print covariance matrices */
optn[5]= -1; /* nrep: use all subsets */
ods exclude EigenRobust Distances DistrRes;
call mve(sc, xmve, dist, optn, mass);
ods select all;
r1 = {"Quantile", "Number of Subsets", "Number of Singular Subsets",
"Number of Nonzero Weights", "Min Objective Function",
"Min Distance", "Chi-Square Cutoff Value"};
RobustCenter = xmve[1,];
RobustCov = xmve[3:4,];
print sc[r=r1],
RobustCenter[c={"X1","X2"}],
RobustCov[r={"X1","X2"} c={"X1","X2"}];
MVEOutliers = loc(dist[3,]=0);
print MVEOutliers;
proc iml;
load module=_all_; /* load graphing modules */
use BrainLog;
read all var _num_ into mass;
close BrainLog;
optn = j(5,1,.); optn[5]= -1;
vnam = {"Log(Body Mass)","Log(Brain Mass)"};
titl = "Estimates of Location and Scale (MVE)";
call MVEScatter(mass, optn, 0.9, vnam, titl);
proc iml;
use BrainLog;
read all var _num_ into mass;
close BrainLog;
/* MCD: Use Random Subsets */
optn = j(5,1,.);
call mcd(sc, xmve, dist, optn, mass);
r1 = {"Quantile", "Number of Subsets", "Number of Singular Subsets",
"Number of Nonzero Weights", "Min Objective Function",
"Min Distance", "Chi-Square Cutoff Value"};
RobustCenter = xmve[1,];
RobustCov = xmve[3:4,];
print sc[r=r1],
RobustCenter[c={"X1","X2"}],
RobustCov[r={"X1","X2"} c={"X1","X2"}];
proc iml;
load module=_all_; /* load graphing modules */
use BrainLog;
read all var _num_ into mass;
close BrainLog;
optn = j(5,1,.); optn[5]= -1;
vnam = { "Log(Body Mass)","Log(Brain Mass)" };
titl = "Estimates of Location and Scale (MCD)";
call MCDScatter(mass, optn, 0.9, vnam, titl);