統計ソフトSTATAの部屋 Ver.2

[0001 ごろ 2010/09/10(金) 10:54:10]

統計ソフト STATAを使って統計解析をするスレです。
日本語の良いソースが少ないので盛り上げましょう。

前スレが順調に!?88まで行ったんですが、落ちてしまったので再開です。

[0002 ごろ 2010/09/10(金) 11:22:57]

前スレで役に立ちそうな内容をまとめておきます。
*リンク集
http://www.ats.ucla.edu/stat/stata/
http://www.stata.com/help.cgi?search
http://www.stata.com/support/statalist/faq

*メーリングリスト
[e-mail]

*ブログ
http://homepage.seesaa.net/category/386689-1.html
http://www.kagitaku.com/
http://statamemo.jugem.jp/
http://d.hatena.ne.jp/sdmt/20061106/1162633144
http://keijisaito.info/arc/stata/japanese_stata.htm Stata日本語化？

*mixi
http://mixi.jp/view_community.pl?id=699582

*PDFなど
http://sitemaker.umich.edu/masahiko.aida/files/outlier.pdf #stataで外れ値の扱い
http://www.stata.com/statalist/archive/2007-06/msg00272.html
http://park1.wakwak.com/~mt_tosiyuki/stata-manual.pdf　日本語PDF

[0003 ごろ 2010/09/13(月) 23:42:18]

前スレのまとめをブログにしようと思っているところ。
出来たらリンク貼ります。
前作ったWikiサイトはサーバーが閉じてしまった。

あと、前スレは88でなくて123で終わっていた。どうでもいいけど。

[0004 １３２人目の素数さん 2010/09/17(金) 20:45:14]

知ってるかい？
このソフトを

[0005 ごろ 2010/09/18(土) 01:30:52]

>>4
確かにマイナーなソフトで日本では知名度は低いですが、この業界では有名だと思います。

[0006 ごろ 2010/09/22(水) 15:52:29]

ちょっと、子供が産まれてバタバタしている。
落ちないように保守。

[0007 １３２人目の素数さん 2010/09/23(木) 02:53:58]

age

[0008 ごろ 2010/09/28(火) 03:04:47]

無事、母子ともに退院しました。
そろそろ、何か有益な情報を書き込めたらといいなぁと思います。

[0009 ごろ 2010/09/29(水) 16:49:54]

何か書き込もうと思いましたが、ネタがなかったのでメーリングリストからのネタです。
FAQらしいですが、個人情報の入ったIDなどから個人情報の入らない新しいIDを作る方法は？というのが質問

回答は、http://www.stata.com/support/faqs/data/group.htmlですが、少し要約

解決案１　１，２，３、・・・、ｎという順のIDを作る
. gen id = _n
解決案２
. egen newid = group(oldid)
これをすると、oldidが文字列でも数字列でも昇順に並び替えて、小さい順に1,2,・・・,nとnewidを作る

別法
. sort oldid
. by oldid: gen newid = 1 if _n==1
. replace newid = sum(newid)
. replace newid = . if missing(oldid)
この処理で、下記のようになっている。
oldid newid (as created) newid (as replaced)
1 1 1
1 . 1
1 . 1
1 . 1
22 1 2
22 . 2
22 . 2
33 1 3
33 . 3
個人的には. replace newid = sum(newid)この使い方が新鮮だった。

[0010 １３２人目の素数さん 2010/09/30(木) 10:37:38]

sage

[0011 １３２人目の素数さん 2010/10/02(土) 04:50:24]

数学板の統計スレここだけになっちゃったね

[0012 ごろ 2010/10/02(土) 11:48:37]

>>11
お、おう、何でも来い！

[0013 １３２人目の素数さん 2010/10/04(月) 21:38:57]

>>ごろ

おめでとう！！
元気でrobustに育ちますように。

[0014 ごろ 2010/10/05(火) 04:49:21]

>>13
ありがとう！何かとってもあったかい気持ちになった。

[0015 １３２人目の素数さん 2010/10/08(金) 18:16:11]

どういう点がいいの？

[0016 １３２人目の素数さん 2010/10/08(金) 18:17:06]

子供の話じゃなく、STATAのことね

[0017 １３２人目の素数さん 2010/10/13(水) 03:50:22]

大学（院）の統計関係の授業でよく使われるSPSSなどの統計パッケージ
に対するアドバンテージってどの辺にあるの？

http://www.lightstone.co.jp/stata/gradplan.htm
にあるGradPlanは魅力的だと思うけど、授業で使ってるのと違うソフト
買う学生は少ないと思うな。

日本の代理店も「初めてのStata 11」くらいなら無料のセミナーとして
潜在ユーザーが触れる機会を多くすればいいのに。

[0018 ごろ 2010/10/14(木) 14:57:02]

>>17
遅くなってゴメンなさい。
大きなアドバンテージは、SPSSなどの高価格統計ソフトと比較すると、
はるかに格安で入手できて、統計処理もそれなりにきっちりできるところでしょうか。

実際、アメリカの大学ではUCLAを初め、統計コースで良く使われているソフトです。

[0019 ごろ 2010/10/14(木) 14:57:55]

まぁ、私自身が、StataとRしかまともに使ったことがないので、比較が難しいというのが本音です。

[0020 １３２人目の素数さん 2010/10/16(土) 10:08:15]

>>17
SPSSとStataを使ってます。最近は少しずつRも勉強してますが、
こちらの方はまだまだ初心者です。

SPSSは多くの企業や教育機関で採用された統計ソフトなので、汎用性があるというか、
規模の経済が利いている点が魅力でしょうか。より内容的な利点としては、
分かりやすいGUIを採用しているので、シンタックスを打ち込まなくても割と直観的に
簡単に分析できる、丁寧で分かりやすい解説書なども出回っていて、初心者に敷居が低い、
機能名などが完全に日本語化されている、日本に支社があるためサポートを得られる、
といったところでしょうか。

それに対してStataは、GUIも採用されていますが、SPSSのようにアイコン化されておらず、
通常はコマンドを打ち込む必要があります。日本ではメジャーなソフトではなく、
日本語の解説書も少ないために初心者にとっては敷居が高いという短所はあります。
しかし、SPSSに比べるとはるかに安く、より複雑な分析やグラフ出力の細かな指定ができ、
英語がある程度読めればサポート、リソース、コマンドはかなり充実しており、
便利なコマンドのアップデートが行われているなどの点が挙げられるでしょう。

[0021 17 2010/10/16(土) 19:47:42]

>>18-20

ご丁寧な説明をありがとうございます。自分でも少し調べてみたの
ですがStataではシンタックスによる操作で説明されていることが多い
ようですね。初心者にはちょっと難しいかも。

でもRでもなくSPSSでもない選択肢としては面白い存在ですね。

[0022 ごろ 2010/10/18(月) 14:44:20]

>>21
私も、Stataを初めて触ったときは、右も左も分からず困りましたが、
適切なインストラクションを受けることが出来て、最初の小さなハードルを超えることが
できれば、とても使いやすいソフトだと思います。

たいていのコマンドは、メニューから選んでできます。
SPSSもそうだと思いますが、メニューで選んで実行したコマンドも、ちゃんと
画面にテキストコマンドとして出力されるので、StataではDoファイルと呼ばれるエディタに
記録しておけば、次回からはメニューから選ばなくても実行できます。

[0023 １３２人目の素数さん 2010/10/19(火) 09:18:48]

たとえば、ここの「データ変換」というところを見るとStata-->SPSSでコレスポンデンス
分析という手順をとっています。
http://hpcgi2.nifty.com/tigerhorse/c-board/c-board.cgi?cmd=all;page=6;id=

ご本人に聞いてみないと分からないことですが、Categories＆SPSSがあるのに、すべてを
SPSSでやるのではなく「できるだけStataで作業したい」と感じさせるほどのことがStataに
はあるのでしょうね。

[0024 ごろ 2010/10/19(火) 15:25:35]

SPSSは知りませんが、Stataは変数名をvar1 var2 var3とか付けておくと
var*で一度に指定できたりとか、そういうところは使いやすいですね。

[0025 ごろ 2010/10/21(木) 03:22:20]

私はRとStataを使っています。まだRは使いこなせないので、本業ではStata、
時間があるときにはRという感じですね。
両方のソフトを使っていて良いことは、統計ソフトの考え方というか、スタンスの違いが分かることです。
その違いを調べることで、統計についてより深く学べている気がします。

今日はその例を紹介します。
間違ったことを書いているかもしれないので、是非つっこんでください。

ANOVAの等分散性について
一般的にANOVA(分散分析)は、多群の分散が等しく、正規分布していることを前提としています。
古典的には、Bartlett、Levene、O'Brien、Brown-Foresythe検定などの手法で当分散性を確認し
てからANOVAが行なわれていました。しかし、最近では、検定の多重性の問題からANOVAを行な
う前に分散が等しいかどうか、正規分布しているかどうかを先に検定することは推奨されなくなっ
てきています。
　そのため等分散性が不明な場合は、等分散を仮定しない検定法（Welch の方法の拡張）を最初か
ら用いるのが良いようです。Stataの標準の設定では、等分散しているとして計算するので、注意が
必要です(もちろんコマンドで等分散性を前提としない計算も出来ます)。Rでは標準で等分散性を仮
定しない検定法を行ないます。

[0026 ごろ 2010/10/21(木) 03:56:00]

＞もちろんコマンドで等分散性を前提としない計算も出来ます
これは簡単でないかも。fstarとかwtestとか導入しないといけないのかな。

[0027 １３２人目の素数さん 2010/10/21(木) 09:55:13]

>>25
>最近では、検定の多重性の問題からANOVAを行な
>う前に分散が等しいかどうか、正規分布しているかどうかを先に検定することは推奨されなくなっ
>てきています。

検定初心者の質問ですが、ある検定の前提となる要件（この場合は等分散）を確認するための検定も、
本来の目的である検定（この場合のANOVA）とともにいわゆる多重検定を構成すると考えるべきなの
でしょうか。群間で同じ検定を繰り返す多重検定とは異なると思うのですが。

[0028 ごろ 2010/10/21(木) 17:10:37]

>>27
ありがとうございます。
私も素人なので、厳密には分かりません。
これまでは、行われてきた方法でもありますし、仰る通り、意味の違う検定なので、
根本的な問題を起こす危険な多重性ではないと思います。
ただ、どちらかの検定が偶然に帰無仮説が棄却されてしまうのは、有意水準を5%ととすると、
1-0.95^2で、9.8%になります。

私の理解では、こういう問題を避けるためにも、最近は古典的な検定から、
モデルを構築する方向へ変わっているのではないかと思っています。
どうなんでしょうか？

[0029 １３２人目の素数さん 2010/10/26(火) 15:28:58]

猫に小判、まで読んだ。

[0030 １３２人目の素数さん 2010/10/27(水) 16:53:43]

猫が寝転んだ。

[0031 ごろ 2010/10/29(金) 05:49:48]

Stataにあって欲しいのにない機能ってありますか？

[0032 １３２人目の素数さん 2010/10/29(金) 22:18:01]

Rとの連携機能？
ってのがJMPにはついたそうで。

[0033 １３２人目の素数さん 2010/10/29(金) 23:15:28]

JMPも面白いソフトだね。

[0034 ごろ 2010/10/29(金) 23:51:10]

>>32
なるほど、最近のRの注目度はスゴイですね。
困るってほどではないけど、Rだと同時にいくつものデータフレームを
扱えるけど、Stataだと一つしか扱えませんね。

[0035 １３２人目の素数さん 2010/10/30(土) 01:21:31]

他の統計ソフトのファイル拡張子での保存。
これだけで随分楽になるんだが。

[0036 １３２人目の素数さん 2010/10/30(土) 05:24:17]

>>34
>Rの注目度はスゴイですね

確かに。ただJMPのような市販パッケージが連携する必要あるのかな。
SPSSとかにお金払いたくない人用の統計ソフトっていうのが自分の捉え方なんだけど。
市販パッケージ使ってて、Rも併用する意味あるの？

[0037 ごろ 2010/10/30(土) 13:18:35]

>>36
私はRも使いますが、Rは統計ソフトというよりは、統計関連のライブラリが充実した
プログラミング言語という感じています。
その点、新しい分野などには逆に強いと思いますが、
（SPSSやJMPは詳しくは知りませんが）市販パッケージの強みは、
洗練された使いやすさと統一されたデザイン（出力結果も含めて）に
あるのではないかと思います。

[0038 ごろ 2010/10/30(土) 14:02:00]

RjpWikiってあるけどStataってWikiないよね。立ち上げるかな。。

[0039 ごろ 2010/11/02(火) 15:48:56]

>>38
サーバーをどこにするかが問題だ。
多分、容量はあまり要らないから長期的に安定して無料か安いとこがいいなぁ。

[0040 １３２人目の素数さん 2010/11/04(木) 11:53:40]

猫が寝転んだ。

[0041 ごろ 2010/11/07(日) 01:15:33]

いろいろあって、急に忙しくなった。自分でレンタルしているサーバーがあるんだが、それを使うか別途契約するか検討中>StataJpWiki

[0042 １３２人目の素数さん 2010/11/13(土) 19:59:28]

Rでデータフレームのヘッダだけコンソールに表示させたい時はどうすれば良い？

[0043 ごろ 2010/11/14(日) 04:42:17]

>>42
Rスレじゃないんですが・・・
データフレームのヘッダって変数名のことかな？
色々方法はあると思うが、例えば
names(データフレーム名)というのはどうでしょう。

[0044 ごろ 2010/11/14(日) 04:43:21]

忙しいはずなのに、朝方までStataJpWikiを立ち上げようと画策していた。
上手くいかないなぁと思ったらPHPに対応していないサーバーだった...orz

[0045 ごろ 2010/11/15(月) 12:47:47]

>>44
忙しいはずなのに、Wiki作っちゃいました。
まだテスト段階なので、テストが終了したら公開します！

[0046 ごろ 2010/11/16(火) 18:19:58]

いかん、いかん
かなりハマってしまった。
URLをもう少しいじりたいけどとりあえず、ここです。
http://statajp.s346.xrea.com/

[0047 １３２人目の素数さん 2010/11/16(火) 22:50:43]

Rじゃだめなんか？

[0048 ごろ 2010/11/16(火) 23:21:03]

>>47
R"も"いいですよ。

[0049 １３２人目の素数さん 2010/11/17(水) 00:09:02]

STATAじゃないと駄目なことって何？
それがないと流行らんでしょ

[0050 ごろ 2010/11/17(水) 13:20:28]

>>49
そんなこと言ったらSPSSじゃないと駄目なことって・・・とも言えそう。
Stataがいい理由はここにどうぞ。
*ttp://statajp.s346.xrea.com/index.php?plugin=related&page=Stata%E3%82%92%E8%96%A6%E3%82%81%E3%82%8B100%2B%E3%81%AE%E7%90%86%E7%94%B1

[0051 ごろ 2010/11/17(水) 13:21:35]

ここだった。
*ttp://statajp.s346.xrea.com/index.php?Stata%E3%82%92%E8%96%A6%E3%82%81%E3%82%8B100%2B%E3%81%AE%E7%90%86%E7%94%B1

[0052 １３２人目の素数さん 2010/11/17(水) 20:33:39]

IP抜かれるの嫌だからここにコピペしろ

[0053 １３２人目の素数さん 2010/11/18(木) 00:57:50]

トライアル版なんかないんすかね

[0054 １３２人目の素数さん 2010/11/18(木) 12:08:32]

いまさら市販ソフトの宣伝スレ？

[0055 １３２人目の素数さん 2010/11/18(木) 14:49:00]

「市販」アレルギーこそ今更だと思うが

[0056 ごろ 2010/11/18(木) 15:32:50]

だんだん盛り上がってきましたね。
Stataがいいのは、やっぱり簡単に使えて、qualityの高いグラフが作図できるところかなと思います。
日本語版がなくて、最初はとっつきにくいんですが、最初のハードルを超えれば使いやすいです。

特に、システマティック・レビュー用のadoファイルは秀逸です。
本当に出版レベルの図が簡単に作成できました。

Rも面白いんですが、簡単に結果がほしい時には、やっぱりStataを使ってしまいます。

[0057 １３２人目の素数さん 2010/11/18(木) 20:15:25]

安いので幾らぐらい？

[0058 ごろ 2010/11/18(木) 23:36:05]

学生向けだとこれくらい。
*ttp://www.lightstone.co.jp/stata/gradplan.htm
私はアメリカで買ったのでもっと安かったです。

普通に買うと、結構高いなぁ。

[0059 １３２人目の素数さん 2010/11/19(金) 02:18:42]

だめだコリャ

[0060 １３２人目の素数さん 2010/11/23(火) 12:39:13]

60で終了
STATAなんてそんなもんだろ

[0061 ごろ 2010/11/23(火) 17:06:15]

>>60
いやいや、これからです。

[0062 １３２人目の素数さん 2010/11/23(火) 17:32:18]

>>5で「この業界」と言っておられるが、何業界？

[0063 １３２人目の素数さん 2010/11/24(水) 01:04:02]

>>62
アカデミックじゃね？

[0064 ごろ 2010/11/24(水) 11:42:00]

>>62
そうですね。あまり深い意味で使ったわけではありませんが、

>Stata（すたーた）は、1985年、StataCorp社により開発された統計分析のソフトである。主に経済学、社会学、政治学、疫学の分野で用いられている。by Wikipedia
という分野を想定しています。

[0065 １３２人目の素数さん 2010/11/24(水) 20:03:23]

>>64
>という分野

だとSPSSなんかとかぶりまくりだな。

[0066 ごろ 2010/11/25(木) 15:38:46]

>>65
まぁそうですよね。
というか、統計ソフト自体がテキストエディタソフトのようにごろごろあるわけじゃないので、
有名な統計ソフトを５本挙げよ、と言われたらStataも入ってくると思います。

[0067 １３２人目の素数さん 2010/11/25(木) 18:02:41]

研究に有用なソフトウェア-統計ソフト
http://www.kenkyuu.net/comp-soft-01.html
医薬系なんで一般論にはならないと思うが、Prismが1位

[0068 １３２人目の素数さん 2010/11/25(木) 20:19:19]

>>67
SASが思ったより少ないのと、Excelが結構多いのが意外。
本郷の生物統計では今でもSASやらされると思う。

[0069 １３２人目の素数さん 2010/11/25(木) 23:04:23]

>>68
Ｏ橋先生のとこ？
あの人ホントに数学分かってる？

[0070 １３２人目の素数さん 2010/11/26(金) 00:28:20]

＞ 有名な統計ソフトを５本挙げよ、

R
Matlab
Mathematica
SAS
SPSS

STATAは絶対上位５位に入りません！
STATAをいいという人は、STATAぐらいしか使った
ことのない人か、金持ち大学の人です。

STATAを導入していない貧乏大学では、誰も知りません。
名前さえ知りません。

東大は予算が余りまくりの、予算を無理しても使い切れ！
というところなので、性能に関係なく金額の高いソフトを
導入するのです。

[0071 ごろ 2010/11/26(金) 15:06:03]

ちょっとググッてみました。

MATLAB
正式版はさすがに個人では買える額ではないけど，
学生版なら10,000円程度で買えるよ． ただし，使える行列サイズに制限があります．
Mathematica
学生優待価格の「学生版」があり，生協に注文すれば30,450円で購入できる．アカデミック版が238,350円であることを考えれば格安といえる
SASは個人では買えない？
SPSS
SPSS Base System 15.0
一般価格　　：176,400
アカデミック：102,900

確かに安いとは言いませんが、Stataは他のソフトと比べると比較的、割安感があると思います。

[0072 １３２人目の素数さん 2010/11/26(金) 19:44:37]

Rが流行るのは高機能（特にグラフィックス）で無料だから
SASが流行るのはお役所公認だから
あとは別にどうでもよい

[0073 １３２人目の素数さん 2010/11/27(土) 10:25:34]

sagemath とかどうなんですか？
http://www.sagemath.org/

[0074 １３２人目の素数さん 2010/11/27(土) 12:17:39]

>>70
>性能に関係なく金額の高いソフトを導入するのです。

そっちの人w？と話した時は性能の点でSASしかないという言い方してた。
上の学年の人の受け売りぽかったけど。

[0075 １３２人目の素数さん 2010/11/27(土) 15:58:43]

コマンドで推定するのでなくて、STATAのプログラミングを基礎から学べるよい本はありませんか？

また、データの加工をSTATAで行いたいのですが、データの行数？が多いとMATAでないと扱えないみたいなのですが、MATAの扱いってSTATA公式のマニュアル以外でなにか分かりよいものってないでしょうか。

社会人になってver10買ってみたものの、ほとんど使ってないので最近仕事で使えるように勉強してみようかと思っています。

[0076 ごろ 2010/11/27(土) 17:03:25]

>>75
質問のレベルが私のレベルを越えているようですが、頑張って回答に挑戦します。
プログラミングについて
やはりこの辺りでしょうか？
http://www.stata.com/bookstore/isp.html
私はStataのちゃんとしたマニュアルを買ったことがありません＿|￣|○
必要なときに、必要な部分をネット検索で済ませてきてしまったからです。

データ加工について
データの量はどれくらいですか？
http://www.stata.com/products/whichstata.html
ここで言うデータ量（変数の数）を超えていますか？

MATAについては知りません。MATAって何ですか？済みません。

[0077 75 2010/11/27(土) 20:02:22]

>>76　ごろさん早速ありがとうございます！

