Emmeans 连续自变量

2024-04-29

我想解释一下Type_f with Type_space实验的内容和速率Exhaustion_product和定量变量Age.

这是我的数据:

res=structure(list(Type_space = structure(c(2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 
    3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 
    3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 
    3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 
    3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 
    3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 
    3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 
    3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 
    3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 
    3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 
    3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 
    3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 
    4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 
    4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 
    4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 
    4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 
    4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 
    4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 5L, 5L, 
    5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 
    5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 
    5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 
    5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 
    5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 
    5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 
    5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 
    5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 
    5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 
    5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L), .Label = c("", 
    "29-v1", "29-v2", "88-v1", "88-v2"), class = "factor"), Id = c(1L, 
    2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 
    16L, 17L, 18L, 19L, 20L, 21L, 22L, 23L, 24L, 25L, 26L, 27L, 28L, 
    29L, 30L, 31L, 32L, 33L, 34L, 35L, 36L, 37L, 38L, 39L, 40L, 41L, 
    42L, 43L, 44L, 45L, 46L, 47L, 48L, 49L, 50L, 51L, 52L, 53L, 54L, 
    55L, 56L, 57L, 58L, 59L, 60L, 61L, 62L, 63L, 64L, 65L, 66L, 67L, 
    68L, 69L, 70L, 71L, 72L, 73L, 74L, 75L, 76L, 77L, 78L, 79L, 80L, 
    81L, 82L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 
    13L, 14L, 15L, 16L, 17L, 18L, 19L, 20L, 21L, 22L, 23L, 24L, 25L, 
    26L, 27L, 28L, 29L, 30L, 31L, 32L, 33L, 34L, 35L, 36L, 37L, 38L, 
    39L, 40L, 41L, 42L, 43L, 44L, 45L, 46L, 47L, 48L, 49L, 50L, 51L, 
    52L, 53L, 54L, 55L, 56L, 57L, 58L, 59L, 60L, 61L, 62L, 63L, 64L, 
    65L, 66L, 67L, 68L, 69L, 70L, 71L, 72L, 73L, 74L, 75L, 76L, 77L, 
    78L, 79L, 80L, 81L, 82L, 83L, 84L, 85L, 86L, 87L, 88L, 89L, 90L, 
    91L, 92L, 93L, 94L, 95L, 96L, 97L, 98L, 99L, 100L, 101L, 102L, 
    103L, 104L, 105L, 106L, 107L, 108L, 109L, 110L, 111L, 112L, 113L, 
    114L, 115L, 116L, 117L, 118L, 119L, 120L, 121L, 122L, 123L, 124L, 
    125L, 126L, 127L, 128L, 129L, 130L, 131L, 132L, 133L, 134L, 135L, 
    136L, 137L, 138L, 139L, 140L, 141L, 142L, 143L, 144L, 145L, 146L, 
    147L, 148L, 149L, 150L, 151L, 152L, 153L, 154L, 155L, 156L, 157L, 
    158L, 159L, 160L, 161L, 162L, 163L, 164L, 165L, 166L, 167L, 1L, 
    2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 
    16L, 17L, 18L, 19L, 20L, 21L, 22L, 23L, 24L, 25L, 26L, 27L, 28L, 
    29L, 30L, 31L, 32L, 33L, 34L, 35L, 36L, 37L, 38L, 39L, 40L, 41L, 
    42L, 43L, 44L, 45L, 46L, 47L, 48L, 49L, 50L, 51L, 52L, 53L, 54L, 
    55L, 56L, 57L, 58L, 59L, 60L, 61L, 62L, 63L, 64L, 65L, 66L, 67L, 
    68L, 69L, 70L, 71L, 72L, 73L, 74L, 75L, 76L, 77L, 78L, 79L, 80L, 
    81L, 82L, 83L, 84L, 85L, 86L, 87L, 88L, 89L, 90L, 91L, 92L, 93L, 
    94L, 95L, 96L, 97L, 98L, 99L, 100L, 101L, 102L, 103L, 1L, 2L, 
    3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 
    17L, 18L, 19L, 20L, 21L, 22L, 23L, 24L, 25L, 26L, 27L, 28L, 29L, 
    30L, 31L, 32L, 33L, 34L, 35L, 36L, 37L, 38L, 39L, 40L, 41L, 42L, 
    43L, 44L, 45L, 46L, 47L, 48L, 49L, 50L, 51L, 52L, 53L, 54L, 55L, 
    56L, 57L, 58L, 59L, 60L, 61L, 62L, 63L, 64L, 65L, 66L, 67L, 68L, 
    69L, 70L, 71L, 72L, 73L, 74L, 75L, 76L, 77L, 78L, 79L, 80L, 81L, 
    82L, 83L, 84L, 85L, 86L, 87L, 88L, 89L, 90L, 91L, 92L, 93L, 94L, 
    95L, 96L, 97L, 98L, 99L, 100L, 101L, 102L, 103L, 104L, 105L, 
    106L, 107L, 108L, 109L, 110L, 111L, 112L, 113L, 114L, 115L, 116L, 
    117L, 118L, 119L, 120L, 121L, 122L, 123L, 124L, 125L, 126L, 127L, 
    128L, 129L, 130L, 131L, 132L, 133L, 134L, 135L, 136L, 137L, 138L, 
    139L, 140L, 141L, 142L, 143L, 144L, 145L, 146L, 147L, 148L, 149L, 
    150L, 151L, 152L, 153L, 154L, 155L, 156L, 157L, 158L, 159L, 160L, 
    161L, 162L, 163L, 164L), Age = c(3, 10, 1, 5, 4, 2, 1, 8, 2, 
    13, 1, 6, 3, 5, 2, 1, 3, 8, 3, 6, 1, 3, 7, 1, 2, 2, 2, 1, 2, 
    5, 4, 1, 6, 3, 6, 8, 2, 3, 4, 7, 3, 2, 6, 2, 3, 7, 1, 5, 4, 1, 
    4, 3, 2, 3, 5, 5, 2, 1, 1, 5, 8, 7, 2, 2, 4, 3, 4, 4, 2, 2, 10, 
    7, 5, 3, 3, 5, 7, 5, 3, 4, 5, 4, 1, 8, 6, 1, 12, 1, 6, 3, 4, 
    4, 13, 5, 2, 7, 7, 20, 1, 1, 1, 7, 1, 4, 3, 8, 2, 2, 4, 1, 1, 
    2, 3, 2, 2, 6, 11, 2, 5, 5, 9, 4, 4, 2, 7, 2, 7, 10, 6, 9, 2, 
    2, 5, 11, 1, 8, 8, 4, 1, 2, 14, 11, 13, 20, 3, 3, 4, 16, 2, 6, 
    11, 9, 11, 4, 5, 6, 19, 5, 2, 6, 1, 7, 11, 3, 9, 2, 3, 6, 20, 
    8, 6, 2, 11, 18, 9, 3, 7, 3, 2, 1, 8, 3, 5, 6, 2, 5, 8, 11, 4, 
    9, 7, 2, 12, 8, 2, 9, 5, 4, 15, 5, 13, 5, 10, 13, 7, 6, 1, 12, 
    12, 10, 4, 2, 16, 7, 17, 11, 18, 4, 3, 12, 1, 3, 7, 3, 6, 5, 
    11, 10, 12, 6, 14, 8, 6, 7, 8, 5, 10, 12, 6, 13, 3, 11, 14, 7, 
    9, 9, 4, 13, 4, 2, 1, 2, 2, 1, 7, 9, 3, 10, 3, 2, 1, 3, 1, 4, 
    2, 4, 5, 4, 2, 13, 4, 1, 3, 1, 11, 4, 1, 3, 3, 7, 5, 4, 5, 6, 
    1, 2, 1, 2, 1, 6, 1, 7, 6, 9, 5, 1, 6, 3, 2, 3, 3, 8, 8, 3, 2, 
    2, 4, 2, 5, 2, 6, 8, 11, 1, 6, 3, 3, 4, 5, 5, 7, 4, 2, 7, 3, 
    3, 1, 3, 9, 5, 2, 4, 12, 1, 4, 5, 2, 7, 6, 1, 2, 6, 4, 2, 7, 
    3, 5, 5, 3, 7, 1, 5, 2, 1, 15, 3, 5, 2, 5, 13, 6, 2, 3, 5, 2, 
    8, 4, 2, 6, 7, 2, 4, 1, 13, 8, 2, 1, 2, 1, 1, 5, 2, 1, 6, 11, 
    4, 1, 7, 7, 4, 3, 5, 1, 4, 10, 1, 2, 6, 1, 11, 3, 8, 9, 2, 6, 
    8, 11, 14, 16, 4, 1, 4, 2, 1, 10, 4, 9, 3, 12, 8, 11, 8, 8, 5, 
    1, 4, 13, 3, 8, 5, 14, 3, 5, 5, 12, 1, 3, 4, 5, 2, 7, 6, 9, 6, 
    10, 5, 2, 3, 2, 10, 10, 10, 10, 10, 1, 14, 3, 5, 9, 6, 2, 2, 
    2, 4, 4, 11, 14, 2, 2, 2, 8, 7, 2, 10, 12, 1, 6, 10, 2, 3, 5, 
    10, 6, 1, 8, 4, 11, 5, 4, 3, 6, 2, 4, 6, 9, 3, 9, 11, 7, 3, 15, 
    3, 7, 3, 5, 4, 6, 9, 13, 8, 5, 7, 8, 8, 5, 10), Type_product = c("f", 
    "s", "f", "f", "f", "f", "s", "c", "s", "f", "c", "f", "f", "f", 
    "s", "s", "f", "f", "c", "f", "s", "f", "f", "s", "f", "c", "f", 
    "f", "s", "f", "f", "c", "f", "c", "f", "f", "f", "f", "f", "c", 
    "c", "c", "f", "f", "c", "c", "f", "c", "c", "c", "c", "c", "s", 
    "f", "c", "c", "c", "s", "f", "c", "f", "f", "c", "c", "f", "c", 
    "c", "c", "f", "c", "c", "c", "c", "c", "c", "c", "f", "c", "c", 
    "c", "c", "f", "c", "f", "f", "s", "f", "c", "f", "f", "f", "c", 
    "f", "f", "f", "f", "f", "s", "c", "c", "f", "f", "c", "c", "f", 
    "f", "c", "c", "f", "f", "s", "f", "c", "c", "f", "f", "f", "c", 
    "f", "f", "f", "c", "f", "f", "f", "f", "f", "f", "c", "f", "f", 
    "f", "f", "c", "s", "f", "c", "f", "f", "c", "f", "f", "f", "c", 
    "f", "c", "c", "c", "f", "f", "f", "f", "c", "c", "c", "f", "f", 
    "c", "c", "f", "c", "f", "f", "c", "c", "c", "c", "f", "f", "f", 
    "c", "c", "c", "f", "c", "f", "c", "f", "f", "f", "c", "f", "c", 
    "c", "c", "c", "c", "f", "c", "c", "c", "c", "c", "c", "c", "f", 
    "f", "f", "c", "f", "c", "f", "f", "c", "c", "f", "f", "f", "c", 
    "c", "c", "f", "c", "c", "c", "c", "c", "f", "c", "f", "f", "c", 
    "c", "f", "c", "f", "c", "f", "c", "c", "c", "f", "c", "c", "c", 
    "c", "c", "c", "c", "f", "c", "c", "f", "c", "c", "f", "f", "c", 
    "f", "f", "s", "c", "s", "c", "f", "c", "c", "s", "c", "c", "s", 
    "c", "m", "c", "c", "f", "f", "f", "f", "f", "f", "s", "f", "f", 
    "c", "c", "f", "c", "f", "f", "f", "c", "f", "f", "f", "s", "f", 
    "f", "c", "f", "c", "f", "m", "c", "c", "c", "f", "s", "f", "f", 
    "f", "c", "s", "c", "m", "f", "c", "m", "c", "f", "c", "f", "f", 
    "f", "c", "m", "f", "c", "c", "f", "c", "f", "c", "c", "c", "c", 
    "c", "f", "f", "f", "c", "m", "f", "m", "m", "c", "c", "c", "c", 
    "m", "m", "c", "f", "m", "m", "m", "m", "m", "m", "m", "m", "m", 
    "c", "c", "f", "f", "f", "f", "c", "f", "m", "f", "f", "f", "c", 
    "f", "f", "f", "c", "f", "f", "c", "c", "f", "c", "f", "c", "m", 
    "f", "c", "f", "c", "f", "f", "f", "f", "c", "c", "f", "f", "c", 
    "c", "f", "f", "f", "f", "f", "f", "c", "f", "c", "c", "f", "c", 
    "f", "f", "f", "f", "f", "f", "f", "c", "f", "c", "f", "c", "f", 
    "c", "f", "c", "f", "f", "c", "c", "c", "c", "c", "f", "f", "f", 
    "c", "f", "c", "f", "f", "c", "c", "f", "f", "c", "f", "c", "f", 
    "c", "c", "c", "f", "f", "c", "f", "c", "c", "f", "c", "f", "c", 
    "f", "c", "f", "c", "m", "c", "c", "m", "c", "c", "f", "c", "c", 
    "f", "c", "c", "c", "f", "c", "c", "m", "c", "m", "m", "c", "c", 
    "f", "c", "c", "c", "c", "m", "c", "c", "c", "m", "m", "m", "c", 
    "c", "c", "c", "m", "m", "f", "m", "m", "m", "m", "m", "m", "m", 
    "m", "m", "m", "m", "m", "m", "m", "m"), Exhaustion_product = structure(c(1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 
    3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 5L, 
    5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 
    7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 
    9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 10L, 10L, 10L, 10L, 10L, 10L, 
    10L, 10L, 10L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 
    3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 
    4L, 4L, 4L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 
    5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 
    6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 
    7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 
    8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 9L, 9L, 9L, 9L, 9L, 9L, 
    9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 10L, 10L, 10L, 10L, 
    10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 
    4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 
    5L, 5L, 5L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 
    7L, 7L, 7L, 7L, 7L, 7L, 7L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 
    8L, 8L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 10L, 10L, 10L, 
    10L, 10L, 10L, 10L, 10L, 10L, 10L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 
    3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 
    4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 5L, 5L, 
    5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 6L, 
    6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 7L, 
    7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 8L, 8L, 8L, 8L, 8L, 
    8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 9L, 9L, 9L, 9L, 
    9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 10L, 10L, 10L, 
    10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 
    10L), .Label = c("(0,10]", "(10,20]", "(20,30]", "(30,40]", "(40,50]", 
    "(50,60]", "(60,70]", "(70,80]", "(80,90]", "(90,100]"), class = "factor"), 
        Type_f = c(1, 0, 1, 1, 1, 1, 0, 0, 0, 1, 0, 1, 1, 1, 0, 0, 
        1, 1, 0, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 1, 0, 1, 0, 1, 
        1, 1, 1, 1, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 
        0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 
        0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 1, 1, 1, 0, 
        1, 1, 1, 1, 1, 0, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 
        1, 0, 0, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0, 1, 
        1, 1, 1, 0, 0, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 0, 0, 0, 1, 
        1, 1, 1, 0, 0, 0, 1, 1, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 1, 
        1, 1, 0, 0, 0, 1, 0, 1, 0, 1, 1, 1, 0, 1, 0, 0, 0, 0, 0, 
        1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 1, 0, 1, 1, 0, 0, 1, 
        1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 1, 0, 
        1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 
        0, 1, 1, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 
        0, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 0, 0, 1, 0, 1, 1, 1, 0, 
        1, 1, 1, 0, 1, 1, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 1, 1, 1, 
        0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 1, 1, 0, 0, 1, 0, 0, 1, 
        0, 1, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 
        0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 
        0, 1, 0, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 0, 0, 1, 0, 1, 0, 
        0, 1, 0, 1, 0, 1, 1, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1, 1, 1, 
        1, 1, 0, 1, 0, 0, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 1, 0, 1, 
        0, 1, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1, 1, 0, 1, 0, 1, 
        1, 0, 0, 1, 1, 0, 1, 0, 1, 0, 0, 0, 1, 1, 0, 1, 0, 0, 1, 
        0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 
        0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 
        0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
        0, 0, 0, 0, 0, 0)), .Names = c("Type_space", "Id", "Age", 
    "Type_product", "Exhaustion_product", "Type_f"), row.names = c(1L, 
    2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 
    16L, 17L, 18L, 19L, 20L, 21L, 22L, 23L, 24L, 25L, 26L, 27L, 28L, 
    29L, 30L, 31L, 32L, 33L, 34L, 35L, 36L, 37L, 38L, 39L, 40L, 41L, 
    42L, 43L, 44L, 45L, 46L, 47L, 48L, 49L, 50L, 51L, 52L, 53L, 54L, 
    55L, 56L, 57L, 58L, 59L, 60L, 61L, 62L, 63L, 64L, 65L, 66L, 67L, 
    68L, 69L, 70L, 71L, 73L, 75L, 76L, 79L, 80L, 81L, 82L, 84L, 85L, 
    86L, 91L, 102L, 103L, 104L, 105L, 106L, 107L, 108L, 109L, 110L, 
    111L, 112L, 113L, 114L, 115L, 116L, 117L, 118L, 119L, 120L, 121L, 
    122L, 123L, 124L, 125L, 126L, 127L, 128L, 129L, 130L, 131L, 132L, 
    133L, 134L, 135L, 136L, 137L, 138L, 139L, 140L, 141L, 142L, 143L, 
    144L, 145L, 146L, 147L, 148L, 149L, 150L, 151L, 152L, 153L, 154L, 
    155L, 156L, 157L, 158L, 159L, 160L, 161L, 162L, 163L, 164L, 165L, 
    166L, 167L, 168L, 169L, 170L, 171L, 172L, 173L, 174L, 175L, 176L, 
    177L, 178L, 179L, 180L, 181L, 182L, 183L, 184L, 185L, 186L, 187L, 
    188L, 189L, 190L, 191L, 192L, 193L, 194L, 195L, 197L, 198L, 199L, 
    201L, 202L, 203L, 204L, 206L, 207L, 208L, 209L, 210L, 212L, 213L, 
    214L, 215L, 217L, 218L, 219L, 220L, 221L, 222L, 223L, 225L, 227L, 
    229L, 230L, 231L, 232L, 233L, 234L, 235L, 236L, 237L, 238L, 239L, 
    242L, 243L, 244L, 246L, 247L, 248L, 249L, 250L, 251L, 253L, 254L, 
    256L, 259L, 260L, 261L, 262L, 263L, 264L, 265L, 266L, 269L, 270L, 
    272L, 273L, 274L, 276L, 277L, 278L, 279L, 280L, 281L, 282L, 283L, 
    284L, 285L, 287L, 289L, 290L, 291L, 292L, 293L, 294L, 295L, 296L, 
    297L, 298L, 300L, 301L, 302L, 303L, 306L, 308L, 309L, 311L, 312L, 
    313L, 314L, 315L, 316L, 317L, 318L, 319L, 320L, 322L, 323L, 325L, 
    326L, 327L, 328L, 329L, 331L, 332L, 334L, 335L, 336L, 338L, 339L, 
    340L, 341L, 342L, 343L, 344L, 345L, 346L, 347L, 348L, 349L, 350L, 
    352L, 353L, 354L, 356L, 357L, 358L, 359L, 360L, 361L, 363L, 364L, 
    365L, 366L, 367L, 368L, 369L, 370L, 372L, 373L, 374L, 375L, 376L, 
    377L, 378L, 379L, 380L, 381L, 382L, 384L, 385L, 387L, 388L, 389L, 
    391L, 393L, 394L, 395L, 396L, 397L, 398L, 399L, 400L, 401L, 402L, 
    404L, 407L, 408L, 409L, 411L, 412L, 413L, 414L, 415L, 416L, 417L, 
    418L, 419L, 420L, 421L, 422L, 423L, 424L, 425L, 426L, 427L, 428L, 
    429L, 430L, 431L, 432L, 433L, 434L, 435L, 436L, 437L, 438L, 439L, 
    440L, 442L, 443L, 444L, 445L, 446L, 447L, 448L, 449L, 450L, 451L, 
    452L, 453L, 454L, 455L, 456L, 457L, 458L, 459L, 460L, 461L, 462L, 
    463L, 464L, 465L, 466L, 467L, 468L, 469L, 470L, 471L, 472L, 473L, 
    474L, 476L, 477L, 478L, 479L, 480L, 481L, 482L, 483L, 484L, 486L, 
    487L, 488L, 489L, 490L, 491L, 492L, 493L, 494L, 495L, 496L, 497L, 
    498L, 500L, 501L, 502L, 503L, 504L, 505L, 506L, 507L, 508L, 509L, 
    510L, 511L, 512L, 513L, 514L, 515L, 516L, 517L, 518L, 519L, 520L, 
    521L, 522L, 523L, 524L, 525L, 526L, 527L, 528L, 529L, 530L, 531L, 
    532L, 534L, 535L, 536L, 537L, 538L, 539L, 540L, 541L, 542L, 543L, 
    547L, 548L, 550L, 551L, 552L, 553L, 554L, 555L, 556L, 557L, 558L, 
    559L, 560L, 561L, 562L, 563L, 565L, 566L, 567L, 568L, 569L, 570L, 
    571L, 572L, 573L, 575L, 577L, 579L, 580L, 581L, 582L, 583L, 585L, 
    586L, 587L, 590L, 592L, 599L, 606L, 608L), class = "data.frame")

    an=Anova(glm(Type_f ~  Type_space  + Exhaustion_product + Age , family=binomial,data=res))
    gl=glm(Type_f ~  Type_space  + Exhaustion_product + Age  , family=binomial,data=res)
    library("emmeans")
    emmp <- emmeans( gl, pairwise ~ Exhaustion_product + Age)
    summary( emmp, infer=TRUE)

