我想为R中的有限混合模型写一个最大似然估计，但模型估计的系数与包flexmix估计的系数不同。我不知道你能否指出我的错误。我的代码如下：slope1 <- -.3;slope2 <- .3;slope3 <- 1.8; slope4 <- 0.5;intercept1 <- 1.5
agefit <- optim(vi,fmm , method = "
浏览 16提问于2019-07-11得票数 1
我正在尝试理解flexmix，特别是我在尝试拟合最简单的二项式混合模型(两个仅截取的模型的混合)时做错了什么。library(flexmix)mo2=FLXMRglm(offset=rep(log(.8/.2), nrow(data)),family="
浏览 4提问于2019-10-01得票数 2
2回答
我将一些数据拟合到flexmix中两个高斯分布的混合分布mod <- flexmix(yn ~ x, data = NPreg,FLXMRglm(yn ~ x, family= "gaussian"),模型</em
浏览 3提问于2015-11-11得票数 3
1回答
我一直试图用R编写一个有限的混合跨栏模型，使用pscl的hurdle实现，我特别感兴趣的是对计数数据有负二项式的障碍。我遇到了flexmix包，不幸的是它没有为pscl的hurdle提供内置的驱动程序，尽管多次尝试，我还是无法正确地编写驱动程序。有没有人在pscl的hurdle和flexmix一起使用过，您能共享驱动程序的工作实现吗？
提前谢谢。
浏览 5提问于2015-02-20得票数 0
回答已采纳
paramHankel.scaled()是mixComp包中的一个函数，用于确定有限混合模型的分量。dat)) * sd(dat) MLE.norm.list <- list("MLE.norm.mean" = MLE.norm.mean, "MLE.norm.sd" = MLE.norm.sdnormLoc.dM)
# plot the results for b
浏览 1提问于2022-05-14得票数 1
2回答
我正在尝试使用将一些数据拟合到混合模型中。在Matlab中，代码如下mixtureModel = ...pguess/180 + (1-pguess)*exp(kappa*cos(2*x/180*pi))/(180*besseli(0,kappa));
options.MaxFunEvals = 20000;
% fit the model using maximum li
浏览 2提问于2014-10-20得票数 2
1回答
我试图使用来自包的函数来拟合一个动态混合模型(散装的Weibull，尾部的Pareto )：library(evmix)hist(data,probability=TRUE)
lines(xeval
浏览 0提问于2015-05-26得票数 0
回答已采纳
我有一个mle2模型，我在这里开发就是为了演示这个问题。我从两个独立的高斯分布x1和x2生成值，将它们组合在一起形成x=c(x1,x2)，然后创建一个MLE，它试图通过xsplit参数将x值重新分类为属于特定x值的左边或特定x值的右侧。问题是所发现的参数并不理想。具体来说，xsplit总是以其起始值的形式返回。如果我改变它的起始值(例如，4或9)，结果的日志可能性就会有很大的差异。mo =
浏览 1提问于2014-02-07得票数 4
回答已采纳
1回答
我想比较以下模型的估计值，但我不知道是否可以用plm软件包来估计模型(5)：fit.pooling <- plm(form, data = pdata, model = "pooling我能用plm软件包来估计那些模型吗？
浏览 1提问于2018-02-23得票数 0
回答已采纳
如何在R中使用多项式响应变量(使用柔性混合包)运行潜在类增长模型(LCGM)？
以及如何用二进制/分类因变量对每个类进行分层？FLXMRmultinom在挠性混合中的二项式等效值可能是FLXMRglm (与family=“二项式”)相反的FLXMRglmfix，因此我怀疑LCGM的限制(例如。这些结果在其他方面是合理的，但模型未能将轨迹相似的男性和女性置于相同的类别中(以下是按性别分列的每一种关
浏览 13提问于2022-11-01得票数 0
1回答
我想在这样的模型中找到参数epsilon和mu的MLE：library(Renext{ -sum(R)fit_norm <- mle2
浏览 4提问于2021-01-24得票数 0
回答已采纳
我正在用nlmer从lme4库测试一些非线性模型。但是这个函数假设正常，而我的数据显然遵循Gamma分布。(b0|id), start = start.df, REML=T)在数据不是高斯的情况下，是否有添加家庭组的方法，如lme4，或者在非线性模型中测试组的任何其他提示？
浏览 13提问于2022-01-10得票数 0
1回答
我想从MLE/渐变下降的角度来理解正规化/收缩。我知道这两个概念，但我不知道/不明白两者是否都用于确定线性模型的系数。如果是的话，所遵循的步骤是什么？为了进一步阐明，采用正则化方法来降低方差，这是通过线性模型的惩罚系数来实现的。调整参数lambda是通过交叉验证确定的。一旦，lambda被确定，系数就会自动确定，对吗？因此，为什么我们需要最小化(RSS +正则化项)来找到系数？这些步骤如
浏览 0提问于2015-09-02得票数 5
回答已采纳
1回答
我想用R中的optim()包找到估计参数，并将我的结果与R中的GLM模型进行比较。-param[1]+param[2]*x1+param[3]*x2 -sum(y*log(p)+(1-y)*log(1-p),na.rm=TRUE)MLE_estimates<-optim(c(Intercept=0.1,age=0.1,sector2=0.1),nlldbin,hessian=TRUE)
MLE_estimates
浏览 18提问于2021-04-04得票数 1
我正在使用API中的Backbone Collections和Model。我的id属性是mongo (_id)的默认id属性，但在我的应用编程接口中，我使用一个不同的唯一键来获取单个模型。像鼻塞一样的东西总是唯一的。现在的问题是，当我想要保存()时，它执行PUT请求而不是POST请求，因为它假设已经有一个具有该键( slug )的模型，因为用户能够自己输入slug(当然要经过适当的
浏览 2提问于2015-04-01得票数 0
1回答
然后给出一个移动和桌面混合的例子：所以.col-md-8和.col-md-4相加为12，因此它们应该保持在一个水平块上。但是，.col-xs-12和.col-xs-6的加起来并不等于12，我在这里的假设是，我们混合了xs和md，因为我们希望它们都出现在各自设备的一个水平块上。如果没有，为什么？
浏览 1提问于2015-05-04得票数 0
回答已采纳
1回答
我想写R代码来建立dirichlet混合模型。用于beta发行版的对数可能性如下所示： }
MLE_estimates# custom Inputs
浏览 6提问于2022-11-23得票数 0
1回答
我正在从学习混合模型，这个问题是关于模块5中的一个练习。nullmodel <- lmer(score ~ (1|schoolid), data = mydata, REML = FALSE)我们通过运行nullmodel的摘要得到
浏览 1提问于2017-12-02得票数 0
1回答
当我在Prophet模型上调用predict时，平均预测是yhat，但当我将平均趋势和季节性列相加时，这并不等于yhat中的值。为什么列中的术语不能与yhat、yhat_lower和yhat_upper相加
浏览 21提问于2019-06-21得票数 1
云服务器
ICP备案
云直播
对象存储
实时音视频
腾讯云开发者
