ThierryO

library(ggplot2)
library(lme4)
library(mvtnorm)

fixed.intercept <- 1
n.subject <- 10
n.replicate <- 10
mean.subject <- c(1, 2, 3, 3, 2, 1)
sd.subject <- 0.1 * diag(length(mean.subject))
sd.noise <- 0.5

ThierryO

library(sp)
library(plyr)
library(reshape2)
library(rgdal)

n <- 1e4
dataset <- data.frame(
  ID = seq_len(n),
  X = runif(n),
  Y = runif(n),

ThierryO

abuff <- function(spgeom, ratio = 1.1, precision = 0.00001, maxsteps = 20){
    require(rgeos)
    spgeom.area <- gArea(spgeom)
    spgeom.perimeter <- gLength(spgeom)
    buffer.width <- (ratio - 1) * spgeom.area / spgeom.perimeter
    buffered.spgeom <- gBuffer(spgeom, width = buffer.width)
    achieved.precision <- gArea(buffered.spgeom) / spgeom.area - ratio
    
    steps <- 1