トップページ⇒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/
0655132人目の素数さん2006/11/19(日) 01:07:36
>>653
ぜひ初心者向けの統計学の本を一冊買いましょう
大抵の場合は載っていますし、統計関連の他の勉強にも役立ちます
0656132人目の素数さん2006/11/19(日) 08:00:06
>>654
中央値を取るのかどうかとは別に、まず、Bartlett検定で等分散か
どうかを検定して、等分散とみなせるのだったら迷わず2way repeated
measures ANOVA。みなせないのだったら二元配置Friedmanもやってみて
都合のよい結果がでるほうを採択(うは


0657132人目の素数さん2006/11/19(日) 10:21:51
>>656
顕著に外れ値が出るってことだったから生データで等分散性は
ないでしょうね。(データが少ないと棄却はされないという可能性はあるが…)

ネットで見ても不等分散ならノンパラって見かけるけどノンパラの
前提は未知だけど同じ分布なんだけどね。まあパラよりロバストだとは
思うけど。中央値を使うという方法はロバスト的にはよさそうだね。
0658132人目の素数さん2006/11/19(日) 17:40:40
無知で文系の私にどなたか教えてください。
3つの尺度を使って分析を行いたいのですが、
2つの尺度の関連性を弱める媒介として残りの1つの尺度が
影響を与えているかどうかを調べたいとき、
どの分析方法を用いればいいのでしょうか?
重回帰⇒パス解析でいいのでしょうか・・?
0659132人目の素数さん2006/11/19(日) 19:08:04
>>658
たった3変数なんだよね?媒介の変数は決まっているんだよね?
単に相関と偏相関を見ればいいんじゃね?
0660132人目の素数さん2006/11/20(月) 07:16:40
>>656
なんで2wayなんですか?3wayじゃないの?
0661132人目の素数さん2006/11/20(月) 13:21:17
なんで?
0662132人目の素数さん2006/11/21(火) 18:40:45
すみません、おしえてください。

2群比較の問題で、処理群,コントロール群が存在する。
処理前観測値をB、処理後観測値をFとする。
B〜Poisson(μ)とし、
FにはBとの加法関係
F=B+Poisson(Q*B)
を想定する。ここに、Q〜Gamma(0.5,2.0)とする。
さらに、処理群には処理効果k(0<k<1)を掛け,kFとする。
このとき、処理効果の検出力を検討する。

英文論文からの転記なんですが、上記の仮定の意味が良くわかりません。
まずPoisson分布のパラメータに確率変数をもってくることが意味不明です。
こういうモデルってあるんでしょうか。
どなたかよろしくお願いいたします。
0663132人目の素数さん2006/11/21(火) 22:37:20
>>662
この分布を知っている訳ではありませんが、類似のものとしては
ベータ・二項分布がありますね。二項分布のパラメータpが
ベータ分布に従うというものです。これも同じような意味ですね。

理解は乱数で考えると分かりやすいですよ。
BとQをそれぞれポワソン乱数とガンマ乱数で作り、
その値でQ*Bをパラメータとしたポワソン乱数を作って
Fを構成するのです。ただし、Bは処理前値として対象ごとに
違うことは自明ですが、Qは対象ごとか処理ごとかが
不明ですね。

06646622006/11/21(火) 23:33:13
>>663
早速にご丁寧な回答をありがとうございます。
成程、パラメータが確率変数とはそういう解釈をすればいいのですね。
一歩前進できました。ベータ二項分布も少し調べてみます。

因みに、処理群とコントロール群の差はkをかけるか否かであると思われますので、
Qは個体毎に与えられると思います。

しかし、BとFが加法関係構造にあるというのは未だ不明瞭です。
たしかにB+Po(QB)で和をとっているから加法というのもわからないでもないですが・・・

もう少し調べてみます。ありがとうございました。
これだ!というのがありましたらまたご教示いただきますと幸いです。
0665132人目の素数さん2006/11/22(水) 07:40:24
国勢調査の細かいデータって電話とかすれば提供してもらえる??
0666132人目の素数さん2006/11/23(木) 08:27:56
Barndorff-Nielsenはバルンドルフ・ニールセンで読み方あってますか?
0667132人目の素数さん2006/11/23(木) 09:40:04
AICで京都賞
ttp://www.inamori-f.or.jp/ja_kp_lau_thi.html
0668132人目の素数さん2006/11/23(木) 17:00:06
>>667
あ、それの公演(12日)いってきた。よかったよ。
0669132人目の素数さん2006/11/24(金) 23:44:29
はじめまして、初歩的な質問させてください。
実験を行い、その施行前と後でその変化に有意差があるかどうかを調べたいのですが、
t-検定:一対の標本による平均の検定
を使用していいのでしょうか?
これだけの情報だとわかりかねる場合はすみませんがレスお願いします
よろしくお願いします
06706692006/11/25(土) 07:34:00
>>669です。
何の知識もなしに質問してすみませんでした。
この件に関しては自己解決しました。
今度は解析をしてみて、結果の解釈に困っています。

ウィルコクソン符号付順位和検定でデータ数は10です。
危険率を5%としたときに、P>0.05となるのですが、
T値は帰無仮説の棄却域に入らないとき、
有意差があると判断してよろしいのでしょうか?
何度もすみませんが、よろしくお願いします。
0671132人目の素数さん2006/11/26(日) 22:13:54
変数同士を合わせて新変数をつくり、新変数間でt検定を行いたい
のですが、変量の単位が%です。各変数は独立です。素データに戻
って新変数の%を算出できない事情があり、なんとか変量同士の合
成を行いたいのですが、標準化したら、%のデータ間の加減乗除で
合成変量をつくり、合成変数間の差の検定をして大丈夫でしょうか。

例えば視覚実験でターゲットの性質により検出の違いを見たいとして、
変数xが、1つの記号をターゲットとした場合の検出率、
変数yが、2つの記号をターゲットとした場合の検出率、
変数uが、1つの漢字をターゲットとした場合の検出率、
変数vが、2つの漢字をターゲットとした場合の検出率
で、それぞれの検出率だけがデータとして手元にある場合です。
ターゲットが1つか2つか、の差を見たい。あるいは記号か
漢字か、で差を見たいので、変数同士の合成が必要なのですが、
そもそも刺激中にそれぞれのターゲットが何個あったのか、
不明なので、検出率の単純合計は危険かと思われます。

変量を標準化したら、加減乗除して新変量をつくり、比較でき
ますか?もっと、ほかの操作が必要でしょうか。
すみませんが、よろしくお願いします。
0672132人目の素数さん2006/11/26(日) 22:35:25
>>671
x, y, u, vを標準化(平均0分散1)にしたら変数を足したものの平均も0になるので
どの組み合わせでも検定できません。それとも分散だけ揃えるのでしょうか?
各変数の分散に差があるかどうかの検定をしてなければ単純合計でいいように
思いますが。(と言っても対応があるなら足せますが、なければ分散分析するので
しょう。)
なお、もしターゲット数が検出率の平均にも影響すると考えているならそもそも
どうやっても比較できないでしょう。(たとえば、ターゲット数が多いと検出率が下がって
くると思っているなら検定で差が出たとしてもそれは個数や種類の差なのかターゲット数
自体のせいなのか区別する方法がありません。)
0673132人目の素数さん2006/11/27(月) 07:50:22
>>672
ご回答ありがとうございます。勉強が足りなくてすみません。

扱っているデータは、作業時間が一定で、被験者ごとに遂行
量が異なり、遂行分で遭遇した異なる種類のターゲットを、
それぞれ正しく検出した割合が変量として自動算出された
ものです。
その変量を合成して、別の角度から検討したいと思っている
のですが、割合は割合にすぎないので、単純合計してはいけ
ないように思ったのです。
変量xもyも同じ被験者から得られたもので、例えば全被験
者が10人だったとして、10人分のxとyの値を合計して
20で割ったものを、新しい変数の平均値としたり、xの平
均とyの平均を単純合計したものを新しい変数の値としては
ダメですよね。と、不安になりまして…。



0674132人目の素数さん2006/11/27(月) 09:38:01
>>673
前回書いたようにターゲット数によって検出率が変化すると思っているのかいないのかが
まず焦点です。思っているならお手上げです。
思っていないのなら対応がついているようなので平均を取ったり単純合計したりして
比較してもなんら問題ありません。
個数の違いならx+u-y-vの1標本の検定、種類の方はx+y-u-vの1標本の検定でしょう。
ただし、元が比率なので変化が大きい場合は線形ではないでしょう。
ロジット変換などの後に上記の計算をした方がより適切でしょう。
0675132人目の素数さん2006/11/27(月) 23:59:20
>>674
何度もありがとうございます。よく分かりました。
ターゲット数と検出率の関係は、遂行量データがあるので
別に考察できます。
とりあえず、無視してよし、として分析していきます。
分散分析をしてみて大体の傾向はつかめましたが、属性の
違いの作用をみたい場合は、どうしてもデータの合成に
直面しますので、単純に足してOKとうかがって、安心し
ました。
本当にありがとうございました。
0676sage2006/11/28(火) 00:27:27
質問ですが、正定置対象行列Σがあったとして、

|∂Σ/∂(Σ^-1)|  ←偏微分して行列式をとる。

を解きたいんですがどうしても解けません。
余因子を使って色々やったんですが、やり方がおかしいんでしょうか?
行列の微分に詳しい方いたら、ご教示下さい。
ちなみに、Σは統計学の分散共分散行列です。
0677132人目の素数さん2006/11/28(火) 09:57:22
>>675
分散分析してるのに属性の違いが見られないなんて…
0678132人目の素数さん2006/11/28(火) 14:25:31

{192 234 12 395}

という数値群があった場合、
このときの平均はもちろん
(192+234+12+395)/4 =208.25
です。
それでは、「この数値群の分散をもとめろ」
という場合、どういう公式・数式を使いますか?
0679132人目の素数さん2006/11/28(火) 14:30:13
分散:それぞれの値と平均との距離(偏差)の二乗の平均
なのだから、
{(192-208.25)^2 + (234-208.25)^2 + (12-208.25)^2 + (395-208.25)}/4
=74316.75

これでいいんですか?




0680132人目の素数さん2006/11/28(火) 14:34:56
質問させてください。
1個あたり物を10で仕入れ、50で販売し、余りは5で製造元に引き取っている店があるとする。
平均400分散120^2の正規分布に従っている。
いくつ仕入れると、最も利益が大きくなるか?
という問題で、正規分布から求めるらしいのですが、excelで正規分布表を作るときに
=(1/(120*SQRT(2*PI())))*(EXP(-(((B2-800)/120)^2)/2))
と入力したのですが、どこかが違っているみたいで正しい値が出ません。
この式のどこがおかしいのでしょうか?
0681132人目の素数さん2006/11/28(火) 15:00:57
>>680
>=(1/(120*SQRT(2*PI())))*(EXP(-(((B2-800)/120)^2)/2))

=(1/(120*SQRT(2*PI())))*(EXP(-(((B2-400)/120)^2)/2))
では?
正規分布の密度関数は normdist という関数で表せる。
0682132人目の素数さん2006/11/28(火) 15:06:20
分散がマイナスになることなんてありますか?
0683132人目の素数さん2006/11/28(火) 15:06:57
>>682
ない。
0684132人目の素数さん2006/11/28(火) 15:11:14
>>679
{(192-208.25)^2 + (234-208.25)^2 + (12-208.25)^2 + (395-208.25)^2}/4
=18579.1875
じゃないかな?
0685132人目の素数さん2006/11/28(火) 15:34:55
>>684
そのとおりです。
計算ミスでした。ご指摘ありがとうございます。
06866802006/11/28(火) 15:47:04
>>681
ありがとうございます。この問題の解き方が分からないのですが、どなたか教えてくれませんか?
06876802006/11/28(火) 15:48:02
>>681
ありがとうございます。この問題の解き方が分からないのですが、どなたか教えてくれませんか?
0688132人目の素数さん2006/11/29(水) 11:55:06
ある報告書に「信用区間による検定とは,ゼロが特定の区間(この場合95%,90%区間)
に含まれるかどうかの検定である」と書かれているのを見たのだけど、
1.信用区間と信頼区間は同じ?
2.95%信頼区間に0が含まれたら有意ではないとするのはよくあること?
教えてください。
06896882006/11/29(水) 12:03:06
1.は自己解決しました。ウキペディアによると
信頼区間(confidence interval)と信用区間(credible interval)は違うもののようです。
ベイズ統計の報告書を読んでいたので、信用区間という言葉が出てきたようです。
0690132人目の素数さん2006/11/29(水) 23:52:34
すごい初心者な質問をさせてください。学校ででた問題で、
X〜N(30,4^2)とするとき、
(1)P(X<1)
(2)P(20<X<24)
(3)P(22<X<26)
これらの確率を求めるものです。(2)と(3)はわかったのですが(1)だけ
どうしてもわかりません。標準化すると-29/4もとい7.25になってしまい、
数表見ても載ってないんです。なんとなく0に限りなく近いかなとは思うのですが。
よろしければ教えてください。
0691132人目の素数さん2006/11/30(木) 01:09:32
>>690
それはいくらなんでも出題ミスじゃないの?
Rで計算したら一応2.083858e-13って答が出るけど
そんなものを求めさせて意味があるとは思えない。
(学校の宿題なら教科書の数表の範囲外を出すのは
まずいでしょう。統計ソフトを普段から使わせる授業なら
ともかく。)
0692132人目の素数さん2006/12/02(土) 00:42:54
質問させていただいてもよろしいでしょうか?

現在、マウスのある細胞に薬剤をかけて、
増殖速度がどのように変わるか、
という実験をやっています。

これをマウスを変えて三回ほど行いました。
(マウスの育った環境は同じ)

この場合、有意差検定というのは、
�@、一つの対照群-試験群間で行うものなのか?

�A、三回の試験群、対照群のそれぞれの平均を取って、
その統計量を元に2群間比較を行うべきなのか?

どちらか分からず立ち往生しております。

�@の場合だと、あるものは有意差が出て、あるものは有意差が出ず、
という感じになり・・・。
また、それをどのようにまとめれば良いのか分からず・・・。

�Aだと、有意差は出るのですが、
それが統計学的に妥当なことかどうかがわかりません。

少しでもいいので、どなたかアドバイスいただけると幸いです。
0693132人目の素数さん2006/12/02(土) 01:03:57
>>692
以下のような構造で分散分析して因子(種類)が棄却されるかどうか見ればいいのでは?

種類 回数 結果
対照群 第1回 xxxx
. . . . . . . .
対照群 第1回 xxxx
試験群 第1回 xxxx
. . . . . . . .
試験群 第1回 xxxx
. . . . . . . . . . . . . . .
. . . . . . . . . . . . . . .
対照群 第3回 xxxx
. . . . . . . .
対照群 第3回 xxxx
試験群 第3回 xxxx
. . . . . . . .
試験群 第3回 xxxx
0694132人目の素数さん2006/12/02(土) 01:38:24
>>693
ありがとうございます。
やってみます!
0695132人目の素数さん2006/12/02(土) 13:32:13
(別スレで質問していましたが、こっちに移動するよう提案されました)

「umbrella sampling」の日本語訳と定義を教えて下さい。
統計学の用語らしいですが、数学の英和辞典を見たり、ググったり
しましたが分かりませんでした。
0696132人目の素数さん2006/12/02(土) 16:46:39
MT法と重回帰分析の違いを教えてください。
0697132人目の素数さん2006/12/02(土) 22:09:06
ワロチのt検定wwwwwwwwおkwwwww
0698132人目の素数さん2006/12/02(土) 22:16:44
>>695
アンブレラサンプリングでいいんじゃないの?
 ググったらいっぱい出てくるから・・・。
0699132人目の素数さん2006/12/02(土) 22:43:33
統計専門の方、突然ですが一つ伺いたいことがあります。
「ある群から標本を抽出した結果あり、
それぞれ偏りなく同様に等しい確率で抽出された(単純無作為抽出)ということを検定したい」
のですが、何か良い処理法はありますか?
具体的には70の中から重複可能で50回抽出した結果があり、
そのうち6回抽出されたのが1つ、4回抽出されたのが1つ、
3回抽出されたのが2つ、2回抽出されたのが8つ、
1回が19個、残りは抽出されなかった…といった感じです。
(海外にいるので適当な本がみつからず、
ネットでもうまくヒットしないのでここに書くことにしました)
07006952006/12/02(土) 23:26:21
>>698
Thx
0701132人目の素数さん2006/12/03(日) 10:56:14
t検定のことでお伺いしたいしたいのですが。
3群比較は偶然誤差が大きくなるので、t検定ではなく分散分析を一般的に言われています。
しかし、次のように考えられないでしょうか。
今3群、ケースA群、ケースB群、コントロール群があるとして、
A対B、A対コントロール、B対コントロールでそれぞれt検定します。
このとき有意水準を5%に設定したならば、実質的には1-(0.95)^3=0.142625、
つまり約14%となってしまいます。しかし、有意水準を1%に設定したならば、
1-(0.99)^3=0.029701、つまり約3%におさまります。
ならば、3群比較に1%有意水準を使うのは差し支えないといえるのでしょうか?
どなたかご教授していただいてよろしいでしょうか。お願いします。
0702132人目の素数さん2006/12/03(日) 10:59:07
>>699
合計抽出回数が51回となっていてどこかにミスがありますね。
概算だけど、重なるのが最大3回までで6割5分ぐらいはいきそうです。
4回までで9割ぐらいと言うところでしょうか。5回までだと9割7分ぐらいに
なってしまうようで6回重なったというのが本当なら棄却されるレベルかも。
ただし、正確な計算はとても大変であくまでも目安です。
0703132人目の素数さん2006/12/03(日) 12:42:08
>>701
普通に多重比較すればいいんじゃね?
独立じゃないからその計算も正確じゃないし。
0704132人目の素数さん2006/12/03(日) 22:13:06
>>702
ありがとうございます!サンプル総数については抽出なしが29個で、全50個です。
6回選択されたサンプルについては、ややバイアスがかかっていた可能性があるので、
棄却する必要があるかと考えていました。具体的な算式はどうすればよいのですか?
当初は正規曲線(高校数学で習う曲線)に回帰させれば良いのかと思いましたが、サンプル間に連続性がないので、
適用できるかわかりませんでした。他にローレンツ曲線がありますが、これはどうも検定とは関係なさそうに思えたのですが。
>>702さんは、どのような算式あるいは解析ソフト(?)を使われたのですか?
0705132人目の素数さん2006/12/04(月) 00:16:52
>>704
6+4+3*2+2*8+19=51なんですけど。1回なのが18でしょうか?
まあそれはさておき計算は単に純粋に確率を求めただけです。
ただそれだと場合の数が多くなるので適当に端折って大体求めてます。
計算はRでプログラムを書いて行いました。計算に間違いがないなら最も確率の
大きな組み合わせは0回が33、1回が26、2回が9、3回が2のときの0.04324218のときで
それ以後それより確率の小さなパターンが並びます。そのようにせずに何らかの
近似をする方法があるのかは知りません。
もし18,8,2,1,0,1が正しい50回の結果ならその確率は1.835095e-06ぐらいと
思われ、パターンの上から373個目ぐらいに当たります。
それより小さな確率となるパターンの合計確率はだいたい0.022です。
(なぜ「ぐらい」かというと枝刈りをして計算しているので刈った枝が計算結果に
影響している可能性があるためです。)
0706132人目の素数さん2006/12/04(月) 08:50:22
∞×0は、∞でしょうか。0でしょうか。
0707132人目の素数さん2006/12/04(月) 11:35:32
デルタ関数の積分は1だよな。
0708132人目の素数さん2006/12/04(月) 15:35:54
質問です。
複数の系列の散布図をエクセルに重ねて表示させています。
重なっている部分が見にくいので、
信頼楕円を換わりに描こうと思っています。

信頼楕円を描くアルゴリズムが知りたいです。
また、研究室のパソコンにオリジンがあり、
これでできそうな気がするのですが、
やり方がわかる人がいたら教えてください。
0709大学生2006/12/05(火) 10:50:12
確率変数Xが平均μ、分散σ^2をもつときE(X^2ー4)をμとσ^2で表せという問題なのですが、
誰か教えてください。よろしくお願いします。
0710132人目の素数さん2006/12/05(火) 17:30:22
お手元の統計の教科書をごらんください。
たとえそれが高校の数学の教科書であっても書いてある筈だから。
0711132人目の素数さん2006/12/05(火) 22:17:13
正規分布を書いて尖度と歪度を表したいのですが、縦軸と横軸には何の数値を書けばいいですか?
0712132人目の素数さん2006/12/05(火) 23:43:41
>>702、705さん
ありがとうございます、わかりました。
またRをインストールしました(でもプログラムがうまくかけない)
データの方ですが再チェックしました。計算法を知りたかったので、ちゃんとみていませんでした。
すみません。0*40+1*19+2*8+3*2+4*1+5*1=50です。
プログラムの書き方を紹介している良いサイトやこの検定(確率の算出)にかかわる用語やキーワード
があればおしえてください!
0713132人目の素数さん2006/12/06(水) 00:12:53
>>712
(40,19,8,2,1,1)のパターンなら確率0.000286で174番目、それより低い確率の合計は
0.0327でまだ棄却されるレベルですが、前よりはましですね。
計算の元になるサイトはないです。自分で考えただけですから。
(70,0,0,0,0,0)から始めて1回の抽出ずつ計算します。次は確実に(69,1,0,0,0,0)ですね。
その次は確率69/70で(68,2,0,0,0,0)で確率1/70で(69,0,1,0,0,0)に分かれます。
あとはこれを丁寧に計算していくだけです。ただし、すべてのパターンを計算するのは
メモリーと時間の関係上無理だと思います。私は6回以上の抽出はひとまとめにして
パターンの保存も200までで打ち切ってそれより小さなものは合計確率のみで計算
していきました。
0714132人目の素数さん2006/12/06(水) 13:10:18
>>702(705,713)さん
ずいぶんと手間をとらせてしまったようで。。。ありがとうございます!
わかました。データはやはり無作為というには無理があるようですね。
当初、データを(40,19,8,2,1,1)と並べて、正規分布(の右半分)と使って比較するのは
どうかと考えたのですが、これは数学的に意味のある処理でしょうか。
他にもデータがあるもので。。。いずれにせよ、お聞きする前はまったくわからなかったので
検定の結果がわかったのはずいぶんと参考になりました。ありがとうございました。
0715132人目の素数さん2006/12/06(水) 18:00:04
すみません、初歩的なのですが質問させてください。。

エクセルを使って発生させた正規分布に従う乱数をもとに、
自由度を設定してカイ2乗分布に従う乱数を作成するには
どうすればよいのでしょうか?
0716132人目の素数さん2006/12/06(水) 18:33:23
>711
正規分布だったらskewnessもkurtosisも定数なような・・・
0717132人目の素数さん2006/12/06(水) 20:01:34
いろいろ考えましたが、そんなに単純には算出できないことがわかりました。
確かに計算は煩雑ですね。>>714より。
0718132人目の素数さん2006/12/06(水) 20:20:52
>>702(705,713)さん
たびたびすみません。モンテカルロ法のように、50回のランダム抽出を繰り返して、
一番多い組み合わせを取るのはどうでしょうか。。。組み合わせが多すぎて駄目かな。。。
0719132人目の素数さん2006/12/07(木) 00:36:14
>>718
私も初めは乱数でできるだろうと思ったのですが、重複が5回以上の所は滅多に
起きないため、正確な頻度を見積もるには何回行えばいいのか分からず結局
正確に計算してみようと思ったわけです。

さて、枝刈りして得られた375のパターンから各重複回数に対する期待値を近似的に
求めてみました。
(確率0.0225に対する残りの多くのパターンは無視してますので合計が70になるように
補正をかけています。)
0回:34.02379, 1回:24.76589, 2回:8.81109, 3回:2.02631, 4回:0.33286, 5回:0.03884,
6回以上:0.00123
下がり方は緩やかなため正規確率の当てはめはよくないですね。
何がいいのか分かる猛者がいたら教えてほしいです。

これで、上記があれば検定は可能かと思ったのですが、カイ二乗近似では4回以上の
少ないところは当てはまりがよくなくて使えませんし、正確に多項分布で行うにはこれまた
大変のようです。
たとえば、6回以上が出てくるのはおかしいと言えるのですが、5回が出てくるともうだめ
なのかというと(35,24,9,1,0,1)だと確率0.14ぐらいで起こります。一律に基準を設けるのは
とても大変そうですね。
また、データがたくさんあるとのことですが、多くあればあるほどまれなパターンも生じます
ので先のデータが30個のうちの1つだとすると確率0.03ですので普通のことかもしれません。
0720132人目の素数さん2006/12/08(金) 17:15:10
真面目な質問なのですが、
∞×0は、∞でしょうか。0でしょうか。
0721132人目の素数さん2006/12/08(金) 17:51:17
>>720
∞と0と×の定義による
0722132人目の素数さん2006/12/08(金) 18:34:17
>>720
0と無限大の濃さによる
結果はより高位なほうをとる
0723132人目の素数さん2006/12/08(金) 21:51:43
確率統計って主観なしの学問なの?
データからある数値(例えば平均値)を抽出して何かを予測してみたり
散布図に回帰直線を書いて相関関係があるな、とか
ようするに帰納法でしょ。
帰納法って数学じゃないんじゃないの?
演繹こそが数学だと思うんだけど。
つまり確率統計って数学じゃないんじゃないの?
0724132人目の素数さん2006/12/08(金) 23:23:54
>>723
二元論にこだわるなら
http://www.nn.iij4u.or.jp/~hsat/misc/math/deducinduc.html
でも読んでみたら?
0725132人目の素数さん2006/12/09(土) 00:04:01
初心者です。臨床研究をしています。
ある疾患群で特定の成分が血液中に多い、というねたでリサーチしてましたら
ある成分の血中濃度に男女差があることがわかりました。この場合
コントロールとある疾患群で、その成分濃度に差があるというためには
男女別々に解析しないとだめなのですか?
教えていただければ幸いです。
0726132人目の素数さん2006/12/09(土) 00:32:56
>>725
「交絡因子」って言っても分からないよね…。
男女差を考慮しながら分析する方法があるのだが、
その用語でぐぐってあとはがんばってくれ。
(検索後もちんぷんかんぷんならまじで素養が低すぎる。
ここではなく誰か生身の統計の分かる人に相談した方がいいよ…。)
0727132人目の素数さん2006/12/09(土) 00:36:42
男女差が交絡因子ですか?

0728132人目の素数さん2006/12/09(土) 00:55:49

わかりました!さっそく明日頑張ります
0729132人目の素数さん2006/12/09(土) 13:02:41
>>719
やはり難しいのですね。でも>>719の説明でかなりすっきりしてきました。
とりあえず明確な解析結果は避けて、無作為抽出にしてはやや偏りがみられるという
判断でいこうと思います!!!!
本当にありがとうございました!
0730132人目の素数さん2006/12/09(土) 23:14:07
>>723
統計科学は数学じゃないよ。経験則的なところもあるし。
帰納的なのが統計科学の特徴じゃん。まずデータありき。
主観なしではないけどね。結果を自分なりに解釈するのが統計家だし。
0731132人目の素数さん2006/12/10(日) 07:02:01
>>730
データがあって、主観的な統計モデルがあって、適当に当てはめて悦になるのが統計科学です。
0732132人目の素数さん2006/12/10(日) 07:43:10
>>731
随分主観的だなw
0733132人目の素数さん2006/12/10(日) 07:53:36
統計って科学なの?
統計に反証可能性ってあるっけ?
0734132人目の素数さん2006/12/10(日) 09:06:14
>>733
教えて君退場!
0735おしえてください君2006/12/10(日) 11:26:37
私も臨床研究をしているのですが、
ある血液検査値が基準値より上昇している「群」についてまとめたいのですが、

�@一般的には、過去のデータと比較しようが無いときに
群のうち何%以上の個体が上昇しているのが「多い」と考えるのでしょうか?

�Aもし基準が無ければ、「多い」という基準はどのように示すのが妥当
でしょうか?

おねげします
0736132人目の素数さん2006/12/10(日) 12:14:40
だいたいものの判断を完全に科学的に出来るなんていうのはありえないし、
そんな科学は社会的に必要とされない
0737おしえてください君2006/12/10(日) 18:10:10
すみません、よくわかりません
0738132人目の素数さん2006/12/10(日) 19:36:07
>>735
上昇してるって事は過去のデータがあるんじゃないの?
0739132人目の素数さん2006/12/10(日) 21:14:16
質問よろしいでしょうか。

年齢、スポーツ歴などランダムに選出した被検者にそれぞれ100m走や
潜水などの運動をしてもらい、平常時の心拍数と比較して最も変化の
あったものを求めたいのですが、この場合は対応のあるt検定を使用して
よろしいでしょうか?

宜しくお願い致します。
07407302006/12/10(日) 22:45:30
>>731
まぁある意味ではそうだねw
0741132人目の素数さん2006/12/11(月) 00:23:38
>>739
最も変化のあったものというのなら個人での平常時との差を種目ごとに取って
その平均値を比べれば分かるよな。それが有意かどうかということになると
多重性を考慮するかどうかって問題になるが。
0742132人目の素数さん2006/12/11(月) 00:26:10
>>704

>6回選択されたサンプルについては、ややバイアスがかかっていた可能性があるので、棄却する必要があるかと考えていました。具体的な算式はどうすればよいのですか?

であれば、

>ある群から標本を抽出した結果あり、
それぞれ偏りなく同様に等しい確率で抽出された(単純無作為抽出)

ではなく、

p1=(6回発生した問題の事象Aが起こる確率)
p2=(のこり69個の事象ひとつひとつについて起こる等しい確率)
(日本語がおかしいですが、適宜解釈してください)

として、「p1=p2 を検定する」が正しいのではないでしょうか?
単純な問題なので、自分で考えても分かると思いますが、

p1=p2=1/70 と仮定し、50回中事象A の起こりうる回数をN とすると、
P(N=0) = 50C0  (1/70)^0 (69/70)^50 = 0.4876
P(N=1)=  50C1  (1/70)^1 (69/70)^49 = 0.3529
P(N=2) = 50C2  (1/70)^2 (69/70)^48 = 0.1253
P(N=3) = 50C3 (1/70)^3 (69/70)^47 = 0.0290

P(N<=3) = 0.9943 なので、1%棄却粋は、N>=4 であり、
6回起きているので、p1 と p2 は等しくない。

こちらのほうがまだいいのではないですか?
0743132人目の素数さん2006/12/11(月) 00:31:41
>>739

t検定は正規分布を仮定しての検定ですので、t検定を行う前に、まず正規分布に適合するかどうかを検定しないと、信頼性は無くなると思うのですが・・・。

あまり、古典理論を勉強していないので、嘘かもしれません。
0744132人目の素数さん2006/12/11(月) 01:04:19
>>742
もし特定のサンプルがあやしいと事前に思っているならそれでいいでしょうね。
どれか分からないが偏っているのでは?という問題として719は解きました。
0745132人目の素数さん2006/12/11(月) 02:18:31
実務的な統計学の本によくある言葉・・・
「著者の経験から・・・・が妥当なのではないかと
考えられます。」
この学問(?)終わってる・・・
0746132人目の素数さん2006/12/11(月) 02:25:21
統計学を人文学に認定します
これにはだれも反論できない
0747132人目の素数さん2006/12/11(月) 08:14:08
>>746
おつかれさんですっ
0748おしえてください君2006/12/11(月) 10:29:42
なるへそ
0749132人目の素数さん2006/12/11(月) 17:07:41
rをデータ(Xi, Yi), i=1, 2, …, n の相関係数とすると、データ(aXi+b, cYi+d)の相関係数は

ac>0のときr
ac<0のとき-r

となりますか?教えてください。
0750132人目の素数さん2006/12/11(月) 21:26:47
てか
ぶっちゃけこのスレ板違いじゃね?
削除依頼ヨロ
0751132人目の素数さん2006/12/11(月) 21:56:42
http://yasai.2ch.net/psycho/kako/1013/10137/1013732790.html

>これは確率空間の数学的構造が我々が住んでいる世界の構造や関係に
>マッチしていないからではないか、マッチする数学的構造は
>別にあるのではないか、

>確率論的なモデルが、 ・・・
>精度がじょじょに上がって来て、そのうちに実用に供されそうだ、
>というレベルでは全くない
>絶望的レベルであり、
>「学習率ゼロ」の惨憺たる状態


絶望的レベルだそうです。。。
マッチしてないそうです。。。
0752132人目の素数さん2006/12/11(月) 21:58:28
>カオス的変動の予測にことごとく、
>しかも必然的に失敗し、


「ことごと」く、
「必然的」に、
失敗するそうです。。。
0753おしえてください君2006/12/12(火) 18:53:12
>>735です。
ピアソンかスピアマンの選択で悩んでおります

ある測定値Aは測定上限がありませんが、
他方の測定値Bは、ある限界値以上は一定の測定限界値として表示されます。
二つの測定値について論じるために臨床研究しているのですが、
相関係数を出すにはどうしたらよかでしょうか
0754132人目の素数さん2006/12/12(火) 21:06:34
>>753
あんまそんなんで悩むな必要無いぞ。
両方計算して、結果をここにさらしてみ。
■ このスレッドは過去ログ倉庫に格納されています