lme4 error: boundary (singular) fit: see ?isSingular

Your model did fit, but it generated that warning because your random effects are very small. You can read more about this in this post or the help page

Let us look at your data:

ggplot(Data,aes(x=PLANT,y=Weight,col=Rep)) + geom_jitter() + geom_boxplot(alpha=0.2) + facet_wrap(~Rep)

The effects of PLANT and in combination with Rep is extremely small. Let’s look at the fitted model:

fit = lmer(Weight ~ 1 + (1|PLANT:Rep),data=Data)
boundary (singular) fit: see ?isSingular
ranef(fit)
$`PLANT:Rep`
    (Intercept)
1:1           0
1:2           0
2:1           0
2:2           0
3:1           0
3:2           0
4:1           0
4:2           0

This is exactly what happened. So we can try to account for some other effects and we still see very small coefficients:

fit = lmer(Weight ~ Line + (1|Rep:PLANT),data=Data)
ranef(fit)
$`Rep:PLANT`
      (Intercept)
1:1  1.397563e-19
1:2  2.811371e-19
1:3  8.112169e-20
1:4  1.813251e-19
2:1 -1.725964e-19
2:2 -2.463986e-20
2:3 -2.027357e-19
2:4 -2.833681e-19

The takehome message is, there’s no really systematic effect coming from PLANT, so you don’t need to specify a highly complicated model, do something like:

fit = lmer(Weight ~ Line + (1|Rep),data=Data)

The data in case anyone is interested:

