1

我正在运行一些测试来尝试确定我的数据遵循什么分布。从我的数据密度来看，我认为它看起来有点像逻辑分布。然后我使用包 MASS 来估计分布的参数。但是，当我将它们绘制在一起时，虽然比正常情况要好，但逻辑仍然不是很好。有没有办法找到更好的分布？感谢您的帮助 ！ 

library(quantmod)
getSymbols("^NDX",src="yahoo", from='1997-6-01', to='2012-6-01')
daily<- allReturns(NDX) [,c('daily')]
dailySerieTemporel<-ts(data=daily)
x<-na.omit(dailySerieTemporel)

library(MASS)
(xFit<-fitdistr(x,"logistic"))
#      location        scale    
#   0.0005210570   0.0106366354 
#  (0.0002941922) (0.0001444678)
xFitEst<-coef(xFit)

plot(density(x))
set.seed(125)
lines(density(rlogis(length(x), xFitEst['location'], xFitEst['scale'])), col=3)
lines(density(rnorm(length(x), mean(x), sd(x))), col=2)
4

1 回答 1

3

这是基本的 R：默认情况下plot()会创建一个新的绘图画布，您应该使用诸如lines()添加到现有绘图之类的命令。

这适用于您的示例：

plot(density(x))
lines(density(rlogis(length(x), location = 0.0005210570,
                     scale = 0.0106366354)), col="blue")

因为它将估计的逻辑拟合以蓝色添加到您现有的情节中。

于 2012-08-06T14:33:11.847 回答