(1) 对于分类变量,结果很明确。但对于年龄在 GLM 中很重要的情况,在 GLM 中生成的值是多少?emmeans ?5.455426.这就是手段吗?我该如何解释这一点?

 (0,10]             5.455426  0.36901411 0.2935894 Inf -0.20641061  0.94443883   1.257  0.2088

(2)我想生成交互的图形表示age and Exhaustion_product。这也是没有意义的。

emmip(gl, Exhaustion_product ~ Age)

编辑1 对比结果

$contrasts
 contrast                                                estimate        SE  df   asymp.LCL asymp.UCL z.ratio p.value
 (0,10],5.45542635658915 - (10,20],5.45542635658915    0.33231353 0.4078967 Inf -0.95814279 1.6227698   0.815  0.9984
 (0,10],5.45542635658915 - (20,30],5.45542635658915   -0.53694399 0.4194460 Inf -1.86393835 0.7900504  -1.280  0.9582
 (0,10],5.45542635658915 - (30,40],5.45542635658915   -0.16100309 0.4139472 Inf -1.47060101 1.1485948  -0.389  1.0000
 (0,10],5.45542635658915 - (40,50],5.45542635658915    0.40113723 0.4021403 Inf -0.87110757 1.6733820   0.998  0.9925
 (0,10],5.45542635658915 - (50,60],5.45542635658915    0.60576562 0.4106536 Inf -0.69341247 1.9049437   1.475  0.9022
 (0,10],5.45542635658915 - (60,70],5.45542635658915    1.38800301 0.4319258 Inf  0.02152631 2.7544797   3.214  0.0430
 (0,10],5.45542635658915 - (70,80],5.45542635658915    1.01677522 0.4147441 Inf -0.29534399 2.3288944   2.452  0.2952
 (0,10],5.45542635658915 - (80,90],5.45542635658915    1.99085692 0.4747929 Inf  0.48876247 3.4929514   4.193  0.0011
 (0,10],5.45542635658915 - (90,100],5.45542635658915   2.03923289 0.4745872 Inf  0.53778910 3.5406767   4.297  0.0007