Data = structure(list(Line = structure(c(1L, 1L, 1L, 1L, 12L, 12L, 12L, 
12L, 23L, 23L, 23L, 23L, 34L, 34L, 34L, 34L, 45L, 45L, 45L, 45L, 
56L, 56L, 56L, 56L, 65L, 65L, 65L, 65L, 66L, 66L, 66L, 66L, 67L, 
67L, 67L, 67L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 
5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 8L, 8L, 8L, 8L, 
9L, 9L, 9L, 9L, 10L, 10L, 10L, 10L, 11L, 11L, 11L, 11L, 13L, 
13L, 13L, 13L, 14L, 14L, 14L, 14L, 15L, 15L, 15L, 15L, 16L, 16L, 
16L, 16L, 8L, 8L, 8L, 8L, 66L, 66L, 66L, 66L, 17L, 17L, 17L, 
17L, 18L, 18L, 18L, 18L, 9L, 9L, 9L, 9L, 19L, 19L, 19L, 19L, 
20L, 20L, 20L, 20L, 21L, 21L, 21L, 21L, 22L, 22L, 22L, 22L, 24L, 
24L, 24L, 24L, 25L, 25L, 25L, 25L, 2L, 2L, 2L, 2L, 26L, 26L, 
26L, 26L, 27L, 27L, 27L, 27L, 10L, 10L, 10L, 10L, 28L, 28L, 28L, 
28L, 29L, 29L, 29L, 29L, 30L, 30L, 30L, 30L, 31L, 31L, 31L, 31L, 
67L, 67L, 67L, 67L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 33L, 
33L, 33L, 33L, 35L, 35L, 35L, 35L, 36L, 36L, 36L, 36L, 37L, 37L, 
37L, 37L, 38L, 38L, 38L, 38L, 39L, 39L, 39L, 39L, 40L, 40L, 40L, 
40L, 25L, 25L, 25L, 25L, 19L, 19L, 19L, 19L, 24L, 24L, 24L, 24L, 
41L, 41L, 41L, 41L, 42L, 42L, 42L, 42L, 30L, 30L, 30L, 30L, 43L, 
43L, 43L, 43L, 44L, 44L, 44L, 44L, 22L, 22L, 22L, 22L, 46L, 46L, 
46L, 46L, 47L, 47L, 47L, 47L, 17L, 17L, 17L, 17L, 48L, 48L, 48L, 
48L, 49L, 49L, 49L, 49L, 27L, 27L, 27L, 27L, 23L, 23L, 23L, 23L, 
50L, 50L, 50L, 50L, 51L, 51L, 51L, 51L, 52L, 52L, 52L, 52L, 41L, 
41L, 41L, 41L, 7L, 7L, 7L, 7L, 46L, 46L, 46L, 46L, 11L, 11L, 
11L, 11L, 33L, 33L, 33L, 33L, 53L, 53L, 53L, 53L, 54L, 54L, 54L, 
54L, 13L, 13L, 13L, 13L, 38L, 38L, 38L, 38L, 4L, 4L, 4L, 4L, 
37L, 37L, 37L, 37L, 55L, 55L, 55L, 55L, 57L, 57L, 57L, 57L, 44L, 
44L, 44L, 44L, 58L, 58L, 58L, 58L, 59L, 59L, 59L, 59L, 12L, 12L, 
12L, 12L, 47L, 47L, 47L, 47L, 48L, 48L, 48L, 48L, 60L, 60L, 60L, 
60L, 21L, 21L, 21L, 21L, 18L, 18L, 18L, 18L, 28L, 28L, 28L, 28L, 
26L, 26L, 26L, 26L, 61L, 61L, 61L, 61L, 31L, 31L, 31L, 31L, 59L, 
59L, 59L, 59L, 52L, 52L, 52L, 52L, 29L, 29L, 29L, 29L, 62L, 62L, 
62L, 62L, 63L, 63L, 63L, 63L, 54L, 54L, 54L, 54L, 55L, 55L, 55L, 
55L, 53L, 53L, 53L, 53L, 51L, 51L, 51L, 51L, 50L, 50L, 50L, 50L, 
64L, 64L, 64L, 64L, 20L, 20L, 20L, 20L, 58L, 58L, 58L, 58L, 16L, 
16L, 16L, 16L, 57L, 57L, 57L, 57L, 14L, 14L, 14L, 14L, 63L, 63L,

Related Posts:

  1. R Hex to RGB converter
  2. How do I select the first row in an R data frame that meets certain criteria?
  3. Error in file(file, “rt”) : cannot open the connection [duplicate]
  4. Error in : object of type ‘closure’ is not subsettable
  5. How to join (merge) data frames (inner, outer, left, right)
  6. Interpreting “condition has length > 1” warning from `if` function
  7. rbind error: “names do not match previous names”
  8. rbind error: “names do not match previous names”
  9. How to coerce a list object to type ‘double’
  10. $ operator is invalid for atomic vectors for dataframe R
  11. Plotting with ggplot2: “Error: Discrete value supplied to continuous scale” on categorical y-axis
  12. ggplot2 error : Discrete value supplied to continuous scale
  13. ggplot2 line chart gives “geom_path: Each group consist of only one observation. Do you need to adjust the group aesthetic?”
  14. R programming: How do I get Euler’s number?
  15. adding x and y axis labels in ggplot2
  16. Set NA to 0 in R
  17. ‘x’ and ‘y’ lengths differ ERROR when plotting
  18. What are the “standard unambiguous date” formats for string-to-date conversion in R?
  19. ‘x’ and ‘y’ lengths differ ERROR when plotting
  20. What are the “standard unambiguous date” formats for string-to-date conversion in R?
  21. Difference between paste() and paste0()
  22. Remove duplicated rows
  23. Error in colMeans(x, na.rm = TRUE) : ‘x’ must be numeric in KNN classification
  24. Error in ggplot.data.frame : Mapping should be created with aes or aes_string
  25. Error in plot.window(…) : need finite ‘xlim’ values
  26. Editing legend (text) labels in ggplot
  27. case_when in mutate pipe
  28. Logistic Regression on factor: Error in eval(family$initialize) : y values must be 0 <= y <= 1
  29. dim(X) must have a positive length when applying function in data frame
  30. Extract year from date
  31. Error: attempt to use zero-length variable name
  32. incorrect number of dimensions and incorrect number of subscripts in array
  33. How to debug “contrasts can be applied only to factors with 2 or more levels” error?
  34. R error: cannot coerce type ‘closure’ to vector of type ‘double’
  35. Error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, …) : NA/NaN/Inf in ‘y’, tried every possible way
  36. R t-test Grouping factor must have exactly 2 levels error
  37. What does c do in R?
  38. Understanding the order() function
  39. Error in lis[[i]] : attempt to select less than one element
  40. Plot a line graph, error in xy.coords(x, y, xlabel, ylabel, log) : ‘x’ and ‘y’ lengths differ
  41. How to compute weighted mean in R?
  42. Calculate correlation for more than two variables?
  43. How to count the number of observations in R like Stata command count
  44. Remove all of x axis labels in ggplot
  45. Sum rows in data.frame or matrix
  46. Create an empty data.frame
  47. Can’t draw Histogram, ‘x’ must be numeric
  48. MLE error in R: initial value in ‘vmmin’ is not finite
  49. Rename multiple columns by names
  50. What is the difference between rm() and rm(list=ls())?
  51. Remove rows with all or some NAs (missing values) in data.frame
  52. Convert date-time string to class Date
  53. How do you create vectors with specific intervals in R?
  54. Linear model function lm() error: NA/NaN/Inf in foreign function call (arg 1)
  55. Remove rows with all or some NAs (missing values) in data.frame
  56. Converting data frame column from character to numeric
  57. Singularity in backsolve at level 0, block 1 in LME model
  58. Nested ifelse statement
  59. How do I get RSS from a linear model output
  60. Warning message: In `…` : invalid factor level, NA generated
  61. Plot two graphs in same plot in R
  62. Mean of a column in a data frame, given the column’s name
  63. How to overcome error:”attempt to set ‘colnames’ on an object with less than two dimension” in xts object
  64. remove all variables except functions
  65. How to update a package in R?
  66. Error: Invalid input: date_trans works with objects of class Date only
  67. Run R script from command line
  68. Matrix expression causes error “requires numeric/complex matrix/vector arguments”?
  69. NAs are not allowed in subscripted assignments
  70. What’s the difference between facet_wrap() and facet_grid() in ggplot2?
  71. ggplot2, facet_grid, free scales?
  72. How do I install an R package from source?
  73. How can I convert Json to data frame in R
  74. lmer error: grouping factor must be < number of observations
  75. R Error (from NA’s to 0): duplicate subscripts for column in Data Frame
  76. Error in rbind(deparse.level, …) : numbers of columns of arguments do not match R
  77. Creating box plot on exercise
  78. What is the meaning of the dollar sign “$” in R function()?
  79. How to sort a data frame by alphabetic order of a character variable in R?
  80. Error in na.fail.default(as.ts(x)) : missing values in object in time series forecasting
  81. Argument “No” is missing, with no default
  82. R apply function with multiple parameters
  83. Geometric Mean: is there a built-in?
  84. Error in eval(predvars, data, env) : invalid ‘envir’ argument of type ‘character’
  85. Error in na.fail.default: missing values in object – but no missing values
  86. Error: Invalid number of ‘breaks’ in R
  87. Efficiently sum across multiple columns in R
  88. Convert the values in a column into row names in an existing data frame
  89. Deleting rows that are duplicated in one column based on the conditions of another column
  90. Replace all 0 values to NA
  91. The difference between bracket [ ] and double bracket [[ ]] for accessing the elements of a list or dataframe
  92. Merge error : negative length vectors are not allowed
  93. non-numeric argument to binary operator
  94. Error Error in storage.mode(x) <- "double" : 'list' object cannot be coerced to type 'double'
  95. Column standard deviation R [duplicate]
  96. Meaning of objects being masked by the global environment
  97. Subset and ggplot2
  98. Create a matrix of scatterplots (pairs() equivalent) in ggplot2
  99. How do I print a hexadecimal number with leading 0 to have width 2 using sprintf?
  100. In R, dealing with Error: ggplot2 doesn’t know how to deal with data of class numeric

Leave a Comment