fitdistr 函数如何使用？

(1)不同的初始值可能会得到不同的极大值，实际中应该多试几次，能稳定最好。
(2)

1. uui$estimate[1]和 uui$estimate[2]

复制代码

1. uui$estimate[1]和uui$estimate[2]

复制代码

maths_hjxk 发表于 2020-6-20 12:56
我想起来你说的CASS包，其实是MASS包，里面有个函数fitdistr一样可以估计
start=list(shape1 =4, shape2  ...

你好！
c<-c(0.1 0.1 0.2 0.2 0.2 0.2 0.3 0.6)
uui<- fitdist(c, distr="beta", method = "mle",start=NULL, fix.arg=NULL)
方便你用，
我用的例子在上面，期待你的见解！
多谢！

maths_hjxk 发表于 2020-6-20 12:56
我想起来你说的CASS包，其实是MASS包，里面有个函数fitdistr一样可以估计
start=list(shape1 =4, shape2  ...

你好！
c<-c(0.1 0.1 0.2 0.2 0.2 0.2 0.3 0.6)
uui<- fitdist(c, distr="beta", method = "mle",start=NULL, fix.arg=NULL)
方便你用，
我用的例子在上面，期待你的见解！
多谢！

maths_hjxk 发表于 2020-6-20 13:08

你好！ 我尝试了
> uui$estimate[1]
  shape1
2.265435
还是不能取出 2.265 这个数字，
多谢！

其实beta密度是单峰函数，理论上无论什么初值始结果都是一样的。

maths_hjxk 发表于 2020-6-20 13:08

我也试了
> uui$estimate[1]
  shape1
2.265435


> uui$estimate[1][1]
  shape1
2.265435

都是不能对别的变量赋值 ，因为都有“shape1”字样

maths_hjxk 发表于 2020-6-20 13:15
你似乎应该先入门一些R的基本常识

多谢指点！
> c<-c(0.1,0.1,0.2,0.2,0.2,0.2,0.3,0.6)
> uui<- fitdist(c, distr="beta", method = "mle",start=NULL, fix.arg=NULL)
> z<-c()
> z[1]<-uui$estimate[1]
> z
[1] 2.265435
>
也就是说虽然uui$estimate[1]  带有shape1字样，但是其实还是一个数字是吗？多谢！

江夏雁 发表于 2020-6-20 13:19
多谢指点！
> c uui z z[1] z
[1] 2.265435

是的，一样可以进行计算。