因为这个问题看起来像是一个自学题,所以我将做一个类似的例子,而不是相同的数据。但结构是相同的,只有一个因子和一个协变量作为预测变量。

例子是emmeans::fiber数据集。其响应变量是纤维强度,连续预测变量是直径,因子是制造它的机器。

Model:

> mod = glm(log(strength) ~ machine + diameter, data = fiber)
> summary(mod)
... (output has been abbreviated) ...
Coefficients:
             Estimate Std. Error t value Pr(>|t|)
(Intercept)  3.124387   0.068374  45.695 6.74e-14
machineB     0.026025   0.023388   1.113    0.290
machineC    -0.044593   0.025564  -1.744    0.109
diameter     0.023557   0.002633   8.946 2.22e-06

(Dispersion parameter for gaussian family taken to be 0.001356412)

emmeans 分析基于参考网格,默认情况下由所有级别的因子和协变量的平均值组成:

> ref_grid(mod)
'emmGrid' object with variables:
    machine = A, B, C
    diameter = 24.133
Transformation: “log” 

您可以在 R 中确认mean(fiber$diameter)是 24.133。我强调这就是diameter值,而不是模型中的任何内容。

> summary(.Last.value)
 machine diameter prediction         SE  df
 A       24.13333   3.692901 0.01670845 Inf
 B       24.13333   3.718925 0.01718853 Inf
 C       24.13333   3.648307 0.01819206 Inf

