Паттерны, возникающие из случайной выборки

Question

Я пытаюсь загрузить некоторые данные о радио-ошейнике, чтобы выявить внезапное, значительное увеличение в домашних диапазонах. Я взял данные ошейника и рандомизировал порядок точек, используя sample(), а затем создал 95% MCP на каждые 15 совокупных точек (используя mcp.area() из пакета adehabitatHR) и сохранил их в data.frame. На нормальном животном это будет область MCP для первых собранных 15 точек, 30 точек, 45 точек и т. Д. Я повторил процесс 1000 раз, поэтому у меня было 1000 поддельных животных (рядов) из рандомизированных реальных точек, и у каждого есть список областей (столбцы).

Затем я создал аналогичный data.frame с процентным изменением от одного расчета площади к следующему, то есть с одного временного шага на следующий, и взял среднее изменение для всех итераций на каждом временном шаге. Это то, что я хотел использовать, чтобы сравнить с реальными данными, чтобы я мог определить сдвиг в домашнем диапазоне, но когда я его получил, я получал странный всплеск изменений при каждом 4-м шаге. Вот график среднего изменения площади за шаг по времени.

Вот как выглядит мой код:

cat1 <- read.csv("cat1.csv")
  coordns <- data.frame(x=cat1$Longitude,y=cat1$Latitude)
  coordns <- coordns[sample(1:nrow(coordns),replace=FALSE, prob = NULL),]    
  coordns <- SpatialPoints(coordns,
                           proj4string = CRS("+proj=longlat +datum=WGS84"))    
  coordns <- spTransform(coordns,CRSobj = "+proj=utm +zone=55 +datum=WGS84")    
  e <- data.frame()  
    t <- seq(15, nrow(cat1), by = 15)
while(nrow(e)<1000){
  ar <- data.frame()
  for(i in t){
    ar <- rbind(ar,data.frame(mcp.area(coordns[1:i],percent=100,unout = "ha", plotit = FALSE)))
  } 
    coordns <- data.frame(x=cat1$Longitude,y=cat1$Latitude)
    coordns <- coordns[sample(1:nrow(coordns),replace=FALSE,prob=NULL),]
    coordns <- SpatialPoints(coordns,
                             proj4string = CRS("+proj=longlat +datum=WGS84"))
    coordns <- spTransform(coordns,CRSobj = "+proj=utm +zone=55 +datum=WGS84")
    ar <- t(data.frame(area=ar$a))
e <- rbind(e,data.frame(ar))
}
f <- (e[,2:95]/e[,1:94])

Я запускал его несколько раз с разными настройками MCP и длительностью временных шагов, но всегда есть шаблон. Это не всегда каждый четвертый или как произносится, но это есть.

У кого-нибудь есть идеи, почему это будет происходить?