在R中找到fitdistrplus包中的残差

时间:2016-10-27 12:57:33

标签: r fitdistrplus

我有代码,我使用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对象.....)

有什么想法吗?非常感谢

0 个答案:

没有答案