如果随机变量和均值正态分布和标准差,然后定义随机变量. 问题是要证明遵循带参数的伽马分布(即卡方分布)。
我有:
但我不确定从这里去哪里,因为随机变量不被认为是独立的。我将不胜感激任何提示!
如果随机变量和均值正态分布和标准差,然后定义随机变量. 问题是要证明遵循带参数的伽马分布(即卡方分布)。
我有:
结果不正确。
作为反例,让具有 Clayton copula 的标准 Normal 边距,如https://stats.stackexchange.com/a/30205所示。如图左下角所示,生成 10,000 个独立实现的这种二元分布,产生 10,000 个实现显然不遵循分布(在卡方拟合检验中,)。图顶部的 qq 图确认边缘看起来是标准的 Normal,而右下角的 qq 图表示边缘的上尾太短了。
该结果可以在假设分布为是中心对称的:也就是说,当它同时否定两者时它是不变的和. 这包括所有双变量法线(均值, 当然)。
关键思想是,对于任何, 事件是事件的区别和. (第一个是最小值不小于,而第二个将排除最小值超过.) 这些事件又可以细分如下:
中心对称假设确保最后两个概率相消。 前两个概率由标准正态 CDF 给出, 产生
那个展品作为半正态分布,因此它的平方将具有与标准正态的平方相同的分布,根据定义,它是分配。
这个演示可以倒过来显示有个当且仅当分布对所有人.
这是R产生这些数字的代码。
library(copula)
n <- 1e4
set.seed(17)
xy <- qnorm(rCopula(n, claytonCopula(1)))
colnames(xy) <- c("X", "Y")
z2 <- pmin(xy[,1], xy[,2])^2
cutpoints <- c(0:10, Inf)
z2.obs <- table(cut(z2, cutpoints))
z2.exp <- diff(pgamma(cutpoints, 1/2, 1/2))
rbind(Observed=z2.obs, Expected=z2.exp * length(z2))
chisq.test(z2.obs, p=z2.exp)
par(mfrow=c(2,2))
qqnorm(xy[,1], ylab="X"); abline(c(0,1), col="Red", lwd=2)
qqnorm(xy[,2], ylab="Y"); abline(c(0,1), col="Red", lwd=2)
plot(xy, pch=19, cex=0.75, col="#00000003", main="Data")
qqplot(qgamma(seq(0, 1, length.out=length(z)), 1/2, 1/2), z2,
xlab="Theoretical Quantiles", ylab="Z2",
main="Gamma(1/2,1/2) Q-Q Plot")
abline(c(0,1), col="Red", lwd=2)
让,