Stata同志コミュのあしあと〜メッセージ

１５０人越えて喜んでいたのですが、又１５０人に戻りましたね。
忘れられないよう、気合をいれて頑張ります。
是非お引き立てのほどを。

現在１５８人。２００人目指してじっくり頑張ります！

はじめまして。大学院の指導教官がかわり、急遽SPSSからSTATAに移行しなければならなくなりました。基本事項をマスターするのでさえアップアップなのに、STATAを使って"メタアナリシス"、"モンテカルロシュミレーション"などをしなければならなくなりました。現在、初学者にも分かりやすい参考書の情報を集めながら、実際の解析の準備をしています。（分かりやすい参考書にはまだ出会えていませんが...）こちらで色々と情報交換させていただければ励みになります。どうぞよろしくお願いいたします。

Tefeさん、こちらもよろしくおねがいしますね。StataにMetaAnalysisのコマンドはありますが、ご存知ですか？英語でよろしかったら情報が載っているサイトをひっぱってきますよ。

StataのMetaAnalysisコマンドは秀逸ですね。
グラフなど出版レベルで非常にきれいにできますね。

世捨て人さん、ののさんコメントありがとうございます。Meta-Analysis in Stata:An Updated Collection from the Stata Journalという参考書を買ってみたのですが、初心者向けではないようで、混乱してしまいました...。英語でも頑張って解読しますので、情報が載っているサイトをおしえていただけたらと思います。よろしくお願いいたします。

Tefeさん、了解です。色々探しますね。

こんな情報がありました。
INSTALL META-ANALYSIS COMMANDS IN STATA
The commands are not included with Stata, but can be easily downloaded directly within Stata. Meta-analysis commands which you will use as part of the Systematic Review course include:

• metan: performs meta-analysis
• metabias: performs tests for funnel plot asymmetry and also produces funnel plots
• metafunnel: produces funnel plots
• metareg: performs meta-regression

To download these commands, open Stata and type the following in the command window:

ssc install metan, replace
ssc install metabias, replace
ssc install metafunnel, replace
ssc install metareg, replace
net install sbe16_2, from (http://www.stata.com/stb/stb43) replace

These packages should remain on the system once installed (you need not reinstall the packages each time you intend to use them).

To calculate the summary estimate, you want to use the metan command. The metan command performs a meta-analysis of studies with two groups and calculates a pooled absolute risk reduction, relative risk, or odds ratio.

There are a number of ways to evaluate for publication bias in Stata. The simplest is to use the straightforward Kendall’s Tau test:

ktau _ES _selogES

The p-value for this test is shown as “Prob > z = (p-value) (continuity corrected). If the p-value is less than 0.5, the test is positive for publication bias.
Remember _ES and _seES are the RR and standard error of the RR generated by Stata. Stata also has some more sophisticated ways to look for publication bias in a sample of studies. The metabias command computes two tests. The first is the Begg adjusted rank correlation test, which is a variant of the Kendall’s Tau, also looking for significant correlation between the effect estimates and their standard errors. This test tends to have a low power for detecting publication bias.

The second test computed by the metabias command is the Egger regression asymmetry test, which uses regression to looks for funnel plot asymmetry. The p-value for the bias of the Egger regression is the p-value we are interested in to assess publication bias.

metabias _ES _selogES, graph(begg)

The graph(begg) option produces a funnel plot of the effect size versus the standard error of the effect size. You can produce a similar funnel plot using the metafunnel command. Note that the standard error of the effect size should be proportional to the sample size, because the smaller the study is the larger the variance.

metafunnel _ES _selogES

A funnel plot like the one illustrated in the syllabus, graphing the relative risk of each study against its sample size, can be produced by giving Stata the following command:

funnel _SS _selogES, xlog sample

The sample option tells Stata that you are graphing against the sample size, and the xlog option tells Stata to create a symmetrical plot (required since our effect size is a relative risk, which is on a log scale).

ののさん、詳しい情報を添付してくださりありがとうございます。

お久しぶりです。
少し、困っているので助けてください。
twoway を使って、横軸が日付、縦軸が測定値の折れ線グラフを書こうと思っています。
データ形式は、reshape long ・・・で変更済みです。

ところが、データ全てをつないだりとかlowess を使った方法ではグラフが作成できるのですが、エクセルなんかで言う普通の折れ線グラフが書けません。
いちいち平均値を出してやるしかないんでしょうか？
お教え下さい。

上記例
twoway line value day
twoway lowess value day

結局、１０の時にはデータの平均値を出してconnectedでグラフがかけました。
ただ、sortコマンドを入れないとぐちゃぐちゃのグラフになってしまいました。

もっと簡単に折れ線グラフをかく方法があればまた教えて下さい。

ののさん、仕事が忙しくてしばらく書き込めずにいてごめんなさい。なんとかグラフがかけてよかったです。私はグラフには詳しくないのですが、周りにもきいてみますね。皆さんの中でわかる人がいらしたらよろしくご回答願います。

グラフですが、こういうサイトがありました。
ご参考までに。
http://www.ats.ucla.edu/stat/stata/topics/graphics.htm