Список функций Statistics Toolbox

В оглавление \ К следующему разделу \ К предыдущему разделу

Синтаксис

means = grpstats(X,group)
[means,sem,counts,name] = grpstats(X,group)
[means,sem,counts,name] = grpstats(x,group,alpha)

Описание

means = grpstats(X,group) функция предназначена для расчета среднего арифметического значения means категоризованной переменной Х. Среднее арифметическое значение рассчитывается для каждой категории. Деление на категории выполняется при помощи входного аргумента group. Выборка негруппированных значений Х может быть задана как вектор или матрица. Значения вектора или матрицы Х принадлежат к одной категории, если равны соответствующие значения переменной group. Входной аргумент group может быть представлен как вектор, массив строк или массив ячеек строковых переменных. Также group может быть массивом ячеек, содержащим несколько сгруппированных переменных, например {G1 G2 G3}. В последнем случае наблюдения принадлежат к одной группе, если равны между собой все сгруппированные переменные. Размерность векторов X и group должна совпадать. Если выборка Х задана как матрица, категоризующая переменная group должна быть представлена как вектор. Количество строк матрицы X и элементов вектора group должно быть одинаковым.

[means,sem,counts,name] = grpstats(x,group) функция позволяет рассчитать точечную means и интервальную sem оценки математического ожидания категоризованной переменной Х, количество элементов в каждой категории counts, и отобразить список названий категорий name.

Переменная name полезна для установления соответствия между полученными результатами расчета means, sem и заданными категориями, если элементы group не являются целыми числами.

[means,sem,counts,name] = grpstats(x,group,alpha) этот вариант синтаксиса функции кроме результатов расчета значений means, sem, counts и вывода названий name строит график средних арифметических значений и их 100*(1-alpha)% доверительных интервалов по категориям в порядке возрастания последних.

Примеры использования функции расчета точечной и интервальной оценок математического ожидания категоризованной переменной

Расчет точечной оценки математического ожидания для категоризованной выборки заданной как вектор. Выборка х из 100 наблюдений делится на 4 категории. Категории кодируются целыми числами в диапазоне от 1 до 4.

>> group = unidrnd(4,100,1);
>> x = normrnd(0,1,100,1);
>> means = grpstats(x,group)
means =
    0.0369
    0.0317
    0.0759
    0.0454

Расчет точечной оценки математического ожидания для категоризованной выборки заданной в виде матрицы. Переменная Х задается как матрица с размерностью 100x5 и представляет собой 5 выборок по 100 элементов в каждой. Математическое ожидание при генерации выборок изменяется от 1 до 5. Значения в каждой выборке делятся на 4 группы. Категоризующая переменная group задается как вектор целых чисел с размерностью 100x1.

>> group = unidrnd(4,100,1);
>> true_mean = 1:5;
>> true_mean = true_mean(ones(100,1),:);
>> x = normrnd(true_mean,1);
>> means = grpstats(x,group)
means =
    1.0369    1.5801    2.8850    4.3252    5.1890
    1.0317    1.9425    3.0184    3.8486    4.9424
    1.0759    1.9433    2.8528    4.0344    4.6057
    1.0454    1.9528    2.9190    3.8409    4.7507

Расчет точечных means и интервальных sem оценок математического ожидания категоризованной переменной Х, количества элементов в каждой категории counts, и отображения списка названий категорий name. Переменная Х задается как матрица с размерностью 100x5. Категории кодируются вектором целых чисел на 4 группы.

>> group = unidrnd(4,100,1);
>> x = normrnd(0,1,100,5);
>> [means,sem,counts,name] = grpstats(x,group)
means =
   -0.2741    0.0263   -0.1199   -0.1949    0.0658
    0.0205    0.1609   -0.0889    0.1441   -0.1134
   -0.2905    0.0505   -0.0386    0.1086   -0.2904
    0.5349   -0.0313    0.0819    0.5119   -0.3288
sem =
    0.2056    0.1608    0.1764    0.1989    0.1852
    0.1506    0.1960    0.2535    0.2000    0.2050
    0.1678    0.1209    0.1544    0.1434    0.1537
    0.2332    0.2603    0.2468    0.1457    0.1845
