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

[表示 : 全て最新50 1-99 101- 201- 2chのread.cgiへ]
Update time : 04/09 19:49 / Filesize : 116 KB / Number-of Response : 243
[このスレッドの書き込みを削除する]
[＋板最近立ったスレ＆熱いスレ一覧 : ＋板最近立ったスレ／記者別一覧] [類似スレッド一覧]


↑キャッシュ検索、類似スレ動作を修正しました、ご迷惑をお掛けしました

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

1 名前：ごろ mailto:sage [2010/09/10(金) 10:54:10 ]: 統計ソフト STATAを使って統計解析をするスレです。
日本語の良いソースが少ないので盛り上げましょう。

前スレが順調に!?88まで行ったんですが、落ちてしまったので再開です。
141 名前：ごろ mailto:sage [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
142 名前：１３２人目の素数さん mailto:sage [2011/03/25(金) 20:00:53.28 ]: ＠保守

ごろさん大丈夫かぁー？？
143 名前：ごろ mailto:sage [2011/03/27(日) 03:08:58.46 ]: >>142
ありがとうございます。私と私の周りの人は皆無事です。
しかし、あまりの事態に、のんきにStataしている気になれませんでした。
地震・津波そのものの被害も大変なものでしたが、原発事故の影響が恐ろしいです。
犠牲になられた方のご冥福をお祈り申し上げます。
144 名前：１３２人目の素数さん [2011/04/10(日) 23:28:41.47 ]: ちなみに、stataのオンライン・レクチャーはDVDでみれたり、あるいは、過去のものをみれたりするのでしょうか。社会人の私としては、決まった時期にみるというより、休みの時に一気に見たいのですが。
145 名前：ごろ mailto:sage [2011/04/12(火) 05:52:11.98 ]: >>144
オンラインレクチャーというよりは、オンラインで置いてあるテキストを読みながら
自分で勉強するという感じですね。
そういう意味では本を購入するのと変わらないので１万以上するのは高い気もしますが、
質問もできるというインタラクティブ性が良いのでしょうね。

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

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

Advanceコースは10月なのでまだ先ですが、一応今のところは受けるつもりです。
他にも興味深いコースはありますが、２万円以上するので、さすがにちょっと躊躇しています。
146 名前：１３２人目の素数さん [2011/04/13(水) 00:04:53.57 ]: >>145

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

ちなみに、どれくらいの期間視聴が可能なのでしょうか。
147 名前：１３２人目の素数さん [2011/04/13(水) 00:32:24.18 ]: >>146

ちなみに、MATAについてはオンライン・セミナーで議論されるでしょうか。ネット上でかかれた.mlibファイルが開くことができず、stata9で開けず苦労しています。このあたりの勉強を講義ベースで聞けたらと思うのですが。
148 名前：ごろ mailto:sage [2011/04/13(水) 02:09:41.64 ]: >>147
視聴期間は基本的に講義開催期間のみですが、先にも書いたとおり、最終的には
サイト自体がダウンロードできるようになるので、講義終了後に繰り返し復習することは可能です。

MATAについてはよく分かりませんが、少なくとも私が受講したコースには有りませんでした。
しかし、講義から離れた内容の質問でも答えてもらえる可能性は高いと思います。
149 名前：ごろ mailto:sage [2011/04/24(日) 23:29:58.92 ]: ノートパソコンのOSをクリーンインストールし直した。
前は、Stataの10と11を両方使っていたけど、今回は11のみでいくことにした。
150 名前：ごろ mailto:sage [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に変数が沢山つくと挙動がおかしくなる。
なんでだろー
151 名前：１３２人目の素数さん mailto:sage [2011/05/05(木) 22:21:25.50 ]: MAC版で日本語を表示させると、詰まって表示されるのですが、
何か対策はありますか。

教科書体を利用すると、多少ましになりますが、他に手段は
ないのでしょうか。
152 名前：ごろ mailto:sage [2011/05/06(金) 02:46:33.01 ]: >>151
私もMac版 v11.2を使っています。

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

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

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

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

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

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

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

あと、Font関連では、ギリシャ文字や±などの特殊文字もStataで表示させるのは
難しいようですね。
既存のフォントを入れ替えて表示させる方法がどこかにありましたが、
また見つけたらリンクを貼ります。
153 名前：１３２人目の素数さん mailto:sage [2011/05/06(金) 22:39:02.59 ]: >>152
早速のお返事ありがとうございます。
うまくいかないものなのですね。

どうしても、日本語を使いたいときは、
教科書体でがんばってみたいと思います。
154 名前：１３２人目の素数さん [2011/05/11(水) 00:21:35.86 ]: stataのマニュアルは購入しなければならないのでしょうか。
time serisとpanelのマニュアルが取得したいのですが。。
155 名前：ごろ mailto:sage [2011/05/12(木) 11:41:33.04 ]: >>154
はっきり覚えていませんが、紙のマニュアルは無かったと思います。
しかし、オンラインマニュアルや、ソフト付属のマニュアルも
英語ですが、なかなか詳しいです。
プリントしてあるマニュアルは買わないといけないと思います。
156 名前：ごろ mailto:sage [2011/05/12(木) 11:43:03.43 ]: >>155
あ、私も紙のマニュアルがあるときっと便利なんだろうなぁと思いながら、
ついつい、Googleで済ませてしまってここまで来ています。
あと、本当に困ったら、Stataの英語メーリングリストがあるので、そこで聞いていますね。
157 名前：ごろ mailto:sage [2011/05/24(火) 09:54:32.84 ]: 最近、Sleepsortという面白い記事を見つけたので、Stataでも出来ないかと
試し中。
しかし、バックグラウンドで関数を動かす方法が分からないのでつまづいています。
d.hatena.ne.jp/shinsuku/20110520/p1
158 名前：あんでぃはアホ ◆AdkZFxa49I mailto:sage [2011/05/29(日) 23:38:17.38 ]: あんでぃ
159 名前：ごろ mailto:sage [2011/05/30(月) 00:21:16.77 ]: どうやら、Stataには二つの関数を同時にうごかす方法はないようです。
残念！
160 名前：ごろ mailto:sage [2011/06/29(水) 16:49:52.89 ]: pwcorrコマンドの結果をreturn listで見るとmatrixがあるんですが、
このmatrixの扱いが良くわからない。
matrixの要素だけを取り出してdisplayする方法はあるのでしょうか？
161 名前：１３２人目の素数さん mailto:sage [2011/07/03(日) 14:22:31.00 ]: ヴァージョン１２が25日にでますね。
共分散構造分析ができるようになりますね。
アップデートしたいが、先立つものがない・・・

www.stata.com/stata12/
162 名前：１３２人目の素数さん mailto:sage [2011/07/08(金) 10:53:39.06 ]: このソフトいれてみたけどなんでまともにデータファイルが読めないんだか
ふつーにエクセルぐらい読めよバカ
163 名前：１３２人目の素数さん mailto:sage [2011/07/10(日) 08:46:36.61 ]: >>162の率直な意見を取り入れてstata12ではExcelファイルを読めるようになりますた
164 名前：１３２人目の素数さん mailto:sage [2011/07/11(月) 11:23:03.75 ]: まったくとろい会社だな。
165 名前：１３２人目の素数さん mailto:sage [2011/07/15(金) 11:24:37.13 ]: STATAって貧乏人のSASですか?
166 名前：ごろ mailto:sage [2011/07/17(日) 09:34:11.52 ]: >>162
まぁ、Excel側でCSVファイルで出力すれば言い訳だから、そんなに困んないけどね。
167 名前：１３２人目の素数さん mailto:sage [2011/07/19(火) 15:04:48.01 ]: 急いでる時にわざわざExcelでセーブしなおすのは結構面倒だぜ。なにせ
俺様売れっ子だから
168 名前：１３２人目の素数さん mailto:sage [2011/07/19(火) 23:13:24.38 ]: >>165
プログラムの書けない奴のRだと思う
169 名前：１３２人目の素数さん mailto:sage [2011/07/19(火) 23:49:35.33 ]: >>168
Rって書ける、書けない、言うほどのものか？
170 名前：１３２人目の素数さん mailto:sage [2011/07/31(日) 08:11:01.67 ]: >>169は優秀なスレストッパー(ﾆﾔﾘ
171 名前：１３２人目の素数さん mailto:sage [2011/08/02(火) 12:16:30.64 ]: >>161
ようやくSEMができるようになったか。
Stata 11を入れたばっかりだけど、買おうかなー。
172 名前： ◆YeNleXj6MwwY mailto:sage [2011/08/15(月) 01:19:56.07 ]: 　|>､＿_ ／|ﾌ
　〉｡,　o　く　/＼
　 b ～　O ヽ〉 /
　（>　＜ノ ∨
　　（､＿（、ノ
　　￣　￣
173 名前：ごろ mailto:sage [2011/08/16(火) 23:39:50.80 ]: おー。昨日アクセスできなくて落ちてしまったかと思った。
174 名前：ごろ mailto:sage [2011/09/12(月) 00:37:57.76 ]: 保守
175 名前：ごろ mailto:sage [2011/09/13(火) 02:31:02.76 ]: そろそろまたやりますよ。
176 名前：１３２人目の素数さん mailto:sage [2011/09/14(水) 19:53:32.28 ]: 　|>､＿_ ／|ﾌ
　〉｡,　o　く　/＼
　 b ～　O ヽ〉 /
　（>　＜ノ ∨
　　（､＿（、ノ
　　￣　￣
177 名前：１３２人目の素数さん mailto:sage [2011/09/25(日) 09:19:35.91 ]: 試用はできないんですか？
178 名前：ごろ mailto:sage [2011/09/25(日) 14:16:20.84 ]: >>177
試用できないのかな？オンラインではないし、聞いたことないなぁ。大学とか知り合いで使っている人がいれば一番いいんですけどね。
179 名前：１３２人目の素数さん mailto:sage [2011/09/25(日) 22:22:05.78 ]: 日本の代理店も本家も探してみましたが、ありませんでした。
SEMはどんな感じかなと思いまして。
AMOSはお絵描きでできると言えば簡単そうなんですが、誤差変数を忘れずに入れろとか、
制約を適切に入れないといかん、とか実は結構面倒なので。
毎年のお布施も結構かかりますし。
Stataはぱっと見た感じでは、核となる部分さえ入力すればいいように見えたものですから。
180 名前：猫は熾猫 ◆s9Thv0XHxQ mailto:age [2011/10/01(土) 11:42:22.56 ]: 猫
181 名前：１３２人目の素数さん mailto:sage [2011/10/01(土) 18:21:49.01 ]: >>179
12買ったけど、その機能はまだ使ってないっす。
いじくる暇を見つけたらレポートします。
182 名前：ごろ mailto:sage [2011/10/04(火) 23:06:51.89 ]: >>181
おお、１２買ったんですね。いいなぁ。
レポートよろしくお願いいたします。
183 名前：ごろ mailto:sage [2011/10/07(金) 12:53:48.80 ]: NetCourse 152始まりました。
また日記形式で綴っていきたいと思います。
184 名前：177=179 mailto:sage [2011/10/10(月) 20:48:18.79 ]: Stata/ICという一番下のクラスのを買いました。
とりあえずauto.dtaというファイルで、mpgとweightとforeignで回帰分析してみました。
次に、SEMのお絵描きで、mpg<-weightとmpt<-foreignの2本の矢印のみを描いてEstimateしてみると、
同じ結果が得られました。誤差変数の○を書き足したり、制約を入れる必要もない、という手軽さ。

どの程度使えるものか、少し試してみます。
185 名前：１３２人目の素数さん mailto:sage [2011/10/11(火) 00:54:36.90 ]: 新参初心者ですが、STATAっていくらから購入可能なんですか？
186 名前：１３２人目の素数さん mailto:sage [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を使いましょう。
最近は入門書も充実してるみたいですしね。
187 名前：１３２人目の素数さん mailto:sage [2011/10/12(水) 00:51:08.07 ]: >>186
ありがとうございます。
188 名前：１３２人目の素数さん mailto:sage [2011/10/13(木) 12:02:44.59 ]: >>184
どーも、なんだか、先を越されてしまいましたね。
買ったはいいが、まだまだまともに使えてないもんで。。。

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

まだ使ってないのですが、何か分かったらここに書きます。
189 名前：184 mailto:sage [2011/10/13(木) 19:02:19.15 ]: >>188
いやいや、色々教えてください。Stataは初めてです。

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

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

ただ、探し方が悪いのか、適合度の指標について、GFIやAGFIという表記がありません。
マニュアルを読まないといかんようです。
190 名前：184 mailto:sage [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ビルダーというお絵描きの方に反映できればいいのですが。
ちょっと難しいか。
191 名前：184 mailto:sage [2011/11/05(土) 11:06:50.82 ]: 遅々として進みませんが、気付いたことを。
regとかlogitとかの基本的な分析コマンドも
adoファイルとして提供されてるんですね。
読んでもちっともわかりませんが。
192 名前：１３２人目の素数さん mailto:sage [2011/11/08(火) 11:03:18.08 ]: >>191

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

まぁ、分析する上ではあまり気にされなくてもよいと思います。
193 名前：ごろ mailto:sage [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"
が代入されるということ
194 名前：ごろ mailto:sage [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
195 名前：ごろ mailto:sage [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が代入される。それ以外は""が代入される。
ここはややこしいところ。
196 名前：ごろ mailto:sage [2011/11/13(日) 00:26:51.92 ]: 3/6続き
オプションは更にオプション名( )でオプションの性質を指定できる。
help syntaxの中に詳細な説明があるが、レクチャーでは次の８つの解説
integer 整数を指定　初期値設定可
real 実数を指定　初期値設定可
string 文字列を指定　
varlist １から複数の変数名を指定。ない変数名を入力するとエラー
numlist 数列を指定
passthru
namelist を指定
cilevel を指定
197 名前：ごろ mailto:sage [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コマンドは、次に来る変数が数字かどうかを判定するコマンド)
198 名前：ごろ mailto:sage [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終わった。
199 名前：ごろ mailto:sage [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 論理式
論理式が真なら何も返さず、偽ならエラーコードを返す。
200 名前：ごろ mailto:sage [2011/11/13(日) 11:30:44.15 ]: 5. r()について
プログラム/関数を実行したときに返す値はr(N)などに格納される。
このリストを調べるにはreturn listを使う
201 名前：ごろ mailto:sage [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は練習問題
202 名前：ごろ mailto:sage [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ではなく、二次式で回帰分析をやる手法の解説
203 名前：ごろ mailto:sage [2011/11/18(金) 14:25:35.30 ]: 5. は長いが、結局二次式による回帰分析の結果を示す方法について書かれている
新しい話は無し。
6. エラー表示のさせ方
エラーコード301はlast estimates not foundという意味らしい
[P]errorに解説がるらしい。
7.if文で条件分けして結果を表示させる問題
8. optionにnoConsを使って結果の一部を表示させないオプションを作る問題。
noConsによってconsマクロに入るものが何かを考えること以外は平易。
Lecture2終了
204 名前：ごろ mailto:sage [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と表示*/
205 名前：１３２人目の素数さん [2011/11/19(土) 13:19:11.95 ]: 魂は幾何学です
誰か(アメリカ)気づいた
ソウルコピー機器
無差別で猥褻、日本は危険知ったかブッタの日本人が失敗作

↓テロ資料です。
206 名前：１３２人目の素数さん [2011/11/19(土) 13:19:49.31 ]: 電波テロ装置の戦争(始)
エンジニアと参加願います公安はサリンオウム信者の子供を40歳まで社会から隔離している
オウム信者が地方で現在も潜伏している
それは新興宗教を配下としている公安の仕事だ
発案で盗聴器を開発したら霊魂が寄って呼ぶ来た
<電波憑依>
スピリチャル全否定なら江原三輪氏、高橋佳子大川隆法氏は、幻聴で強制入院矛盾する日本宗教と精神科
<コードレス盗聴>
2004既に国民20%被害250～700台数中国工作員3～7000万円2005ソウルコピー2010ソウルイン医者アカギ絡む<盗聴証拠>
今年５月に日本の警視庁防課は被害者SDカード15分を保持した有る国民に出せ!!<創価幹部>
キタオカ1962年東北生は二十代で2人の女性をレイプ殺害して入信した創価本尊はこれだけで潰せる<<<韓国工作員鸛<<<創価公明党 <テロ装置>>東芝部品)>>ヤクザ<宗教<同和<<公安<<魂複<<官憲>日本終Googl検索
207 名前：１３２人目の素数さん mailto:sage [2011/11/22(火) 00:52:44.75 ]: >>205-206
完全にスレ違いだろ。。。
ガイドラインは無いけど、ちゃんと空気を読んでから投稿してね☆
208 名前：１３２人目の素数さん mailto:sage [2011/11/22(火) 22:32:56.09 ]: >>202の4.と>>203の5.はどんな話？
209 名前：ごろ mailto:sage [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と表示*/
210 名前：ごろ mailto:sage [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として入力することが解説されている。
211 名前：ごろ mailto:sage [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を使う方法の解説
212 名前：ごろ mailto:sage [2011/11/25(金) 15:00:09.14 ]: >>211のLecture2-3続きはLecture3-2続きの間違い
Lecture3-4
マクロとscalarとどちらが正しいか？
マクロは１２桁まで、Scalarは１６桁までなのでScalarがより正確
迷うならScalarを使っておけということ！？
213 名前：ごろ mailto:sage [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を使ってもプログラム終了後は順番が元に戻る。
しかし、必要な時のみ使うこと。
214 名前：ごろ mailto:sage [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の後ろの数は文字列の数
欠損値は.で入力
215 名前：ごろ mailto:sage [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
として、変数を決めていく手法が詳細に説明されていました。
216 名前：ごろ mailto:sage [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した時に不完全な新しい変数が出来るのを防ぐ方法
一時的な変数名を作成して、全てが終わった後、新しい変数名に変更する。
217 名前：ごろ mailto:sage [2011/11/26(土) 00:06:54.14 ]: Lecture3-12
練習問題
大部分が自作ｔテストのプログラムの書き換え問題。
218 名前：ごろ mailto:sage [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)などはプログラム途中でも
代入される。
219 名前：１３２人目の素数さん mailto:sage [2011/11/26(土) 19:42:35.40 ]: >>215
わかりました。
ただ、それは、weight*weightではなく、c.weight#c.weightという項を入れるといいんですよね。
220 名前：ごろ mailto:sage [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) 変数マトリックス
221 名前：ごろ mailto:sage [2011/11/27(日) 01:08:16.72 ]: >>219
ありがとうございます。
ちなみに
>>215
> regress weight2 weight foreign
> として、変数を決めていく手法が詳細に説明されていました。
regress mpg weight2 weight foreign
mpgが抜けていました。
222 名前：SEMやってます mailto:sage [2011/11/28(月) 23:16:18.42 ]: 何レスか前にSEMを始めた者です。
大石展緒・都竹浩生『Amosで学ぶ調査系データ解析』の1章にあった、
食物摂取と、大腸ガン、直腸ガンの関連のデータで、同じようにして推定してみようと思ったのですが、推定できません。
このデータ自体は豊田・前田・柳井『原因をさぐる統計学』(1992)のデータだそうなのですが。

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

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

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

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

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

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

もう少し調べます。
223 名前：１３２人目の素数さん mailto:sage [2011/11/29(火) 17:24:16.59 ]: >>222

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

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

下部消化管のガン傾向（潜在変数）を構成している変数の指定はあっていますか？
224 名前：１３２人目の素数さん mailto:sage [2011/11/29(火) 22:18:55.99 ]: >>223
レスどうもです。
中部大学の小塩さんのページ
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ビルダー上では攪乱項を外すことができない、という状況です。
何か設定があるのか、コマンドでやればいけるのか、というのを調べるつもりなのですが、今日はそろそろ寝ます。
225 名前：SEM mailto:sage [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
ですが、推定方法が違うためか、若干数値に差が出ますね。
226 名前：１３２人目の素数さん mailto:sage [2011/12/13(火) 15:14:09.22 ]: ワイは日本人やが聞いた情報によると、もうじき中国はバブルがはじけて昔の貧乏な中国に戻るらしいで
みんなも知っての通りもう経済は破綻してて、取り戻すのは無理なんだそうや

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

まあお前らほどの頭の良い連中には、今さらなくらいのネタやな、かえって失礼かもしれん
お前らからすればもう常識的なくらいの知識やろ！！
227 名前：１３２人目の素数さん [2012/01/02(月) 04:55:30.74 ]: Hello,King!
228 名前：ごろ mailto:sage [2012/01/05(木) 11:05:07.39 ]: 明けましておめでとうございます。
今年は皆様にとって素晴らしい一年になることをお祈りいたします。

Stataオンラインレクチャー途中でほったらかしているので時間を見つけて終わらせます。
229 名前：あんでぃ ◆AfX9iez2fA mailto:sage [2012/02/04(土) 23:51:25.79 ]: あんでぃ
230 名前：１３２人目の素数さん [2012/02/05(日) 21:56:05.19 ]: 　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　／:::::; -‐''"　　　　　　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　　　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　　　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　私は只の数ヲタなんかとは付き合わないわ。
.　　　　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　頭が良くて数学が出来てかっこいい人。それが必要条件よ。
　　　　 | 　　　　l＾,人| 　` `-' 　　ゝ　 |　　　　　　　　さらに Ann.of Math に論文書けば十分条件にもなるわよ。
　　　 | 　　　　 ` -'＼　　　　　ｰ' 　人　　　　　　　　　　一番嫌いなのは論文数を増やすためにくだらない論文を書いて
　　　　|　　　　　　　 /(l 　　　＿＿／　ヽ、　　　　　　　　　　　良い論文の出版を遅らせるお馬鹿な人。
　　　 |　　　　　　　（:::::`‐-､__ 　|::::`､　　ﾋﾆﾆヽ、　　　　　　　　　あなたの論文が Ann of Math に accept される確率は?
　　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　それとも最近は Inv. Math. の方が上かしら？
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼　 |　'､　　＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/
231 名前：１３２人目の素数さん mailto:sage [2012/02/09(木) 18:25:45.89 ]: 静まりかえっていますが、Stata/Transferってソフトの案内が来た。
ただのデータフォーマット変換ソフトに金とるのね。
おそれいった。
232 名前：１３２人目の素数さん mailto:sage [2012/02/25(土) 16:52:29.31 ]: 　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　／:::::; -‐''"　　　　　　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　　　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　　　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　私は只の数ヲタなんかとは付き合わないわ。
.　　　　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　頭が良くて数学が出来てかっこいい人。それが必要条件よ。
　　　　 | 　　　　l＾,人| 　` `-' 　　ゝ　 |　　　　　　　　さらに Ann.of Math に論文書けば十分条件にもなるわよ。
　　　 | 　　　　 ` -'＼　　　　　ｰ' 　人　　　　　　　　　　一番嫌いなのは論文数を増やすためにくだらない論文を書いて
　　　　|　　　　　　　 /(l 　　　＿＿／　ヽ、　　　　　　　　　　　良い論文の出版を遅らせるお馬鹿な人。
　　　 |　　　　　　　（:::::`‐-､__ 　|::::`､　　ﾋﾆﾆヽ、　　　　　　　　　あなたの論文が Ann of Math に accept される確率は?
　　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　それとも最近は Inv. Math. の方が上かしら？
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼　 |　'､　　＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/
233 名前：１３２人目の素数さん [2012/03/08(木) 19:21:43.89 ]: STATAを今使用して重回帰分析を行なっているんですが、非標準化回帰係数の95%信頼区間は表示されるのですが、"標準化回帰係数"の95%信頼区間はどうやったら表示されるのでしょうか？
ご存知の方、よろしくお願いします。
234 名前：１３２人目の素数さん mailto:sage [2012/03/08(木) 19:24:17.96 ]: マルチポストに答える気はなくなりました。
235 名前：１３２人目の素数さん mailto:sage [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

と貼り付けておけば、あとで見たときに区別がつくわけです。
236 名前：１３２人目の素数さん mailto:sage [2012/03/10(土) 17:38:11.84 ]: 知らなかった
237 名前：１３２人目の素数さん mailto:sage [2012/03/23(金) 02:35:06.57 ]: 　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　／:::::; -‐''"　　　　　　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　　　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　　　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　私は只の数ヲタなんかとは付き合わないわ。
.　　　　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　頭が良くて数学が出来てかっこいい人。それが必要条件よ。
　　　　 | 　　　　l＾,人| 　` `-' 　　ゝ　 |　　　　　　　　さらに Ann.of Math に論文書けば十分条件にもなるわよ。
　　　 | 　　　　 ` -'＼　　　　　ｰ' 　人　　　　　　　　　　一番嫌いなのは論文数を増やすためにくだらない論文を書いて
　　　　|　　　　　　　 /(l 　　　＿＿／　ヽ、　　　　　　　　　　　良い論文の出版を遅らせるお馬鹿な人。
　　　 |　　　　　　　（:::::`‐-､__ 　|::::`､　　ﾋﾆﾆヽ、　　　　　　　　　あなたの論文が Ann of Math に accept される確率は?
　　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　それとも最近は Inv. Math. の方が上かしら？
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼　 |　'､　　＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/
238 名前：１３２人目の素数さん mailto:sage [2012/03/25(日) 23:15:22.24 ]: 　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　／:::::; -‐''"　　　　　　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　　　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　　　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　私は只の数ヲタなんかとは付き合わないわ。
.　　　　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　頭が良くて数学が出来てかっこいい人。それが必要条件よ。
　　　　 | 　　　　l＾,人| 　` `-' 　　ゝ　 |　　　　　　　　さらに Ann.of Math に論文書けば十分条件にもなるわよ。
　　　 | 　　　　 ` -'＼　　　　　ｰ' 　人　　　　　　　　　　一番嫌いなのは論文数を増やすためにくだらない論文を書いて
　　　　|　　　　　　　 /(l 　　　＿＿／　ヽ、　　　　　　　　　　　良い論文の出版を遅らせるお馬鹿な人。
　　　 |　　　　　　　（:::::`‐-､__ 　|::::`､　　ﾋﾆﾆヽ、　　　　　　　　　あなたの論文が Ann of Math に accept される確率は?
　　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　それとも最近は Inv. Math. の方が上かしら？
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼　 |　'､　　＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/
239 名前：１３２人目の素数さん mailto:sage [2012/03/26(月) 12:40:49.52 ]: 　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　／:::::; -‐''"　　　　　　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　　　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　　　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　私は只の数ヲタなんかとは付き合わないわ。
.　　　　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　頭が良くて数学が出来てかっこいい人。それが必要条件よ。
　　　　 | 　　　　l＾,人| 　` `-' 　　ゝ　 |　　　　　　　　さらに Ann.of Math に論文書けば十分条件にもなるわよ。
　　　 | 　　　　 ` -'＼　　　　　ｰ' 　人　　　　　　　　　　一番嫌いなのは論文数を増やすためにくだらない論文を書いて
　　　　|　　　　　　　 /(l 　　　＿＿／　ヽ、　　　　　　　　　　　良い論文の出版を遅らせるお馬鹿な人。
　　　 |　　　　　　　（:::::`‐-､__ 　|::::`､　　ﾋﾆﾆヽ、　　　　　　　　　あなたの論文が Ann of Math に accept される確率は?
　　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　それとも最近は Inv. Math. の方が上かしら？
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼　 |　'､　　＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/
240 名前：１３２人目の素数さん mailto:sage [2012/03/26(月) 16:41:13.45 ]: 　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　／:::::; -‐''"　　　　　　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　　　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　　　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　私は只の数ヲタなんかとは付き合わないわ。
.　　　　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　頭が良くて数学が出来てかっこいい人。それが必要条件よ。
　　　　 | 　　　　l＾,人| 　` `-' 　　ゝ　 |　　　　　　　　さらに Ann.of Math に論文書けば十分条件にもなるわよ。
　　　 | 　　　　 ` -'＼　　　　　ｰ' 　人　　　　　　　　　　一番嫌いなのは論文数を増やすためにくだらない論文を書いて
　　　　|　　　　　　　 /(l 　　　＿＿／　ヽ、　　　　　　　　　　　良い論文の出版を遅らせるお馬鹿な人。
　　　 |　　　　　　　（:::::`‐-､__ 　|::::`､　　ﾋﾆﾆヽ、　　　　　　　　　あなたの論文が Ann of Math に accept される確率は?
　　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　それとも最近は Inv. Math. の方が上かしら？
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼　 |　'､　　＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/
241 名前：１３２人目の素数さん mailto:sage [2012/03/26(月) 17:16:15.84 ]: 　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　／:::::; -‐''"　　　　　　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　　　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　　　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　私は只の数ヲタなんかとは付き合わないわ。
.　　　　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　頭が良くて数学が出来てかっこいい人。それが必要条件よ。
　　　　 | 　　　　l＾,人| 　` `-' 　　ゝ　 |　　　　　　　　さらに Ann.of Math に論文書けば十分条件にもなるわよ。
　　　 | 　　　　 ` -'＼　　　　　ｰ' 　人　　　　　　　　　　一番嫌いなのは論文数を増やすためにくだらない論文を書いて
　　　　|　　　　　　　 /(l 　　　＿＿／　ヽ、　　　　　　　　　　　良い論文の出版を遅らせるお馬鹿な人。
　　　 |　　　　　　　（:::::`‐-､__ 　|::::`､　　ﾋﾆﾆヽ、　　　　　　　　　あなたの論文が Ann of Math に accept される確率は?
　　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　それとも最近は Inv. Math. の方が上かしら？
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼　 |　'､　　＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/
242 名前：１３２人目の素数さん [2012/03/26(月) 18:15:41.40 ]: 　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　／:::::; -‐''"　　　　　　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　　　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　　　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　私は只の数ヲタなんかとは付き合わないわ。
.　　　　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　頭が良くて数学が出来てかっこいい人。それが必要条件よ。
　　　　 | 　　　　l＾,人| 　` `-' 　　ゝ　 |　　　　　　　　さらに Ann.of Math に論文書けば十分条件にもなるわよ。
　　　 | 　　　　 ` -'＼　　　　　ｰ' 　人　　　　　　　　　　一番嫌いなのは論文数を増やすためにくだらない論文を書いて
　　　　|　　　　　　　 /(l 　　　＿＿／　ヽ、　　　　　　　　　　　良い論文の出版を遅らせるお馬鹿な人。
　　　 |　　　　　　　（:::::`‐-､__ 　|::::`､　　ﾋﾆﾆヽ、　　　　　　　　　あなたの論文が Ann of Math に accept される確率は?
　　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　それとも最近は Inv. Math. の方が上かしら？
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼　 |　'､　　＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/

[ 新着レスの取得/表示 (agate) ] / [ 携帯版 ]

read.cgi ver5.27 [feat.BBS2 +1.6] / e.0.2 (02/09/03) / eucaly.net products.
担当:undef