Questions tagged «marginal»

边际分布是指联合分布中包含的变量子集的概率分布。

2
随机效应模型,固定效应模型和边际模型之间有什么区别?
我正在尝试扩展我的统计知识。我来自物理科学背景,采用“基于配方”的方法进行统计测试,我们说它是连续的,是否呈正态分布-OLS回归。 在阅读中,我遇到了以下术语:随机效应模型,固定效应模型,边际模型。我的问题是: 简单来说,它们是什么? 它们之间有什么区别? 他们有同义词吗? 传统测试(例如OLS回归,ANOVA和ANCOVA)在哪里分类? 只是尝试决定自学的下一步。

5
关于Copulas的入门阅读
一段时间以来,我一直在为我的研讨会寻找有关Copulas的良好介绍性阅读。我发现有很多关于理论方面的材料,这是很好的,但是在我将其介绍之前,我希望对这一主题建立良好的直观理解。 谁能提出建议为初学者打好基础的好论文(我在合理的程度上开设了1-2门统计学课程,并了解边际,多元分布,逆变换等)?

1
lmer模型使用哪种多重比较方法:lsmeans或glht?
我正在使用具有一个固定效果(条件)和两个随机效果(由于主题设计和配对而导致的参与者)的混合效果模型分析数据集。该模型是使用lme4包生成的exp.model<-lmer(outcome~condition+(1|participant)+(1|pair),data=exp)。 接下来,我针对没有固定效果(条件)的模型对该模型进行了似然比检验,结果有显着差异。我的数据集中有3个条件,因此我想进行多重比较,但不确定使用哪种方法。我在CrossValidated和其他论坛上发现了许多类似的问题,但我仍然很困惑。 据我所见,人们建议使用 1.该lsmeans包- lsmeans(exp.model,pairwise~condition)这给了我下面的输出: condition lsmean SE df lower.CL upper.CL Condition1 0.6538060 0.03272705 47.98 0.5880030 0.7196089 Condition2 0.7027413 0.03272705 47.98 0.6369384 0.7685443 Condition3 0.7580522 0.03272705 47.98 0.6922493 0.8238552 Confidence level used: 0.95 $contrasts contrast estimate SE df t.ratio p.value Condition1 - Condition2 -0.04893538 0.03813262 62.07 -1.283 0.4099 Condition1 - …

3
仅给出边际计数的联合分布的最大似然估计
令是两个类别变量的联合分布,其中。说从该分布中抽取了样本,但仅给出了边际计数,即: X ,ÿ X ,ÿ ∈ { 1 ,... ,ķ } Ñ Ĵ = 1 ,... ,ķpx,ypx,yp_{x,y}X,YX,YX,Yx,y∈{1,…,K}x,y∈{1,…,K}x,y\in\{1,\ldots,K\}nnnj=1,…,Kj=1,…,Kj=1,\ldots,K Sj=∑i=1nδ(Xi=l),Tj=∑i=1nδ(Yi=j),Sj=∑i=1nδ(Xi=l),Tj=∑i=1nδ(Yi=j), S_j = \sum_{i=1}^{n}{\delta(X_i=l)}, T_j = \sum_{i=1}^{n}{\delta(Y_i=j)}, 给定,的最大似然估计是?这是已知的吗?计算上可行吗?除了机器学习之外,还有其他合理的方法来解决这个问题吗?小号Ĵ,Ť Ĵpx,ypx,yp_{x,y}Sj,TjSj,TjS_j,T_j


