我将使用R中的Hmisc包在估算缺失值后生成置信区间。 使用Hmisc进行插补过程后,是否有任何特定的代码来生成置信区间? 还是t.test()函数适合用于生成置信区间? 这是我当前的代码:
library(Hmisc)
PM <-c(1,4,10,22,15,34,35,55,100,102,103,110,115,116,118,130,135,150,145,200)
data(faithful)
fData <- faithful
XNanData = DataFull = fData
pointsMissing = PM
missIdx = length(pointsMissing)
XNanData[pointsMissing,2] = NA
imputePmm <- aregImpute(~V2 + V1 , data=XNanData, type='pmm',pmmtype=1)