У меня есть список широт и долгот, которые представляют средний центр аварий и местоположений больниц, рассчитанный с использованием гравитационного метода, то есть для различных аварий у меня есть 1 средняя широта и долгота для года и 1 средняя широта и долгота для больницы. Теперь я хочу создать диаграмму роз для несчастного случая и больницы, рассчитать их корреляцию, значимость и направление движения по годам, а также среднее направление и другую сводную статистику с использованием круговой статистики. Как мне поступить? Некоторые идеи в коде R и R были бы полезны. Я просмотрел пакет циркуляр , но не понял, как преобразовать координаты x, y в один столбец для создания диаграммы роз и другого анализа. Любая помощь приветствуется.
Мои данные таковы: LA - это долгота несчастных случаев, LTA - это широта несчастных случаев, LH - это долгота больниц, а LTH - это долгота больниц.
df <- read.table(text='Year LA LTA LH LTH
Year1 -84.3213 33.8488 -84.3281 33.8779
Year2 -84.322 33.8470 -84.3284 33.8782
year3 -84.323 33.84461 -84.3293 33.8791
year4 -84.3165 33.8359 -84.3452 33.8404
year5 -84.3257 33.8330 -84.3413 33.8340', header=TRUE)
Это единственные данные, которые я имею на данный момент, я конвертировал эту широту и долготу в азимут, который дал мне угол от источника, я верю. Затем я использовал этот угол, чтобы нарисовать диаграмму розы. Я не уверен, должен ли я преобразовать эти отрицательные координаты в положительный угол. Любые вклады от статистического лица будут оценены. Так как у меня есть два угла от несчастных случаев и больниц, я считаю, что могу сопоставить и сделать тест Уотсона. Ниже приведен мой код, который до сих пор работает над ним: диаграмма Роуз предназначена только для несчастных случаев, а не для больниц, хорошо ли рисовать их вместе? что еще нужно сделать, пожалуйста, предложения руды
library(geosphere)
library(circular)
library(CircStats)
data <- read.table("trial_circ.csv", header = T, sep = ",")
dfr<-data.frame(data[3:4])
dist<- trackAzimuth(as.matrix(dfr))
d<-rose.diag(dist)
Заранее спасибо