トップページ⇒math
1001コメント421KB

統計学なんでもスレッド5

■ このスレッドは過去ログ倉庫に格納されています
0001132人目の素数さん2006/04/18(火) 21:18:41
前スレ
統計学なんでもスレッド4
http://science4.2ch.net/test/read.cgi/math/1123896809/
統計学なんでもスレッド3
http://science3.2ch.net/test/read.cgi/math/1097491056/
統計学なんでもスレッド2
http://science3.2ch.net/test/read.cgi/math/1068288283/
統計学なんでもスレッド
http://science.2ch.net/test/read.cgi/math/1012782106/

関連スレ
= 統計解析フリーソフト R =
http://science3.2ch.net/test/read.cgi/math/1062650510/
0002132人目の素数さん2006/04/18(火) 21:20:30
= 統計解析フリーソフト R =
http://science4.2ch.net/test/read.cgi/math/1062650510/
0003132人目の素数さん2006/04/18(火) 21:22:20
多変量で、多元グラフ(3次元以上)をかけるやつない?
0004132人目の素数さん2006/04/18(火) 23:12:50
いきなり質問なんですけど
ヒズミが極端に大きい場合の分布はどのようになるのでしょうか?
0005132人目の素数さん2006/04/19(水) 08:20:01
帰無仮説ってなんでこういう名前なのか教えて欲しい。
英語の直訳ってのは分かるんだけど、そもそも直感的に分かりにくい。

対立仮説は口語っぽいし直感的だし。

「おまえとは対立したくない」とは言うけど
「おまえとは帰無したくない」とはいわないし。。。

0006132人目の素数さん2006/04/19(水) 13:10:28
無に帰すから金仮説(なぜかマジで変換できない)
0007前スレ9912006/04/19(水) 15:35:47
レスつく前に埋まってしまったので、こちらでお願いします。

カイ二乗分布に従っているデータを正規分布に従うように
変換する方法ってありますか?
あればその方法を教えてください。
なければ、この変換方法は今後も開発されることはないものなのか
(理論的に無理なのか)教えてください。
0008132人目の素数さん2006/04/19(水) 22:39:42
統計解析初心者です。統計解析の内容そのものの質問では無いですが、予算15万以下で統計解析ソフトを購入しようと思います。比較的高機能で、初心者が使いやすいソフトで、お勧めのものはあるでしょうか?
0009132人目の素数さん2006/04/20(木) 02:43:13
とりあえずR
タダ

物足りなくなったらS
0010132人目の素数さん2006/04/20(木) 09:52:59
Rは初心者には向かないやろ。
金出すんならJMPかな。
0011132人目の素数さん2006/04/20(木) 12:54:44
用途にもよるけど
・gretl
0012132人目の素数さん2006/04/20(木) 20:32:51
>>8
excel/calc
0013132人目の素数さん2006/04/21(金) 03:19:31
日本人しか使ってるのみたことないけど
GAUSS
0014132人目の素数さん2006/04/21(金) 09:17:05
アゼルバイジャン人キタ━━━━━━(゚∀゚)━━━━━━ !!
0015132人目の素数さん2006/04/22(土) 22:58:10
>>5

今はあんまり使われなくなりましたが、

「ゼロ仮説」

って言い方もありますよ。


元々学問用語なんて英語の直訳でヘンな訳が多いし、統計自体も日本の学者の世界的な貢献なんて
ほとんど無いからしょーがないかもね。
0016132人目の素数さん2006/04/22(土) 23:02:42
>>7

これ試してみたら?

T得点:
http://aoki2.si.gunma-u.ac.jp/lecture/Bunpu/normdist/hensati.html

これが貴方が望んでるモノかどうかは正直分からんけど。
0017132人目の素数さん2006/04/23(日) 13:37:25
>>7
カイ自乗からガウスって非線形変換だよね。
それやって意味あるのかどうか?