1
使用条件分布从边际分布采样?
我想从单变量密度采样,但我只知道这种关系:FXfXf_X FX(x )= ∫FX| ÿ(x | y)fÿ(y)dÿ。fX(x)=∫fX|Y(x|y)fY(y)dy.f_X(x) = \int f_{X\vert Y}(x\vert y)f_Y(y) dy. 我想避免使用MCMC(直接在整数表示上),并且由于和易于采样,因此我在考虑使用以下采样器:f Y(y )FX| ÿ(x | y)fX|Y(x|y)f_{X\vert Y}(x\vert y)Fÿ(y)Fÿ(ÿ)f_Y(y) 对于。j = 1 ,… ,NĴ=1个,…,ñj=1,\dots, N 样本。ÿĴ〜˚FÿÿĴ〜Fÿy_j \sim f_Y 样本。XĴ〜˚FX| ÿ(⋅ | yĴ)XĴ〜FX|ÿ(⋅|ÿĴ)x_j \sim f_{X\vert Y}(\cdot\vert y_j) 然后,我将得到对,仅获取边际样本。 它是否正确?(x 1,… ,x N)(x1个,ÿ1个),。。。,(xñ,ÿñ)(X1个,ÿ1个),。。。,(Xñ,ÿñ)(x_1,y_1),...,(x_N,y_N)(x1个,… ,xñ)(X1个,…,Xñ)(x_1,\dots,x_N)