Results are given on the log (not the response) scale.

这些汇总值是来自的预测mod在每个组合machine and diameter。现在看看 EMMmachine

> emmeans(mod, "machine")
 machine   emmean         SE  df asymp.LCL asymp.UCL
 A       3.692901 0.01670845 Inf  3.660153  3.725649
 B       3.718925 0.01718853 Inf  3.685237  3.752614
 C       3.648307 0.01819206 Inf  3.612652  3.683963

Results are given on the log (not the response) scale. 
Confidence level used: 0.95

...我们得到了完全相同的三个预测。但如果我们看看diameter:

> emmeans(mod, "diameter")
 diameter   emmean          SE  df asymp.LCL asymp.UCL
 24.13333 3.686711 0.009509334 Inf  3.668073  3.705349

Results are averaged over the levels of: machine 
Results are given on the log (not the response) scale. 
Confidence level used: 0.95

...我们得到 EMM 等于参考网格中三个预测值的平均值。请注意,它在注释中表示结果的平均值为machine,所以值得一读。

为了获得模型结果的图形表示,我们可以这样做

> emmip(mod, machine ~ diameter, cov.reduce = range)

论点cov.reduce = range添加以使参考网格使用最小和最大直径,而不是其平均值。如果没有这个,我们就会得到三个点而不是三条线。该图仍然显示模型预测,只是在更详细的值网格上。请注意,所有三条线都有相同的斜率。那是因为模型是这样指定的:diameter效果是added to the machine影响。因此,每条线的共同斜率为 0.023557(参见输出summary(mod).

没有post hoc需要测试diameter,由于其one效果已经测试过summary(mod).

最后一件事。使用的型号log(strength)作为回应。如果我们希望 EMM 的规模与strength,只需添加type = "response":

> emmeans(mod, "machine", type = "response")
 machine response        SE  df asymp.LCL asymp.UCL
 A       40.16118 0.6710311 Inf  38.86728  41.49815
 B       41.22008 0.7085126 Inf  39.85455  42.63239
 C       38.40960 0.6987496 Inf  37.06421  39.80384

Confidence level used: 0.95 
Intervals are back-transformed from the log scale

同样,结果下方的注释有助于解释输出。

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

Emmeans 连续自变量 的相关文章

随机推荐

  • PHP内部hashCode函数

    我正在寻找与 JAVA 等效的 PHP SomeString hashCode 功能 我正在寻找的 hashCode 应该与 PHP 中用于索引 Hashmap 的 hashCode 相同 我希望你可以帮助我 EDIT 好的 找到了我正在搜
  • WPF 重复按钮 MouseUp

    有没有办法让重复按钮上的 MouseUpevent 在不再按下按钮时触发 我试图使用 MouseMove 事件来跟踪按下按钮时鼠标的位置 但 MouseDown 和 MouseUp 都不会触发鼠标左键的事件 关于可以做什么的任何想法或建议
  • 如何将图像(位图)调整为给定尺寸? [复制]

    这个问题在这里已经有答案了 如何以编程方式将图像 位图 大小调整为例如 800 480 我在我的应用程序中检索了一张约 1MB 的图片 我需要将其缩小到 800 480 我已经加载了该图片并对其进行了压缩 但是如何使其变小 ByteArra
  • 选择组合框文本

    我有一个组合框声明如下
  • 将Pycharm项目移动到另一个目录

    我的主目录空间不足 我想将 PyCharm 项目移动到另一个目录 我最终复制了它 因为重构不起作用 我删除了 pycache 和 zip 异常 出现的异常消失了 现在一切都按预期工作 不过 从 settings python interpr
  • 如何使用 GitHub README 中的 Markdown 渲染多列?

    为了呈现三列中的项目 我尝试将以下 CSS3 指令添加到我的项目中README md文件 但样式被删除 div div style display inline block div div div div
  • 从我的应用程序创建日历事件,无需默认提醒

    我正在开发一个应用程序 可以创建 更新和删除原生 Google 日历中的事件 我正在通过以下代码创建一个事件 ContentValues cvEvent new ContentValues cvEvent put Events DTSTAR
  • `geom_line()` 连接映射到不同组的点

    我想group我的数据基于两个变量的相互作用 但仅将美学映射到其中一个变量 另一个变量代表重复 理论上 它们应该彼此相等 我可以找到不优雅的方法来做到这一点 但似乎应该有更优雅的方法来做到这一点 例如 Data frame with two
  • 如何使用 url_for() 将路径和查询数据传递到使用 FastAPI 和 Jinja2 的路由?

    我想从 Jinja2 模板中调用 FastAPI 路由 并传递两者path and query该路由的数据 参数 我在 Jinja2 模板中尝试过如下所示 url for function1 uustr data uustr interva
  • 从 R 中的 Sys.Date() 获取当前时间(包括秒)

    我确信这是非常基本的 但我无法从系统日期中提取各种时间 简而言之 我不确定为什么我无法使用 M 或 S 准确提取当前的分钟和秒数 有什么想法吗 我的机器上当前时间为中午 12 38 gt format Sys Date c 1 12 16
  • “冲突可序列化”和“冲突等效”有什么区别?

    在数据库理论中 冲突可序列化 和 冲突等效 有什么区别 我的教科书有一个关于可序列化冲突的部分 但掩盖了冲突等价性 这些可能都是我熟悉的概念 但我不熟悉术语 所以我正在寻找解释 DBMS 中的冲突可以定义为两个或多个不同的事务访问同一变量
  • 将带有变量的循环转换为 BigQuery SQL

    我有数千个脚本 其中包括循环数据集并使用变量进行累积 例如 assuming that ids is populated from some BQ table ids 1 2 3 4 5 var1 v1 initialize variabl
  • 如何使用 RealmSwift 存储字典?

    考虑以下模型 class Person Object dynamic var name let hobbies Dictionary
  • 清漆配置文件在哪里?

    我刚刚安装了清漆并准备好提供网页服务 问题是我们达到了打开文件限制 我们已经设定了ulimit给清漆用户 我们想要找到一个可以写入此配置的文件 Maximum number of open files NFILES 131072 Locke
  • ios7 UITableViewCell SelectionStyle 不会变回蓝色

    Xcode 5 0 iOS 7 并更新现有应用程序 UITableView所选行现在是灰色的 而不是蓝色的 据我所知 他们更改了默认值selectionStyle至灰色 但 蓝色 仍然是 IB 的一个选项UITableViewCellSel
  • 在内核模块中执行shell命令

    是否可以在内核模块中执行shell命令 我知道我们可以在用户空间 C 代码中使用system子程序 我正在调试一个存在内存泄漏问题的内核模块 在无限循环中执行 insmod 和 rmmod module ko 后 8G RAM 的系统在几分
  • boost Deadline_timer 最小示例:我应该替换“睡眠”吗?

    我有一个线程 我需要每 10 毫秒执行一些操作 所以我有非常简单的代码 如下所示 while work do something Sleep 10000 boost sleep can be also used 我听说Sleep一般不推荐
  • 避免在react + typescript中使用通用函数类型

    我正在寻找一条规则来阻止使用 函数 作为类型 myMethod Function 我没有找到任何东西 所以我愿意接受建议 您可以使用 typescript eslint ban types rule 规则链接 https github co
  • Caffeine Expiry 中如何设置多个过期标准?

    我正在使用 Caffeine v2 8 5 我想创建一个具有可变到期时间的缓存 基于 值的创建 更新以及 该值的最后一次访问 读取 无论先发生什么都应该触发该条目的删除 缓存将成为三层值解析的一部分 The key is present i
  • Emmeans 连续自变量

    我想解释一下Type f with Type space实验的内容和速率Exhaustion product和定量变量Age 这是我的数据 res structure list Type space structure c 2L 2L 2L