BeginPackage["MyPackages`MyFunctions`"]; FlatNumericList[alist_List] := Module[{l, n, newlist, templist}, newlist = {}; templist = Flatten[alist]; n = 0; Do[n++; If[NumericQ[templist\[LeftDoubleBracket]n\[RightDoubleBracket]], newlist = Join[newlist, {templist\[LeftDoubleBracket] n\[RightDoubleBracket]}]], {Length[templist]}]; Return[newlist]] ComputeStatistics[alist_List] := Module[{numericlist}, Needs["Statistics`DescriptiveStatistics`"]; numericlist = FlatNumericList[alist]; Print["Mean: ", Mean[numericlist]]; Print["Standard Deviation: ", StandardDeviation[numericlist]]; Print["Median: ", Median[numericlist]]; Print["Variance: ", Variance[numericlist]]; If[Length[Flatten[alist]] > Length[numericlist], Print["Excluding non-numeric values in input."], Print["Including all input values."]]] EndPackage[]