コンピュータで乱数を作るとき
一様乱数から他の乱数作るにはよくやるけど。
0018132人目の素数さん2006/04/23(日) 17:04:45
>>7の問題について
>>16
順位を通して正規化しているので元の分布は問わないようだね。
あくまでもカイ2乗からにこだわるなら少し違うかな。

>>17
確かに意味があるのかどうかが重要だね。
一様分布を経由すればどんな分布も正規分布に従うように
変換できるが、なぜ元のままのカイ2乗ではいけないのか
なぜ変換にこだわるかが分からないね。

>>7の質問が元々平均0の正規分布に従うデータを2乗したら
カイ2乗分布に従うようにできるが、それを元に戻せるかという
意味だったら無理というのが答だな。
0019132人目の素数さん2006/04/23(日) 20:35:11
>>16-18ありがとうございます。

正規分布だと
1σ以内に68%
2σ以内に95%
3σ以内に99%
と言えますよね。
ひとつは、これが言いたいです。
自由度nのときσ=√2nを使って
1σ以内にxx%入っている
という風に。

もうひとつは、左右対称になる分布にしたいです。
今は、明らかに対称でない分布(自由度n=3〜5程度のカイ二乗)を
正規分布と仮定して話を進めているので。

上記二つが目的なのですが>>16のサイトの方法でいいのでしょうか?

対数正規分布のようにも見えるのでlogをとって無理矢理正規分布に
しようかとも考えているのですが。

>>18
元に戻せるかという意味ではないです。

002019=72006/04/23(日) 20:36:08
>>19=>>7です
0021132人目の素数さん2006/04/23(日) 21:42:20
>>19

> 自由度nのときσ=√2nを使って
1σ以内にxx%入っているという風に。

いや、これは言えるんじゃないかな?ただ必要ないから誰も言わないだけで(笑)。
例えば、自由度3のχ^2分布だと、平均3、標準偏差2.449490として、

χ^2=3+2.449490の点の上側確率⇒0.1416911
χ^2=3−2.449490の点の下側確率⇒0.09233764
∴平均±σには全体の77%が含まれる

とかね。
基本的には理論分布は端の確率しか使わないから言わないだけでしょ。
(あと負の値取れない例があるから言わないんじゃないでしょうか?)

>明らかに対称でない分布(自由度n=3〜5程度のカイ二乗)を
正規分布と仮定して話を進めているので。

その仮定が成り立つか成り立たないか以前にそれは現実のデータに合わせたんですか?
適合度検定してみての結論ですか?
明らかに正規分布に合わないのに、正規分布仮定して話を進めると無理が生じると思いますよ。

>>16は試してみないと何とも言えないです。
今のところ、何がやりたいのかイマイチ把握していないので。
0022192006/04/24(月) 00:32:05

>>21

>∴平均±σには全体の77%が含まれる
よく分かりました。
ありがとうございます。

(ある特定の統計解析をしたいというよりも)
正規分布を仮定した上で進める統計処理が多く、
目安として自由度nが30以上のカイ二乗分布は正規分布で近似できる
ことは分かったのですが、
n=3〜5程度のカイ二乗分布に従うデータでも、正規に従うように
変換できれば、色々な統計処理が進められるので、変換したい
というのが目的です。
(少なくとも左右対称な分布にはしたいです)

>明らかに正規分布に合わないのに、正規分布仮定して話を進めると無理が生じると思いますよ。
正規性の検定はしてないのですが、間違いなく正規ではないデータを正規分布と
仮定してしまっています。
正規性を仮定しないと話が進められない場合には、正規変換もできないと、
無理を承知で正規を仮定せざるを得ないと思うのですが。
0023132人目の素数さん2006/04/24(月) 01:26:14
>>22

確かに正規性仮定した方法論は多いですよね。
取り合えずT得点に変換してやってみればイイとは思いますよ。
統計解析はどの道現実のデータ相手に必ずしも成功するわけではないですし、
不具合が見つかるのも大事な分析過程だとは思いますんで。
0024132人目の素数さん2006/04/24(月) 01:46:18
( Σ )
0025182006/04/24(月) 11:24:51
>>22
>>23
T得点で正規分布を仮定した分析をするというのは結局ノンパラと同じようなものなので
素直にノンパラの方法ではだめなのかと思いますね。
相手が1群ならともかく多群ならT得点では無理ですしね。
0026132人目の素数さん2006/04/24(月) 12:23:15
>>25

