R中的Bootstrap-错误:下标超出范围

时间:2019-02-01 04:41:51

标签: r overlap bootstrapping

我正在使用代码进行引导以发现两个物种之间活动模式的差异。我收到一个错误消息,说“下标超出范围”。我在下面粘贴了我正在使用的代码。关于我哪里出错了的任何想法?

#Load packages
library(data.table)
library(activity)
library(overlap)

#Load dataset 
<-fread(file="data.csv")
names(data)  #these are the column names
unique(data$Subproject)  #these are the counties in the dataset    

##### Turn the time into Radians 
setDT(data)[,paste0("Begin.Time",1:2):= tstrsplit(Begin," ")] #Splits the time stamp into Date and Time Columns
setnames(data, old = c('Begin.Time1','Begin.Time2'), new = c('Date','Time')) #Renames the new columns 

setDT(data)[,paste0("Times",1:3):= tstrsplit(Time,":")]
data$hours<-as.numeric(data$Times1)
data$mins<-as.numeric(data$Times2)
data$seconds<-as.numeric(data$Times3)

data$hours<-data$hours * 60
data$seconds<-data$seconds / 60
data$totalminutes<-data$hours + data$mins + data$seconds

data[,!c("Times1","Times2","Times3"), with=F]
data$totalminutes<-apply(data, MARGIN=2, FUN=function(x) (data$totalminutes-min(data$totalminutes))/diff(range(data$totalminutes)))

summary(data$totalminutes)

data$radians <- data$totalminutes * 2 * pi  #This turns total minutes into radians

######################## Bootstrapping ###################################
Deer.bootstrap <- resample (Deer$radians, 1000)
dim(Deer.bootstrap)
coyote.bootstrap <- resample (Coyote$radians, 1000)
dim(coyote.bootstrap)

coyote.Deer.boot <- bootEst(coyote.bootstrap, Deer.bootstrap, adjust = c (NA,1,NA))
dim(coyote.Deer.boot)
Boot.mean <- colMeans (coyote.Deer.boot)
Boot.mean

Boot.Bias <- Boot.mean -D.hat.estimates 
Boot.Bias    

####################### confidence intervals #############################
#use bootCI to calculate confidence intervals for D.hat / boot data
Conf.deer.coyote <- coyote.Deer.boot [, 2]
bootCI(D.hat.estimates[2], Conf.deer.coyote, conf = 0.95)

######## Calculating difference between 2 species in relation to time #####
#to compare activity across certain desired parts of the day per species: reps 
# represents bootstrapping amount
fdeer <- fitact (Deer$radians, reps = 500)
compareTimes(fdeer, c(5.5,6, 0.5, 1))

这是我从最后一行代码中得到的错误(其他所有方法都正常):

Error in fit@pdf[, 3] : subscript out of bounds

0 个答案:

没有答案