http://www.stata.com/bookstore/isp.html
　↑タイトルはまさしくぴったりです。早速購入してみます。
　STATAだけでなく、プログラミング自体にまったく素養がないので、これを機会に・・・

データ加工ですが、言葉足らずでした・・・
データの行や列を分解したりくっつけたりして、コンディショナル・ロジット分析用のデータをつくる必要がありまして、
いままではエクセルでごにょごにょやってたのですが・・・データが多いとエクセルではきつくて・・・
ある方に聞いたところGAUSSという行列扱うソフトで加工している、との話でしたが私は持っていないので、

どうせ推定にSTATAつかうなら、いっそ加工もSTATAの行列計算の機能をつかってできないか！と思い調べていたら、
どうもSTATAの通常の行列コマンドでは扱える行や列のサイズに制限があり、
（STATAのオンラインヘルプ　「help matsize」　に書いてあります）
それ以上はMATAというSTATAの新しい行列をあつかうやり方だとできそう。。
とういうところまで発見していま止まっています。(;´д｀)

そのやり方考えるのついでに、プログラミングの知識がないので、
STATAでプログラミング学ぼうかなと・・・無謀かな・・

お返事ありがとうございました。

[0078 １３２人目の素数さん 2010/11/28(日) 13:19:15]

>>76

clogitなんてしばらく使ってないから忘れてしまったけど、
最初はcase2altコマンドでデータをスタックしてたんだけど、
これじゃぁうまくいかないの？基本的に、

. case2alt, choice(y) casevars(X1 X2 X3), gen(Y)

だとオレのあんちょこは言っている。

[0079 ごろ 2010/11/28(日) 13:38:31]

>>77
Stataには、Mataなんてモードがあったんですね。
私の勉強不足でお恥ずかしい限りです。
Stataからhelpで参照できるMataのマニュアルも結構詳しそうですね。

Mataとmatsizeについては
http://www.stata.com/support/faqs/mata/matsize.html
このリンクにあるとおり，matsizeのリミットを超えて扱えるようですね。

私は普段、matsizeを超えるようなデータを扱わないのですが、
もし宜しければ一緒に勉強したいと思います。
Mataで大きな行列を扱う方法を見つけたらまた書き込みます。

[0080 75 2010/11/30(火) 01:35:36]

>>78

case2altコマンド あ、こんなものがあったとは・・・
J. Scott Longという人の
Regression Models for Categorical Dependent Variables Using Stata
http://www.stata.com/bookstore/regmodcdvs.html

のP295に書いてありますね。本棚でほこりかぶってたの今発見しました　orz
ありがとうございます。

>>79　ごろさん
　私はまったくの初心者＆統計も専門ではない文系人間なので、ぼちぼち
としか進めませんが、いろいろ勉強させていただきたいです。
ではまたのぞきにきます

[0081 ごろ 2010/12/04(土) 06:34:01]

>>75
http://www.stata.com/netcourse/nc151.html
これもいいかも

[0082 ごろ 2010/12/10(金) 03:00:33]

>>81
値段も手頃ですし、受けて感想を報告してみます。

[0083 75 2010/12/11(土) 00:53:36]

>>82

ver11じゃないとだめみないなことが書いてありますね。


Prerequisites:
Stata 11, installed and working
や
http://www.stata.com/netcourse/enrollment.html
みると。

ネットのレクチャー、受けてみたいですが・・・
まあ、英語というところで少し厳しい

[0084 ごろ 2010/12/11(土) 08:45:55]

>>83
まだ申し込んでいませんが、人柱として受講して、可能な範囲で情報を載せたいと思います。

[0085 １３２人目の素数さん 2010/12/11(土) 18:01:51]

いきなりSTATAの代理店から売り込みの電話が来るようになったね
このスレも営利目的臭がプンプンする
ウザイ代理店営業の仕返しに２４しておこっと

[0086 １３２人目の素数さん 2010/12/12(日) 20:09:20]

このスレ見ると電話番号が代理店にバレるわけですね

[0087 １３２人目の素数さん 2010/12/12(日) 21:33:09]

>>85-86

SPSS業者おつ。

[0088 １３２人目の素数さん 2010/12/13(月) 09:39:19]

？

SPSSって売ってもらうの面倒だった。
パッケージで即購入できるのかと思ったら違うのね。

[0089 75 2010/12/13(月) 20:35:03]

以外とみてる人いるのですね(；・∀・)

>>84 ごろさん　無理ない範囲で人柱お願いします（矛盾・・・

そういえばおとといSTATA　NEWSがエアメールで来たが、スウェーデンからになってました。
なんでだろ・・・

[0090 ごろ 2010/12/14(火) 07:13:55]

>>89
見てくれている人が多いのに、あまり有益な情報が提供できずに残念です。
ネットレクチャーのプログラミング初級を申し込みました。
Stata１１のシリアルナンバーがないと申し込めないみたいですね。
11でなくても普通に使えるのに、こういうところで、制限をかけてしまうのは、
ユーザーに反感を買うので、ちょっと残念ですね。

コンディショナル・ロジット分析は経済学では有名な分析手法のようですね。
ちょっと勉強してみましたが、まだ雰囲気しか分かっていません。
何か、小さなデータを手に入れて実際にStataに分析させれると身につくと思うのですが、
そこまで手が回っていません。

MATAの方は、もともとのマニュアルもしっかりしているので自分で勉強できそうですが、
まだMATAを使わないといけない状況がないので、学習意欲が湧いてきていません。

[0091 75 2010/12/15(水) 23:32:10]

私も素人なので話題の提供なかなかできないですが、話題としてはSTATA中心がよいでしょうが、
もっと統計学や多変量解析、データマイニングなんかの
話もあってもよい気が個人的にはします。ほかにスレあれば別ですが・・

コンディショナルロジットは経済学では企業の立地の分析に使いますね。
あと、嗜好や中毒、行動経済学の実証にも使います。
マーケティングだと、コンジョイント分析という、商品の望ましい仕様の研究で定番です。

私もなんであの式でいいのか証明は追ってないのですが、こんなに簡単な尤度関数になるのか！
と思いました。市場シェアのシミュレーションみたいなことすると面白いです。
100円価格上げるとシェアがどんなふうに変化するのかなどみられるので、そのあたりを仕事で
ちょこちょこ自主的な実験レベルですがやったことがあります。

[0092 ごろ 2010/12/16(木) 14:07:57]

>>91
ありがとうございます。
もうちょっと勉強してみます。

[0093 １３２人目の素数さん 2010/12/18(土) 20:23:25]

>>88
そう？
AMOSだけしか持ってないが、即購入したような記憶があるけど。

[0094 ごろ 2010/12/23(木) 12:39:09]

本来のスレの趣旨に戻って、Stataの基本的な使い方を書こう。

最初、私が何も知らないでStataを触った時には、やはり英語ということと、
どのコマンドがどこにあるのかがさっぱりで、苦労しました。

実際は数時間レベルのレクチャーを受けるか、周りによく知っている人がいれば、
すぐに使えるようになると思いますが、何も知らないで使い始めると戸惑うと思います。

最初は、Statistics>Summaries, tables, and tests>この中だけで、古典的な手法はほとんど
カバーされているので、ここをいじってみると良いと思います。
あとは、メニューから実行するときに、吐き出されるコマンドを覚えたり、Doファイルに
メモしていくうちに自然と覚えていくと思います。

[0095 １３２人目の素数さん 2010/12/25(土) 01:03:47]

明石家サンタでも見るかな。

[0096 ごろ 2010/12/26(日) 10:52:05]

>>95
ちょっと遅いけど、メリークリスマスです。

[0097 １３２人目の素数さん 2010/12/27(月) 16:08:44]

ちょっと早いけどあけましておめでとう。
STATAはHLM、いわゆるマルチレベルモデルは扱えますか

[0098 ごろ 2010/12/29(水) 15:00:49]

http://www.stata.com/statalist/archive/2007-01/msg00260.html
ここに参考になる質問がありました。
gllammかxtmixedというコマンドでできるようです。
gllammはstata10の標準パッケージには入っていませんが、無料でインストールできます。

gllamm estimates Generalized Linear Latent And Mixed Models.
These models include multilevel (hierarchical) regression models
with an arbitrary number of levels, generalized linear mixed models,
multilevel factor models and some types of latent class models.
We refer to the random effects (random intercepts, slopes or coefficients),
factors, etc. as latent variables or random effects.

If the latent variables are assumed to be multivariate normal,
gllamm uses Gauss-Hermite quadrature, or adaptive quadrature
if the adapt option is also specified. Adaptive quadrature
can be considerably more accurate than ordinary quadrature,
see first reference at the bottom of this help file.

With the ip(f) option, the latent variables are specified
as discrete with freely estimated probabilities (masses) and locations.

More information on the models is available from
http://www.gllamm.org

私の興味ある分野なのでもっと勉強したいと思っています。

[0099 １３２人目の素数さん 2011/01/01(土) 15:36:05]

>>97
できますよ。私はよく使っています。

>>98でごろさんが既に答えているけど、gllammやxt-あるいはxtme-コマンドを使います。
テキストとしては、

Rabe-Hesketh, Sophia and Anders Skrondal. 2008. Multilevel and Longitudinal Modeling Using Stata. College Station, TX: Stata Press.

があります。

xt-やxtme-はそこそこ動いてくれるのですが、gllammはシャレにならないほど
結果の出力に時間がかかるので注意が必要です。
Nやモデルの内容、マシンパワーにもよりますが、3レベル以上になると、
さすがにgllammで何度も結果を出力させるのは無理があると思います。

[0100 ごろ 2011/01/01(土) 16:46:19]

>>99
新年おめでとうございます。今年も宜しくお願いいたします。
補足ありがとうございます。
もっと勉強したいと思います。

[0101 １３２人目の素数さん 2011/01/02(日) 12:49:45]

>>99
googleブックスでその本が読めるようですね。
勉強します。

[0102 １３２人目の素数さん 2011/01/09(日) 17:38:30]

【R言語】統計解析フリーソフトＲ 第４章【GNU R】
http://kamome.2ch.net/test/read.cgi/math/1294561909/

[0103 １３２人目の素数さん 2011/01/10(月) 15:21:05]

遅くなったけど、新年おめでとう！
今年もよろしく！

[0104 ごろ 2011/01/11(火) 00:28:21]

>>103
今年も宜しくお願いいたします。
>>102
ついにRスレも復活したんですね。

今、R用のレクチャーノートを作っている。
Stataにはlabelという概念があって、変数の意味をlabelに書くことが出来る。
表を作った時などに、labelがあると、変数名の変わりにlabelが表示される。

Rにはこれと全く同じ概念のものはないが、comment関数で変数を含むオブジェクト全てに
commentを書きこむことが出来る。Stataみたいに、作表時に自動的にタイトルになるという
ことはないが、きっと上手く指定してやればStataと同様なことは出来るはず。

[0105 ごろ 2011/01/20(木) 06:00:08]

Thank you for enrolling in NetCourse 151.

The course will begin January 21, 2011. You will receive the preliminary
information on the course, including the URL, password, and username,
via email on the starting day of the course.

なんか１月もあっという間だな。もうすぐNet course始まる。

[0106 ごろ 2011/01/22(土) 09:39:11]

ネットコース始まった。

さっそくいくつか学んだ。
Stata内でも簡単なUnixコマンドが使えること。
cdとかmkdirとか

レクチャー１は全１７ステップ
２まで終わったところ。
最初のプログラムはやっぱりHello Worldを表示させることだった。

[0107 ごろ 2011/01/22(土) 15:49:01]

全体17の４まで読んだ
出てきたコマンドは
program
which
set trace on
set more off
set more on
set trace off
など

結構基礎から丁寧に教えてくれていて分かりやすい。
doファイルでコマンドの最後に改行を入れ忘れると駄目など、
結構落とし穴も注意してくれる。

楽しみ。

[0108 ごろ 2011/01/22(土) 19:00:16]

今ステップ７のadoファイルを作る所まで来た。ここまでは簡単。
でも次から急にややこしくなってきたので休憩。

レクチャーは全部で４回みたいだ。火曜日までに宿題や講義への質問かコメントを書かないと
いけないようだ。

[0109 ごろ 2011/01/25(火) 07:23:03]

あれからあまり進んでいないが、
ステップ８−１０は、結局沢山のdoファイルを混乱しないようにオーガナイズ
する方法の紹介だった。
例えば、データセットを作るためのdoファイルとデータ解析をするための
doファイルの名前を先頭の数文字で区別できるようにするなど。
ここで勉強になったのは、doファイルのなかに、do doファイル名と書きこむ
ことで、doファイルから別のdoファイルを呼び出すことができること。
ステップ11は、infileコマンドの解説だが、ややこしそう。

ちなみにLecture１はWordにコピペしてみたら５４ページ１３０００単語を超えていた。
残り２９ページ７０００単語。

[0110 ごろ 2011/01/26(水) 03:55:57]

ステップ１１はinfileの使い方、dictionary formatの説明、階層データの
読み込み方。
巨大データを非力なコンピュータで扱う方法
（日中はデータサイズを縮小して仕事をして、寝ている間に全データで作業する）
新しいコマンド
sample 10 /*データから１０％のデータに間引きする*/
ステップ１２は再現性について
Stataの乱数について
コマンド　
set seed 数字
sort 変数名, stable <-このstableオプション
ステップ１３は添字
変数名[_n]
[_n+1],[_N]などの使い方
gen y=x[_N-_n+1]　/*これはxの変数を逆順にｙに入力*/
ステップ１４
Stataは条件分岐とループが遅いから出来るだけ避けるべきとのこと
assertの進んだ使い方。このassertで避けられるそうだが、理解不十分
ステップ１５
summarize(x)の後 r(mean)で結果を取り出す方法など
ステップ１６
まとめ
ステップ１７は確認のための練習問題

[0111 ごろ 2011/01/26(水) 06:06:36]

実は、Stata11は持っているんですが、実習はStata10でやっていました。
これまでほとんど問題なく出来ていましたが、
mergeコマンドで問題発生

stata11からは
merge m:1 ・・・
というのが使えるがstata10では使えなかった。

[0112 ごろ 2011/01/26(水) 13:40:40]

英語力の問題もあるんだろうけど、Lec1の演習問題はあんまり面白くないなぁ。
今日正解発表があった。

一応、提出義務はないようだ。というか、送ってくるなと書いてある。

[0113 ごろ 2011/01/28(金) 09:12:40]

Lecture1は講義内容は８−９割以上理解できた。
演習問題はまだ６−７割しか理解できていないので後でもう一度。

すでにLecture2の講義がUpされている。

[0114 １３２人目の素数さん 2011/01/28(金) 18:35:09]

いくら何でも日記はいらんよ
自分が立てたからといっても、スレはあんたのものではない

[0115 75 2011/01/30(日) 22:34:48]

STATAレクチャーよさそうですね。やはり、講義のよさってありますよね（強制的に進んでいくとか・・・）

いま、ちょっとわけありで直交計画がSTATAで出せないかと探していましたが、
可能かどうかがいまいちわからず・・・

Rで、oa.design()という関数でできるということが分かりましたが、できればSTATAで
やりたいので、ご存じの方いましたら教えていただきたく。

[0116 ごろ 2011/02/01(火) 18:53:28]

>>115
直交計画という分析方法もあるんですね。勉強になります。
調べてみたらorthogという関数がStataにあるようです。
help orthogonal
詳しくは、で見てみてください。

Rに、oa.design()という関数はありますか？
RjWikiで調べたところS-Plusにはこの関数はあるが・・・という記載は見つけました。
Rでは、パッケージは何を使うのでしょうか？

Stataレクチャーの方は、ざっと2回目までの講義を見たところ、プログラミングの基本というだけ
あって、かなり基礎的な内容の印象です。
2回目のレクチャーの内容の詳細はまたおいおいUpしたいと思いますが、とりあえず、１週間の休み
があるので、のんびりしています。

しかし、知っていると思っていたことでも新しい気付きがあったりして、受講して良かったと思っています。
内容的にはStata１１でなくても十分なので、Stata１１所有者という縛りはない方がいいのにと思います。

Stataプログラミング関連で何か質問があれば、講義に直接関係なくても聞くことができるかもしれません。
もし、代わりに聞いて欲しい質問があれば教えて下さい。

[0117 75 2011/02/03(木) 00:03:24]

>>116
Rの関数　oa.design　は、DoE.baseというパッケージに含まれています。
意外と知られてないんでしょうか・・・といっても私もたまたま見つけたのですが。

STATA　orthog　ありがとうございます。早速しらべます。

[0118 Stataの為にパソコンを買い替えようとしてるアホンダラ 2011/02/04(金) 10:49:44]

最近、古くなってきたパソコンを買い替えようと思っているんですが、
折角なので、Stataの演算速度やらデータの容量やらについてちょっと考えてみたいと思い、
StataのFaqページなどを調べつつ、カキコすることにしました。

需要があるかどうかは分かりませんが、もともと過疎ってるスレですし（失礼ｗ）
分析方法に比べると、この手のトピックは割と手薄だと思うので、
>>114さんは嫌がるでしょうが、自分の日記代わりに書かせていただきます。
そうです。私は、人の嫌がることでも自ら進んでしますｗｗｗ

ちなみに、パソコンに関してはほぼド素人に近い人間ですので、
ここを見てる多くの方には当たり前の情報ばかりでかもしれませんが、
自分も勉強するつもりで頑張りますんで、よろしくお願いします。
主に、「Stataが遅い！」　「Stataを効率よく使うのにはどんなパソコンがいいのか知りたい！」
と思ってる希有な人々、あるいはパソコンを買い替えるときには役に立つかもしれません。
Windows OSを利用しています。MacやLinuxは分からんのですが、
多分内容的にはほとんど同じだと思われます。

[0119 Stataの為にパソコンを買い替えようとしてるアホンダラ 2011/02/04(金) 11:45:49]

さて、最初の疑問に思ったのは、Stataの処理を素早く行うには、
パソコンのどの部分をどうしたらよいかとということです。
当然ですが、コンピューターのデータ処理には、
CPUとメモリが大きく関わってきます。

まずCPUですが、これはコンピューターの頭のよさに該当する場所で、
これが早ければ当然、演算速度も速くなります。
Stataからは、「CPUクロックとStataのパフォーマンスの間には線的な関係がある」
というレポートが出ています。

ここで「CPUクロック」値というのは、Intel Core 2 DuoとかAMD Phenom IIとかの、
CPUの製品名の後に書いてある、「2.3GHz」やら、「3.0GHz」やらの数値のことを指しています。
つまり、他の条件を一定にすれば、この値が大きければ大きいほど、
Stataの計算・処理速度も速くなる、ということだそうです。
（詳細：http://www.stata.com/support/faqs/win/specs.html）

ところで、3-4年前のCPUと、現在販売されている最新型のCPUの
クロック値だけ比べてみると、同じような値であることに気付きます。
例えば、2006年頃に発売されたIntelのCore 2 DuoというCPUの、
一般ユーザー向けCPUクロック値は、およそ1.3GHz〜2.6GHz
といったところでした。これに対して、現在売られている
最新のIntel Core i5あるいはCore i7シリーズの標準的なCPUクロック値は、
下は1.8GHzくらいから上はせいぜい3.4GHzで、
その値を見る限りは、代わり映えがしないように見えます。

しかし、クロック値自体は変わらなくとも、一般的な処理スピード、
とりわけマルチタスク（複数のプログラムを動かすこと）
への対処速度などは、間違いなく上がっています。
これは、�@搭載されているCPUの数が2つから3つ、4つ、あるいはそれ以上に増えてきていること、
�ACPUの中身の構造（アーキテクチャ）や複数のCPU同士で仕事を効率よく
振り分ける機能（スレッディング等）が改善されたことが原因のようです。

[0120 Stataの為にパソコンを買い替えようとしてるアホンダラ 2011/02/04(金) 11:48:29]

「そんなパソコン話はどうでもいいから、�@と�Aのお陰でStataは速くなるんかい！？」
というのが知りたいところでしょうが、Stata MPを使っていれば、これらの技術的進化の
恩恵を受けることができるようです。Stata SEでは、1つのCPUのみを用いる設定
になっているそうですが、Stata MPはCPUの数に対応して処理を効率的に行えるようで、
テストでは、Stata MPを1つのCPUで走らせたときと比べると、2つのCPUのときには
およそ1.6倍ほどの速さになることもあると報告されています。（データの内容、分析方法などによる）
（詳細：http://www.stata.com/statamp/）
（より詳細：http://stata.com/statamp/report.pdf）

更に、これは自己責任になりますが、規定のCPUクロック値を上回る働きをさせる機能
（オーバークロック）を利用すれば、実装されたクロック値以上の処理速度を期待できるかもしれません。
最近では、そのような機能を行わせることを前提としたCPUも売り出されているようです。
（つい最近不良欠陥が見つかって、回収騒ぎになってしまいますが。）
分析データが大きい場合や、分析に時間がかかるモデルを利用している人は、
このような機能を使ったりすることで処理速度を速くできるかもしれません。

少し書きすぎました。
次は気が向いたらメモリの話を書きます。

[0121 ごろ 2011/02/05(土) 18:36:10]

>>118-120
ありがとうございます。私はあまり巨大なデータを扱わないので、
CPU数による処理速度の違いとかを気にしたことはありませんが、
確かにStataのパッケージの違いで、CPU数がいくつまで対応しているかという
のは見た記憶があります。