1
如何从具有多变量相关性的联合分布中找到边际分布?
我的教科书中的一个问题如下。二维随机连续向量具有以下密度函数: fX,Y(x,y)={15xy20if 0 &lt; x &lt; 1 and 0 &lt; y &lt; xotherwisefX,Y(x,y)={15xy2if 0 &lt; x &lt; 1 and 0 &lt; y &lt; x0otherwise f_{X,Y}(x,y)= \begin{cases} 15xy^2 & \text{if 0 < x < 1 and 0 < y < x}\\ 0 & \text{otherwise}\\ \end{cases} 证明边际密度函数和为:fXfXf_XfYfYf_Y fX(x)={5x40if 0 &lt; x &lt; …

1
更新贝叶斯因子
在贝叶斯假设检验和贝叶斯模型选择中,贝叶斯因子是通过两个边际可能性的比率来定义的:给定iid样本以及各自的采样密度和,具有相应的先验和,用于比较两个模型的贝叶斯因子为 一书,我目前正在审查有奇怪的声明,上面的贝叶斯因子(x1,…,xn)(x1,…,xn)(x_1,\ldots,x_n)f1(x|θ)f1(x|θ)f_1(x|\theta)f2(x|η)f2(x|η)f_2(x|\eta)π1π1\pi_1π2π2\pi_2B12(x1,…,xn)=defm1(x1,…,xn)m2(x1,…,xn)=def∫∏ni=1f1(xi|θ)π1(dθ)∫∏ni=1f2(xi|η)π2(dη)B12(x1,…,xn)=defm1(x1,…,xn)m2(x1,…,xn)=def∫∏i=1nf1(xi|θ)π1(dθ)∫∏i=1nf2(xi|η)π2(dη)\mathfrak{B}_{12}(x_1,\ldots,x_n)\stackrel{\text{def}}{=}\frac{m_1(x_1,\ldots,x_n)}{m_2(x_1,\ldots,x_n)}\stackrel{\text{def}}{=}\frac{\int \prod_{i=1}^n f_1(x_i|\theta)\pi_1(\text{d}\theta)}{\int \prod_{i=1}^n f_2(x_i|\eta)\pi_2(\text{d}\eta)}B12(x1,…,xn)B12(x1,…,xn)\mathfrak{B}_{12}(x_1,\ldots,x_n)是“通过将各个[贝叶斯因子]相乘而形成的”(第118页)。如果使用分解 但我看不到此分解的计算优势,因为需要与的原始计算相同的计算量B12(x1,…,xn)=m1(x1,…,xn)m2(x1,…,xn)=m1(xn|x1,…,xn−1)m2(xn|x1,…,xn−1)×m1(xn−1|xn−2,…,x1)m2(xn−1|xn−2,…,x1)×⋯⋯×m1(x1)m2(x1)B12(x1,…,xn)=m1(x1,…,xn)m2(x1,…,xn)=m1(xn|x1,…,xn−1)m2(xn|x1,…,xn−1)×m1(xn−1|xn−2,…,x1)m2(xn−1|xn−2,…,x1)×⋯⋯×m1(x1)m2(x1)\begin{align*}\mathfrak{B}_{12}(x_1,\ldots,x_n)&=\frac{m_1(x_1,\ldots,x_n)}{m_2(x_1,\ldots,x_n)}\\&=\frac{m_1(x_n|x_1,\ldots,x_{n-1})}{m_2(x_n|x_1,\ldots,x_{n-1})}\times \frac{m_1(x_{n-1}|x_{n-2},\ldots,x_1)}{m_2(x_{n-1}|x_{n-2},\ldots,x_1)}\times\cdots\\&\qquad\cdots\times\frac{m_1(x_1)}{m_2(x_1)}\end{align*}m1(xn|x1,…,xn−1)m2(xn|x1,…,xn−1)m1(xn|x1,…,xn−1)m2(xn|x1,…,xn−1)\frac{m_1(x_n|x_1,\ldots,x_{n-1})}{m_2(x_n|x_1,\ldots,x_{n-1})}m1(x1,…,xn)m2(x1,…,xn)m1(x1,…,xn)m2(x1,…,xn)\frac{m_1(x_1,\ldots,x_n)}{m_2(x_1,\ldots,x_n)}外面的人造玩具的例子。 问题:是否存在将Bayes因子从B12(x1,…,xn)B12(x1,…,xn)\mathfrak{B}_{12}(x_1,\ldots,x_n)为 B12(x1,…,xn+1)B12(x1,…,xn+1)\mathfrak{B}_{12}(x_1,\ldots,x_{n+1})的通用且计算有效的方法{n + 1})不需要重新计算整个边际m1(x1,…,xn)m1(x1,…,xn)m_1(x_1,\ldots,x_n)和 m2(x1,…,xn)m2(x1,…,xn)m_2(x_1,\ldots,x_n)? 我的直觉是,除了粒子滤波器实际上确实是在估计贝叶斯因子B12(x1,…,xn)B12(x1,…,xn)\mathfrak{B}_{12}(x_1,\ldots,x_n)中进行一次,没有一个自然的方法可以回答这个问题。

2
稳健的边际可能性MCMC估计器?
我正在尝试通过蒙特卡洛方法来计算统计模型的边际可能性: f(x)=∫f(x∣θ)π(θ)dθf(x)=∫f(x∣θ)π(θ)dθf(x) = \int f(x\mid\theta) \pi(\theta)\, d\theta 可能性表现良好-平滑,对数凹入-但维数高。我已经尝试过重要性抽样,但是结果很奇怪,并且在很大程度上取决于我使用的建议。我简要地考虑了假设哈密顿量在前一个统一的基础上,进行哈密顿量计算θθ\theta并以谐波均值,直到我看到了。经验教训,谐波均值可以具有无限方差。是否存在替代MCMC估算器,该估算器几乎一样简单,但具有良好的方差?

1
寻找的边际密度
如标题所示,我正在寻找的边际密度F(x ,y)= c1 -X2-ÿ2---------√,X2+ÿ2≤1 。F(X,ÿ)=C1个-X2-ÿ2,X2+ÿ2≤1。f (x,y) = c \sqrt{1 - x^2 - y^2}, x^2 + y^2 \leq 1. 到目前为止,我发现为。我通过将转换为极坐标并在积分来弄清楚,这就是为什么我被困在边际密度部分上的原因。我知道,但是我不确定如何在没有大的混乱积分的情况下解决这个问题,我知道答案是“应该是一个很大的混乱积分。是否有可能找到,然后采用来找到CCc32个π32π\frac{3}{2 \pi}F(x ,y)F(X,ÿ)f(x,y)d[R dθd[Rdθdrd\thetaFX(x )=∫∞- ∞F(x ,y)dÿFX(X)=∫-∞∞F(X,ÿ)dÿf_x(x) = \int_{-\infty}^\infty f(x,y)dyF(x ,y)F(X,ÿ)F(x,y)dFdXdFdX\frac{dF}{dx}FX(x )FX(X)f_x(x)?这似乎是一种直观的方法,但我似乎无法在教科书中找到说明这些关系的任何内容,因此我不想做出错误的假设。
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.