Загрузка магистральных данных смертельного случая

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 и нажмите кнопку Section Break. Код в разделе может быть запущен независимо, который дает возможность исследовать идеи, когда вы написали свою программу.

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

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, кликните по Кнопке управления и выберите из доступных параметров.

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

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 = макс. (смертельные случаи. (dataToPlot));
p = polyfit (смертельные случаи. (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')