Двунаправленная таблица распределения частот для количественных данных в R - PullRequest
Новая
       41

Двунаправленная таблица распределения частот для количественных данных в R

0 голосов
/ 22 мая 2018

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

dbh(cm) tht(m)
3   53.35
19  13.37
27  16.53
22  17.8
9   8.33
10  8.76
24  15.62
44  30.3
17  14.91
10  8.93

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

для диаметра: 

8 - 11.9
12 - 15.9
16 - 19.9 and so.

для высоты: 

3 - 4.9
5 - 6.9
7 - 8.9 and so.

Итак, у меня есть тысячи строк данных и я суммирую каждую частоту всама по себе полная боль.Я написал следующие строки (dbh представляет диаметр, а tht - высоту): 

> data <- read.csv('data.csv')
> diameter <- data$dbh
> range(diameter)
[1]  6.0 60.5
> breaks <- seq(6, 61, by=4)
> diameter.cut <- cut(diameter, breaks, right = FALSE)
> diameter.frq <- table(diameter.cut)
> cbind(diameter.frq)
        diameter.frq
[6,10)            35
[10,14)           77
[14,18)           59
[18,22)           25
[22,26)           25
[26,30)           51
[30,34)           38
[34,38)           28
[38,42)           21
[42,46)           22
[46,50)           14
[50,54)            2
[54,58)            6
>

Я могу сделать то же самое с высотой дерева (tht).Но проблема в том, как мне создать «таблицу перекрестных частот»

Например: в классе диаметров [6,10) есть 35 деревьев.Но мне нужно распределить эти 35 деревьев в каждом классе высоты.Допустим, 12 из этих деревьев относятся к классу [3,5) высоты, 8 из них в [5,7) ... и т. Д.

PS: я довольно новичок в R. Так что моя проблема может показатьсядовольно глупо, но я действительно копался в интернете, прежде чем отправлять здесь.Я прошу прощения за это.

1 Ответ

0 голосов
/ 22 мая 2018

Вы имели в виду это?После настройки вашего кода я получил ... 

library(dplyr)

df %>%
  mutate(diameter.cm = cut(diameter.cm, seq(range(diameter.cm)[1], range(diameter.cm)[2]+4, by=4), right = F),
         height.m = cut(height.m, seq(range(height.m)[1], range(height.m)[2]+2, by=2), right = F)) %>%
  group_by(diameter.cm, height.m) %>%
  tally()

Вывод: 

  diameter.cm height.m        n
1 [3,7)       [52.3,54.3)     1
2 [7,11)      [8.33,10.3)     3
3 [15,19)     [14.3,16.3)     1
4 [19,23)     [12.3,14.3)     1
5 [19,23)     [16.3,18.3)     1
6 [23,27)     [14.3,16.3)     1
7 [27,31)     [16.3,18.3)     1
8 [43,47)     [28.3,30.3)     1

Пример данных: 

df <- structure(list(diameter.cm = c(3L, 19L, 27L, 22L, 9L, 10L, 24L, 
44L, 17L, 10L), height.m = c(53.35, 13.37, 16.53, 17.8, 8.33, 
8.76, 15.62, 30.3, 14.91, 8.93)), class = "data.frame", row.names = c(NA, 
-10L))
...