また、NetCourseのレクチャー１でも、巨大データを扱う方法で、昼間は、
間引いたデータで行い、最終的な処理は、夜中寝ている間にやらせるという
Tipsがあり、興味深く読みました。

さて、日記の続きを書きます。
今日は、NetCourseのレクチャー２です。全部で１５ステップあります。
ステップ１
レクチャー２用のパッケージ導入
StataとRしか知りませんが、最近の統計ソフトは、コマンドを打つだけで、
勝手にネット上から、拡張パッケージをインストールできるのが良いですね。

2はマクロについて
C言語のマクロと非常に似ているなと思う。
意識せずに使っていたが、
local i=1 という命令はマクロを定義していたのだと、気がついた。
C言語ライクに、=を除いて、
local i 1　と書いても同じ。

3はマクロの具体的な使い方について
とりあえず今日の日記はここまで。

[0122 ごろ 2011/02/08(火) 05:59:36]

ステップ3/15からもう一度
マクロ定義で=を使ってもいいし、使わなくてもいいという違いがあるが、
これには、重要な違いがあるということ。まだ理由は説明されていない。

#delimit ;命令を使って、文の終わりを改行ではなくて、;で終わりにするように変更する方法。

keep/drop関数の使い方

マクロで間違えやすいことの例
local x = -1
di `x'^2
di (`x')^2

summarize 変数名, detail
の結果をr(p5),r(p95)で参照する方法

ステップ4/15
これは知らなかった。
do doファイル名 マクロ名１ マクロ名２
の命令で、doファイル内に書いた`1',`2'にそれぞれマクロ名１、マクロ名２を代入できる。
今日はここまで。

[0123 Stataの為にパソコンを買い替えようとしてるアホンダラ 2011/02/08(火) 15:32:14]

CPUが頭の良さを表すなら、メモリは作業をする机の広さを表すそうなんですが、
これもStataの処理速度と機能に関わる重要な要素です。メモリは「システムメモリ」とか、「メインメモリ」、
あるいは「RAM (Random Access Memory)」と呼ばれている部分で、1GBとか、2GBとかの数値がその大きさを表しています。

さて、作業机は大きい方が良いので、1GBよりも2GB、2GBよりも3GBのメモリを持つコンピューターの方が
もちろん良いわけですが、Stataがコンピューターにあるメモリのうち、どのくらいを利用するかは、
Stata自身で決めることになります。比較的大きなデータを扱ったことがある人はご存じでしょうが、

. query

と打ち込むと、最初の方に

Memory settings (also see query memory for a more complete report)
set memory 100M

と表示され、上の例では、コンピューターが元々持っているメモリ領域のうち、
100MBがStataに割り当てられ、Stataはそのスペースにデータを広げていることになります。

基本的に、仮想メモリを使わない限り、この値以上のデータを広げることはできないので、
例えば、上の例で言えば、100MBが割り当てられているときに200MBのデータを開けようとすると、

"no room to add more observations"

云々というエラーが返されます。そこで、

. set memory 300M

と割り当てる了を300MBに増やしてあげれば、それ以下の量のデータを読むことができますし、
その割り当てを常に保持しておきたければ、後に", permanently"のオプションを打ち込めば、
次回起動時にその設定が保持されます。

[0124 Stataの為にパソコンを買い替えようとしてるアホンダラ 2011/02/08(火) 15:35:06]

さて、Stataにどのくらいのメモリを割り当てられるかは、コンピューターによって制限がかかっており、
例えば、もともと1GBのメモリーしかないコンピューターに、

. set memory 2GB

と指定しても、

"op. sys. refuses to provide memory"
r(909)

のようなエラーが返ってきます。
以前のバージョンでは、Stata起動時に"/k"というオプションを使うと、
仮想メモリー（机の上に臨時の机を用意するようなもの）にもデータを割り当てられることができたため、
このオプションを使えば大きなデータも読み込めたようですが、
現在のバージョン（Stata 11）で同様の方法は確認できませんでした。
（参考：http://www.stata.com/support/faqs/data/large.html
　　　　http://www.stata.com/support/faqs/win/memory.html）

[0125 Stataの為にパソコンを買い替えようとしてるアホンダラ 2011/02/08(火) 15:37:56]

32ビットのWindowsシステムでは、仮にコンピューターが2GBのメモリを積んでいても、
他のプログラムが利用しているメモリの割り当てなどによって、
Stataには最高でも1GB程度のメモリしか割り当てられないことがあるそうです。
（参考：http://www.stata.com/support/faqs/win/winmemory.html）
Stataでは、Windowsで1GB以上の大きなデータを扱う場合は、
64ビットのWindowsシステムを使うことを推奨しています。
（参考：http://www.stata.com/support/faqs/win/misc.html）

また、これは仕方がないことなのですが、Stataが計算のために物理メモリだけでなく、
仮想メモリを必要とするとき、（ハードディスクに多くアクセスしている場合は、たいていそう）
Stataの処理速度がかなり遅くなることがあります。これを回避するためには、
コンピュータにより大きなメモリを導入するか、増設する必要があるようです。
（http://www.stata.com/support/faqs/win/misc.html）

[0126 ごろ 2011/02/09(水) 04:25:08]

レクチャー２
ステップ5/15
レクチャー１で習ったdictionary format形式のデータベースからの
データの取り込みをマクロを使ってもっと華麗に行うという実例。
一応、眺めて理解したつもりだが、自分で作れと言われたら出来ないだろうな。
ステップ6/15
マクロを使ってデータの取り込みをする別の例
この例を学んで、Doファイルの中に
> do doファイル名　マクロ名
というコマンドを書きこむメリットが良く理解できた。
簡単に説明すると、
同じようなデータ構造をした複数のテキストデータからデータを取り込む際、
一つのdoファイルに一つのテキストデータからの具体的なデータの取り込み手順
を記載しておく。ただし、テキストデータのファイル名は`1'としておく。
すると、別のdoファイルに
> do 今作成したdoファイル名 最初のテキストデータファイル名
> do 今作成したdoファイル名 ２番目のテキストデータファイル名
とdo コマンドを列記するだけで、複数のテキストデータからのデータ取り込み
が完了する。

次の例は、上記のように別のdoファイルに記述した内容を
programコマンドを使って一つのコマンドにして、
一つのdoファイル内で完結する方法

[0127 ごろ 2011/02/09(水) 04:27:32]

ステップ7/15
ステップ8/15
ここはちょっとややこしいから、飛ばして後回し
ステップ9/15
if, forvalues, while関数の紹介
> local i = 1
> while `i' <= 10 {
> display `i'
> local i = `i' + 1
> }
など
ステップ10/15
プログラムの書き方のスタイルについての話し
Stataは{}の使い方を注意しないといけない。
例えば
if `i'==1 {
display "yes"
} else {
display "no"
}
これだと} else {ここでエラーが起きる
今日はここまで。

[0128 ごろ 2011/02/15(火) 14:21:23]

また日記だけど、
のんびりしているうちに既に３回目のレクチャーのDiscussion期限が来てしまった。
追いつかないと・・・

[0129 ごろ 2011/02/18(金) 07:25:20]

だいぶ遅れている。
レクチャー２の残りはステップ７・８と１１−１５
理解しやすいところから終わらせる。
ステップ１２
forvaluesコマンドの使い方
ここはそれほど新しい話はなし。
local N = _Nとする使い方と
displayコマンドで　_column() と _skip()　コマンドを使う方法がためになった。

続きはまたあとで。
lecture3は日付を扱う関数がでてくる。
stataは02/14/2011という文字列をdate()関数で日付と認識できるが、
02/14/11という文字列からはそのままでは、日付として認識できない。
これに関して質問したら、非常に有用な回答をもらった。
解決方法
date("02/14/10","MD20Y")　あるいは、date("02/14/10","MDY",2011)
特に後者の方法は上手くつかえば、1999とか2010とか異なる年代が混在する場合でも
上手くはたらく。
やはり、このコースのメリットはその道の専門家に質問が出来、回答を得られることだと認識した。

[0130 ごろ 2011/02/20(日) 10:51:09.96]

時間ができたのでちょっと頑張った。

レクチャー２　ステップ7/15
ややこしくて飛ばしたステップ７に戻る。
ここは、localで定義したマクロが及ぶ範囲についての話。
例えば、doファイルの最初に定義したlocalマクロは、その後に出てくるプログラム内には影響を及ぼさないなど。
globalマクロの使い方も少し出てくるが、危険なので、NetCourse152までは出てこないということ。

ステップ8/15
local マクロ名 1 とlocal マクロ名 = 1 の違いなどについての詳細な説明。
=を使うとコンパイラーが起動されるなど動作が多少遅くなる可能性があるようだが、通常はあまり気にしなくて良さそう。

ステップ11/15
foreachの使い方
>do doファイル名　マクロ名１　マクロ名２
このときlocal 0でマクロ名が全て読み込める。
だからforeach var of local 0 {で、マクロ名リストを読み込んで一つずつ実行できるようになる。

[0131 ごろ 2011/02/20(日) 15:21:17.19]

ちょっと頑張っておいつかないと。
ステップ13/15
プログラムとしてのif文の使い方。
特に変わった使い方はなさそう。

レクチャー３
ステップ1/15　レクチャー３の準備
ステップ2/15　複雑なデータの紹介とこれまでの方法でのデータ取り込み方法
ステップ3/15　clock関数の使い方
ステップ4/15
以前にも出てきたassertを使ってデータを確認する方法。
mod()関数を使ったテクニックの紹介
by id: drop if ...というのは結構使えそうな気がする。
ステップ5/15
summarize関数などを使ったときに自動的に保存される値の扱い方。
return list, ereturn list
またプログラムの中で使うと少しややこしい。
キーワードはprogram プログラム名, rclass
ステップ5.1/15
プログラムの中で
return scalar
return local　/* マクロを返す*/
return matrix
の使い方
ステップ5.2/15
他のプログラムから保存される値を呼び出す方法
`r()'の使い方がややこしい。

[0132 ごろ 2011/02/23(水) 13:44:37.37]

Stataの条件設定って結構融通が効かない
replace var1=0 if no==1 | no==8
replace var1=0 if no== (1 | 8)
上の式はOKだけど、下の式では8は無視される

replace var=1 in 1/3 /* change from var1[1] to var1[3].
in を使うとこういう表現ができる。ifでは使えない。使い分けが必要。

[0133 ごろ 2011/03/02(水) 15:06:55.87]

今週いっぱいでネットコースが終わる。
一応全てのレクチャーにひと通り目を通した。

まだ改善点はあると思うが、全４回のレクチャーを終了すると、
adoファイルに書かれている内容が　"ちょっとだけ"　読めるようになる。
値段の割にはお得だった気がする。

ステップ5.3/15
return()とr()の違い return()は自分で定義した変数の値を返す。r()はStata関数の実行した値を取り出す。

ステップ6/15
ブートストラップ法
bootstrap r(mean), reps(1000) nodots: summarize mpg
ここでreps()は繰り返す回数
nodotsは＃を表示しないオプション
プログラムを利用したブートストラップの方法

ステップ7/15
trace on/offでプログラムのエラーを見つける方法

ステップ8/15
ブートストラップの前に欠損データを削除する方法
keep if 変数名 < .　この表現重要！
Stataのhelpファイルより
Numeric missing values are represented by large positive values. The ordering is
all nonmissing numbers < . < .a < .b < ... < .z
ステップ9/15
ブートストラップでスピードをUpする方法
ブートストラップを走らせる前に不必要な変数を削除する

[0134 ごろ 2011/03/02(水) 15:08:08.80]

ステップ10-13/15
プログラムを使ってモンテカルロ法をする方法
モンテカルロ法は、計算された95%CIが本当に正しいかどうかを確かめる時などにも使える手法

ちょっとややこしいが、モンテカルロ法を行う時などに使う新規データベースの作成方法
program プログラム名
postfile 新規データベースを作るための名前 新変数名１　新変数名２　新変数名３　 using 作成するデータベース名,replace
forvalues i = 1/作成するデータ数 {
drop _all
　　/* ここに計算式 */
post mysim (新変数名１に代入する値)　(新変数名２に代入する値)　(新変数名３に代入する値)
}
postclose 新規データベースを作るための名前
end
quietlyコマンドの使い方
quietlyを適応する部分を別のプログラムにして
プログラムの中から
quietly 別プログラム名
で呼び出す方法

14/15
モンテカルロ法をスピードアップする方法
which コマンド名
これでbuilt-in command:と表示されるコマンドは実行速度が速いが、
/Applications/o_Stata/ado/base/c/ci.ado
*! version 3.3.15 17mar2005
このように表示される.adoコマンドは実行速度が遅いのでモンテカルロ法などで使うのを避けるべき

[0135 ごろ 2011/03/03(木) 05:09:56.70]

Lecture3-15/15
練習問題
date関数は以外と融通が利かない。以下例
display %td date("Feb02","MD") /* OK, but Year is 1960. */
display %td date("Feb02","MD1993") /* This doesn't work. */
display %td date("Feb02"+"1993","MDY") /* OK */
display %td date("Feb02"+ string(1993),"MDY") /* OK */
display %td date("Feb02","MD",1959) /* Strange, Year is 1960. */
display %td date("2001/02/03","YMD") /* OK: 03feb2001 */
display %tc date("2001/02/03","YMD") /* NO: 01jan1960 00:00:15 */

display string(1993) + "Feb02" /* OK: 1993Feb02 */
display "Feb02" + string(1993) /* Error: Feb02unknown function +string() */

ブートストラップの結果表示のための関数
. matrix define 変数名 = もともとのデータから計算　/*計算された値はMatrixに直す必要がある。*/
. bstat, stat(変数名) n(対象数)

[0136 ごろ 2011/03/03(木) 05:11:59.82]

Lectrure4
1/24 pwd /* 現在の作業ディレクトリを表示 */
2-6/24 Adoファイルを用いて、capture log closeなどのコマンドをサポートする。
discard
adoファイルはメモリに読み込まれるため、元ファイルを変更してもすぐには反映されない。
このdiscardコマンドで、adoファイルを再読込する。
capture noisily ・・・
このコマンドで、途中でエラーが起きてもlogファイルが作成されるようになる。
7/24 sysdirコマンドでStataが使う７つのディレクトリが表示される。
通常は、自作のadoファイルはPERSONAL:で表示される以下に格納する。Macなら例えば
PERSONAL: ~/Library/Application Support/Stata/ado/personal/
となる。
8/24 Stata helpファイルの作成方法
9/24 Doファイル、プログラム、adoファイルの使い分け

[0137 ごろ 2011/03/03(木) 05:13:20.91]

