Вычислите средние высоты по полу для групп пациентов и отобразите результаты.

Загрузите высоты пациентов и полы из файла данных patients.mat.

load patients
whos Gender Height

  Name          Size            Bytes  Class     Attributes

  Gender      100x1             11412  cell                
  Height      100x1               800  double              

Задайте группы по полу с findgroups.

G = findgroups(Gender);

Разделение Height в группы, заданные G. Вычислите среднюю высоту по полу. Первая строка выходного аргумента является средним ростом пациентов женского пола, а вторая строка - средним ростом пациентов мужского пола.

splitapply(@mean,Height,G)

ans = 2×1

   65.1509
   69.2340

Вычислите отклонения различий в показаниях артериального давления для групп пациентов и отобразите результаты. Показания артериального давления содержатся в двух переменных данных. Чтобы вычислить различия, используйте функцию, которая принимает два входных параметров.

Загрузите показания артериального давления и данные о курении для 100 пациентов из файла данных patients.mat.

load patients
whos Systolic Diastolic Smoker

  Name             Size            Bytes  Class      Attributes

  Diastolic      100x1               800  double               
  Smoker         100x1               100  logical              
  Systolic       100x1               800  double               

Определите func как функция, которая вычисляет отклонения различий между стадолическими и диастолическими показаниями артериального давления для курильщиков и некурящих. func требуется два входных параметров.

func = @(x,y) var(x-y);

Использование findgroups и splitapply чтобы разделить данные о пациентах на группы и вычислить отклонения различий. findgroups также возвращает идентификаторы групп в smokers. The splitapply вызовы функций func один раз в группу, с Systolic и Diastolic как два входных параметров.

[G,smokers] = findgroups(Smoker);
varBP = splitapply(func,Systolic,Diastolic,G)

varBP = 2×1

   44.4459
   48.6783

Составьте таблицу, которая содержит отклонения различий с количеством пациентов в каждой группе.

numPatients = splitapply(@numel,Smoker,G);
T = table(smokers,numPatients,varBP)

T=2×3 table
    smokers    numPatients    varBP 
    _______    ___________    ______

     false         66         44.446
     true          34         48.678

Вычислите минимальный, медианный и максимальный веса для групп пациентов и верните эти результаты в виде массивов для каждой группы. splitapply объединяет выходные аргументы так, чтобы можно было отличить выход для каждой группы от выхода для других групп.

Задайте функцию, которая возвращает минимум, медиану и максимум как вектор-строка.

mystats = @(x)[min(x) median(x) max(x)];

Загрузите веса пациентов, полы и статус курильщиков из patients.mat.

load patients
whos Weight Gender Smoker

  Name          Size            Bytes  Class      Attributes

  Gender      100x1             11412  cell                 
  Smoker      100x1               100  logical              
  Weight      100x1               800  double               

Использование findgroups и splitapply чтобы разделить веса пациентов на группы и вычислить статистику для каждой группы.

G = findgroups(Gender,Smoker);
Y = splitapply(mystats,Weight,G)

Y = 4×3

  111.0000  131.0000  147.0000
  115.0000  131.0000  146.0000
  158.0000  181.5000  194.0000
  164.0000  181.0000  202.0000

В этом примере можно вернуть нескаляр выхода как векторы-строки, потому что данные и сгруппированные переменные векторов-столбцов. Каждая строка Y содержит статистику для другой группы пациентов.

Вычислите средний индекс массы тела (ИМТ) из таблиц данных о пациентах. Сгруппировать пациентов по полу и статусу курильщиков или некурящих.

Загрузите данные о пациентах и сгруппируйте переменные в таблицы.

load patients
DT = table(Height,Weight);
GT = table(Gender,Smoker);

Определите функцию, которая вычисляет средний ИМТ из весов и высот групп или пациентов.

meanBMIFcn = @(h,w)mean((w ./ (h.^2)) * 703);

Составьте таблицу, которая содержит средний ИМТ для каждой группы.

[G,results] = findgroups(GT);
meanBMI = splitapply(meanBMIFcn,DT,G);
results.meanBMI = meanBMI

results=4×3 table
      Gender      Smoker    meanBMI
    __________    ______    _______

    {'Female'}    false     21.672 
    {'Female'}    true      21.669 
    {'Male'  }    false     26.578 
    {'Male'  }    true      26.458 

Рассчитать минимальные, средние и максимальные высоты для групп пациентов и вернуть результаты в таблицу.

% Copyright 2015 The MathWorks, Inc.

function [lo,avg,hi] = multiStats(x)
lo = min(x);
avg = mean(x);
hi = max(x);
end

load patients
T = table(Gender,Height);
summary(T)

Variables:

    Gender: 100x1 cell array of character vectors

    Height: 100x1 double

        Values:

            Min          60   
            Median       67   
            Max          72   

[G,gender] = findgroups(T.Gender);
[minHeight,meanHeight,maxHeight] = splitapply(@multiStats,T.Height,G);
result = table(gender,minHeight,meanHeight,maxHeight)

result =

  2x4 table

      gender      minHeight    meanHeight    maxHeight
    __________    _________    __________    _________

    {'Female'}       60          65.151         70    
    {'Male'  }       66          69.234         72