#r #loops #ggplot2 #formatting
Вопрос:
Я провел анализ видов индикаторов в трех наборах данных (multipatt {indicspecies}) и хочу распечатать их сводки в виде таблицы.
Я хотел бы получить совет о том, как объединить или связать эти сводки, и как использовать ggplot (или аналогичный) для настройки вывода (изменение меток, полужирное/курсивное/форматирование границ строк и т. Д.). Я понимаю, что приведенный ниже пример может быть полезен, если его написать в виде цикла — я еще не освоил это, поэтому я бы также оценил помощь в этом отделе.
Спасибо всем!
Пример типа выходных данных, которые у меня есть:
library(indicspecies) # Data for Year 2000 site lt;- c(1:10) Year lt;- 2000 sp1 lt;- c(0,1,4,4,3,3,2,1,0,2) sp2 lt;- c(4,3,0,0,2,2,3,2,1,3) sp3 lt;- c(3,2,1,1,2,2,3,2,1,3) sp4 lt;- c(2,4,1,0,1,2,3,4,3,1) df_1 lt;- data.frame(site, Year, sp1, sp2, sp3, sp4) clusters_1 lt;- df_1 %gt;% mutate(Association = case_when(site %in% c(3, 4)~"A", site %in% c(2, 8, 9)~"B", site %in% c(1, 5, 6, 7, 10)~"C")) abund_1 = clusters_1[3:6] assoc_1 = clusters_1$Association indic_1 = multipatt(abund_1, assoc_1, func = "r.g", control = how(nperm=99)) summary(indic_1) #~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ # Data for Year 2005 site lt;- c(1:10) Year lt;- 2005 sp1 lt;- c(0,1,4,4,3,3,2,1,0,2) sp3 lt;- c(4,3,0,0,2,2,3,2,1,3) sp5 lt;- c(3,2,1,1,2,2,3,2,1,3) sp7 lt;- c(2,4,1,0,1,2,3,4,3,1) df_2 lt;- data.frame(site, Year, sp1, sp3, sp5, sp7) clusters_2 lt;- df_2 %gt;% mutate(Association = case_when(site %in% c(3, 4)~"A", site %in% c(2, 8, 9)~"B", site %in% c(1, 5, 6, 7, 10)~"C")) abund_2 = clusters_2[3:6] assoc_2 = clusters_2$Association indic_2 = multipatt(abund_2, assoc_2, func = "r.g", control = how(nperm=99)) summary(indic_2) #~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ # Data for Year 2010 site lt;- c(1:10) Year lt;- 2010 sp2 lt;- c(0,1,4,4,3,3,2,1,0,2) sp4 lt;- c(4,3,0,0,2,2,3,2,1,3) sp6 lt;- c(3,2,1,1,2,2,3,2,1,3) sp8 lt;- c(2,4,1,0,1,2,3,4,3,1) df_3 lt;- data.frame(site, Year, sp2, sp4, sp6, sp8) clusters_3 lt;- df_3 %gt;% mutate(Association = case_when(site %in% c(3, 4)~"A", site %in% c(2, 8, 9)~"B", site %in% c(1, 5, 6, 7, 10)~"C")) abund_3 = clusters_3[3:6] assoc_3 = clusters_3$Association indic_3 = multipatt(abund_3, assoc_3, func = "r.g", control = how(nperm=99)) summary(indic_3) Я действительно ломаюсь, пытаясь создать пример кода нужной мне таблицы, поэтому я создал желаемый вывод в Excel и включил изображение желаемой выходной таблицы. Примечание: «sp» — это аббревиатура для «видов», поэтому «sp1» = «Вид 1». 