いずれにせよ、その「χ^2に従うと思われる」データを見てないんで、何とも言えないんですよね(苦笑)。
ヒストグラムとか作ったんでしょうかね?オーソドックスには取り合えず多項分布化してみてからあれこれ仮定を
始める、ってのがテだとは思うんですけど・・・・・・。
この状況じゃ良く分かりませんね(笑)。
0027222006/04/24(月) 13:03:26
>>23 >>25-26
色々とアドバイスありがとうございます。

データは1群なのでT得点でやってみようと思います。
どうしても「データが正規分布する」という前提があるので
ノンパラは今までは考えてなかったです。

データはヒストグラムにしてみて(検定はしていないのですが)
カイ二乗のように見えました。
生のデータなので諸事情で公開できなくて申し訳ないです。
0028132人目の素数さん2006/04/24(月) 13:38:52
>>27

>カイ二乗のように見えました。

確かに「見えた/見えない」ってのは大事ですよね。
ところで離散型ですか?ポアソン分布は試してみました?
多分自由度3〜5辺りの「χ^2分布に見える」って事はポアソンにも見える、
とも言えると思うのですが・・・・・・。
0029272006/04/24(月) 14:06:13
>>28
データは連続型です。
形はパラメータλ=2〜3程度のポアソン分布(確率変数x)にも見えます。
ポアソン分布だとλが小さくても正規近似できるみたいですね。
(実際のデータはx=0をとらない場合もあるので、その場合にはポアソンは
まずいかもしれません。)
問題は横軸のxのとる値が整数とは限らないことです。
xに実数を入れても問題なければ解決しそうな気もするのですが。
確率関数[e^(-λ)*λ^x]/x! はx≧0が実数でも計算はできるようなので。
0030132人目の素数さん2006/04/24(月) 22:26:20
コーシー分布のモーメントが存在しないのは何故ですか?
級数1+(1/2)+(1/3)+(1/4)+...が収束しないのと同じ理由といいますが、
良くわかりません。
0031132人目の素数さん2006/04/25(火) 21:27:10
二項分布を正規分布で近似するときに0.5の調整をするのはなぜですか?

二項分布の確率は高さ、正規分布の確率は面積だとすると
プラスマイナス0.5の調整をすることで正規分布の確率(面積)が二項分布の高さ×1になり
二項分布の確率とちかくなるという理解でいいのでしょうか?
0032132人目の素数さん2006/04/26(水) 01:11:20
>>30
モーメントの定義に従って計算すると発散するからとしか言えないのでは?
それは積分だが、その級数と計算が似ているのでそう説明されたのでしょう。

>>31
そうですね。柱は中心から0.5はみ出て立つことになるのでそれを補正した方が
面積が近くなるという理解で正しいです。
0033132人目の素数さん2006/04/27(木) 11:21:24
すみません教えてください。
対数正規分布のパーセント点ってどう求めるのでしょうか?
μ+zα σの指数変換でよいのでしょうか?
お願いします
0034132人目の素数さん2006/04/27(木) 12:38:38
すごく初心者的な、おそらく的外れな質問で申し訳ないのですが、
中心極限定理について質問させてください。
最尤推定量の漸近有効性は
√n(\hat{θ}_n−θ) → N(0,I(θ)^{-1}),
where \hat{θ}_n:最尤推定量,θ:真のパラメータ(ベクトル),
と記述されますが、これは最尤推定量に√nをかけたものの
漸近分散(共分散行列)がI(θ)^{-1}だということですよね?
すると、最尤推定量そのものの漸近分散はどうなるのでしょう?
0になるんでしょうか。そんなはずもなさそうなんですが・・・。
どなたかよろしくお願いいたします。
0035132人目の素数さん2006/04/27(木) 22:12:13
便乗質問
漸近分散ってどういう意味をなすんですか?
無限回推定したときの分散って意味?
でもサンプルサイズのnを無限大に飛ばしてるんだよね?
無限大のサンプルをつかって推定したとするなら、
何回も推定するって意味じゃないですよね。
無限大のサンプルつかって一回だけ推定するなら,
分散の意味自体がわからない・・・。
誰か教えてください。
0036132人目の素数さん2006/04/27(木) 22:51:13
>>33
いいんじゃないかな。