10/24 adoファイルで有効な一時的な変数を作成するtempvarの使い方
11/24 preserveコマンドの使い方。
プログラムが終わったり、エラーが出るとこのpreserveコマンドを実行した時点のデータセットに戻る。
12/24 一時的にファイルを作成するtempfileコマンドの説明
13/24 adoファイルでもassertコマンドを使いなさいという教え
14/24 変数名をマクロで読み込むadoファイルの紹介
15-16/24 summarize,detailコマンドで得られる結果のうち、interquartile rangeのみを表示させる新しいプログラムを
作成しながら、どのように上手く結果を表示させるかなどの説明
使われている関数は、
_col()関数、Stataは初期設定で%9.0gで結果を返すので、９文字以上の数列にはならないことを利用。
17-18/24
syntax [varlist]
foreach i of local varlist {
この組み合わせで、Stataの通常のコマンドのような操作が出来るようになる。
19/24 syntax [varlist]の詳しい説明 help syntaxで表示される内容とほぼ同じ。
20/24 自作プログラムで、if,in, =exp, usingf, w aw pw iwコマンドが使えるようにする

[0138 ごろ 2011/03/03(木) 05:50:41.84]

21/24
syntax [varlist]
syntax varlist の違い。[ ]が無いときには、変数名が指定されていないとエラーになる。
特定の処理をさせる場合にはエラーを避けるためになるべくifオプションなどは使わない。
使う場合には次のコマンドでデータを守る。
if "`if'" != "" {
preserve
quietly keep `if'
}
22/24
" ",`', `" "'の使い方。ここは非常に参考になった。
マクロのに""を含む文字列を代入すると、エラーが起きることがある。
それを避けるために`" "'を使う。
ここを学んで、何故Stataが`'などというおかしな引用符を使うのかが少し分かった気がする。
23/24
version 10
というようなバージョンコントロールの説明
*! （adoファイル内）で始まるコメントは、which コマンド名　を実行したときにコメントが表示される。

24/24
練習問題
error エラー番号　エラーの説明が表示される
abbrev(変数名,表示させる長さ)　変数名が長いとき、省略してくれる。
unab マクロ名 : _all マクロ名に、データに含まれる省略されていない変数名が全て格納される

これで一応全てのレクチャー内容の報告を終わります。
参考にして頂ければ幸いです。

[0139 ごろ 2011/03/04(金) 03:00:03.77]

Stataのバグ？
. display string(1) + "A"
1A
. display "A" + string(1)
Aunknown function +string()
r(133);
. display = "A" + string(1)
A1

string()と文字列を足し算するときに順番を入れ替えるとエラーになる
意味の分からない=だが、これを付けるとOK

[0140 ごろ 2011/03/04(金) 03:03:37.69]

今、気がついた
A1だと変数名と勘違いされてしまうからエラーになるのかな？

[0141 ごろ 2011/03/05(土) 08:33:06.57]

Advancedも受けようかなぁ

NetCourse® 152: Advanced Stata Programming
Content:
This course teaches you how to create and debug new commands that
are indistinguishable from the commands in Stata.
It is assumed that you know why and when to program and to some extent
how.
You will learn how to parse both standard and nonstandard Stata syntax
using the intuitive syntax command, how to manage and process saved
results, how to process by-groups, and more.
Course leaders:
Kerry Kammire, technical services analyst at StataCorp
Jennifer Rolfes, technical services representative at StataCorp
Course length:
7 weeks (5 lectures)
Dates:
October 7–November 25, 2011
Prerequisites:
* Stata 11, installed and working
* Course content of NetCourse 151 or equivalent knowledge
* Internet web browser, installed and working
(course is platform independent)
Price:
$150.00

[0142 １３２人目の素数さん 2011/03/25(金) 20:00:53.28]

＠保守

ごろさん大丈夫かぁー？？

[0143 ごろ 2011/03/27(日) 03:08:58.46]

>>142
ありがとうございます。私と私の周りの人は皆無事です。
しかし、あまりの事態に、のんきにStataしている気になれませんでした。
地震・津波そのものの被害も大変なものでしたが、原発事故の影響が恐ろしいです。
犠牲になられた方のご冥福をお祈り申し上げます。

[0144 １３２人目の素数さん 2011/04/10(日) 23:28:41.47]

ちなみに、stataのオンライン・レクチャーはDVDでみれたり、あるいは、過去のものをみれたりするのでしょうか。社会人の私としては、決まった時期にみるというより、休みの時に一気に見たいのですが。

[0145 ごろ 2011/04/12(火) 05:52:11.98]

>>144
オンラインレクチャーというよりは、オンラインで置いてあるテキストを読みながら
自分で勉強するという感じですね。
そういう意味では本を購入するのと変わらないので１万以上するのは高い気もしますが、
質問もできるというインタラクティブ性が良いのでしょうね。

テキストファイルなので、時間があるときに、自分のペースで勉強出来ます。
質問も、24時間OKです。返信は２４時間いつでも対応ではありませんが、数時間以内には
全て答えが返ってきていましたし、とても参考になる答えもありました。

講義や、質疑応答の資料は後日、ダウンロードできるようになるので、復習も出来ます。
本と比較すると高いですが、巷の講演会を受講するよりは、質問も何度でも好きなだけできますし、
ある程度期間があるので、自分を強制的に勉強させることが出来たし、結果的には良かったのかと
思っています。


Advanceコースは10月なのでまだ先ですが、一応今のところは受けるつもりです。
他にも興味深いコースはありますが、２万円以上するので、さすがにちょっと躊躇しています。

[0146 １３２人目の素数さん 2011/04/13(水) 00:04:53.57]

>>145

せっかくオンラインなのだから、DVDなどで販売してくれたほうがうれしいですよね。
何度も見返すことができるし、第一、好きなときに勉強できます。
次回も、上級レベルと通常レベルがかぶっていて、両方受講するにはちょっとタイミングが悪いです。
といういうニーズの対応はTACやLECは普通やっていることですので、臨機応変にしてほしいところです。

ちなみに、どれくらいの期間視聴が可能なのでしょうか。

[0147 １３２人目の素数さん 2011/04/13(水) 00:32:24.18]

>>146

ちなみに、MATAについてはオンライン・セミナーで議論されるでしょうか。ネット上でかかれた.mlibファイルが開くことができず、stata9で開けず苦労しています。このあたりの勉強を講義ベースで聞けたらと思うのですが。

[0148 ごろ 2011/04/13(水) 02:09:41.64]

>>147
視聴期間は基本的に講義開催期間のみですが、先にも書いたとおり、最終的には
サイト自体がダウンロードできるようになるので、講義終了後に繰り返し復習することは可能です。

MATAについてはよく分かりませんが、少なくとも私が受講したコースには有りませんでした。
しかし、講義から離れた内容の質問でも答えてもらえる可能性は高いと思います。

[0149 ごろ 2011/04/24(日) 23:29:58.92]

ノートパソコンのOSをクリーンインストールし直した。
前は、Stataの10と11を両方使っていたけど、今回は11のみでいくことにした。

[0150 ごろ 2011/05/04(水) 18:13:09.80]

今byableなプログラムを書いているけど結構難しいなぁ

capture program drop myprog
program myprog , byable(recall)
syntax varlist
marksample touse
foreach i of local varlist {
quietly:sum `i' if `touse'
di abbrev("`i'",9) _col(10) "N ="r(N) _col(25) "Mean +/- SD="%4.1f r(mean) " +/- " %4.1f r(sd) _col(55) "Min =" %4.1f r(min) _col(70) "Max =" %4.1f r(max)
}
end

こんな感じになるんだけど、myprogに変数が沢山つくと挙動がおかしくなる。
なんでだろー

[0151 １３２人目の素数さん 2011/05/05(木) 22:21:25.50]

MAC版で日本語を表示させると、詰まって表示されるのですが、
何か対策はありますか。

教科書体を利用すると、多少ましになりますが、他に手段は
ないのでしょうか。

[0152 ごろ 2011/05/06(金) 02:46:33.01]

>>151
私もMac版 v11.2を使っています。

Stataで日本語表示をさせようと考えたことが無かったので
気が付きませんでしたが、同じ状況は再現できました。

試したコマンド
di "日本語を表示させる"

OsakaとかもPreferenceから選べますが、ダメでしたね。
やっぱりDF教科書体が一番まともでした。

コピーしようと選択すると文字化けみたいになりますね。
ペーストはちゃんとできますが、日本語表示のために何か複雑な仕組みでも
あるのでしょうかね。

詰まってしまうのは、等幅フォント設定の為ではないか、と思いますが、
それを解除する初期設定項目もないので難しそうですね。

明確な解決案は思いつきませんが、
例えば、EmacsからStataを操作することが出来れば
うまくいきそうな気がします。

以前試そうと頑張ってみましたが途中で諦めました。

あと、Font関連では、ギリシャ文字や±などの特殊文字もStataで表示させるのは
難しいようですね。
既存のフォントを入れ替えて表示させる方法がどこかにありましたが、
また見つけたらリンクを貼ります。

[0153 １３２人目の素数さん 2011/05/06(金) 22:39:02.59]

>>152
早速のお返事ありがとうございます。
うまくいかないものなのですね。

どうしても、日本語を使いたいときは、
教科書体でがんばってみたいと思います。

[0154 １３２人目の素数さん 2011/05/11(水) 00:21:35.86]

stataのマニュアルは購入しなければならないのでしょうか。
time serisとpanelのマニュアルが取得したいのですが。。

[0155 ごろ 2011/05/12(木) 11:41:33.04]

>>154
はっきり覚えていませんが、紙のマニュアルは無かったと思います。
しかし、オンラインマニュアルや、ソフト付属のマニュアルも
英語ですが、なかなか詳しいです。
プリントしてあるマニュアルは買わないといけないと思います。

[0156 ごろ 2011/05/12(木) 11:43:03.43]

>>155
あ、私も紙のマニュアルがあるときっと便利なんだろうなぁと思いながら、
ついつい、Googleで済ませてしまってここまで来ています。
あと、本当に困ったら、Stataの英語メーリングリストがあるので、そこで聞いていますね。

[0157 ごろ 2011/05/24(火) 09:54:32.84]

最近、Sleepsortという面白い記事を見つけたので、Stataでも出来ないかと
試し中。
しかし、バックグラウンドで関数を動かす方法が分からないのでつまづいています。
http://d.hatena.ne.jp/shinsuku/20110520/p1

[0158 あんでぃはアホ ◆AdkZFxa49I 2011/05/29(日) 23:38:17.38]

あんでぃ

[0159 ごろ 2011/05/30(月) 00:21:16.77]

どうやら、Stataには二つの関数を同時にうごかす方法はないようです。
残念！

[0160 ごろ 2011/06/29(水) 16:49:52.89]

pwcorrコマンドの結果をreturn listで見るとmatrixがあるんですが、
このmatrixの扱いが良くわからない。
matrixの要素だけを取り出してdisplayする方法はあるのでしょうか？

[0161 １３２人目の素数さん 2011/07/03(日) 14:22:31.00]

ヴァージョン１２が25日にでますね。
共分散構造分析ができるようになりますね。
アップデートしたいが、先立つものがない・・・

http://www.stata.com/stata12/

[0162 １３２人目の素数さん 2011/07/08(金) 10:53:39.06]

このソフトいれてみたけどなんでまともにデータファイルが読めないんだか
ふつーにエクセルぐらい読めよバカ

[0163 １３２人目の素数さん 2011/07/10(日) 08:46:36.61]

>>162の率直な意見を取り入れてstata12ではExcelファイルを読めるようになりますた

[0164 １３２人目の素数さん 2011/07/11(月) 11:23:03.75]

まったくとろい会社だな。

[0165 １３２人目の素数さん 2011/07/15(金) 11:24:37.13]

STATAって貧乏人のSASですか?

[0166 ごろ 2011/07/17(日) 09:34:11.52]

>>162
まぁ、Excel側でCSVファイルで出力すれば言い訳だから、そんなに困んないけどね。

[0167 １３２人目の素数さん 2011/07/19(火) 15:04:48.01]

急いでる時にわざわざExcelでセーブしなおすのは結構面倒だぜ。なにせ
俺様売れっ子だから

[0168 １３２人目の素数さん 2011/07/19(火) 23:13:24.38]

>>165
プログラムの書けない奴のRだと思う

[0169 １３２人目の素数さん 2011/07/19(火) 23:49:35.33]

>>168
Rって書ける、書けない、言うほどのものか？

[0170 １３２人目の素数さん 2011/07/31(日) 08:11:01.67]

>>169は優秀なスレストッパー(ﾆﾔﾘ

[0171 １３２人目の素数さん 2011/08/02(火) 12:16:30.64]

>>161
ようやくSEMができるようになったか。
Stata 11を入れたばっかりだけど、買おうかなー。

[0172 ◆YeNleXj6MwwY 2011/08/15(月) 01:19:56.07]

　|>､＿_ ／|ﾌ
　 〉｡,　o　く　/＼
　 b 〜　O ヽ〉 /
　 （>　 ＜ノ ∨
　　（､＿（、 ノ
　　￣ 　 ￣

[0173 ごろ 2011/08/16(火) 23:39:50.80]

おー。昨日アクセスできなくて落ちてしまったかと思った。

[0174 ごろ 2011/09/12(月) 00:37:57.76]

保守

[0175 ごろ 2011/09/13(火) 02:31:02.76]

そろそろまたやりますよ。

[0176 １３２人目の素数さん 2011/09/14(水) 19:53:32.28]

　|>､＿_ ／|ﾌ
　 〉｡,　o　く　/＼
　 b 〜　O ヽ〉 /
　 （>　 ＜ノ ∨
　　（､＿（、 ノ
　　￣ 　 ￣

[0177 １３２人目の素数さん 2011/09/25(日) 09:19:35.91]

試用はできないんですか？

[0178 ごろ 2011/09/25(日) 14:16:20.84]

>>177
試用できないのかな？オンラインではないし、聞いたことないなぁ。大学とか知り合いで使っている人がいれば一番いいんですけどね。

[0179 １３２人目の素数さん 2011/09/25(日) 22:22:05.78]

日本の代理店も本家も探してみましたが、ありませんでした。
SEMはどんな感じかなと思いまして。
AMOSはお絵描きでできると言えば簡単そうなんですが、誤差変数を忘れずに入れろとか、
制約を適切に入れないといかん、とか実は結構面倒なので。
毎年のお布施も結構かかりますし。
Stataはぱっと見た感じでは、核となる部分さえ入力すればいいように見えたものですから。

[0180 猫は熾猫 ◆s9Thv0XHxQ 2011/10/01(土) 11:42:22.56]

猫

[0181 １３２人目の素数さん 2011/10/01(土) 18:21:49.01]

>>179
12買ったけど、その機能はまだ使ってないっす。
いじくる暇を見つけたらレポートします。

[0182 ごろ 2011/10/04(火) 23:06:51.89]

>>181
おお、１２買ったんですね。いいなぁ。
レポートよろしくお願いいたします。

[0183 ごろ 2011/10/07(金) 12:53:48.80]

NetCourse 152始まりました。
また日記形式で綴っていきたいと思います。

[0184 177=179 2011/10/10(月) 20:48:18.79]

Stata/ICという一番下のクラスのを買いました。
とりあえずauto.dtaというファイルで、mpgとweightとforeignで回帰分析してみました。
次に、SEMのお絵描きで、mpg<-weightとmpt<-foreignの2本の矢印のみを描いてEstimateしてみると、
同じ結果が得られました。誤差変数の○を書き足したり、制約を入れる必要もない、という手軽さ。

どの程度使えるものか、少し試してみます。

[0185 １３２人目の素数さん 2011/10/11(火) 00:54:36.90]

新参初心者ですが、STATAっていくらから購入可能なんですか？

[0186 １３２人目の素数さん 2011/10/11(火) 23:41:01.11]

>>185

ステータス（学生なのか、一般なのか）、住んでいる国、
パッケージプランなどによって大きく違いますよ。
日本であれば、ライトストーンという会社が正規ディーラー（？）として
ライセンス委託され、Stataの代わりに販売をしているようなので、価格表は以下を参照してください。
（私は日本に住んでいないため、この辺の事情は詳しくありません。）

一般：ttp://www.lightstone.co.jp/stata/price.htm
学生：ttp://www.lightstone.co.jp/stata/gradplan.htm

一般の方であれば、最新のStata 12 IC（一番下のクラス）が52,500円、
学生の方であれば、同製品が21,000円です。

一番安いのは、おそらく学生用のStata 12 smallで6,825円ですが、
1年間しか使えませんし、扱えるサンプルの数にも大きな制約があるようです。

「お金が無くてそんなに出せない！」って言う人は頑張ってRを使いましょう。
最近は入門書も充実してるみたいですしね。

[0187 １３２人目の素数さん 2011/10/12(水) 00:51:08.07]

>>186
ありがとうございます。

[0188 １３２人目の素数さん 2011/10/13(木) 12:02:44.59]

>>184
どーも、なんだか、先を越されてしまいましたね。
買ったはいいが、まだまだまともに使えてないもんで。。。

Stata 12では、marginsplotという限界効用のグラフを書くのに便利なコマンドが搭載されたようです。
11までは、marginsで算出した結果を一度マトリクスに書きだして、
その情報の一部を取り出したり移し替えたりしてグラフを作っていたので、
今度は大分便利になってるはずです。

まだ使ってないのですが、何か分かったらここに書きます。

[0189 184 2011/10/13(木) 19:02:19.15]

>>188
いやいや、色々教えてください。Stataは初めてです。

で、次に、mpg<-weightとmpt<-foreignと直接矢線を入れず、潜在変数を作って、
その潜在変数に対してweightとforeignから矢線を引いたところ、これも攪乱項を自動的につけてくれました。
で、潜在変数からmpgに矢線を引いてお絵描きは完了。

これで実行すると、名前のない変数があるよ、との指摘。
さすがにそこまで手抜きは許してくれないようです。
潜在変数をf1という名前にして、もう一度。今度はIterationが止まりません。
f1についた攪乱項のVarianceを1にとめて実行すると、ようやく推定完了です。

ただ、探し方が悪いのか、適合度の指標について、GFIやAGFIという表記がありません。
マニュアルを読まないといかんようです。

[0190 184 2011/10/29(土) 15:19:56.64]

その後、いろいろと探索中。
適合度の指標としては、RMSEAとCFIが算出されます。AIC、BICも出ますね。
RMSEAについては、Probability RMSEA <= 0.05が出てます。これはいいかも。
TLIというのはあまり聞かない。
SRMS(Standardized root mean squared residual)というのと、CD(Coefficient of determination)というのも出てます。
return listやereturn listで返されるものからGFIを計算する方法がないか、調べ中です。
finditしてもちょっとわかりませんね。
あと、コマンドで指定する方法もわかりましたが、これをSEMビルダーというお絵描きの方に反映できればいいのですが。
ちょっと難しいか。

[0191 184 2011/11/05(土) 11:06:50.82]

遅々として進みませんが、気付いたことを。
regとかlogitとかの基本的な分析コマンドも
adoファイルとして提供されてるんですね。
読んでもちっともわかりませんが。

[0192 １３２人目の素数さん 2011/11/08(火) 11:03:18.08]

>>191

私もStataプログラムの構造的な話はよく分からないのですが、
確かに、分析コマンドは全部adoファイル経由だったと思います。
ただ、後から追加されたものと、基本的な分析用のファイルとでは格納場所が違うようです。

まぁ、分析する上ではあまり気にされなくてもよいと思います。

[0193 ごろ 2011/11/10(木) 13:12:01.34]

NetCourse 152記録
だいぶ遅れてしまいましたが、ぼちぼちやります。
Lecture１　1/6 前回レクチャー(NC-151)のレビュー
基本的な内容
ポイントは、例えばmytestというプログラムに
syntax [varlist] [if]
と書いてある時
mytest 変数１ 変数２ if 変数２==0
と入力すると
このコマンドを実行した瞬間に
マクロ`0'には”変数１ 変数２ if 変数２==0”がスペースの数も含めて入力されたまま
マクロ`1'には"変数１"
マクロ`2'には"変数２"
マクロ`3'には"if"
マクロ`3'には"変数２==0"
が代入されて、
プログラムの中の
syntax [varlist] [if]
が実行された後初めて
マクロ`varlist'には"変数１ 変数２"
マクロ`if'には"if 変数２==0"
が代入されるということ

[0194 ごろ 2011/11/10(木) 13:12:16.49]

2/6 Parsing(構文解析）
結局ここから延々とSyntaxコマンドの解説が始まる。
マクロ`0'はスペースを含めて、入力したまま
マクロ`*'はスペースは一つに調整される。
Sample code
program sample0
syntax [varlist] [if]
di "line1 " "`varlist'"
di "line2 " "`if'"
di "line3 " "`0'"
di "line4 " "`*'"
end

実行結果(bbb,ccc,dddの変数を定義しておく必要あり）
sample0 bbb ccc ddd if ddd==5
line1 bbb ccc ddd
line2 if ddd==5
line3 bbb ccc ddd if ddd==5
line4 bbb ccc ddd if ddd==5

[0195 ごろ 2011/11/10(木) 16:15:23.29]

2/6続き
syntax コマンドの後、要素を書かないと使えない。[　]に入れるとどちらでも良い、直接書くと、必須になる
ex. syntax varlist [if]
と書くとvarlistは必須、ifはどちらでもOK、その他の例えばbyとかは使えない
syntax [varlist]でプログラムコマンドの後変数名を省略すると
Stataはデータセットの全ての変数を入力したものとして取り扱う。
これは、syntax[varlist(default=none)]で変更出来る。
細かいことは全てhelp syntaxに書いてある！！

3/6 Parsing options
syntaxコマンドの,（カンマ）以降はオプションになる。
オプションにはどんな名前も指定できて例えば
syntax ・・・,[nandemo]
とすると、
プログラム名　変数名, nandemo
と入力すると
マクロ`nandemo'にnandemoが入力される
nandemo以外の文字を入れるとエラーになる
オプションを指定しないと
マクロ`nandemo'は""が入力される
また、オプションのマクロ名はNANdemoなどとsyntaxで指定することにより
プログラム実行時にnandemoと全て入力しなくてもNANで入力OKになる。
また、noからはじまる(大文字NOならOK)オプション名は使えない
プログラム実行時にオプション名を大文字にするとエラーになる。
nonandemoとオプションを指定すると
マクロは`nandemo'が作られて、nonandemoと入力された時のみ
マクロ`nandemo'にnonandemoが代入される。それ以外は""が代入される。
ここはややこしいところ。

[0196 ごろ 2011/11/13(日) 00:26:51.92]

3/6続き
オプションは更にオプション名( )でオプションの性質を指定できる。
help syntaxの中に詳細な説明があるが、レクチャーでは次の８つの解説
integer 整数を指定　初期値設定可
real 実数を指定　初期値設定可
string 文字列を指定　
varlist １から複数の変数名を指定。ない変数名を入力するとエラー
numlist 数列を指定
passthru
namelist を指定
cilevel を指定

[0197 ごろ 2011/11/13(日) 00:27:31.09]

>>196は書きかけだった

3/6続き
オプションは更にオプション名( )でオプションの性質を指定できる。
help syntaxの中に詳細な説明があるが、レクチャーでは次の８つの解説
integer 整数を指定 初期値設定可
real 実数を指定 初期値設定可
string 文字列を指定 初期値設定不可
varlist １から複数の変数名を指定。無い変数名を入力するとエラー
numlist 数列を指定
passthru stringに似ているがオプション名と(　)も含めてユーザーが入力したものを返す。
namelist １から複数の変数名を指定。新しい変数名や未定義の変数名、マトリックスも扱える
cilevel 信頼区間の範囲を指定　初期値は95となっている　指定できる値は10-99.99

オプションに*をつけると、optionsというマクロに使われなかった残りの文字列が代入される。

4.複雑なsyntax構文解析
ttestの３つのパターンttest varname==#(数字), ttest varname1==varnmae2, ttest varname, by(groupvar)の
３種類を例にして解説

今は不必要になったtokenizeコマンドの解説　varlistを`1',`2'・・・というマクロに代入し直す。
syntax varname [=/exp]
/をつけると=なしで`exp'に論理式が代入される。/がないと先頭に=とスペースが入る

ttestの３つのパターン分けにifと次の構文を使う
syntax rvarnmae [=/exp] ・・・
capture confirm number `exp'
`exp'が数字のみならcaptureで得られたエラーメッセージがないため_rc==0となる
エラーが出るとエラー番号の７が代入されて_rc==7となる
(confirm numberコマンドは、次に来る変数が数字かどうかを判定するコマンド)

[0198 ごろ 2011/11/13(日) 03:32:43.61]

5.subprogramの話
.adoファイルの中でいくつもプログラムを定義することが出来るが、その中で定義されたサブプログラムは、
.adoファイルの中でしか使われない。
調べるためにはprogram dirを実行する。
例えば、mytest.adoファイルの中にsubprogramなどがあると、
mytest.suprogram
mytest
と表示される。

6.Exercises
忘れていたが、Stataには予め使えるデータセットがある。help datasetなどでチェック可。
使う場合は、例えばauto.dtaデータセットならば
sysuse auto.dta
で使える。
Exerciseは４問あり、最初の２問は簡単。
３問目は、95%CIを自分で計算するもの。
４問目は、stataのif文は融通が効かないのでしっかり書かないとエラーを起こす。
ex.
if `var'=="" { di "nothing" } else { di `var'}
ではダメで
if `var'=="" {
di "nothing"
}
else {
di `var'
}
と書かないとダメ。
やっとLecture1終わった。

[0199 ごろ 2011/11/13(日) 11:17:27.43]

Lecture2 1/8
1. Quotes
set trace on でエラーの箇所を調べる
netCourse 151にもあったと思うが、
代入される文字列に" "が入ると混乱するため`" "'を使うという内容

2. 続き
Lecture1からの続きで３つのタイプのttestを行うコマンドを自分で作成する。

3. Temporary variables
プログラム内で計算をするために新しい変数を作成する場合は、
gen double 新しい変数名
を使う。これだと最も高精度で計算をする。

プログラム内のみで使う変数には
tempvar 変数名
を使う。
変数名には__000001のような名前が付けられるが、
プログラム内ではユーザーは上記で指定した変数名が使える。

4.続き
エラーの出し方
if "`by'" == {
display as error "by() required"
exit 198
}
error 198はsyntax error を意味する。しかし0以外ならOK

assert 論理式
論理式が真なら何も返さず、偽ならエラーコードを返す。

[0200 ごろ 2011/11/13(日) 11:30:44.15]

5. r()について
プログラム/関数を実行したときに返す値はr(N)などに格納される。
このリストを調べるにはreturn listを使う

[0201 ごろ 2011/11/15(火) 16:39:40.53]

Lecture2-6　式のプログラミング
[if/]でif がマクロ`if'に入力されなくなる。
markとmarksample
mark 変数名
を実行すると全て１の変数が作成される。
mark 変数名 if 論理式
論理式に合うものは１，合わないものは０にコードされた変数名が出来る。
marksampleは新しく変数を作成しない方法。
markで変数名を作成したくなければtempvarを前もって使う。
marksampleはifが使えないので使いたい時には次の行でreplaceコマンドを使う。
program sampleone
syntax [if][in]
marksample tempvarname
di `tempvarname' in 1/50
end
と作成し
sysuse auto.dta
sampleone if mpg <20 in 10/20
とすると、10-20の中でmpgが<20以下のものはtempvarnameに１、そうでないものは0がコードされる。

markはsyntaxコマンドと一緒に使っても上記のようにならない。
ちょっと違いについて混乱中
7はコードのまとめ
8は練習問題

[0202 ごろ 2011/11/18(金) 12:29:38.95]

Lecture2-8練習問題
1. マクロ`by'でなく`"`by'"'を使う理由
2. program define プログラム名
通常はdefineは省略できるが、例えばdefineという名前のプログラムを作る場合などには必要。
varlist(min=2 max=2)とすることでvarlistの数を２つに特定できる。
predictコマンド
直前にregressなどを実行すると、その結果に基づいて新しい変数名に予測値が保存される。
tokenizeの有用性
varlistを使う時、tokenizeを使わないと何故かvarlist名の後ろに,(カンマ)も入る。
これはプログラムにエラーを起こす原因になりうる。
3. オプションsaving(ファイル名) or saving("directory/フアイル名",replace)で
ファイルが保存される。
4. linear regressionではなく、二次式で回帰分析をやる手法の解説

[0203 ごろ 2011/11/18(金) 14:25:35.30]

5. は長いが、結局二次式による回帰分析の結果を示す方法について書かれている
新しい話は無し。
6. エラー表示のさせ方
エラーコード301はlast estimates not foundという意味らしい
[P]errorに解説がるらしい。
7.if文で条件分けして結果を表示させる問題
8. optionにnoConsを使って結果の一部を表示させないオプションを作る問題。
noConsによってconsマクロに入るものが何かを考えること以外は平易。
Lecture2終了

[0204 ごろ 2011/11/19(土) 00:14:17.52]

Lecture3　1/12
1. 何を学ぶべきか
syntaxコマンド, localマクロ(`',"`'"の使い分け),
サブルーチンをマスターする。
マクロに分かりやすい名前を付ける。
単純なプログラムを長く重ねることでちゃんと動くプログラムを作る。

2. Scalars
Scalarsとマクロの違い
scalar x=2　　 　local x=2
di x di `x'
左と右は同じ結果が出る
マクロはどこでも使えるのでtt_`x'というようにも使える。

注意！ Stataでは同じ文字列に３種類の可能性があり、それぞれ別のものを代入出来る。
1. 変数名 2. Scalar 3. マクロ
マクロは`'で区別がつくが、1と2は区別がつきにくい。
例
scalar x=1
di x /* 1 と表示される */
input x /*変数xの定義開始*/
2
3
end /*変数xに一つ目に2,２つ目に3を代入*/
di x /*変数xがある場合には変数xの最初の値2を返す*/
di scalar(x) /*しかしscalarのxも残っていて、これで呼び出せし1と表示*/

[0205 １３２人目の素数さん 2011/11/19(土) 13:19:11.95]

魂は幾何学です
誰か(アメリカ)気づいた
ソウルコピー機器
無差別で猥褻、日本は危険知ったかブッタの日本人が失敗作

↓テロ資料です。

[0206 １３２人目の素数さん 2011/11/19(土) 13:19:49.31]

電波テロ装置の戦争(始)
エンジニアと参加願います公安はサリンオウム信者の子供を40歳まで社会から隔離している
オウム信者が地方で現在も潜伏している
それは新興宗教を配下としている公安の仕事だ
発案で盗聴器を開発したら霊魂が寄って呼ぶ来た
<電波憑依>
スピリチャル全否定なら江原三輪氏、高橋佳子大川隆法氏は、幻聴で強制入院矛盾する日本宗教と精神科
<コードレス盗聴>
2004既に国民20%被害250〜700台数中国工作員3〜7000万円2005ソウルコピー2010ソウルイン医者アカギ絡む<盗聴証拠>
今年５月に日本の警視庁防課は被害者SDカード15分を保持した有る国民に出せ!!<創価幹部>
キタオカ1962年東北生は二十代で2人の女性をレイプ殺害して入信した創価本尊はこれだけで潰せる<<<韓国工作員鸛<<<創価公明党 <テロ装置>>東芝部品)>>ヤクザ<宗教<同和<<公安<<魂複<<官憲>日本終Googl検索

[0207 １３２人目の素数さん 2011/11/22(火) 00:52:44.75]

>>205-206
完全にスレ違いだろ。。。
ガイドラインは無いけど、ちゃんと空気を読んでから投稿してね☆

[0208 １３２人目の素数さん 2011/11/22(火) 22:32:56.09]

>>202の4.と>>203の5.はどんな話？

[0209 ごろ 2011/11/23(水) 03:11:07.57]

Lecture3　1/12
1. 何を学ぶべきか
syntaxコマンド, localマクロ(`',"`'"の使い分け),
サブルーチンをマスターする。
マクロに分かりやすい名前を付ける。
単純なプログラムを長く重ねることでちゃんと動くプログラムを作る。

2. Scalars
Scalarsとマクロの違い
scalar x=2　　 　local x=2
di x di `x'
左と右は同じ結果が出る
マクロはどこでも使えるのでtt_`x'というようにも使える。

注意！ Stataでは同じ文字列に３種類の可能性があり、それぞれ別のものを代入出来る。
1. 変数名 2. Scalar 3. マクロ
マクロは`'で区別がつくが、1と2は区別がつきにくい。
例
scalar x=1
di x /* 1 と表示される */
input x /*変数xの定義開始*/
2
3
end /*変数xに一つ目に2,２つ目に3を代入*/
di x /*変数xがある場合には変数xの最初の値2を返す*/
di scalar(x) /*しかしscalarのxも残っていて、これで呼び出せて1と表示*/

[0210 ごろ 2011/11/23(水) 03:24:28.13]

Lecture3
2. Scalars続き
programの中で
args マクロ名1 マクロ名2 ・・・と書いて
program名 インプット１　インプット２　・・・と実行すると
それぞれローカルなマクロ名１にインプット１が入力される。

次は、backslashの話
日本語環境のStataではこの辺りはややこしい話題
backslashとYen記号問題も絡む。
実際backslashでDoファイルに入力して保存すると、
次開くとYen記号になったり、？マークになったり、
なぜかIに置き換わったりと挙動不審。

ここの解説ではエスケープ記号としてbackslashを使う場合、
マクロにちゃんと入力するにはmacval()コマンドを使う。
そうでなければローカルScalarとして入力することが解説されている。

[0211 ごろ 2011/11/24(木) 06:21:57.28]

Lecture2-3続き
macval()コマンドの使い方
local mymacro=`"Nakamiha Nandemo OK.For example \$ sign"'
di "`macval(mymacro)'"
でOK

2-3 Binary accuracy
コンピュータの2進法で計算するために起きてくる問題の話。
具体例が挙げられている。
それを避ける方法としてscalrを使う方法の解説

[0212 ごろ 2011/11/25(金) 15:00:09.14]

>>211のLecture2-3続きはLecture3-2続きの間違い
Lecture3-4
マクロとscalarとどちらが正しいか？
マクロは１２桁まで、Scalarは１６桁までなのでScalarがより正確
迷うならScalarを使っておけということ！？

[0213 ごろ 2011/11/25(金) 20:15:31.35]

Lecture3-5は、マクロで書いていたプログラムをscalarで書きなおす説明
例
quietly sum `varlist'
local xbar=r(mean)
なら
tempname xbar
quietly sum `varlist'
scalar `xbar`=r(mean)
となる。

3-6. byオプションの使い方
変数xの数を数えるテクニック

tempvar n_variable
sort x
quietly by x: gen `n_variable' = 1 if _n==1
quietly replace `n_variable' = sum(`n_variable')
di `n_variable'[_N]
欠損値もひとつの変数と数えるので注意

Lecture3-7
program プログラム名, sortpreserve
プログラムの記載
end
これでこのプログラム内でsortを使ってもプログラム終了後は順番が元に戻る。
しかし、必要な時のみ使うこと。

[0214 ごろ 2011/11/25(金) 20:16:44.10]

Lecture3-8
8. Low-level parsing
tokenize "文字列", parse("区切り文字のリスト")
これを実行すると区切り文字で文字列を分割し、前からマクロ`1'`2'へ代入していく
例
tokenize "abc+d efg+3*8", parse(" +*")
di "`1'" -> abc
di "`2'" -> +
di "`3'" -> d
di "`4'" -> efg
di "`5'" -> +
di "`6'" -> 3
di "`7'" -> *
di "`8'" -> 8
di "`9'" ->
di "`10'"->
tokenize "文字列"
parse以下を省略した場合は,parse(" ")と同じでスペースで分割となる。
スペースで分割してもスペース自体はマクロには代入されない。
スペースは重複してもひとつのスペースと同じく扱われるが、スペース以外の区切り文字では、その回数分マクロに代入される。

レクチャーと関係ないが
inputコマンドについて
文字列を入力する時は、str3 文字列変数　と指定する必要がある。strの後ろの数は文字列の数
欠損値は.で入力

[0215 ごろ 2011/11/25(金) 20:24:24.96]

>>208
見逃していました！
このレクチャーはプログラミングが主体のレクチャーなので回帰分析の細かい手法の正否は分かりません。
概要は次の通りです。

例えばsysuse autoのデータで行うとすると
mpg= a*weight^2+b*weight+c*foreign+d
という二次式の回帰曲線モデルを考えます。
この変数a,b,c,dを
gen weight2=weight*weight
regress weight2 weight foreign
として、変数を決めていく手法が詳細に説明されていました。

[0216 ごろ 2011/11/25(金) 23:43:18.19]

Lecture3-9
9. Programming immediate commands
気が付かなかった。プログラム名は大文字小文字を区別する。
macro shift 数
`1' `2' `3' `4'とあると
macro shift 2
とすると`1' `2'が削除されて`3' `4'がそれぞれ`1' `2'になる`0'には影響しない。

3-10
10. Rewriting mytt in terms of mytti
これまで学んだことの応用問題
プログラムの中に小さなプログラムを組み込むことでわかりやすくすることと、
プログラムを変更する時に、同じ動作部分はひとつのプログラム(サブルーチン）に
まとめることで、一つの修正で全てが修正できるようになるという内容

3-11
11. Parsing new variables
syntax newvarnameコマンドで新しい変数を作成する方法
途中でbreakした時に不完全な新しい変数が出来るのを防ぐ方法
一時的な変数名を作成して、全てが終わった後、新しい変数名に変更する。

[0217 ごろ 2011/11/26(土) 00:06:54.14]

Lecture3-12
練習問題
大部分が自作ｔテストのプログラムの書き換え問題。

[0218 ごろ 2011/11/26(土) 01:25:15.02]

やっとLecture4 1/14
4-1. Where are we?
このレクチャー３までのまとめ

4-2. Saved results
program プログラム名, rclass
とすることでこのプログラムがreturn listで
r(N)・・・などの値を返す様になる。
programの中にreturn listを返すようなプログラム（例sumなど）を使うと
return addを加えることでプログラム中に最初に実行した結果がr(N)・・・に
入力される。
自分で出力させたい値を作る場合は
プログラム中に
return scalar 名前 = 値で
r(名前)=値と表示されるように成る。
returnコマンドはプログラム中のどこに置いても良い

4-3. What can be returned in r()?
scalar, macro, matrixが返せる。

mac osxの場合
matrix A = (1,2\3,4)　ではダメで
matrix A = (1,2\3,4)　でないといけない。

4-4. Referring to returned results in other programs
di _nで改行
mat lsit A /*Aはmatrix行列名 */
でmatrixそのものを表示

4-5. Referring to returned results in the program that sets them
r()はプログラムの最後に代入される。return(sum)などはプログラム途中でも
代入される。

[0219 １３２人目の素数さん 2011/11/26(土) 19:42:35.40]

>>215
わかりました。
ただ、それは、weight*weightではなく、c.weight#c.weightという項を入れるといいんですよね。

[0220 ごろ 2011/11/27(日) 01:02:59.86]

4-6. Other types of returned values: s() and e()
r() rclassとs() sclass e() eclassは同時には使えない

4-7. s-class returned values
プログラムの引数にr(mean)などr-classのコマンドが入ってくる
可能性のあるプログラムを作成する際、r-classでプログラムを作成すると
不具合が生じることがある。そのためにs-classは存在する。
しかし、使用頻度は稀。s-classはマクロにしか値を返せない。
substr(文字列，左から何文字目か（-をつけると右から）,何文字取るか)
substrコマンドで文字列から指定した文字数のみ取り出せる。

4-8. e-class returned results
e-classはregressコマンドなどの結果を表示するclass
e(cmd) コマンドそのもの
e(ll) the log-likelihood value
e(b) the coefficient value
e(V) 変数マトリックス

[0221 ごろ 2011/11/27(日) 01:08:16.72]

>>219
ありがとうございます。
ちなみに
>>215
> regress weight2 weight foreign
> として、変数を決めていく手法が詳細に説明されていました。
regress mpg weight2 weight foreign
mpgが抜けていました。

[0222 SEMやってます 2011/11/28(月) 23:16:18.42]

何レスか前にSEMを始めた者です。
大石展緒・都竹浩生『Amosで学ぶ調査系データ解析』の1章にあった、
食物摂取と、大腸ガン、直腸ガンの関連のデータで、同じようにして推定してみようと思ったのですが、推定できません。
このデータ自体は豊田・前田・柳井『原因をさぐる統計学』(1992)のデータだそうなのですが。

総熱量→洋食傾向(潜在変数)
肉類→洋食傾向
乳製品→洋食傾向
酒類→洋食傾向

総熱量・肉類・乳製品・酒類にすべて共分散

洋食傾向→下部消化管のガン傾向(潜在変数)

下部消化管のガン傾向(潜在変数)→大腸ガン
下部消化管のガン傾向(潜在変数)→直腸ガン

で図を描くと、攪乱項を自動的につけてくれるので、その状態ですべてデフォルトでEstimationするものの、エラーが出ます。

「洋食傾向」を削除し、潜在変数を「下部消化管のガン傾向」一つにしてやってみると、推定完了しますが、文献とは数値が合いません。

もう少し調べます。

[0223 １３２人目の素数さん 2011/11/29(火) 17:24:16.59]

>>222

同書を見たことが無いので、詳しいことは分かりませんが、エラーは何と出ているのでしょうか？

洋食傾向（潜在変数）は実際に計量された変数から導き出されているようですが、
下部消化管のガン傾向（潜在変数）は、洋食傾向と同様、計量された変数によって
作られている潜在変数でしょうか？

下部消化管のガン傾向（潜在変数）を構成している変数の指定はあっていますか？

[0224 １３２人目の素数さん 2011/11/29(火) 22:18:55.99]

>>223
レスどうもです。
中部大学の小塩さんのページ
http://psy.isc.chubu.ac.jp/~oshiolab/teaching_folder/datakaiseki_folder/add_folder/daad_02.html
の下から3分の1くらいのところにある「PLSモデル」と同じ形です。
このページでいうところの「知的能力」にあたる部分が「洋食傾向」で、「対人関係能力」にあたるところが「下部消化管のガン傾向」になります。
「知性」「信頼性」「素直さ」のところが、「総熱量」「肉類」「乳製品」「酒類」で、「外向性」「社交性」「積極性」のところが「大腸ガン」と「直腸ガン」になります。

エラーメッセージは、ダイアログボックスに、Estimation failedと出、メイン画面上には、
_sem_iv(): 3301 subscript invalid
_sem_build__params(): - function returned error
_sem_build(): - function returned error
_sem_parse(): - function returned error
st_sem_parse_and_build(): - function returned error
<istmt>: - function returned error
r(3301);
と出ています。

AMOSでは、「洋食傾向」(小塩さんの図では「知的能力」にあたる)にささる攪乱項を外すことができるのに対し、StataのSEMビルダー上では攪乱項を外すことができない、という状況です。
何か設定があるのか、コマンドでやればいけるのか、というのを調べるつもりなのですが、今日はそろそろ寝ます。

[0225 SEM 2011/12/04(日) 19:57:51.01]

どうもうまくいかないので、ちょっと違うことをして遊んでいます。
sysuse auto.dta
reg mpg weight foreign
で重回帰分析ができるわけですが、同様のことをSEMを使ってやってみますと、
clear
set more off
sysuse auto.dta
sem (mpg <- weight foreign) , nocapslatent
estat eqgof
disp "R2=" r(CD)
matrix eqf=r(eqfit)
disp "F=" (eqf[1,2]/(e(df_m)-2))/(eqf[1,3]/(e(N)-e(df_m)+2-1))
disp "p=" Ftail((e(df_m)-2),(e(N)-e(df_m)+2-1),(eqf[1,2]/(e(df_m)-2))/(eqf[1,3]/(e(N)-e(df_m)+2-1)))
set more on
ですが、推定方法が違うためか、若干数値に差が出ますね。

[0226 １３２人目の素数さん 2011/12/13(火) 15:14:09.22]

ワイは日本人やが聞いた情報によると、もうじき中国はバブルがはじけて昔の貧乏な中国に戻るらしいで
みんなも知っての通りもう経済は破綻してて、取り戻すのは無理なんだそうや


その世界ではごっつい有名な政府関係者筋から聞いた確かな情報やで

まあお前らほどの頭の良い連中には、今さらなくらいのネタやな、かえって失礼かもしれん
お前らからすればもう常識的なくらいの知識やろ！！

[0227 １３２人目の素数さん 2012/01/02(月) 04:55:30.74]

Hello,King!

[0228 ごろ 2012/01/05(木) 11:05:07.39]

明けましておめでとうございます。
今年は皆様にとって素晴らしい一年になることをお祈りいたします。

Stataオンラインレクチャー途中でほったらかしているので時間を見つけて終わらせます。

[0229 あんでぃ ◆AfX9iez2fA 2012/02/04(土) 23:51:25.79]

あんでぃ

[0230 １３２人目の素数さん 2012/02/05(日) 21:56:05.19]

　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　 ／:::::; -‐''"　　　　　 　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　 　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　 　 　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　 　 　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　私は只の数ヲタなんかとは付き合わないわ。
.　　　 　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　頭が良くて数学が出来てかっこいい人。それが必要条件よ。
　　　　 | 　　　　l＾,人| 　` `-' 　 　 ゝ　 |　　　　　　　　さらに Ann.of Math に論文書けば十分条件にもなるわよ。
　 　 　 | 　　　　 ` -'＼　　　 　 　ｰ' 　人　　　　　　　　　　一番嫌いなのは論文数を増やすためにくだらない論文を書いて
　　　　|　　　　　　　 /(l 　　 　＿＿／　 ヽ、　　　　　　　　　　　良い論文の出版を遅らせるお馬鹿な人。
　　 　 |　　　　　　　（:::::`‐-､__ 　|::::`､ 　 　 ﾋﾆﾆヽ、　　　　　　　　　あなたの論文が Ann of Math に accept される確率は?
　 　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　それとも最近は Inv. Math. の方が上かしら？
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　 ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ 　 |　'､　　 ＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ 　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/

[0231 １３２人目の素数さん 2012/02/09(木) 18:25:45.89]

静まりかえっていますが、Stata/Transferってソフトの案内が来た。
ただのデータフォーマット変換ソフトに金とるのね。
おそれいった。

[0232 １３２人目の素数さん 2012/02/25(土) 16:52:29.31]

　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　 ／:::::; -‐''"　　　　　 　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　 　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　 　 　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　 　 　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　私は只の数ヲタなんかとは付き合わないわ。
.　　　 　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　頭が良くて数学が出来てかっこいい人。それが必要条件よ。
　　　　 | 　　　　l＾,人| 　` `-' 　 　 ゝ　 |　　　　　　　　さらに Ann.of Math に論文書けば十分条件にもなるわよ。
　 　 　 | 　　　　 ` -'＼　　　 　 　ｰ' 　人　　　　　　　　　　一番嫌いなのは論文数を増やすためにくだらない論文を書いて
　　　　|　　　　　　　 /(l 　　 　＿＿／　 ヽ、　　　　　　　　　　　良い論文の出版を遅らせるお馬鹿な人。
　　 　 |　　　　　　　（:::::`‐-､__ 　|::::`､ 　 　 ﾋﾆﾆヽ、　　　　　　　　　あなたの論文が Ann of Math に accept される確率は?
　 　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　それとも最近は Inv. Math. の方が上かしら？
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　 ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ 　 |　'､　　 ＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ 　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/

[0233 １３２人目の素数さん 2012/03/08(木) 19:21:43.89]

STATAを今使用して重回帰分析を行なっているんですが、非標準化回帰係数の95%信頼区間は表示されるのですが、"標準化回帰係数"の95%信頼区間はどうやったら表示されるのでしょうか？
ご存知の方、よろしくお願いします。

[0234 １３２人目の素数さん 2012/03/08(木) 19:24:17.96]

マルチポストに答える気はなくなりました。

[0235 １３２人目の素数さん 2012/03/09(金) 21:47:03.78]

久々に。
よく、調査データで、「わからない」を88、「無回答」を99なんてコーディングして、ともに欠損値にしたいなんてことがありますが。
Stataで欠損値は「 . 」ということはわかったのですが、「わからない」と「無回答」は区別しておきたいわけです。
解説本を見ても書いてないし、どうするのかと思ってましたが、helpであっさりわかりました。
Stataの欠損値は「 . 」だけではなく、「 .a 」、「 .b 」、・・・と使えるのだと。
recode var1 99=.a 88=.b
として、無回答99を.a、わからない88を.bとしまして、
ラベルとして、
label define example .a "NA" .b "DK"
label values var1 example

と貼り付けておけば、あとで見たときに区別がつくわけです。

[0236 １３２人目の素数さん 2012/03/10(土) 17:38:11.84]

知 ら な か っ た

[0237 １３２人目の素数さん 2012/03/23(金) 02:35:06.57]

　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　 ／:::::; -‐''"　　　　　 　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　 　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　 　 　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　 　 　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　私は只の数ヲタなんかとは付き合わないわ。
.　　　 　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　頭が良くて数学が出来てかっこいい人。それが必要条件よ。
　　　　 | 　　　　l＾,人| 　` `-' 　 　 ゝ　 |　　　　　　　　さらに Ann.of Math に論文書けば十分条件にもなるわよ。
　 　 　 | 　　　　 ` -'＼　　　 　 　ｰ' 　人　　　　　　　　　　一番嫌いなのは論文数を増やすためにくだらない論文を書いて
　　　　|　　　　　　　 /(l 　　 　＿＿／　 ヽ、　　　　　　　　　　　良い論文の出版を遅らせるお馬鹿な人。
　　 　 |　　　　　　　（:::::`‐-､__ 　|::::`､ 　 　 ﾋﾆﾆヽ、　　　　　　　　　あなたの論文が Ann of Math に accept される確率は?
　 　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　それとも最近は Inv. Math. の方が上かしら？
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　 ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ 　 |　'､　　 ＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ 　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/

[0238 １３２人目の素数さん 2012/03/25(日) 23:15:22.24]

　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　 ／:::::; -‐''"　　　　　 　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　 　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　 　 　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　 　 　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　私は只の数ヲタなんかとは付き合わないわ。
.　　　 　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　頭が良くて数学が出来てかっこいい人。それが必要条件よ。
　　　　 | 　　　　l＾,人| 　` `-' 　 　 ゝ　 |　　　　　　　　さらに Ann.of Math に論文書けば十分条件にもなるわよ。
　 　 　 | 　　　　 ` -'＼　　　 　 　ｰ' 　人　　　　　　　　　　一番嫌いなのは論文数を増やすためにくだらない論文を書いて
　　　　|　　　　　　　 /(l 　　 　＿＿／　 ヽ、　　　　　　　　　　　良い論文の出版を遅らせるお馬鹿な人。
　　 　 |　　　　　　　（:::::`‐-､__ 　|::::`､ 　 　 ﾋﾆﾆヽ、　　　　　　　　　あなたの論文が Ann of Math に accept される確率は?
　 　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　それとも最近は Inv. Math. の方が上かしら？
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　 ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ 　 |　'､　　 ＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ 　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/

[0239 １３２人目の素数さん 2012/03/26(月) 12:40:49.52]

　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　 ／:::::; -‐''"　　　　　 　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　 　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　 　 　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　 　 　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　私は只の数ヲタなんかとは付き合わないわ。
.　　　 　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　頭が良くて数学が出来てかっこいい人。それが必要条件よ。
　　　　 | 　　　　l＾,人| 　` `-' 　 　 ゝ　 |　　　　　　　　さらに Ann.of Math に論文書けば十分条件にもなるわよ。
　 　 　 | 　　　　 ` -'＼　　　 　 　ｰ' 　人　　　　　　　　　　一番嫌いなのは論文数を増やすためにくだらない論文を書いて
　　　　|　　　　　　　 /(l 　　 　＿＿／　 ヽ、　　　　　　　　　　　良い論文の出版を遅らせるお馬鹿な人。
　　 　 |　　　　　　　（:::::`‐-､__ 　|::::`､ 　 　 ﾋﾆﾆヽ、　　　　　　　　　あなたの論文が Ann of Math に accept される確率は?
　 　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　それとも最近は Inv. Math. の方が上かしら？
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　 ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ 　 |　'､　　 ＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ 　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/

[0240 １３２人目の素数さん 2012/03/26(月) 16:41:13.45]

　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　 ／:::::; -‐''"　　　　　 　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　 　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　 　 　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　 　 　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　私は只の数ヲタなんかとは付き合わないわ。
.　　　 　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　頭が良くて数学が出来てかっこいい人。それが必要条件よ。
　　　　 | 　　　　l＾,人| 　` `-' 　 　 ゝ　 |　　　　　　　　さらに Ann.of Math に論文書けば十分条件にもなるわよ。
　 　 　 | 　　　　 ` -'＼　　　 　 　ｰ' 　人　　　　　　　　　　一番嫌いなのは論文数を増やすためにくだらない論文を書いて
　　　　|　　　　　　　 /(l 　　 　＿＿／　 ヽ、　　　　　　　　　　　良い論文の出版を遅らせるお馬鹿な人。
　　 　 |　　　　　　　（:::::`‐-､__ 　|::::`､ 　 　 ﾋﾆﾆヽ、　　　　　　　　　あなたの論文が Ann of Math に accept される確率は?
　 　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　それとも最近は Inv. Math. の方が上かしら？
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　 ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ 　 |　'､　　 ＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ 　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/

[0241 １３２人目の素数さん 2012/03/26(月) 17:16:15.84]

　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　 ／:::::; -‐''"　　　　　 　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　 　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　 　 　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　 　 　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　私は只の数ヲタなんかとは付き合わないわ。
.　　　 　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　頭が良くて数学が出来てかっこいい人。それが必要条件よ。
　　　　 | 　　　　l＾,人| 　` `-' 　 　 ゝ　 |　　　　　　　　さらに Ann.of Math に論文書けば十分条件にもなるわよ。
　 　 　 | 　　　　 ` -'＼　　　 　 　ｰ' 　人　　　　　　　　　　一番嫌いなのは論文数を増やすためにくだらない論文を書いて
　　　　|　　　　　　　 /(l 　　 　＿＿／　 ヽ、　　　　　　　　　　　良い論文の出版を遅らせるお馬鹿な人。
　　 　 |　　　　　　　（:::::`‐-､__ 　|::::`､ 　 　 ﾋﾆﾆヽ、　　　　　　　　　あなたの論文が Ann of Math に accept される確率は?
　 　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　それとも最近は Inv. Math. の方が上かしら？
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　 ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ 　 |　'､　　 ＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ 　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/

[0242 １３２人目の素数さん 2012/03/26(月) 18:15:41.40]

　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　 ／:::::; -‐''"　　　　　 　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　 　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　 　 　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　 　 　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　私は只の数ヲタなんかとは付き合わないわ。
.　　　 　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　頭が良くて数学が出来てかっこいい人。それが必要条件よ。
　　　　 | 　　　　l＾,人| 　` `-' 　 　 ゝ　 |　　　　　　　　さらに Ann.of Math に論文書けば十分条件にもなるわよ。
　 　 　 | 　　　　 ` -'＼　　　 　 　ｰ' 　人　　　　　　　　　　一番嫌いなのは論文数を増やすためにくだらない論文を書いて
　　　　|　　　　　　　 /(l 　　 　＿＿／　 ヽ、　　　　　　　　　　　良い論文の出版を遅らせるお馬鹿な人。
　　 　 |　　　　　　　（:::::`‐-､__ 　|::::`､ 　 　 ﾋﾆﾆヽ、　　　　　　　　　あなたの論文が Ann of Math に accept される確率は?
　 　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　それとも最近は Inv. Math. の方が上かしら？
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　 ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ 　 |　'､　　 ＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ 　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/

[0243 １３２人目の素数さん 2012/04/22(日) 12:54:52.37]

　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　 ／:::::; -‐''"　　　　　 　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　 　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　 　 　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　 　 　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　私は只の数ヲタなんかとは付き合わないわ。
.　　　 　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　頭が良くて数学が出来てかっこいい人。それが必要条件よ。
　　　　 | 　　　　l＾,人| 　` `-' 　 　 ゝ　 |　　　　　　　　さらに Ann.of Math に論文書けば十分条件にもなるわよ。
　 　 　 | 　　　　 ` -'＼　　　 　 　ｰ' 　人　　　　　　　　　　一番嫌いなのは論文数を増やすためにくだらない論文を書いて
　　　　|　　　　　　　 /(l 　　 　＿＿／　 ヽ、　　　　　　　　　　　良い論文の出版を遅らせるお馬鹿な人。
　　 　 |　　　　　　　（:::::`‐-､__ 　|::::`､ 　 　 ﾋﾆﾆヽ、　　　　　　　　　あなたの論文が Ann of Math に accept される確率は?
　 　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　それとも最近は Inv. Math. の方が上かしら？
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　 ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ 　 |　'､　　 ＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ 　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/

[0244 １３２人目の素数さん 2012/04/22(日) 18:36:52.52]

　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　 ／:::::; -‐''"　　　　　 　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　 　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　 　 　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　 　 　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　私は只の数ヲタなんかとは付き合わないわ。
.　　　 　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　頭が良くて数学が出来てかっこいい人。それが必要条件よ。
　　　　 | 　　　　l＾,人| 　` `-' 　 　 ゝ　 |　　　　　　　　さらに Ann.of Math に論文書けば十分条件にもなるわよ。
　 　 　 | 　　　　 ` -'＼　　　 　 　ｰ' 　人　　　　　　　　　　一番嫌いなのは論文数を増やすためにくだらない論文を書いて
　　　　|　　　　　　　 /(l 　　 　＿＿／　 ヽ、　　　　　　　　　　　良い論文の出版を遅らせるお馬鹿な人。
　　 　 |　　　　　　　（:::::`‐-､__ 　|::::`､ 　 　 ﾋﾆﾆヽ、　　　　　　　　　あなたの論文が Ann of Math に accept される確率は?
　 　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　それとも最近は Inv. Math. の方が上かしら？
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　 ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ 　 |　'､　　 ＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ 　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/

[0245 １３２人目の素数さん 2012/04/26(木) 01:14:27.15]

　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　 ／:::::; -‐''"　　　　　 　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　 　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　 　 　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　 　 　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　私は只の数ヲタなんかとは付き合わないわ。
.　　　 　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　頭が良くて数学が出来てかっこいい人。それが必要条件よ。
　　　　 | 　　　　l＾,人| 　` `-' 　 　 ゝ　 |　　　　　　　　さらに Ann.of Math に論文書けば十分条件にもなるわよ。
　 　 　 | 　　　　 ` -'＼　　　 　 　ｰ' 　人　　　　　　　　　　一番嫌いなのは論文数を増やすためにくだらない論文を書いて
　　　　|　　　　　　　 /(l 　　 　＿＿／　 ヽ、　　　　　　　　　　　良い論文の出版を遅らせるお馬鹿な人。
　　 　 |　　　　　　　（:::::`‐-､__ 　|::::`､ 　 　 ﾋﾆﾆヽ、　　　　　　　　　あなたの論文が Ann of Math に accept される確率は?
　 　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　それとも最近は Inv. Math. の方が上かしら？
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　 ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ 　 |　'､　　 ＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ 　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/

[0246 １３２人目の素数さん 2012/04/27(金) 01:11:42.78]

　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　 ／:::::; -‐''"　　　　　 　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　 　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　 　 　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　 　 　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　私は只の数ヲタなんかとは付き合わないわ。
.　　　 　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　頭が良くて数学が出来てかっこいい人。それが必要条件よ。
　　　　 | 　　　　l＾,人| 　` `-' 　 　 ゝ　 |　　　　　　　　さらに Ann.of Math に論文書けば十分条件にもなるわよ。
　 　 　 | 　　　　 ` -'＼　　　 　 　ｰ' 　人　　　　　　　　　　一番嫌いなのは論文数を増やすためにくだらない論文を書いて
　　　　|　　　　　　　 /(l 　　 　＿＿／　 ヽ、　　　　　　　　　　　良い論文の出版を遅らせるお馬鹿な人。
　　 　 |　　　　　　　（:::::`‐-､__ 　|::::`､ 　 　 ﾋﾆﾆヽ、　　　　　　　　　あなたの論文が Ann of Math に accept される確率は?
　 　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　それとも最近は Inv. Math. の方が上かしら？
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　 ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ 　 |　'､　　 ＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ 　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/

[0247 １３２人目の素数さん 2012/04/28(土) 15:07:19.94]

　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　 ／:::::; -‐''"　　　　　 　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　 　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　 　 　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　 　 　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　私は只の数ヲタなんかとは付き合わないわ。
.　　　 　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　頭が良くて数学が出来てかっこいい人。それが必要条件よ。
　　　　 | 　　　　l＾,人| 　` `-' 　 　 ゝ　 |　　　　　　　　さらに Ann.of Math に論文書けば十分条件にもなるわよ。
　 　 　 | 　　　　 ` -'＼　　　 　 　ｰ' 　人　　　　　　　　　　一番嫌いなのは論文数を増やすためにくだらない論文を書いて
　　　　|　　　　　　　 /(l 　　 　＿＿／　 ヽ、　　　　　　　　　　　良い論文の出版を遅らせるお馬鹿な人。
　　 　 |　　　　　　　（:::::`‐-､__ 　|::::`､ 　 　 ﾋﾆﾆヽ、　　　　　　　　　あなたの論文が Ann of Math に accept される確率は?
　 　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　それとも最近は Inv. Math. の方が上かしら？
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　 ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ 　 |　'､　　 ＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ 　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/

[0248 １３２人目の素数さん 2012/04/28(土) 23:40:17.01]

　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　 ／:::::; -‐''"　　　　　 　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　 　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　 　 　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　 　 　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　私は只の数ヲタなんかとは付き合わないわ。
.　　　 　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　頭が良くて数学が出来てかっこいい人。それが必要条件よ。
　　　　 | 　　　　l＾,人| 　` `-' 　 　 ゝ　 |　　　　　　　　さらに Ann.of Math に論文書けば十分条件にもなるわよ。
　 　 　 | 　　　　 ` -'＼　　　 　 　ｰ' 　人　　　　　　　　　　一番嫌いなのは論文数を増やすためにくだらない論文を書いて
　　　　|　　　　　　　 /(l 　　 　＿＿／　 ヽ、　　　　　　　　　　　良い論文の出版を遅らせるお馬鹿な人。
　　 　 |　　　　　　　（:::::`‐-､__ 　|::::`､ 　 　 ﾋﾆﾆヽ、　　　　　　　　　あなたの論文が Ann of Math に accept される確率は?
　 　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　それとも最近は Inv. Math. の方が上かしら？
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　 ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ 　 |　'､　　 ＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ 　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/

[0249 １３２人目の素数さん 2012/04/29(日) 15:01:57.99]

　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　 ／:::::; -‐''"　　　　　 　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　 　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　 　 　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　 　 　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　私は只の数ヲタなんかとは付き合わないわ。
.　　　 　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　頭が良くて数学が出来てかっこいい人。それが必要条件よ。
　　　　 | 　　　　l＾,人| 　` `-' 　 　 ゝ　 |　　　　　　　　さらに Ann.of Math に論文書けば十分条件にもなるわよ。
　 　 　 | 　　　　 ` -'＼　　　 　 　ｰ' 　人　　　　　　　　　　一番嫌いなのは論文数を増やすためにくだらない論文を書いて
　　　　|　　　　　　　 /(l 　　 　＿＿／　 ヽ、　　　　　　　　　　　良い論文の出版を遅らせるお馬鹿な人。
　　 　 |　　　　　　　（:::::`‐-､__ 　|::::`､ 　 　 ﾋﾆﾆヽ、　　　　　　　　　あなたの論文が Ann of Math に accept される確率は?
　 　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　それとも最近は Inv. Math. の方が上かしら？
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　 ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ 　 |　'､　　 ＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ 　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/

[0250 ごろ 2012/04/30(月) 04:46:13.30]

あ。最近忘れてた。

[0251 あぼーん NG NG]

あぼーん

[0252 あぼーん NG NG]

あぼーん

[0253 あぼーん NG NG]

あぼーん

[0254 あぼーん NG NG]

あぼーん

[0255 あぼーん NG NG]

あぼーん

[0256 あぼーん NG NG]

あぼーん

[0257 あぼーん NG NG]

あぼーん

[0258 あぼーん NG NG]

あぼーん

[0259 あぼーん NG NG]

あぼーん

[0260 あぼーん NG NG]

あぼーん

[0261 あぼーん NG NG]

あぼーん

[0262 あぼーん NG NG]

あぼーん

[0263 あぼーん NG NG]

あぼーん

[0264 あぼーん NG NG]

あぼーん

[0265 あぼーん NG NG]

あぼーん

[0266 １３２人目の素数さん 2012/06/24(日) 17:02:54.94]

Multilevel and Longitudinal Modeling Using Stata
つう本を買いました。
xtコマンド関連の勉強をしようと思いまして。

[0267 あぼーん NG NG]

あぼーん

[0268 ごろ 2012/07/05(木) 23:17:10.65]

>>266
勉強して教えてください。
教えることは最高の学びともいいますし。

(^O^)／ただで勉強しようと思うなヽ(`Д´)ﾉﾌﾟﾝﾌﾟﾝ

[0269 初心者 2012/07/21(土) 15:38:27.94]

はじめまして。STATAのコマンドを調べていて辿り着きました。
そして早速ですが質問させてください。
グラフの凡例（legend）を右に表示させ、かつラベルを縦1列に並べたいのですが、この縦1列に並べるコマンドが分かりません。
helpを参照するとlegendオプションのstyle( )というのが怪しいのですが上手くいかず。
Graph EditorでlegendのColumnsを1にすれば目的の形にはできるのですが、コマンドが知りたいです。
ご存知の方がいらっしゃいましたらどうぞよろしくお願い致します。

参考までにサンプルデータで作ったものを書いておきます。
・データ
year group freq
2000 1 1
2010 1 2
2000 2 3
2010 2 4
2000 3 5
2010 3 6
2000 4 7
2010 4 8

・コマンド
twoway (line freq year if group==1)/*
*/ ( line freq year if group==2)/*
*/ ( line freq year if group==3)/*
*/ ( line freq year if group==4),/*
*/ legend(position(3))

4本の線グラフが表示される。
今の状態ではグラフ右側に2×2の形でラベルが表示されるがこれを1列にしたい。

[0270 １３２人目の素数さん 2012/07/21(土) 16:19:04.66]

>>269
最後のposition(3)の次にスペース1個入れてc(1)ではいかが？
legend(position(3) c(1))

[0271 初心者 2012/07/21(土) 17:02:20.38]

>>270
269です。
素早いご解答ありがとうございます
そして無事c(1)でできました。感動！
column(1)とか色々打ち込んでみてたのですが、c( )だったんですね。

ちなみにそういう情報はSTATAのhelpなどで検索できますか？
コマンドを調べるのにいつも手間取っていて。。
後学のため、ご存知でしたら教えて頂けると大変助かります。

[0272 １３２人目の素数さん 2012/07/21(土) 18:36:26.76]

>>271
省略しない場合、colsのようです。
検索は、Help→Searchで、適宜どうぞ。今回の場合はlegendで。

[0273 馬鹿を焼く描写 ◆ghclfYsc82 2012/07/22(日) 12:01:58.85]

勉強や努力が足りなくて優秀になれない奴が惨めな思いをするのは当然
なんだよ。それを自分で何もせずに優秀な人間の足を引っ張るとは言語
道断である。他人を貶めるだけで自分は楽をする奴は恥を知れ。今後も
そういう馬鹿者を発見次第、即刻攻撃を掛けて当該スレを焼け野が原に
するので、覚悟をする様に願いたい。こういう考え方が国家を滅ぼす。
無能な馬鹿は自滅するに任せ、優秀な人材こそを選択的に抽出し、それ
を国家が意図して保護しなければならない。そうする事が国家が生き残
る唯一の道である。繰り返す。何の努力もしない馬鹿を無条件に保護す
れば、その結果として誰も努力しなくなるだけである。だから馬鹿を保
護しては絶対にならない。

描

>みんなで優秀な人間の足を引っ張って沈もうよ。
>そうすれば自分だけが馬鹿で惨めな思いをしなくて
>すむから楽チン。
>一億総白痴可で横並びになれば怖くは無い
>

[0274 初心者 2012/07/25(水) 11:28:30.70]

>>272
ありがとうございます。
help確認してみました。確かにcolsですね。
私もlegendのヘルプは何度か見ていたのですが、cols見逃していました。。
お手数をおかけしてすみません。

syntaxをいつも頑張って解読しようとするのですがなかなか力及ばず。
helpやマニュアル（黒い冊子）片手にこれからも頑張ります。

初心者なりに最近驚いたことを一つ。
stataで小数点を含む数値を扱うときはfloat()関数が必要。
count if 変数名==0.1　ではなく　count if 変数名==float(0.1)　としなければならない。

気づくまでは数が合わなくてオロオロしてました。

[0275 描者の品格 ◆ghclfYsc82 2012/07/25(水) 12:15:25.77]

描

>367 ：匿名希望：2012/07/23(月) 17:23:38.69
> >>365
> 頭悪いのはお前の方だろ？
> 「猫」という字を「描」に間違えやがってｗｗｗ
> 小学生並みの頭の悪さだぞ、お前！
> 焼かれるのはお前の方だろ？
> 外国では猫の丸焼きというゲテモノ料理もあるらしいぜｗｗｗ
>

[0276 １３２人目の素数さん 2012/09/08(土) 19:56:43.18]

STATAで、重回帰分析にpropensity scoreって使えるんですかね？
2群比較やロジスティックはわかるんですけど。

[0277 １３２人目の素数さん 2012/09/09(日) 01:44:35.92]

猫は？

[0278 １３２人目の素数さん 2012/09/09(日) 09:35:57.51]

やつは中の人

[0279 ごろ 2012/10/01(月) 18:16:35.89]

お久しぶりです。
>>276
＞propensity score使えるようですよ
SPSSでは一応使えるけどっていうレベルみたい。SPSSの中の人が言っていた。
やっぱりRかな。でもStataいいっす。

[0280 ごろ 2012/10/02(火) 17:21:13.92]

>>274
小数点問題
ちょうど私もハマってしまいました。
ここにまとめておくと調べるのにも便利ですね。
情報ありがとうございます。
この情報の公式な発表とかありますか？

ちょっとググったら、
ttp://www.google.co.jp/url?sa=t&rct=j&q=&esrc=s&source=web&cd=20&ved=0CHkQFjAJOAo&url=http%3A%2F%2Fwww.rpi.edu%2F~simonk%2Fpdf%2FUsefulStataCommands.pdf&ei=H6NqUP3UIe3LmAX0kYDABQ&usg=AFQjCNGWGyuzTiMvneGJkgFjA5Wr38qvYw&sig2=PP6RQkXbjVehu3TUwAvPgQ
のここに書いてありました。
K.Data Types,When 3.3 3.3,andMissingValues

[0281 ごろ 2012/10/09(火) 15:41:22.02]

Stataについて書いているブログを見つけた
http://www.masabu.com/category/stata/
勉強させてもらおう

[0282 ごろ 2012/10/09(火) 17:21:41.65]

foreach n of numlist 1(1)255 {
di `n' , char(`n')
}
charコマンドで特殊な文字列を表現できる。でも日本語フォントはダメみたい。

ずっと文字列を変数名と認識させる方法を探しているが見つからない。
. input testdesu
testdesu
1. 1009
2. 9001
3. end

. di testdesu
1009

. di ("test"+"desu")
testdesu

最後のところでちゃんと1009と表示させる方法を探している。
. di ("test"+"desu")ここは実際はもっと複雑なstring functionコマンドです。
単純な方法がありそうなんだけどなぁ

[0283 あぼーん NG NG]

あぼーん

[0284 １３２人目の素数さん 2012/10/09(火) 18:46:38.77]

　ε⌒ﾍ ⌒ヽﾌ
（ 　　( 　・ω・）　ブヒブヒ
　 しー し─Ｊ

[0285 ごろ 2012/10/10(水) 12:41:57.29]

>>282
解決しました！
local i = "test"+"desu"
di `i'
でOKでした。
単純でしたね。

[0286 ごろ 2012/10/11(木) 11:04:15.04]

*log using forabstract.log, txt append /* easy to read, but Japanese will dissapear*/
do "C:\Users\MSC12-03\Dropbox\HOME\2013JPS\forconfidence.do"

replace confidence=. if confidence==999
save forconfidence2012-10-10.dta,replace

ttest confidence ,by(experience)
ranksum confidence, by(experience)

[0287 ごろ 2012/10/11(木) 11:07:45.12]

Doファイルの中でDoファイルを作成して実行する方法
capture program drop makedofile
program makedofile
syntax anything
qui{ /*これがポイント。これでログファクィルのヘッダーを出力しない。*/
capture log close
log using newdofile.do , text replace
}
di "clear"
di "input str10 newvarname data1 data2"
foreach i of local anything {
local x = _N
foreach j of numlist 1(1)`x' {
di originalvarname " " /*ここにdata1, data2用の命令*/
}
}
di "end"
qui {
log close
}
end

makedofile /*originaldataを呼び出す変数名など*/

[0288 １３２人目の素数さん 2012/10/11(木) 12:39:49.81]

回帰係数がｔ分布で検定できることの説明ができる人、どうぞ

[0289 ごろ 2012/10/11(木) 13:19:56.61]

>>285
di `="test"+"desu"'
でも出来ました〜。

[0290 あぼーん NG NG]

あぼーん

[0291 あぼーん NG NG]

あぼーん

[0292 ごろ 2012/12/26(水) 16:58:41.83]

はやいものですなぁ
もう今年も残す所あとわずか・・

[0293 あぼーん NG NG]

あぼーん

[0294 あぼーん NG NG]

あぼーん

[0295 ごろ 2013/11/11(月) 11:14:41.10]

おお、まだ残っていましたね。
もう1年たつんだ。
>>159
今日散歩していたら急に思いついた。StataでSleepsortする方法。
上手くできたらまた書き込みます。

[0296 １３２人目の素数さん 2013/11/11(月) 14:58:14.93]

>>295
よくよく考えなおしたら、Sleep Sortの理論上、バックグラウンドで複数の関数が実行出来ないと
無理なんですね。あーダメだ。やるとしたらStataから別のプログラムを操作する方法を考えないとダメかな。

[0297 ごろ 2014/03/13(木) 23:06:29.99]

このページ、我ながら本当に役に立つ。
ネットさえ繋がれば、忘れたSyntaxのコツなど検索ですぐに見つかるから助かる。

今日も自分のためにメモ。



/*
あるテキストデータの先頭にある文字列に
特定の文字列(複数可)を含んでいるかどうかを判定し、
新しい変数名にその結果を格納するプログラム

下記はやっと動いた例。
上手く行かなかったコードも残しておきます。
理由は良くわかりません。
分かる方お教え下さい。
*/
長くなったのでプログラム本体は次の投稿で

[0298 ごろ 2014/03/13(木) 23:08:41.03]

capture program drop testprogram
program testprogram
syntax namelist
local newvarname = "`1'"
gen `1' = 0
foreach i of local namelist {
local j= length("`i'")
replace `newvarname' = 1 if substr(textdata,1,`j')== "`i'"
}
list
tab `1'
end
/*
下記の方法では上手く行かなかった。
そのためわざわざnewvarnameというlocal変数を使っています。
replace "`1'" = 1 if substr(textdata,1,`j')== "`i'"
*/
testprogram list_pine pine
testprogram pine_and_apple pine apple

[0299 ごろ 2014/03/13(木) 23:10:05.73]

データはこれを使って下さい。
clear
set more off
input id str10 textdata
1 pine
2 apple
3 pineapple
4 pines
5 "pine juice"
6 "pink juice"
7 orange
8 oranges
9 spine
10 "blue pine"
end

[0300 超初心者 2014/03/15(土) 14:50:47.03]

群間の多重比較に関して質問させて頂きます。
カイ2乗検定で、群間の有意差を出した後にさらに2つの群間において多重比較を行いたいのですが、どうすればよいのでしょうか？
STATAで"tukeyhsd"というコマンドを発見したのですが、どう使っていいのか分かりません。

[0301 ごろ 2014/04/11(金) 05:51:27.61]

おお、放置プレーしてしまいました。
まず、そのパッケージは、標準で無いので、ネットから、インストールする必要があります。
net install tukeyhsd
で、インストール出来ます。
qsturng.adoも必要そうなのでインストールしましょう。
後は、ヘルプ(help tukeyhsd)の通りですが
tukeyhsd 連続変数, dv(カテゴリカル変数)
で出来ます。

例
sysuse citytemp
tukeyhsd region, dv(tempjan)
これで、地域毎に1月の気温の比較が、出来ます。

しかし、多重比較する前に、多重比較が本当に良い選択なのか、研究デザインを考えることをお勧めします。

[0302 １３２人目の素数さん 2014/04/17(木) 19:29:28.10]

stataを使いながらstataの使い方や統計学について知識を得られるサイトと
いうのはないでしょうか？あるいは本はないでしょうか？

[0303 ごろ 2014/04/25(金) 05:55:39.80]

>>302
あまり親切な回答ではないですが、UCLAのサイトとか有名ですね。

[0304 ごろ 2014/04/25(金) 06:13:22.73]

>>302
あ、まじ回答すると、StataのHelpがかなり優秀なんで、
StataのHelpを見ていれば大抵のことは分かります。
サンプルプログラムもあります。

>>2
にもまとめていましたね。

[e-mail]
このメーリングリストの回答はググればネット上でも閲覧できます。
内容が高度なことが多いですが、困ったときにググると大抵回答があります。
それでも見つからない時には、メーリングリストの会員になって、
質問したら良いと思います。

[0305 ごろ 2014/04/25(金) 06:21:03.89]

http://s497044000.onlinehome.us/docs/workshop-2013/Handbook%20of%20Stat%20Analysis%20Using%20Stata%203rd%20Ed.pdf

↑　これはかなり良さ気
Introduction to Statistics using state .pdf
でぐぐったら見つかった

[0306 １３２人目の素数さん 2014/04/25(金) 23:18:51.25]

>>305
>>304
>>303

教えてくれてありがとう

[0307 １３２人目の素数さん 2014/05/05(月) 04:11:39.20]

STATAとRって一般にどっちが優れているんですか？
やっぱり一長一短ですか？

[0308 ごろ 2014/05/13(火) 11:38:51.76]

Rがスラスラと使いこなせる人にはStataは不要かもしれない。
しかし、RはバージョンやOSシステムが変わるだけで挙動が不安定・不確実になるので、
講義で使用する場合など難しいケースが多い。

経済学の最新のパッケージが平易なコマンドで利用できる点など
を考えるとやはりStataは優秀。

[0309 ごろ 2014/07/12(土) 00:39:01.39]

なんか書き込めない。テスト

[0310 ごろ 2014/07/12(土) 00:42:28.91]

URLが変わったせいかな！？

メモ書き
set more on
set more off
実行結果を位置画面毎に止めるか止めないか。

set trace on
set trace off
Stataの処理を一行ずつ丁寧に吐き出してくれるコマンド

Keywords　デバック　モード　一行毎　プログラム　プログラミング　結果をその都度　変数

[0311 ごろ 2014/09/19(金) 10:11:29.60]

どうでもいいけど、面白かったんで。
Stata helpに存在しない関数の説明があった！！
strcat(s1,s2)
Domain s1: strings
Domain s2: strings
Range: strings
Description: There is no strcat() function. Instead the addition operator is used to concatenate strings:
"hello " + "world" = "hello world"
"a" + "b" = "ab"

[0312 ごろ 2014/11/16(日) 09:54:51.58]

Stataの日本語環境どうにかならないかなぁ

[0313 １３２人目の素数さん 2014/11/16(日) 16:03:11.84]

stataを勉強し始めるのに何かいいものないですか？おしえてください

[0314 １３２人目の素数さん 2014/12/29(月) 10:43:19.32 ID:Phz2K+Xn]

>>307
SPSSが学内サーバーを通して使えるので、Rを覚える必要がないと思う学生も多い。
でも一旦卒業するとその環境が継続する人は稀なわけで、Rも併行して覚えておいた方が良いよと言ってる。
STATAも提供されてるけど、稼働率はSPSSよりかなり低い模様。

[0315 イスラム金融系最高指導者遅獄先生−主対性卓上理論より実戦 2014/12/30(火) 16:57:39.87 ID:ISKEAmyr]

統計で時間を早められない。

[0316 イスラム金融系最高指導者遅獄先生−主対性卓上理論より実戦 2014/12/30(火) 16:58:09.54 ID:ISKEAmyr]

それがオジサンおばさんの限界

[0317 １３２人目の素数さん 2015/04/11(土) 19:34:12.48 ID:6j6wzkXy]

Stata14リリース記念あげ

[0318 ごろ 2015/06/24(水) 09:48:34.05 ID:iMJRcMt7]

Stata14最高ですね。日本語対応なんて・・と思っていましたが、Unicodeがちゃんと使えることでデータ解析がこんな居捗るとは思ってもいませんでした。
Stata８から基本コマンドは変わらないので、古くても大丈夫！と思っていましたが、Stata14は早速ポチりました！

[0319 ごろ 2015/07/03(金) 12:30:32.10 ID:VcMQhWuQ]

これはStata分かりにくいだろ！
by class: egen number_of_male=count(gender=="male")

なんと上記では、男性の数を数えるのではなくてクラスの総人数を数えてしまう！
欠損値はカウントしない様子。

クラスごとの男子生徒の人数を出したいなら下記でOKということ。
egen number_of_male=sum(gender=="male"), by(class)

[0320 １３２人目の素数さん 2015/07/03(金) 20:01:47.60 ID:bWkL4jtp]

なぜ有料ソフトをつかいたがるんだろう
R使いなよ

[0321 ごろ 2015/07/17(金) 11:42:38.84 ID:0TOemjko]

>320
それは抜群の安定性があるからです。
Rはバージョンやプラットフォームが変わるととたんに挙動が怪しくなります。

[0322 １３２人目の素数さん 2015/09/12(土) 04:39:58.52 ID:ns6AxouT]

test1

[0323 ごろ 2015/09/13(日) 01:00:54.60 ID:DXV7Mj5O]

>>322
ttest

[0324 ごろ 2015/09/13(日) 11:13:06.46 ID:DXV7Mj5O]

>>321
Stataに不満は少ないんだが、少ない不満のいくつかとして挙げるとすると
1. バージョンが変わると.data形式が変わって読めなくなる。
直近のバージョンだと互換性は保たれているが少し離れると、
複数のバージョンでデータ形式を整えないといけない。
まぁCSVを介したりすればいいんだけどね。

2. 図表の調整はコマンドで全て記録されない。
これは多分開発側も考えていてRecordボタンなど有るんでしょうけどね。

私が使う範囲では、これくらいかな。贅沢な悩みかもしれない。

[0325 ごろ 2015/09/13(日) 11:14:24.47 ID:DXV7Mj5O]

>>324
あ、.dtaだったね。データ形式。

[0326 ごろ 2015/12/04(金) 21:59:25.03 ID:sdVKKjOJ]

お久しぶりです。
Stataの正規表現って最小マッチできないのかな？

[0327 ごろ 2016/01/01(金) 09:36:13.69 ID:xTr3fpeg]

あけましておめでとうございます。今年もよろしくお願い申し上げます。

[0328 １３２人目の素数さん 2016/04/16(土) 22:56:08.77 ID:6RAjlRZE]

公衆衛生の授業で触ってよかったなと思って
大学院に入って使い始めたら周りにStataユーザーが皆無で
困ったなと思ったらHelpが超充実していてなんとかなりそうな気がしてきた。

[0329 １３２人目の素数さん 2016/05/04(水) 07:57:56.62 ID:liMuA23p]

複数の変数について、特定のidだけ表示できる方法
やっと見つけたんで、一応共有を
なんて、みんな知ってるもんなんかな。

browse id syussinti if inlist(id, 1256, 1338)

[0330 ごろ 2016/05/18(水) 11:43:54.97 ID:PzPRpvWP]

書き込みありがとうございます。
> 複数の変数について、特定のidだけ表示できる方法
> やっと見つけたんで、一応共有を なんて、みんな知ってるもんなんかな。
面白そうですが、いまいち問題定義が分かりません。
何か具体例を出して頂けますか？
解決案は少し自分で考えてみたいので後で投稿していただけると面白いです！！

[0331 ￥ ◆2VB8wsVUoo 2016/05/18(水) 14:20:20.89 ID:aEUwscPn]

￥

[0332 ￥ ◆2VB8wsVUoo 2016/05/18(水) 14:20:39.86 ID:aEUwscPn]

￥

[0333 ￥ ◆2VB8wsVUoo 2016/05/18(水) 14:20:58.24 ID:aEUwscPn]

￥

[0334 ￥ ◆2VB8wsVUoo 2016/05/18(水) 14:21:16.85 ID:aEUwscPn]

￥

[0335 ￥ ◆2VB8wsVUoo 2016/05/18(水) 14:21:36.39 ID:aEUwscPn]

￥

[0336 ￥ ◆2VB8wsVUoo 2016/05/18(水) 14:21:53.79 ID:aEUwscPn]

￥

[0337 ￥ ◆2VB8wsVUoo 2016/05/18(水) 14:22:10.09 ID:aEUwscPn]

￥

[0338 ￥ ◆2VB8wsVUoo 2016/05/18(水) 14:22:27.84 ID:aEUwscPn]

￥

[0339 ￥ ◆2VB8wsVUoo 2016/05/18(水) 14:22:48.98 ID:aEUwscPn]

￥

[0340 ￥ ◆2VB8wsVUoo 2016/05/18(水) 14:23:10.30 ID:aEUwscPn]

￥

[0341 ￥ ◆2VB8wsVUoo 2016/05/19(木) 05:49:31.47 ID:FZLtRdgq]

￥

[0342 ￥ ◆2VB8wsVUoo 2016/05/19(木) 12:24:45.04 ID:FZLtRdgq]

￥

[0343 ￥ ◆2VB8wsVUoo 2016/05/19(木) 14:04:42.78 ID:FZLtRdgq]

￥

[0344 ￥ ◆2VB8wsVUoo 2016/05/20(金) 06:49:56.91 ID:9ZaMs54t]

￥

[0345 ￥ ◆2VB8wsVUoo 2016/05/20(金) 07:10:10.33 ID:9ZaMs54t]

￥

[0346 ￥ ◆2VB8wsVUoo 2016/05/20(金) 08:25:37.29 ID:9ZaMs54t]

￥

[0347 ￥ ◆2VB8wsVUoo 2016/05/20(金) 23:08:29.79 ID:9ZaMs54t]

￥

[0348 １３２人目の素数さん 2016/05/22(日) 00:38:50.70 ID:eiV1vWEJ]

マ イ ン ド コ ン ト ロ ー ル の手法

・沢山の人が、偏った意見を一貫して支持する
　偏った意見でも、集団の中でその意見が信じられていれば、自分の考え方は間違っているのか、等と思わせる手法

・不利な質問をさせなくしたり、不利な質問には答えない、スルーする
　誰にも質問や反論をさせないことにより、誰もが皆、疑いなど無いんだと信じ込ませる手法

偏った思想や考え方に染まっていたり、常識が通じない人間は、頭が悪いフリをしているカルト工作員の可能性が高い

靖 国 参 拝、皇 族、国 旗 国 歌、神 社 神 道を嫌う カ ル ト

10人に一人は カ ル ト か 外 国 人

「ガ ス ラ イ テ ィ ン グ」 で 検 索 を ！

[0349 １３２人目の素数さん 2016/06/26(日) 09:03:20.73 ID:vflKx0ln]

統計学は統計学 数学ではない

[0350 １３２人目の素数さん 2016/07/14(木) 07:10:34.15 ID:NLJZUe4v]

Mac良いよ〜

[0351 ￥ ◆2VB8wsVUoo 2016/07/14(木) 07:24:00.46 ID:W/q5qfz4]

日本人の躾けは『大人の都合』、その目的は威厳に屈服させる為：
ある父親：クマが出没する山林に息子を放置、しかも嘘を吐いて保身。
別の父親：勉強の邪魔をして進路を妨害し、学歴を砕く。出世を強要。
ソレでも「親の行為は子供の為」という傲慢な常識を振り回す世間、しかも
「親を尊敬して大切に扱え」という無根拠な思想を押し付ける儒教文化。

お父さん、お母さんを大切にしましょう！！！ソレが世間体というモノ！

ケケケ￥

政治家も、お教授も、権力を振り回すのが大好きな低能人種：
ある男：ボクは都民の為に湯河原で休んでるんだ。知事が信じられんのかっ！
別の男：オレは哲也の為に指導してやってるんだ。父親が信じられんのかっ！
上から目線で強弁すれば、自分の言い分は何でも通る国があるらしい…

ああ、素晴らしき日本文化よ。キミ達も国会議員を見習い給え。何せ多数決で選
ばれた『皆の代表』なので。だからある男も別の男もエラいんだよォ〜〜〜んw

コココ￥

終わり良ければ全てヨシ。途中経過はどうでもヨシ。
大学：学生の知能なんてどうでもヨシ。カネが儲かる教室を巧みに運営シロ。
狸研：研究の詳細なんてどうでもヨシ。世間が驚く大論文を外国に発表シロ。
芳雄：学問の中身なんてどうでもヨシ。安易に教授になれる分野を専攻シロ。

学問なんて所詮は出世の道具。周囲に秀才っぽく見せ掛けられたらソレでヨシ。
社会的に高い地位、そして豪華で贅沢な暮らし。世間が羨む大学教授のポスト。
ソレさえ手に入れば学問そのものなんて洋梨よォ〜〜〜ん。

よよよ、よ〜〜〜しお。そやしノ〜ベル賞が欲しいよォ〜〜〜んんんwww

シシシ￥

[0352 １３２人目の素数さん 2016/08/03(水) 22:34:25.89 ID:RetPXQAc]

ここ数年分を、ど初心者にもわかりやすく書いてくれら、私はあなたを神と呼ぶことにするであろう

[0353 １３２人目の素数さん 2016/08/24(水) 13:14:53.63 ID:w89vX+Fy]

コレ使ってcommerceの売り上げをアプしたいんや
手取り足とり頭とり教えてや

[0354 １３２人目の素数さん 2016/08/24(水) 13:23:49.03 ID:82xMB2d4]

>>353
普通にダミー使っての重回帰で良くね？あまり考えなくてもさ…
流入元とかあるなら色々面倒だけどろ

[0355 ごろ 2016/09/07(水) 16:38:13.01 ID:BSJMfCpj]

久しぶりにメモ
Stataはクオーテーションマーク関連に煩い。
1.

[0356 ごろ 2016/09/07(水) 16:38:48.40 ID:BSJMfCpj]

書き込み失敗した。
久しぶりにメモ
Stataはクオーテーションマーク関連に煩い。
1. "
2. '
3. “
4. ”
5. ‘
6. ’
7. `
上記全て違う扱い1,2,7はStataで使える。
7はマクロの左側で使うやつ　`macro'
全角・半角・日本語・英語を気をつけていてもエディタによっては間違える
ちなみに文字コードは下記の通り
１．print(chr(34))
２．print(chr(39))
３．print(chr(8220))
４．print(chr(8221))
５．print(chr(8216))
６．print(chr(8217))
７．print(chr(96))

[0357 ごろ 2016/09/07(水) 16:52:47.89 ID:BSJMfCpj]

もうひとつ
Stata Programmingでどつぼにはまったのでメモ
プログラムにオプションを付けるときsyntaxのみではカンマ(,)を上手く取り除いてくれないのtokenizeを使う。実は、
>>202
に概出
解決方法
program testprogram
syntax namelist [,Option]
tokenize "`0'" ,parse(" ,")
...
end
この ,parse(" ,")がキモ ,parse(",")は駄目
””の中に半角スペースと,が入っていてこれで分けるというコマンド
この後、`1' , `2'で取り出し可能。
ちなみに[,Option]は省略可能で,o や　,opt や ,optionと記載OK
これ以外、例えば,otと書くとエラー
プログラム内では記載が無ければ`option'というマクロに""が
記載があれば`option'というマクロに"option"が入る。

も一つ
capture hogehogeと記載すると
_rcというシステム変数に入力される
エラーが無ければ０が、エラーがあればそのエラー番号が入力される
だからこれを利用して
capture confirm variable 変数名
if _rc {
変数名が無い場合の処理
}
else {
変数名がある場合の処理
}
裏返しで
capture confirm variable 変数名
if !_rc {
変数名がある場合の処理
}
else {
変数名が無い場合の処理
}
と記載できる。

[0358 ￥ ◆2VB8wsVUoo 2016/09/09(金) 04:40:06.26 ID:c0Am8D5Q]

￥

[0359 ￥ ◆2VB8wsVUoo 2016/09/09(金) 04:40:25.28 ID:c0Am8D5Q]

￥

[0360 ￥ ◆2VB8wsVUoo 2016/09/09(金) 04:40:42.91 ID:c0Am8D5Q]

￥

[0361 ￥ ◆2VB8wsVUoo 2016/09/09(金) 04:41:03.16 ID:c0Am8D5Q]

￥

[0362 ￥ ◆2VB8wsVUoo 2016/09/09(金) 04:41:21.51 ID:c0Am8D5Q]

￥

[0363 ￥ ◆2VB8wsVUoo 2016/09/09(金) 04:41:39.60 ID:c0Am8D5Q]

￥

[0364 ￥ ◆2VB8wsVUoo 2016/09/09(金) 04:41:56.81 ID:c0Am8D5Q]

￥

[0365 ￥ ◆2VB8wsVUoo 2016/09/09(金) 04:42:30.02 ID:c0Am8D5Q]

￥

[0366 ￥ ◆2VB8wsVUoo 2016/09/09(金) 04:43:01.26 ID:c0Am8D5Q]

￥

[0367 ￥ ◆2VB8wsVUoo 2016/09/09(金) 04:43:23.17 ID:c0Am8D5Q]

￥

[0368 ￥ ◆2VB8wsVUoo 2016/09/10(土) 17:07:25.33 ID:9+jIN1kB]

￥

[0369 ￥ ◆2VB8wsVUoo 2016/09/10(土) 17:46:36.63 ID:9+jIN1kB]

￥

[0370 ￥ ◆2VB8wsVUoo 2016/09/11(日) 03:23:25.86 ID:M5fLCYO8]

￥

[0371 ￥ ◆2VB8wsVUoo 2016/09/12(月) 02:47:18.42 ID:wdbNdCQa]

￥

[0372 ￥ ◆2VB8wsVUoo 2016/09/16(金) 01:24:42.78 ID:Lx0ol+eh]

￥

[0373 ￥ ◆2VB8wsVUoo 2016/09/16(金) 05:13:15.28 ID:Lx0ol+eh]

￥

[0374 １３２人目の素数さん 2016/09/28(水) 18:37:14.49 ID:fWNjEo2B]

なんかもうすこし活発な日本のStataコミュニティって無いのかな。

[0375 ￥ ◆2VB8wsVUoo 2016/09/28(水) 18:49:49.75 ID:WIEHiU/H]

￥

[0376 ￥ ◆2VB8wsVUoo 2016/09/28(水) 20:26:36.09 ID:WIEHiU/H]

￥

[0377 １３２人目の素数さん 2016/09/28(水) 20:53:14.60 ID:KZvQj/bA]

>>374
stata使える時点で英語読めないとイマイチなわけだし、日本語で聞くよりも前に英語で聞くか、ラインストーンに電話してしまうんじゃね

[0378 １３２人目の素数さん 2016/10/03(月) 14:09:32.04 ID:mEM7hnKf]

stataのdo fileとか出力結果ってsdkとして読み出したりパイソンあたりと連携できたりするのかな？

[0379 １３２人目の素数さん 2016/10/21(金) 19:17:36.41 ID:wtCNBYoA]

c()やp()の使い方と意味がわからないんだけど、誰か教えてください。

[0380 のの 2017/01/09(月) 00:58:12.07 ID:ei+5V4I3]

今年もよろしくお願い致します。

[0381 １３２人目の素数さん 2017/01/09(月) 08:34:51.36 ID:nI8BECga]

こちらこそ、よろしくお願いします

[0382 ごろ 2017/03/08(水) 14:12:32.89 ID:b/g5dyBn]

今日のメモ

知らなかった。
for num 1/10 \ num 2/11 \ num 3/12: gen aX=X \ gen aY=Y \gen aZ=Z

for numの使い方

[0383 ごろ 2017/03/09(木) 15:09:12.49 ID:IiwnUzov]

これ忘れていた・・・
プログラムで作成すれば余計なコマンドが表示されないんだった！

Doファイルの中でDoファイルを作成して実行する方法
capture program drop makedofile
program makedofile
syntax anything
qui{ /*これがポイント。これでログファクィルのヘッダーを出力しない。*/
capture log close
log using newdofile.do , text replace
}
di "clear"
di "input str10 newvarname data1 data2"
foreach i of local anything {
local x = _N
foreach j of numlist 1(1)`x' {
di originalvarname " " /*ここにdata1, data2用の命令*/
}
}
di "end"
qui {
log close
}
end

makedofile /*originaldataを呼び出す変数名など*/

[0384 ごろ 2017/03/29(水) 17:00:27.17 ID:mBNHDgR4]

tabstat 今まであまり使っていなかったけど便利。
tabstat 複数変数名, stat(n sum mean)で複数変数名の観察数、0 or 1変数なら１の合計, 平均値が出る。
ただ残念なのことに、変数名の表示方法を指定できない（ハズだ）から、変数名を省略しないで表示したい場合には、
adoファイルまで戻って、colwdithを書き換えないといけない（ハズ）。
誰かもっと良い方法があれば教えて欲しい。
Statalistでもsummarizeコマンドで同様のDiscussionがあった。

[0385 ごろ 2017/03/30(木) 15:44:53.75 ID:IwKOh0nf]

私のためのメモ
capture ds
foreach i in `r(varlist)'{
コマンド
}
これ便利。
dsでr(varlist)というローカル変数に変数名リストが全てリストアップされる。
全ての変数に対してあるコマンドを繰り返し実行することが可能。

dsコマンドはヘルプを見ると分かるが、色々と使い勝手が良さそう。

[0386 ごろ 2017/04/18(火) 19:22:49.08 ID:aaD2jafi]

変数名が32文字までしかダメなためプログラム内でエラーを起こさない工夫。
多分、使い捨ての変数を使う方法がこの2ch内の記事にもあったが
以下の方法でも出来たので私のためのメモ。

clear
capture program drop test
program test
syntax varlist
foreach i of local varlist{
local j `=strtoname(abbrev("`i'",29))'
gen new`j' = 1
}
end

input thisvarname_has_32characterslong short_name
1 1
2 2
3 3
end

test thisvarname_has_32characterslong short_name
des ,f

[0387 ￥ ◆2VB8wsVUoo 2017/05/08(月) 09:23:35.71 ID:OR+quqWp]

￥

[0388 ￥ ◆2VB8wsVUoo 2017/05/08(月) 09:23:57.25 ID:OR+quqWp]

￥

[0389 ￥ ◆2VB8wsVUoo 2017/05/08(月) 09:24:18.45 ID:OR+quqWp]

￥

[0390 ￥ ◆2VB8wsVUoo 2017/05/08(月) 09:24:41.23 ID:OR+quqWp]

￥

[0391 ￥ ◆2VB8wsVUoo 2017/05/08(月) 09:25:05.22 ID:OR+quqWp]

￥

[0392 ￥ ◆2VB8wsVUoo 2017/05/08(月) 09:25:27.65 ID:OR+quqWp]

￥

[0393 ￥ ◆2VB8wsVUoo 2017/05/08(月) 09:25:51.39 ID:OR+quqWp]

￥

[0394 ￥ ◆2VB8wsVUoo 2017/05/08(月) 09:26:15.51 ID:OR+quqWp]

￥

[0395 ￥ ◆2VB8wsVUoo 2017/05/08(月) 09:26:38.15 ID:OR+quqWp]

￥

[0396 ￥ ◆2VB8wsVUoo 2017/05/08(月) 09:26:59.71 ID:OR+quqWp]

￥

[0397 １３２人目の素数さん 2017/06/08(木) 23:43:11.92 ID:B17WMjKl]

Stata15リリース記念あげ

[0398 ￥ ◆2VB8wsVUoo 2017/06/09(金) 09:53:17.88 ID:XVb2Gvc9]

￥

[0399 ごろ 2017/07/03(月) 00:28:57.24 ID:SG4Prbcu]

>>386
一時変数（使い捨て変数）の作成方法
>>199
に記載があった。
再掲
プログラム内のみで使う変数には
tempvar 変数名
を使う。
変数名には__000001のような名前が付けられるが、
プログラム内ではユーザーは上記で指定した変数名が使える。

[0400 ごろ 2017/07/06(木) 07:14:36.32 ID:EZM9tbus]

>>300
多重比較の問題ですね。
Stataのコマンドというよりは、多重比較について考えると良いと思います。
下記サイトなど参考になると思います。
http://www.ibaraki-kodomo.com/toukei/posthoc.html

tukeyhsdはユーザが作ったコマンドです。
使い方としては、ネットがつながった環境で、Stataから下記コマンドを打つとリスト出てくるので
net from http://stats.idre.ucla.edu/stat/stata/ado/analysis
リストからtukeyhsdを選択し次の画面でclick install
同様にqsturngもインストールします（tukeyhsdを使うために必要らしい）。
インストールが完了すればヘルプファイルから分かると思いますが、
例えば下記を試してみて下さい。
webuse systolic,clear
anova systolic drug
tukeyhsd drug
です。
ああ、2年以上も前の質問にまた回答してしまった。
進歩がないなぁ。

[0401 １３２人目の素数さん 2017/08/06(日) 17:52:05.11 ID:oDKJI1vJ]

耳栓をしたら世界が変わってワロタ

[0402 ￥ ◆2VB8wsVUoo 2017/08/06(日) 18:08:46.27 ID:+CYdGQny]

☆☆☆馬鹿板は数学徒の脳を腐らせる悪い板であり、そやし廃止してナシにすべき。☆☆☆

￥

[0403 １３２人目の素数さん 2017/08/06(日) 18:08:57.47 ID:oDKJI1vJ]

耳栓をしたら世界が変わってワロタ

[0404 ￥ ◆2VB8wsVUoo 2017/08/06(日) 18:18:34.65 ID:+CYdGQny]

￥

[0405 １３２人目の素数さん 2017/08/07(月) 08:59:12.21 ID:0YzkEl/p]

耳栓をしたら世界が変わってワロタ

[0406 ￥ ◆2VB8wsVUoo 2017/08/07(月) 09:26:30.63 ID:/rspiZFz]

￥

[0407 ￥ ◆2VB8wsVUoo 2017/08/07(月) 10:24:50.38 ID:/rspiZFz]

￥

[0408 ￥ ◆2VB8wsVUoo 2017/08/07(月) 10:25:09.17 ID:/rspiZFz]

￥

[0409 ￥ ◆2VB8wsVUoo 2017/08/07(月) 10:25:28.46 ID:/rspiZFz]

￥

[0410 ￥ ◆2VB8wsVUoo 2017/08/07(月) 10:25:46.37 ID:/rspiZFz]

￥

[0411 ￥ ◆2VB8wsVUoo 2017/08/07(月) 11:18:05.71 ID:/rspiZFz]

￥

[0412 ￥ ◆2VB8wsVUoo 2017/08/07(月) 11:18:23.08 ID:/rspiZFz]

￥

[0413 ￥ ◆2VB8wsVUoo 2017/08/07(月) 11:18:46.65 ID:/rspiZFz]

￥

[0414 ￥ ◆2VB8wsVUoo 2017/08/07(月) 11:19:03.83 ID:/rspiZFz]

￥

[0415 ￥ ◆2VB8wsVUoo 2017/08/07(月) 11:19:20.62 ID:/rspiZFz]

￥

[0416 １３２人目の素数さん 2017/09/02(土) 12:50:17.34 ID:3V8qFOPU]

耳栓をしたら世界が変わってワロタ

[0417 ￥ ◆2VB8wsVUoo 2017/09/02(土) 13:06:56.61 ID:z17/uuYO]

￥

[0418 １３２人目の素数さん 2017/09/05(火) 06:10:57.32 ID:PCfG056b]

耳栓をしたら世界が変わってワロタ

[0419 ￥ ◆2VB8wsVUoo 2017/10/24(火) 11:05:17.27 ID:n4VIAFJn]

￥

[0420 ￥ ◆2VB8wsVUoo 2017/10/24(火) 11:05:35.75 ID:n4VIAFJn]

￥

[0421 ￥ ◆2VB8wsVUoo 2017/10/24(火) 11:05:54.82 ID:n4VIAFJn]

￥

[0422 ￥ ◆2VB8wsVUoo 2017/10/24(火) 11:06:12.72 ID:n4VIAFJn]

￥

[0423 ￥ ◆2VB8wsVUoo 2017/10/24(火) 11:06:30.48 ID:n4VIAFJn]

￥

[0424 ￥ ◆2VB8wsVUoo 2017/10/24(火) 11:06:47.83 ID:n4VIAFJn]

￥

[0425 ￥ ◆2VB8wsVUoo 2017/10/24(火) 11:07:11.67 ID:n4VIAFJn]

￥

[0426 ￥ ◆2VB8wsVUoo 2017/10/24(火) 11:07:28.30 ID:n4VIAFJn]

￥

[0427 ￥ ◆2VB8wsVUoo 2017/10/24(火) 11:07:46.16 ID:n4VIAFJn]

￥

[0428 ￥ ◆2VB8wsVUoo 2017/10/24(火) 11:08:04.01 ID:n4VIAFJn]

￥

[0429 １３２人目の素数さん 2017/11/25(土) 00:57:18.24 ID:vOvLAPzg]

stataを使おうとしている経済系院生です
RCL: Stata module for estimation and simulation of random coefficient logit models
を使って需要関数の推計をして合併シミュレーションを行おうとしているのですが、都内で誰か教えてくれる人はいますか。
謝礼はします。
修士レベルのミクロ経済・産業組織論の話なのですが…

[0430 ￥ ◆2VB8wsVUoo 2017/11/25(土) 10:59:49.23 ID:kb2pRtxG]

￥

[0431 ￥ ◆2VB8wsVUoo 2017/11/25(土) 11:00:03.91 ID:kb2pRtxG]

￥

[0432 ￥ ◆2VB8wsVUoo 2017/11/25(土) 11:00:18.56 ID:kb2pRtxG]

￥

[0433 ￥ ◆2VB8wsVUoo 2017/11/25(土) 11:00:50.06 ID:kb2pRtxG]

￥

[0434 ￥ ◆2VB8wsVUoo 2017/11/25(土) 11:01:06.67 ID:kb2pRtxG]

￥

[0435 ￥ ◆2VB8wsVUoo 2017/11/25(土) 11:01:22.02 ID:kb2pRtxG]

￥

[0436 ￥ ◆2VB8wsVUoo 2017/11/25(土) 11:01:39.33 ID:kb2pRtxG]

￥

[0437 ￥ ◆2VB8wsVUoo 2017/11/25(土) 11:01:55.74 ID:kb2pRtxG]

￥

[0438 ￥ ◆2VB8wsVUoo 2017/11/25(土) 11:02:12.77 ID:kb2pRtxG]

￥

[0439 ￥ ◆2VB8wsVUoo 2017/11/25(土) 11:02:33.93 ID:kb2pRtxG]

￥

[0440 ごろ 2018/03/13(火) 11:47:35.06 ID:yzMoRI6M]

いつの間にか5chになっていて驚いた

[0441 １３２人目の素数さん 2018/03/14(水) 22:41:49.65 ID:AlbLRV+J]

>>440
統計の話をしたいね

[0442 １３２人目の素数さん 2018/03/17(土) 15:23:21.63 ID:/k6W/g7+]

Rとstata比較するとどうなんですか？

[0443 ￥ ◆2VB8wsVUoo 2018/04/06(金) 03:09:40.54 ID:I+Mybrk/]

￥

[0444 ￥ ◆2VB8wsVUoo 2018/04/06(金) 03:09:57.81 ID:I+Mybrk/]

￥

[0445 ￥ ◆2VB8wsVUoo 2018/04/06(金) 03:10:16.66 ID:I+Mybrk/]

￥

[0446 ￥ ◆2VB8wsVUoo 2018/04/06(金) 03:10:33.76 ID:I+Mybrk/]

￥

[0447 ￥ ◆2VB8wsVUoo 2018/04/06(金) 03:10:51.54 ID:I+Mybrk/]

￥

[0448 ￥ ◆2VB8wsVUoo 2018/04/06(金) 03:11:04.17 ID:I+Mybrk/]

￥

[0449 ￥ ◆2VB8wsVUoo 2018/04/06(金) 03:11:23.58 ID:I+Mybrk/]

￥

[0450 ￥ ◆2VB8wsVUoo 2018/04/06(金) 03:11:41.31 ID:I+Mybrk/]

￥

[0451 ￥ ◆2VB8wsVUoo 2018/04/06(金) 03:11:59.07 ID:I+Mybrk/]

￥

[0452 ￥ ◆2VB8wsVUoo 2018/04/06(金) 03:12:19.77 ID:I+Mybrk/]

￥

[0453 １３２人目の素数さん 2018/04/07(土) 19:19:24.03 ID:NNMRscPu]

耳栓をしたら世界が変わってワロタ

[0454 １３２人目の素数さん 2019/03/29(金) 21:05:14.57 ID:g6RZxVSs]

「統計」は「疑似科学」な

[0455 １３２人目の素数さん 2019/06/10(月) 12:06:38.87 ID:UAW6/NkJ]

STATA　なんて使ってる人いまいるの？

[0456 １３２人目の素数さん 2019/06/11(火) 18:20:01.95 ID:2LvhxiIt]

>>455
上級者はいるやろ

[0457 １３２人目の素数さん 2019/06/29(土) 16:49:54.88 ID:DHiuKlHq]

統計ソフトSTATAの部屋 Ver.2
ふうL@Fu_L12345654321
学コン1傑いただきました！
とても嬉しいです！

https://pbs.twimg.com/media/D-IuUuqVUAALnAB.jpg

https://twitter.com/Fu_L12345654321/status/1144528199654633477
https://twitter.com/5chan_nel (5ch newer account)

[0458 １３２人目の素数さん 2019/07/03(水) 19:39:12.04 ID:dqLWAG/2]

3915
ふうL@Fu_L12345654321
学コン1傑いただきました！
とても嬉しいです！

https://pbs.twimg.com/media/D-IuUuqVUAALnAB.jpg

https://twitter.com/Fu_L12345654321/status/1144528199654633477
https://twitter.com/5chan_nel (5ch newer account)

[0459 高添沼田エロ老義父の告発（葛飾区青戸６−２６−６） 2021/02/05(金) 07:43:58.25 ID:WSArCHeX]

491色川高志「井口千明の金属バット集団殴打撲殺を熱望します」2018/10/18(木) 18:33:15.90ID:78662J73
龍神連合五代目総長・井口千明（葛飾区青戸6−２３−１９）の挑発
井口千明「糞関東連合文句があったらいつでも俺様を金属バットで殴り殺しに来やがれっ!! 糞関東連合の見立・石元・伊藤リオンの糞野郎どもは
龍神連合四代目総長の俺様がぶちのめしてやるぜっ!! 賞金をやるからいつでもかかって来いっ!!　糞バエ関東連合どもっ!! 待ってるぜっ!!」(挑戦状)

492盗聴盗撮犯罪者色川高志（青戸６−２３−２１ハイツニュー青戸１０３2021/02/03(水) 13:53:22.55ID:QtP78E4Z
●青戸六丁目被害者住民一同「盗聴盗撮犯罪者の高添沼田ハゲエロ老義父の逮捕を要請します」
長木親父＆長木よしあき（盗聴盗撮犯罪者の高添沼田ハゲエロ老義父を逮捕に追い込む会＆被害者の会会長）住所＝東京都葛飾区青戸6−23−20
●盗聴盗撮つきまとい嫌がらせ犯罪者／アナル挿入食糞愛好家で息子の嫁で自慰行為をしている高添沼田ハゲエロ老義父
高添沼田ハゲエロ老義父の住所＝東京都葛飾区青戸6−26−６
【通報先】亀有警察署＝東京都葛飾区新宿4ー22ー19　�рO３ー３６０７ー０１１０

盗聴盗撮つきまとい嫌がらせ犯罪者／アナル挿入食糞愛好家で息子の嫁で自慰行為をしている高添沼田ハゲエロ老義父の盗聴盗撮つきまとい嫌がらせ犯罪者／愛人変態メス豚家畜清水婆婆（青戸6−23−19）の
五十路後半強制脱糞
http://img.erogazou-pinkline.com/img/2169/scatology_anal_injection-2169-027.jpg




アナル挿入食糞愛好家で息子の嫁で自慰行為をしている高添沼田ハゲエロ老義父によりバスタブで清水婆婆の巨尻の肛門にシャワーのキャップをはずしてずっぽり挿入。
そして、大量浣腸。 勢い良く噴出！腸内洗浄状態です。
http://101.dtiblog.com/b/bodytk9690/file/kan01.jpg




浣腸器と異なりどくどくと直腸内に注入され清水婆婆は激しくあえぎます