在R中绘制对数似然函数
齐格弗里德
嗨，我收到的错误x和y长度在我的代码中不同。datasim是大小为1000的模拟样本。请帮助我。
x <- datasim
loglik <- function(theta){
k<- theta[1]
lambda<- theta[2]
out <- sum(dweibull(x,shape = k, scale=lambda, log = TRUE) )
return(out)
}
theta<- c(0.5,1.5)
plot(theta, loglik(theta), type="l", lwd=3, main="logliklihood_Weibull, n=1000")
笨狼
您可以使用一个包装，因为它很适合您：
library(fitdistrplus)
x = rweibull(1000,20,10)
fit <- fitdist(x, "weibull")
llplot(fitg, expand=5)
或使用您的loglik函数：
library(lattice)
da = expand.grid(k=seq(8,22,length.out=50),
lambda = seq(9,12,length.out=50))
da$LL = apply(da,1,loglik)
wireframe(LL ~ k * lambda, data = da,
scales = list(arrows = FALSE),drape = TRUE, colorkey = TRUE)
本文收集自互联网，转载请注明来源。
如有侵权，请联系[email protected] 删除。
编辑于
- 上一篇：AWS System Manager（SSM）在InvalidPluginName中获取命令调用结果
- 下一篇：将Firebase CLI更新到8.4.0会出现有关“不受支持的引擎”的错误，并显示“ {“ node”：“> = 10”}'
相关文章
Related 相关文章
- 1
用ggplot2绘制对数似然函数的问题
- 2
如何在R中编码多参数对数似然函数
- 3
如何在R中编写多参数对数似然函数
- 4
对数似然成本函数：均值或总和？
- 5
对数似然成本函数：均值或总和？
- 6
R：对数似然优化产生的NaN
- 7
R从plm对象提取对数似然
- 8
statsmodels与pymc中的对数似然
- 9
对数似然在Delphi中的实现
- 10
R中的最大似然
- 11
对数似然性：R中的NA Flexmix软件包
- 12
GDA的对数似然函数（高斯判别分析）
- 13
anova（）函数返回的对数似然的值是多少
- 14
Mathematica中的对数似然和多项式分布
- 15
在R中实现最大似然
- 16
具有“似然”方法的R调查包函数Svyciprop
- 17
马尔可夫网络的对数似然
- 18
从概率后缀树获取对数似然
- 19
scipy.stats.rv_continuous.fit生成的对数似然函数
- 20
lmer（来自R包lme4）如何计算对数似然率？
- 21
使用R计算从两个正态分布的混合物中采样的一组观测值的对数似然
- 22
如何在sklearn GMM中获得每次迭代的对数似然性？
- 23
Python中的似然比测试
- 24
R中的对数函数
- 25
主题模型：对数似然或困惑的交叉验证
- 26
Statsmodels OLS回归：对数似然，用法和解释
- 27
主题模型：对数似然或困惑的交叉验证
- 28
计算对数似然（MATLAB）时避免使用-inf
- 29
为 WAIC（逻辑层次斯坦模型）编写对数似然
我来说两句