fitdistr 函数如何使用？

已知一个样本，并知道总体符合Beta 分布。想求beta分布的参数
查到CASS包下面有一个fitdistr 函数
但是使用时，发现有人的案例如下

res =fitdistr(rate, "beta", start=list(shape1 =4, shape2 =5))
请问这里的4跟5是什么含义，难道是beta分布的2个参数？但是我这里现在就是想求这两个参数啊

这里又要求输入是什么意思？

多谢！

江夏雁 发表于 2020-6-20 13:12
你好！ 我尝试了
> uui$estimate[1]
  shape1

你似乎应该先入门一些R的基本常识

1. unname(uui$estimate[1]);
   
2. as.numeric(uui$estimate[1]);

复制代码

不知道 你这个什么扩展包，下面这个函数可以实现利用最大似然法来估计beta分布的两个参数

1. fitdistrplus::fitdist(x,"beta")

复制代码

maths_hjxk 发表于 2020-6-20 11:42
不知道 你这个什么扩展包，下面这个函数可以实现利用最大似然法来估计beta分布的两个参数

> uui<- fitdist(c, distr="beta", method = "mle",start=NULL, fix.arg=NULL)
>
> uui
Fitting of the distribution ' beta ' by maximum likelihood
Parameters:
       estimate Std. Error
shape1 2.265435   1.062954
shape2 7.039041   3.559200



请问如果把2.265435   这个结果提取出来，并赋值给另外的变量？应该如何写？多谢！

1. uui$estimate

复制代码

其他的请查看

1. str(uui)

复制代码

我想起来你说的CASS包，其实是MASS包，里面有个函数fitdistr一样可以估计
start=list(shape1 =4, shape2 =5)是初始值，可设可不设

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

多谢大牛！
我在网上看还有 fitdistr(x,"beta",start=list(shape1=0,shape2=1))
请问这个start里参数的不同有什么意义吗？多谢！

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

多谢大牛！
我在网上看还有 fitdistr(x,"beta",start=list(shape1=0,shape2=1))
请问这个start里参数的不同有什么意义吗？多谢！

maths_hjxk 发表于 2020-6-20 12:49
其他的请查看

你好！
我试了
> uui$estimate
  shape1   shape2
2.265435 7.039041
还是不能取出数字，请问可以多指点一下吗?多谢！

maths_hjxk 发表于 2020-6-20 12:49
其他的请查看

你好！我试了一下
> uui$estimate
  shape1   shape2
2.265435 7.039041
好像还是不能取出数字？
请多指点，多谢！