我有代码,我使用fitdistrplus将Weibull分布拟合到各种数据集,并将原始数据绘制到CDF。
##Performs Weibull fit using the Maximum GOF Estimations
zFIT<-fitdist(x,"weibull",keepdata = TRUE, method = "mge", gof ="KS") ##Changedfrom CvM
zbFIT <- bootdist(zFIT, niter= 251)
zparam<- c(coef(zFIT))
z_alpha <-zparam[2]
z_beta <-zparam[1]
##Write data output to file
zDat<-c(summary(zFIT))
zDataList <-as.vector(zDat[10])
weibull_Data <-lapply(z, write, "Weibull_Data.txt", append = TRUE)
weibull_Data <-lapply(zDataList, write, "Weibull_Data.txt", append = TRUE)
##Write the parameters to a txt file...
weibull_params <-lapply(z, write, "Weibull_Fits.txt", append = TRUE)
weibull_params <-lapply(zparam, write, "Weibull_Fits.txt", append = TRUE)
我想从此计算中访问拟合的残差,以应用+/- 5伪西格玛屏幕来识别和计算异常值的数量,因为这些将被视为外在缺陷群体。
如何获得残差的向量,列表或矩阵。 fitdistr使用resid()函数,fitdistrplus应该继承这个函数,但它似乎不能通过调用resid(zFIT)(其中zFIT是我的fitdist对象.....)
有什么想法吗?非常感谢