Чтобы создать фрейм данных:
num <- sample(1:25, 20)
x <- data.frame("Day_eclosion" = num, "Developmental" = c("AP", "MA",
"JU", "L"), "Replicate" = 1:5)
model <- glmer(Day_eclosion ~ Developmental + (1 | Replicate), family =
"poisson", data= x)
Я получаю этот возврат от:
a <- lsmeans(model, pairwise~Developmental, adjust = "tukey")
a$contrasts
contrast estimate SE df z.ratio p.value
AP - JU 0.2051 0.0168 Inf 12.172 <.0001
AP - L 0.3009 0.0212 Inf 14.164 <.0001
AP - MA 0.3889 0.0209 Inf 18.631 <.0001
JU - L 0.0958 0.0182 Inf 5.265 <.0001
JU - MA 0.1839 0.0177 Inf 10.387 <.0001
L - MA 0.0881 0.0222 Inf 3.964 0.0004
Я ищу простой способ превратить этот вывод (только значения p) в:
AP MA JU L
AP - <.0001 <.0001 <.0001
MA - - <.0001 0.0004
JU - - - <.0001
L - - -
У меня есть около 20 таких наборов, которые мне нужно превратить в таблицы, так что чем проще и универсальнее, тем лучше.
Бонусные баллы, если выходные данные разделены табуляцией и т. д., чтобы я мог легко вставить их в word/excel.
Спасибо!
pval = as.data.frame(a$contrasts)$p.value
. Это будут фактические значения p, а не отформатированные, как в выводе, поэтому, если вы этого хотите, потребуется дополнительное кодирование. - person Russ Lenth   schedule 02.04.2020