R中的分布拟合

Question

我想分配一个分布。 如果我有一个数据集，我可以很容易做到：

library("fitdistrplus")
data_raw <- c(1018259, 1191258, 1265953, 1278234, 1630327, 1780896, 1831466, 1850446, 1859801, 1928695, 2839345, 2918672, 3058274, 3303089, 3392047, 3581341, 4189346, 5966833, 11451508)
fitdist(data_raw, "lnorm")

这就是我要使对数正态分布适合数据集的方法。
但是，如果我没有一个数据集，只是平均值，标准差和一些分位数，该怎么办？ 例如：

均值：2965042
标准设备：2338555

位数：
0.1：1251014
0.5：1928695
0.8：3467765
0.9：4544483
0.95：6515300
0.999：11352784

您将如何继续对此类数据进行估算？

感谢你并致以真诚的问候
Norbi

Answer 1

只需使用nls拟合模型：

DF <- read.table(text = "0.1: 1251014
                 0.5: 1928695
                 0.8: 3467765
                 0.9: 4544843
                 0.95: 6515300
                 0.999: 11352784 ", sep = ":")
plot(V1 ~ V2, data = DF, 
     xlim = c(0, 1.2e7),ylim = c(0, 1))


fit <- nls(V1 ~ plnorm(V2, meanlog, sdlog), data = DF, 
           start = list(meanlog = 13, sdlog = 2), trace = TRUE, algorithm = "port",
           lower = c(0, 0))

summary(fit)

curve(plnorm(x, coef(fit)[[1]], coef(fit)[[2]]), add = TRUE, col = "blue")