s*****r
发帖数: 11 | 1 t(n)^2 = F(1, n) 如何用R验证?
> pt(1, 10)
[1] 0.8295534
> sqrt(pf(1, 1, 10))
[1] 0.811854
结果不对啊. 是应该用cdf吧。 | m*****s
发帖数: 156 | 2 这两个一样
> pt(1,10)-pt(-1,10)
[1] 0.6591069
> pf(1,1,10)
[1] 0.6591069
【在 s*****r 的大作中提到】
: t(n)^2 = F(1, n) 如何用R验证?
: > pt(1, 10)
: [1] 0.8295534
: > sqrt(pf(1, 1, 10))
: [1] 0.811854
: 结果不对啊. 是应该用cdf吧。
| s*****r
发帖数: 11 | 3 刚发现PDF之间有这个关系： df(x^2, 1, 10) = dt(x,10)/x
> df(1.2^2,1,10)
[1] 0.1547199
> dt(1.2, 10)/1.2
[1] 0.1547199
CDF 之间好像无关联。
但一般统计书（wiki也有）上的说法好像不是这个意思：
"If X is t(df=n), then Y=X^2 is an F(df1=1, df2=n)"。那位大侠能澄清一下。 | a****r
发帖数: 1486 | 4 你这儿瞎比什么呢？？！
t^2 < a 等价于 -a
你那个pt(1, 10)是从负无穷开始的
【在 s*****r 的大作中提到】
: t(n)^2 = F(1, n) 如何用R验证?
: > pt(1, 10)
: [1] 0.8295534
: > sqrt(pf(1, 1, 10))
: [1] 0.811854
: 结果不对啊. 是应该用cdf吧。
| q********i
发帖数: 795 | 5 F is the ratio of two chisq/df, t is z over sqrt of chisq/df, z^2 is chisq 1
. that is why f(1,df2)=t(df2)^2
【在 s*****r 的大作中提到】
: 刚发现PDF之间有这个关系： df(x^2, 1, 10) = dt(x,10)/x
: > df(1.2^2,1,10)
: [1] 0.1547199
: > dt(1.2, 10)/1.2
: [1] 0.1547199
: CDF 之间好像无关联。
: 但一般统计书（wiki也有）上的说法好像不是这个意思：
: "If X is t(df=n), then Y=X^2 is an F(df1=1, df2=n)"。那位大侠能澄清一下。
|
|
