Загрузка данных о смертности на шоссе

Live Editor отображает выход вместе с кодом, который произвел его. Чтобы запустить раздел, перейдите на вкладку Live Editor и нажмите кнопку Run Section. Можно также кликнуть синюю полосу, которая появляется при перемещении мыши к левому краю сечения.

В этом примере мы исследуем некоторые данные о смертности на шоссе. Начните, загрузив данные. Переменные показаны в виде заголовков столбцов таблицы.

load fatalities
fatalities(1:10,:)

ans=10×8 table
                            longitude    latitude    deaths    drivers    vehicles    vehicleMiles    alcoholRelated    urbanPopulation
                            _________    ________    ______    _______    ________    ____________    ______________    _______________

    Wyoming                  -107.56      43.033      164      380.18      671.53         9261              54              65.226     
    District_of_Columbia     -77.027      38.892       43      349.12       240.4         3742              12                 100     
    Vermont                  -72.556      44.043       98      550.46      551.52         7855              20              38.196     
    North_Dakota               -99.5      47.469      100      461.78      721.84         7594              35              55.807     
    South_Dakota             -99.679      44.272      197       563.3      882.77         8784              76              51.923     
    Delaware                 -75.494      39.107      134      533.94      728.52         9301              48              80.021     
    Montana                  -110.58      46.867      229      712.88      1056.7        11207             100              54.031     
    Rhode_Island             -71.434      41.589       83      741.84       834.5         8473              41              90.936     
    New_Hampshire            -71.559      43.908      171      985.77      1244.6        13216              51              59.181     
    Maine                    -69.081      44.886      194      984.83      1106.8        14948              58              40.206     

Вычислите показатели смертности

Live Editor позволяет вам разделить программу на разделы, содержащие текст, код и вывод. Чтобы создать новый раздел, перейдите на вкладку Live Editor и нажмите кнопку «Пропуск раздела». Код в разделе может запускаться независимо, что облегчает изучение идей, когда вы пишете свою программу.

Рассчитать смертность на миллион транспортных миль. Из этих значений мы можем найти состояния с самым низким и самым высоким уровнем смертности.

states = fatalities.Properties.RowNames;
rate = fatalities.deaths./fatalities.vehicleMiles;
[~, minIdx] = min(rate);                  % Minimum accident rate
[~, maxIdx] = max(rate);                  % Maximum accident rate
disp([states{minIdx} ' has the lowest fatality rate at ' num2str(rate(minIdx))])

Massachusetts has the lowest fatality rate at 0.0086907

disp([states{maxIdx} ' has the highest fatality rate at ' num2str(rate(maxIdx))])

Mississippi has the highest fatality rate at 0.022825

Распределение смертельных случаев

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

Мы можем использовать столбчатую диаграмму, чтобы увидеть распределение смертности среди состояний. Существует 11 состояния, смертность которых превышает 0,02 на миллион транспортных миль.

histogram(rate,10)
xlabel('Fatalities per Million Vehicle Miles')
ylabel('Number of States')

Поиск корреляций в данных

Можно быстро исследовать данные в Live Editor, экспериментировав со значениями параметров, чтобы увидеть, как изменятся результаты. Добавьте элементы управления, чтобы изменить значения параметров в интерактивном режиме. Чтобы добавить элементы управления, перейдите на вкладку Live Editor, нажмите кнопку Control и выберите из доступных опций.

Мы можем экспериментировать с данными, чтобы увидеть, коррелирует ли какая-либо из переменных в таблице с смертельными исходами на шоссе. Например, похоже, что уровень смертности на шоссе ниже в штатах с более высоким процентом городского населения.

dataToPlot = "urbanPopulation";
закрыть                                      % Close any open figures
рассеяние (смертельные случаи. (dataToPlot),          скорость)% Plot fatalities vs. selected variable
xlabel (dataToPlot)
ylabel ('Percent Fatalities per Million Vehicle Miles') 

держаться on
xmin = min (смертельные случаи. (dataToPlot));
xmax = max (смертельные случаи. (dataToPlot));
p = полифит (смертельные случаи. (dataToPlot), скорость,   1);% Calculate & plot least squares line
график ([xmin xmax], polyval (p, [xmin xmax]))

Постройте графики смертей и урбанизации на карте США

Результирующие результаты и поделитесь live скриптом с коллегами. Используя ваш live скрипт, они могут воссоздать и расширить ваш анализ. Можно также сохранить анализ как HTML, Microsoft ® Word или PDF для публикации.

Основываясь на этом анализе, мы можем обобщить наши результаты, используя график уровня смертности и городского населения на карте континентальных Соединенных Штатов.

load usastates.mat
figure
geoplot([usastates.Lat], [usastates.Lon], 'black')
geobasemap darkwater
hold on
geoscatter(fatalities.latitude,fatalities.longitude,2000*rate,fatalities.urbanPopulation,'filled')
c = colorbar;
title(c,'Percent Urban')