API tools faq
paste
Advertisement
Guest User

Untitled

a guest
Jun 25th, 2019
559
0
Never
Add comment
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 10.28 KB | None | 0 0
raw download clone embed print report
  1. mod1 <- gls(data = my_data, model = Height ~ Origin*Genotype*Substrate*Date)
  2.  
  3. Error in glsEstimate(glsSt, control = glsEstControl) :
  4. l'ajustement "gls" calculé est singulier, de rang 61
  5.  
  6. mod1 <- gls(data = my_data, model = Height ~ Genotype*Substrate*Date)
  7.  
  8. mod2 <- gls(data = my_data, model = Height ~ Genotype*Substrate*Date,
  9. correlation = corCAR1(form = ~ Date | Sample_ID))
  10.  
  11. > AICc(mod1, mod2)
  12. df AICc
  13. mod1 61 2385.981
  14. mod2 62 2389.311
  15.  
  16. my_data <- structure(list(Sample_ID = structure(c(1L, 2L, 3L, 4L, 5L, 6L,
  17. 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L,
  18. 20L, 21L, 22L, 23L, 24L, 25L, 26L, 27L, 28L, 29L, 30L, 31L, 32L,
  19. 33L, 34L, 35L, 36L, 37L, 38L, 39L, 40L, 41L, 42L, 43L, 44L, 45L,
  20. 46L, 47L, 48L, 49L, 50L, 51L, 52L, 53L, 54L, 55L, 56L, 57L, 58L,
  21. 59L, 60L, 61L, 62L, 63L, 64L, 65L, 66L, 67L, 68L, 69L, 70L, 71L,
  22. 72L, 73L, 74L, 75L, 76L, 77L, 78L, 79L, 80L, 81L, 82L, 83L, 84L,
  23. 85L, 86L, 87L, 88L, 89L, 90L, 91L, 92L, 1L, 2L, 3L, 4L, 5L, 6L,
  24. 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L,
  25. 20L, 21L, 22L, 23L, 24L, 25L, 26L, 27L, 28L, 29L, 30L, 31L, 32L,
  26. 33L, 34L, 35L, 36L, 37L, 38L, 39L, 40L, 41L, 42L, 43L, 44L, 45L,
  27. 46L, 47L, 48L, 49L, 50L, 51L, 52L, 53L, 54L, 55L, 56L, 57L, 58L,
  28. 59L, 60L, 61L, 62L, 63L, 64L, 65L, 66L, 67L, 68L, 69L, 70L, 71L,
  29. 72L, 73L, 74L, 75L, 76L, 77L, 78L, 79L, 80L, 81L, 82L, 83L, 84L,
  30. 85L, 86L, 87L, 88L, 89L, 90L, 91L, 92L, 1L, 2L, 3L, 4L, 5L, 6L,
  31. 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L,
  32. 20L, 21L, 22L, 23L, 24L, 25L, 26L, 27L, 28L, 29L, 30L, 31L, 32L,
  33. 33L, 34L, 35L, 36L, 37L, 38L, 39L, 40L, 41L, 42L, 43L, 44L, 45L,
  34. 46L, 47L, 48L, 49L, 50L, 51L, 52L, 53L, 54L, 55L, 56L, 57L, 58L,
  35. 59L, 60L, 61L, 62L, 63L, 64L, 65L, 66L, 67L, 68L, 69L, 70L, 71L,
  36. 72L, 73L, 74L, 75L, 76L, 77L, 78L, 79L, 80L, 81L, 82L, 83L, 84L,
  37. 85L, 86L, 87L, 88L, 89L, 90L, 91L, 92L), .Label = c("08C1", "08C2",
  38. "08C3", "08M1", "08M2", "08M3", "08W1", "08W2", "08W3", "08W4",
  39. "09C1", "09C2", "09C3", "09M1", "09M2", "09M3", "09W1", "09W2",
  40. "09W3", "09W4", "10C1", "10C2", "10C3", "10M1", "10M2", "10M3",
  41. "10W1", "10W2", "10W3", "13C1", "13C2", "13C3", "13M1", "13M2",
  42. "13M3", "13W1", "13W2", "13W3", "16C1", "16C2", "16M1", "16M2",
  43. "16W1", "16W2", "16W3", "21C1", "21C2", "21C3", "21M1", "21M2",
  44. "21M3", "21W1", "21W2", "21W3", "23C1", "23C2", "23C3", "23M1",
  45. "23M2", "23M3", "23M4", "23W1", "23W2", "23W3", "25C1", "25C2",
  46. "25C3", "25M1", "25M2", "25M3", "25W1", "25W2", "25W3", "29C1",
  47. "29C2", "29C3", "29M1", "29M2", "29M3", "29W1", "29W2", "29W3",
  48. "33C1", "33C2", "33C3", "33M1", "33M2", "33M3", "33W1", "33W2",
  49. "33W3", "33W4"), class = "factor"), Origin = structure(c(3L,
  50. 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
  51. 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
  52. 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L,
  53. 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
  54. 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
  55. 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L,
  56. 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
  57. 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
  58. 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
  59. 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
  60. 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L,
  61. 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
  62. 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
  63. 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L,
  64. 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
  65. 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
  66. 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
  67. 2L, 2L, 2L), .Label = c("Moly", "Road", "West"), class = "factor"),
  68. Genotype = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
  69. 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L,
  70. 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 5L,
  71. 5L, 5L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L,
  72. 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 8L, 8L, 8L, 8L, 8L,
  73. 8L, 8L, 8L, 8L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 10L,
  74. 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 1L, 1L, 1L,
  75. 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
  76. 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L,
  77. 4L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 6L, 6L, 6L,
  78. 6L, 6L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L,
  79. 7L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 9L, 9L, 9L, 9L, 9L,
  80. 9L, 9L, 9L, 9L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L,
  81. 10L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L,
  82. 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
  83. 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 5L,
  84. 5L, 5L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L,
  85. 7L, 7L, 7L, 7L, 7L, 7L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L,
  86. 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 10L, 10L, 10L, 10L, 10L,
  87. 10L, 10L, 10L, 10L, 10L), .Label = c("PB08", "PB09", "PB10",
  88. "PB13", "PB16", "PB21", "PB23", "PB25", "PB29", "PB33"), class = "factor"),
  89. Substrate = structure(c(1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L,
  90. 3L, 1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 1L, 1L, 1L, 2L,
  91. 2L, 2L, 3L, 3L, 3L, 1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L, 1L,
  92. 1L, 2L, 2L, 3L, 3L, 3L, 1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L,
  93. 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 1L, 1L, 1L, 2L, 2L,
  94. 2L, 3L, 3L, 3L, 1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L, 1L, 1L,
  95. 1L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 1L, 1L, 1L, 2L, 2L, 2L, 3L,
  96. 3L, 3L, 3L, 1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 1L, 1L,
  97. 1L, 2L, 2L, 2L, 3L, 3L, 3L, 1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L,
  98. 3L, 1L, 1L, 2L, 2L, 3L, 3L, 3L, 1L, 1L, 1L, 2L, 2L, 2L, 3L,
  99. 3L, 3L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 1L, 1L, 1L,
  100. 2L, 2L, 2L, 3L, 3L, 3L, 1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L,
  101. 1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 1L, 1L, 1L, 2L, 2L,
  102. 2L, 3L, 3L, 3L, 3L, 1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L, 3L,
  103. 1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L, 1L, 1L, 1L, 2L, 2L, 2L,
  104. 3L, 3L, 3L, 1L, 1L, 2L, 2L, 3L, 3L, 3L, 1L, 1L, 1L, 2L, 2L,
  105. 2L, 3L, 3L, 3L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 1L,
  106. 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L, 1L, 1L, 1L, 2L, 2L, 2L, 3L,
  107. 3L, 3L, 1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L, 3L), .Label = c("Control",
  108. "Molybdenum", "Westwood"), class = "factor"), Date = structure(c(17280,
  109. 17280, 17280, 17280, 17280, 17280, 17280, 17280, 17280, 17280,
  110. 17280, 17280, 17280, 17280, 17280, 17280, 17280, 17280, 17280,
  111. 17280, 17280, 17280, 17280, 17280, 17280, 17280, 17280, 17280,
  112. 17280, 17280, 17280, 17280, 17280, 17280, 17280, 17280, 17280,
  113. 17280, 17280, 17280, 17280, 17280, 17280, 17280, 17280, 17280,
  114. 17280, 17280, 17280, 17280, 17280, 17280, 17280, 17280, 17280,
  115. 17280, 17280, 17280, 17280, 17280, 17280, 17280, 17280, 17280,
  116. 17280, 17280, 17280, 17280, 17280, 17280, 17280, 17280, 17280,
  117. 17280, 17280, 17280, 17280, 17280, 17280, 17280, 17280, 17280,
  118. 17280, 17280, 17280, 17280, 17280, 17280, 17280, 17280, 17280,
  119. 17280, 17333, 17333, 17333, 17333, 17333, 17333, 17333, 17333,
  120. 17333, 17333, 17333, 17333, 17333, 17333, 17333, 17333, 17333,
  121. 17333, 17333, 17333, 17333, 17333, 17333, 17333, 17333, 17333,
  122. 17333, 17333, 17333, 17333, 17333, 17333, 17333, 17333, 17333,
  123. 17333, 17333, 17333, 17333, 17333, 17333, 17333, 17333, 17333,
  124. 17333, 17333, 17333, 17333, 17333, 17333, 17333, 17333, 17333,
  125. 17333, 17333, 17333, 17333, 17333, 17333, 17333, 17333, 17333,
  126. 17333, 17333, 17333, 17333, 17333, 17333, 17333, 17333, 17333,
  127. 17333, 17333, 17333, 17333, 17333, 17333, 17333, 17333, 17333,
  128. 17333, 17333, 17333, 17333, 17333, 17333, 17333, 17333, 17333,
  129. 17333, 17333, 17333, 17364, 17364, 17364, 17364, 17364, 17364,
  130. 17364, 17364, 17364, 17364, 17364, 17364, 17364, 17364, 17364,
  131. 17364, 17364, 17364, 17364, 17364, 17364, 17364, 17364, 17364,
  132. 17364, 17364, 17364, 17364, 17364, 17364, 17364, 17364, 17364,
  133. 17364, 17364, 17364, 17364, 17364, 17364, 17364, 17364, 17364,
  134. 17364, 17364, 17364, 17364, 17364, 17364, 17364, 17364, 17364,
  135. 17364, 17364, 17364, 17364, 17364, 17364, 17364, 17364, 17364,
  136. 17364, 17364, 17364, 17364, 17364, 17364, 17364, 17364, 17364,
  137. 17364, 17364, 17364, 17364, 17364, 17364, 17364, 17364, 17364,
  138. 17364, 17364, 17364, 17364, 17364, 17364, 17364, 17364, 17364,
  139. 17364, 17364, 17364, 17364, 17364), class = "Date"), Height = c(129.5,
  140. 140, 122, 112, 109, 127, 109, 134.5, 96.5, 109, 114.5, 104,
  141. 107, 117, 79, 119.5, 101.5, 63.5, 122, 74, 109.5, 112, 91.5,
  142. 104, 109, 104, 117, 114.5, 117, 114, 124.5, 111.5, 122, 99,
  143. 134.5, 84, 142, 117, 112, 71, 74, 76, 117, 48.5, 114.5, 81.5,
  144. 89, 109, 114.5, 78.5, 106.5, 89, 119.5, 119.5, 96.5, 96.5,
  145. 99, 91.5, 106.5, 33, 101.5, 89, 112, 117, 104, 122, 78.5,
  146. 89, 66, 76, 106.5, 104, 104, 94, 127, 132, 122, 122, 109,
  147. 129.5, 139.5, 78.5, 104, 139.5, 122, 122, 81.5, 134.5, 91.5,
  148. 139.5, 122, 114.5, 162.5, 193, 167.5, 137, 160, 172.5, 117,
  149. 167.5, 114.5, 172.5, 127, 109, 134.5, 144.5, 101.5, 137,
  150. 117, 111.5, 147.5, 134.5, 139.5, 165, 129.5, 152.5, 155,
  151. 137, 147.5, 142, 137, 147.5, 157.5, 160, 124.5, 142, 160,
  152. 144.5, 157.5, 127, 180.5, 147.5, 147.5, 142, 172.5, 122,
  153. 157.5, 157.5, 137, 172.5, 160, 111.5, 157.5, 144.5, 175.5,
  154. 157.5, 134.5, 150, 134.5, 137, 157.5, 111.5, 144.5, 122,
  155. 139.5, 132, 147.5, 167.5, 152.5, 142, 134.5, 152.5, 160,
  156. 160, 142, 178, 134.5, 185.5, 178, 160, 160, 160, 162.5, 139.5,
  157. 160, 185.5, 180.5, 170, 142, 183, 129.5, 165, 172.5, 172.5,
  158. 175, 193, 170, 137, 160, 172.5, 122, 167.5, 114.5, 172.5,
  159. 127, 124.5, 134.5, 144.5, 101.5, 139.5, 122, 114, 147.5,
  160. 134.5, 139.5, 165, 132, 155, 155, 142, 147.5, 142, 137, 147.5,
  161. 157.5, 160, 124.5, 142, 160, 144.5, 160, 127, 180.5, 150,
  162. 147.5, 142, 172.5, 124.5, 157.5, 157.5, 137, 175, 160, 114,
  163. 157.5, 144.5, 175.5, 157.5, 134.5, 152, 134.5, 144.5, 157.5,
  164. 119, 147, 124.5, 155, 132, 147.5, 167.5, 152.5, 147, 137,
  165. 155, 160, 160, 142, 178, 134.5, 185.5, 178, 160, 160, 160,
  166. 162.5, 139.5, 162.5, 188, 183, 172.5, 144.5, 183, 129.5,
  167. 165, 172.5, 172.5)), row.names = c(NA, -276L), class = "data.frame")
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement
Public Pastes
Advertisement
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy.  OK, I Understand
Not a member of Pastebin yet?
Sign Up, it unlocks many cool features!