1

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

有任何想法吗？非常感谢

4

0 回答 0