counts =
    29    29    29    29    29
    25    25    25    25    25
    29    29    29    29    29
    17    17    17    17    17
name = 
    '1'
    '2'
    '3'
    '4'
        
Расчет значений means, sem, counts, вывод name и графика средних арифметических значений с 95% доверительными
интервалами по соответствующим категориям.

>> group = unidrnd(4,100,1);
>> x = normrnd(0,1,100,5);
>> [means,sem,counts,name] = grpstats(x,group,0.05)
means =
    0.0135   -0.4496    0.0543    0.2349   -0.1285
    0.1529    0.2906   -0.1331    0.2133    0.2041
    0.1151    0.0340    0.0273    0.1469    0.2279
   -0.3031   -0.4670   -0.2492    0.0589    0.2070
sem =
    0.2215    0.2651    0.1741    0.2816    0.1873
    0.2205    0.1940    0.1795    0.2196    0.1597
    0.2330    0.1574    0.2300    0.2046    0.1756
    0.2227    0.2238    0.1720    0.2088    0.2229
counts =
    19    19    19    19    19
    29    29    29    29    29
    28    28    28    28    28
    24    24    24    24    24
name = 
    '1'
    '2'
    '3'
    '4'

Расчет значений means, sem, counts и вывод списка категорий name. Переменная Х задается как матрица с размерностью 20x5. Категории определяются строковым массивом и выборка делится на 4 группы.

>> group ={'A' 'B' 'C' 'D' 'A' 'B' 'A' 'B' 'C' 'D' 'C' 'D' 'A' 'A' 'B' 'C' 'D' 'B' 'C' 'D'}'
group = 
    'A'
    'B'
    'C'
    'D'
    'A'
    'B'
    'A'
    'B'
    'C'
    'D'
    'C'
    'D'
    'A'
    'A'
    'B'
    'C'
    'D'
    'B'
    'C'
    'D'
>> x = normrnd(0,1,20,5)
x =
   -0.2463   -0.2979   -0.5299    0.8683    0.7193
   -0.1457    1.1543    0.5411   -0.8048   -0.2831
   -1.1690    1.0461    0.6817   -0.7527   -1.4250
   -0.0220    2.1269    0.5386   -0.7458    0.4615
    0.6183   -0.6558   -0.5100   -0.3097    1.0915
    1.8659   -1.1424   -1.3221   -1.5219   -1.0443
    0.0819    0.9490   -0.6107    0.8265   -2.8428
    1.6080   -0.4046   -0.5653   -0.6130    0.9968
   -0.3807   -0.3843    0.0862    0.9597    0.0765
   -1.2996    0.4820    0.6915    1.9730   -1.8667
   -0.7240    0.4438    2.1338    0.2950   -0.6136
   -0.5650    0.3811   -0.0029   -0.3927    1.1694
    0.6217    1.1023   -0.0895    0.5759   -0.5750
   -1.3355    0.8564   -0.2550   -1.1414   -0.2648
   -0.1231   -1.1785   -0.8742    0.0611    0.0047
   -1.1028    0.4020    0.4229    0.0123   -0.0394
   -2.7532   -0.5842   -0.1334   -0.1681   -0.5054
    0.2520   -0.9795    0.5396   -0.6873   -1.1578
   -0.8581    0.1151    0.8752   -0.9907    0.7104
    1.1354    0.0685   -1.2508   -0.0498    0.7282
>> [means,sem,counts,name] = grpstats(x,group)
means =
    0.3502   -0.8565   -0.5624    0.2891   -0.0456
    0.2269   -0.1071   -0.5212   -0.7214    0.4791
   -1.2206    0.6396    0.3454   -0.0956   -0.2041
   -0.2389   -0.2455    0.0123    0.4581    0.4169
sem =
    0.6679    0.3892    0.6913    0.7837    0.3769
    0.2133    0.4039    0.5240    0.3131    0.3355
    0.5199    0.3487    0.6232    0.2142    0.2502
    0.4466    0.2645    0.4143    0.4821    0.3391
counts =
     5     5     5     5     5
     5     5     5     5     5
     5     5     5     5     5
     5     5     5     5     5
name = 
    'A'
    'B'
    'C'
    'D'

В оглавление \ К следующему разделу \ К предыдущему разделу