Видовые богачи с богатой проблемой упаковки с АН

Question

Мне нужно рассчитать видовое богатство для кадра данных, называемого riquesa1, который состоит из 13 переменных (видов) и 59 наблюдений (частота каждого вида в общей выборке). Я использую функцию rich из пакета rich с этим кодом:

rich(matrix = riquesa1, verbose = FALSE, nrandom = NULL)

Однако, когда я запускаю функцию, я получаю эту ошибку:

Error in if (S[is] == 0) next : missing value where TRUE/FALSE needed

Я предполагаю, что это может быть из-за присутствия Na в кадре данных, но я не уверен. Как я могу решить это?

Заранее спасибо !!

Answer 1

reconstruct <- structure(list(`Amphora lybica` = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), `Amphora spp.` = c(0.405, NA, 0.127, NA, NA), `Cocconeis placentula` = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), `Cyclotella ocellata` = c(NA, NA, NA, NA, 0.0606), `C. distinguenda complex` = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), `Diploneis ovalis` = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), `Encyonopsis subminuta` = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), `Encyonopsis caespitosa` = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), `Hantzschia amphioxus` = c(NA, 0.187, 0.08, NA, NA), `Luticula mutica` = c(NA, 0.283, 0.163, NA, NA), `Mastogloia smithii` = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), `Mastogloia spp.` = c(0.092, NA, NA, NA, NA), `Pseudostaurosira brevistriata` = c(0.428, 0.063, 0.331, 0.964, 0.798)), row.names = c(NA, -5L), class = c("tbl_df", "tbl", "data.frame"))