R中相类型密度的对数变换

问题描述

我正在处理相类型分布，但我似乎无法正确获得密度的对数变换。

library(matrixdist)
data <- rexp(500,1.5)
paretoData <- exp(data)-1
hist1 <- hist(paretoData,breaks=30,probability =T)
y <- seq(0.01,116,by=0.1)
ph5 <- ph(structure="general",dimension=5)
fit5 <- fit(ph5,paretoData,stepsEM=1)
fit5 <- fit(fit5,stepsEM=10000,every=5)
lines(y,dens(fit5,y)$dens,col = "red")

从中我得到以下情节

然后我想记录变换，以便更好地了解发生的情况，但我似乎无法让它对密度函数起作用。我使用以下直方图

hist1 <- hist(log(paretoData),probability =T)

但我不知道如何处理密度才能让它发挥作用。

如果我应该提供更多信息，请告诉我

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）