Re: st: ''option ysi not allowed'' when using metan command

[Nick Cox <njcoxstata@gmail.com>]

The reference for that update is
http://www.stata-journal.com/sjpdf.html?articlenum=up0025 and the
version on SSC is more up-to-date.

With -meta- from SSC, once again I can't reproduce your results but
get a different error message "type mismatch".

Nick
njcoxstata@gmail.com


On 10 July 2013 15:47, Amal Khanolkar <Amal.Khanolkar@ki.se> wrote:
> It should be ''sbe24_3.pkg'' which i think is the latest available.
>
> Sorry that was a typo that occurred while I was copying in the syntax from my do file. It's here again, corrected:
>
> *Generate dataset with RR, corresponding CIs, and grouping variables; group (5 nos) and outcome (3 no):
>>
>> input id rr lci uci group outcome
>> 1 1.12 0.78 1.62 2 1
>> 2 0.96 0.83 1.30 2 2
>> 3 0.96 0.85 1.09 2 3
>> 4 0.84 0.65 1.07 3 1
>> 5 0.98 0.90 1.07 3 2
>> 6 0.95 0.88 1.02 3 3
>> 7 1.11 0.69 1.80 4 1
>> 8 1.08 0.90 1.30 4 2
>> 9 0.85 0.73 0.97 4 3
>> 10 1.76 1.24 2.50 5 1
>> 11 1.15 0.98 1.34 5 2
>> 12 0.86 0.75 0.99 5 3
>>  end
>>
>> *Create labels for the group variable:
>>
>> gen str immigrant_group=""
>> replace immigrant_group ="Western europe and North America" if group==2
>> replace immigrant_group ="Finland" if group==3
>> replace immigrant_group ="Eastern Europe" if group==4
>> replace immigrant_group ="Poland" if group==5
>>
>> *Create labels for the Outcome variable:
>>
>> gen str Outcome_type=""
>> replace Outcome_type ="Early preterm" if outcome==1
>> replace Outcome_type ="Late preterm" if outcome==2
>> replace Outcome_type ="Postterm" if outcome==3
>>
>> *Create a sorter variable to instruct stata on the right order:
>>
>> ***defining a sorter variable to bring the graph results into the intended order
>> gen sorter =.
>> replace sorter =1 if immigrant_group=="Western europe and North America" & outcome==1
>> replace sorter =2 if immigrant_group=="Western europe and North America" & outcome==2
>> replace sorter =3 if immigrant_group=="Western europe and North America" & outcome==3
>> replace sorter =4 if immigrant_group=="Finland" & outcome==1
>> replace sorter =5 if immigrant_group=="Finland" & outcome==2
>> replace sorter =6 if immigrant_group=="Finland" & outcome==3
>> replace sorter =7 if immigrant_group=="Eastern Europe" & outcome==1
>> replace sorter =8 if immigrant_group=="Eastern Europe" & outcome==2
>> replace sorter =9 if immigrant_group=="Eastern Europe" & outcome==3
>> replace sorter =10 if immigrant_group=="Poland" & outcome==1
>> replace sorter =11 if immigrant_group=="Poland" & outcome==2
>> replace sorter =12 if immigrant_group=="Poland" & outcome==3
>>
>>
>> *sorting data into the correct order
>> sort sorter immigrant_group outcome
>>
>> *Generate the forestplot:
>>
>> #delimit ;
>> metan rr lci uci,
>> by(immigrant_group) nosubgroup
>>         xlabel( 0.4, 0.8, 1, 1.2, 1.4, 1.8, 2.2, 2.6, 3.0)
>>         xtick(0.4, 0.8, 1, 1.2, 1.4, 1.6, 1.8, 2.0, 2.2, 2.4, 2.6, 3.0)
>>         effect(Risk ratio)
>>         force nowt nooverall nobox
>>         label(namevar=immigrant_group)
>>         sortby(sorter immigrant_group group)
>>         lcols(immigrant_group)
>>         title("Adjusted relative risks for preterm and postterm birth by immigrant groups in Sweden", size(small) color(black))
>>     pointopt( msymbol(T) mcolor(black) msize(medsmall))
>>     ciopt( lcolor(gs8) lwidth(medium) )
>>    favours(Pretermbirth # Postterm birth)
>>     graphregion(fcolor(white) lcolor(gs15) ifcolor(white) ilcolor(white)) ;
>> #delimit cr
>>
>>
>> The output that I get:
>>
>>
>> Study     |     ES    [95% Conf. Interval]
>> ---------------------+---------------------------------------------------
>>      Western europe and N
>> Western europe and N |    1.120     0.780     1.620
>> Western europe and N |    0.960     0.830     1.300
>> Western europe and N |    0.960     0.850     1.090
>> ---------------------+---------------------------------------------------
>>      Finland
>> Finland              |    0.840     0.650     1.070
>> Finland              |    0.980     0.900     1.070
>> Finland              |    0.950     0.880     1.020
>> ---------------------+---------------------------------------------------
>>      Eastern Europe
>> Eastern Europe       |    1.110     0.690     1.800
>> Eastern Europe       |    1.080     0.900     1.300
>> Eastern Europe       |    0.850     0.730     0.970
>> ---------------------+---------------------------------------------------
>>      Poland
>> Poland               |    1.760     1.240     2.500
>> Poland               |    1.150     0.980     1.340
>> Poland               |    0.860     0.750     0.990
>>
>> option ysi not allowed
>> r(198);
>>
>> end of do-file
>>
>> r(198);
*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/faqs/resources/statalist-faq/
*   http://www.ats.ucla.edu/stat/stata/