当采样集中的所有值都相等时，避免在R boot.ci函数中发生错误

Question

我有许多输入到函数的数据集。数据存储在数据表中，我正在计算函数输出的置信区间。但是，在某些情况下，所有输入数据都相同，从而导致错误：“ x的所有值都等于100 \ n无法计算置信区间”如何避免此错误（例如，仅设置置信度如果所有值都相等，则间隔为0或NA之类的任意值）？例如：

library(boot)
library(data.table)

problem=1

data<-data.table(column1=c(1:100),column2=c(rep(100,99),problem))
resample.number=1000
confidence=0.95

sample.mean<-function(indata,x){mean(indata[x])}

boot_obj<-lapply(data,boot,statistic = sample.mean,R = resample.number)

boot.mean.f<-function(x,column){
    x[column][1]
}

means<-data.table(sapply(boot_obj,boot.mean.f))
bootci_obj<-lapply(boot_obj,boot.ci, conf = confidence, type = "perc")
bootci.f<-function(x,column){
    x<-x[column][4]
    x<-unlist(strsplit(as.character(x[1]),","))
    x<-sub("[:punct:].*","",x)
    x<-sub("lis.*","",x)
    x<-sub(").?","",x)
    x<-na.omit(as.numeric(x))
}

cis<-data.table(t(sapply(bootci_obj,bootci.f)))
setnames(means,"V1","stat")

cis[,V1:=NULL]
cis[,V2:=NULL]
setnames(cis,c("V3","V4"),c("lci","uci"))

return(cbind(means,cis))

返回：

stat      lci       uci
1:  50.5 44.96025  56.26797
2: 99.01 97.03000 100.00000

更改

problem=1

返回：“ t的所有值都等于100 \ n无法计算置信区间”，这会导致其他错误。

我希望结果是：

stat      lci       uci
1:  50.5 44.96025  56.26797
2: 100.0 0.0000 0.00000

Answer 1

我堆叠了data.table，因为使用长格式的data.table效率更高。如果所有值都相等，我还希望将置信度限制设置为与平均值相同的值。根据需要进行调整。

library(boot)
library(data.table)

DT <- data.table(column1=1:100,column2=rep(100,100))
DT <- data.table(stack(DT))

resample.number=1000
confidence=0.95

sample.mean <- function(indata,x){mean(indata[x])}
ci.mean <- function(x, resample.number,confidence) {
  if(length(unique(x)) > 1) {
    temp <- boot.ci(boot(x,statistic = sample.mean,R = resample.number), conf = confidence, type = "perc")$percent
    list(mean=mean(x),lwr=temp[,4],upr=temp[,5])
  } else {
    list(mean=mean(x),lwr=mean(x),upr=mean(x)  
  }
}

set.seed(42)
DT[,ci.mean(values,resample.number,confidence),by=ind]

#       ind  mean       lwr       upr
#1: column1  50.5  44.92305  55.93949
#2: column2 100.0 100.00000 100.00000

请注意，如果所有值都相等，boot.ci只会发出警告并返回NA值。没有错误，并且如果您可以使用NA，则不需要if条件。