2

我有一个带有变量和两个因素的数据框

test <- data.frame(variable = rnorm(100,10,5),
                   factor_1 = as.factor( rep(1:2,50)),
                   factor_2 = as.factor( rep(1:5,20)))

test <- test[test$variable>0,]

我将 Weibull 分布拟合到全因子子集并收集拟合参数，如下所示：

DF <- split(test, list(test$factor_1, test$factor_2), drop = TRUE)

library(fitdistrplus)
fit.weibull <- function(DF) {
   scale <- summary(fitdist(DF[,1],'weibull'))[[1]][2]
   scale_sd <- summary(fitdist(DF[,1],'weibull'))[[3]][2]
   rbind(scale, scale_sd)
}
params <- lapply(DF, fit.weibull)

我本质上想要作为最终输出的是一个数据框，其中包括scale,scale_sd和factor_2列。但是我目前的输出很难转换成这样的数据框。

任何帮助深表感谢。

4

1 回答 1

0

好的，我找到了一个有效的丑陋解决方案：

scale <- c()
scale_sd <- c()
name <- c()

for (i in 1:length(params)) {

  scale[[i]] <- params[[i]][[1]]
  scale_sd[[i]] <- params[[i]][[2]]
  name[[i]] <- names(params)[i]

}

DF <- data.frame("name"= as.character(name),
                 "scale"= as.numeric(scale),
                 "scale_sd"= as.numeric(scale_sd))

该name列实际上包含两个因子，因此数据框的每一行都代表一个完整的因子子集。

于 2019-02-20T16:28:55.587 回答