>>34
データが増えると推定量は真値に近づくので
分散が0に行くのは自明では。

>>35
データ数が増えると推定量の分散は0に行くけどその行き方に違いはある。
それがルートn倍すると見えてくるけどその意味は、大きなデータを取ったときに
早く安定するのかどうかの違いということじゃないかな。
0037332006/04/27(木) 23:21:08
>36
ありがとうございます。
単純にこれでいいのかなと不安になりまして・・・
0038352006/04/27(木) 23:46:24
>>36
んー、その推定量が安定しているかの測度という感じですかね?
推定量Aと推定量Bの漸近分散を比較してこちらの方が小さいから良い推定量みたいな?
むずかしいですね。ありがとうございました。
0039132人目の素数さん2006/04/28(金) 10:39:15
「平均」の統計学的な意味って何ですかね?
0040342006/04/28(金) 13:58:54
>>36
そうですか・・・。やはり、分散は0になるんですね。
早く安定するかどうかということなんですね。
ありがとうございました!
0041132人目の素数さん2006/04/28(金) 23:09:21
>>39
統計学的と言われても難しいなあ。
結局、重心という意味しかないのでは?
それともそんなものをなぜ後生大事に使っているのかということかな?
0042132人目の素数さん2006/04/28(金) 23:44:38
質問させてください。確率変数Xが標準正規分布に従い,その密度関数をφ(x)とするとき、
1/Xの期待値

∫1/x φ(x)dx

って存在しますか?よろしくお願いします。
0043132人目の素数さん2006/04/28(金) 23:46:50
>39
まぁ、単に標本平均っていうならば、一般にいわれる「平均」と考え方は
変わりないんじゃね?
0044132人目の素数さん2006/04/29(土) 01:29:18
>>42
発散しますね。
0045422006/04/29(土) 10:23:09
>44
やはりそうですよね・・・.ありがとうございます.

たとえば,正規分布の位置パラメータが
大きく(>0),0以下が生起する確率が非常に低いときでも,0以下になる確率が
0ではないから,0のポイントで発散しますよね.この場合何らかの形
(近似にしろ何らかで抑えるにしろ)で期待値を記述することはできないですかね?
0046442006/04/29(土) 11:39:08
そもそもどういう場合にそんなことが必要なのでしょうか?
xが負なら0近くで−∞に行きますから0を挟んで極端な
挙動を示すのですが。
0047132人目の素数さん2006/04/30(日) 05:14:52
分からない問題があるので良かったら回答お願いします
英語の文で申し訳ないのですが・・・
Supervisores were asked to rate the productivity of employees on their jobs. A researcher wishes to see whether married men receive higher ratings
than single men. A rating scale of 1 to 50 yielded the data shown here. At アルファー=0.01 is there evidence to support this claim?
この問題は、The Wilcoxon Rank Sum Test を使って求めなければならないのですが、
Single ... 48, 46, 42, 50, 38, 36, 40, 31, 28, 24, 49, 34
Married ... 44, 35, 41, 37, 42, 43, 29, 31, 37, 32, 36
がデータです。
自分の回答としては、R=122.5
ミュウ=132
ロウ=16.2
Z=-0.586
が出ています。この問題はtwo-tailed ですか?それとも、right-tailed testですか?
もしお分かりであれば、その時のHypothesisも教えてください。よろしくお願いします。
■ このスレッドは過去ログ倉庫に格納されています