1

我有一个包含数千个独特测试 (TestNum) 及其相关响应 (Response) 的大型数据集，响应的长度在测试中有所不同。测试基于某些标准被丢弃，因此 TestNum 值缺乏序列。这里有一个简化的例子：

dat <- data.frame(Response=c(rlnorm(10, 2.9, 0.3), rlnorm(14, 2.88, 0.38), rlnorm(19, 2.44, 0.08)),TestNum=rep(c(1,4,9), times=c(10,14,19)))
dat$TestNum<-factor(dat$TestNum)
dat

我正在为每个 TestNum 拟合 lnorm 分布并提取系数

dat_fit1 <- with(dat,
        by(dat[,1], TestNum, fitdist, "lnorm"))
dat_fit2 <-t(sapply(dat_fit1, coef))

我想测试其他分布，但需要每个分布的拟合优度统计（gofstat；例如“chi”chisqpvalue”、“cvm”、“ad”、“ks”、“aic”、“bic”） TestNum 的拟合曲线。我可以用下面的代码得到“aic”和“bic”，但不能得到其余的统计数据。

gof_f<-do.call(rbind, dat_fit1) 
gof_f<-gof_f[,7:8]

任何建议，将不胜感激！

4

1 回答 1

0 
dat.lnorm <- with(dat,
              by(dat[,1], dat[,2], 
        function(x){
        fit<-fitdist(x,"lnorm", method="mme")                                   
        coef_meanlog <-fit[[1]][[1]]
        coef_sdlog   <-fit[[1]][[2]]
        ks <-ks.test(jitter(x),"plnorm", meanlog=coef_meanlog, sdlog=coef_sdlog)$p.value
        ad <-ad.test(plnorm(x, meanlog=coef_meanlog, sdlog=coef_sdlog))$p.value
        return(list(cbind(rbind(fit[7:8]), ks, ad)))
}))
于 2016-10-13T15:53:09.753 回答