Announcement

Collapse
No announcement yet.
X
  • Filter
  • Time
  • Show
Clear All
new posts

  • Confidence interval by groups

    Hello everyone,
    I have a dataset, where I looked for confidence interval for groups of "a test" coded from 0-10. I have around 1000 data points, and want the confidence interval for group 0 (110 have a test result of 0) for group 2 (50 have a test result 2) and so on. How is this possible in stata. I tried to use the frequencies but it still won't work, when I make a dataset with the amount of data points per each group and the test result (0-10), so:

    Code:
    input frequency test
    110 0
    50 2
    and so on...
    then I tried ci prop frequency by(test) but got:
    "factor variables and time-series operators not allowed"

    The best would be if I could get a command how to do it in my original dataset where only the id and the result is written, but when it ís only possible in my new dataset with frequencies it is also ok

    I manged to do it with the cii command:
    cii prop 1000 110
    but would be happy to find an easier way

    Best regards
    Henrik

  • #2
    I'm not sure that I follow you, but it seems as if you're ignoring the dependence between the score categories, given that they tally to 1000. If the test results are ordered categories, then you can try something like ologit and if unordered categories, then something like mlogit. Both give confidence intervals for k - 1 groups relative to the reference group (test-score category).

    .ÿversionÿ14.1

    .ÿ
    .ÿclearÿ*

    .ÿsetÿmoreÿoff

    .ÿsetÿseedÿ`=date("2015-11-04",ÿ"YMD")'

    .ÿquietlyÿsetÿobsÿ1000

    .ÿranintÿtest_result,ÿb(10)

    .ÿcontractÿtest_result,ÿfreq(count)

    .ÿlist,ÿnoobsÿseparator(0)ÿabbreviate(20)

    ÿÿ+---------------------+
    ÿÿ|ÿtest_resultÿÿÿcountÿ|
    ÿÿ|---------------------|
    ÿÿ|ÿÿÿÿÿÿÿÿÿÿÿ0ÿÿÿÿÿÿ82ÿ|
    ÿÿ|ÿÿÿÿÿÿÿÿÿÿÿ1ÿÿÿÿÿÿ92ÿ|
    ÿÿ|ÿÿÿÿÿÿÿÿÿÿÿ2ÿÿÿÿÿÿ79ÿ|
    ÿÿ|ÿÿÿÿÿÿÿÿÿÿÿ3ÿÿÿÿÿÿ98ÿ|
    ÿÿ|ÿÿÿÿÿÿÿÿÿÿÿ4ÿÿÿÿÿÿ91ÿ|
    ÿÿ|ÿÿÿÿÿÿÿÿÿÿÿ5ÿÿÿÿÿÿ93ÿ|
    ÿÿ|ÿÿÿÿÿÿÿÿÿÿÿ6ÿÿÿÿÿ102ÿ|
    ÿÿ|ÿÿÿÿÿÿÿÿÿÿÿ7ÿÿÿÿÿÿ94ÿ|
    ÿÿ|ÿÿÿÿÿÿÿÿÿÿÿ8ÿÿÿÿÿÿ94ÿ|
    ÿÿ|ÿÿÿÿÿÿÿÿÿÿÿ9ÿÿÿÿÿÿ78ÿ|
    ÿÿ|ÿÿÿÿÿÿÿÿÿÿ10ÿÿÿÿÿÿ97ÿ|
    ÿÿ+---------------------+

    .ÿologitÿtest_resultÿ[fweight=count],ÿorÿnolog

    OrderedÿlogisticÿregressionÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿNumberÿofÿobsÿÿÿÿÿ=ÿÿÿÿÿÿ1,000
    ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿLRÿchi2(0)ÿÿÿÿÿÿÿÿ=ÿÿÿÿÿÿÿ0.00
    ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿProbÿ>ÿchi2ÿÿÿÿÿÿÿ=ÿÿÿÿÿÿÿÿÿÿ.
    Logÿlikelihoodÿ=ÿ-2394.4097ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿPseudoÿR2ÿÿÿÿÿÿÿÿÿ=ÿÿÿÿÿ0.0000

    ------------------------------------------------------------------------------
    ÿtest_resultÿ|ÿOddsÿRatioÿÿÿStd.ÿErr.ÿÿÿÿÿÿzÿÿÿÿP>|z|ÿÿÿÿÿ[95%ÿConf.ÿInterval]
    -------------+----------------------------------------------------------------
    ÿÿÿÿÿÿÿ/cut1ÿ|ÿÿÿ.0893246ÿÿÿ.0102954ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.0712629ÿÿÿÿ.1119642
    ÿÿÿÿÿÿÿ/cut2ÿ|ÿÿ-1.557539ÿÿÿ.0834133ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ-1.721027ÿÿÿ-1.394052
    ÿÿÿÿÿÿÿ/cut3ÿ|ÿÿ-1.082676ÿÿÿÿ.072741ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ-1.225245ÿÿÿ-.9401059
    ÿÿÿÿÿÿÿ/cut4ÿ|ÿÿ-.6146465ÿÿÿ.0662558ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ-.7445055ÿÿÿ-.4847875
    ÿÿÿÿÿÿÿ/cut5ÿ|ÿÿ-.2330491ÿÿÿ.0636754ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ-.3578506ÿÿÿ-.1082476
    ÿÿÿÿÿÿÿ/cut6ÿ|ÿÿÿ.1402293ÿÿÿ.0634011ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.0159655ÿÿÿÿ.2644932
    ÿÿÿÿÿÿÿ/cut7ÿ|ÿÿÿ.5623668ÿÿÿ.0657623ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.4334751ÿÿÿÿ.6912586
    ÿÿÿÿÿÿÿ/cut8ÿ|ÿÿÿ.9997021ÿÿÿ.0713124ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.8599323ÿÿÿÿ1.139472
    ÿÿÿÿÿÿÿ/cut9ÿ|ÿÿÿ1.550597ÿÿÿÿ.083225ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ1.387479ÿÿÿÿ1.713715
    ÿÿÿÿÿÿ/cut10ÿ|ÿÿÿ2.231012ÿÿÿ.1068489ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ2.021591ÿÿÿÿ2.440432
    ------------------------------------------------------------------------------

    .ÿ
    .ÿmlogitÿtest_resultÿ[fweight=count],ÿbaseoutcome(0)ÿrrrÿnolog

    MultinomialÿlogisticÿregressionÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿNumberÿofÿobsÿÿÿÿÿ=ÿÿÿÿÿÿ1,000
    ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿLRÿchi2(0)ÿÿÿÿÿÿÿÿ=ÿÿÿÿÿÿÿ0.00
    ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿProbÿ>ÿchi2ÿÿÿÿÿÿÿ=ÿÿÿÿÿÿÿÿÿÿ.
    Logÿlikelihoodÿ=ÿ-2394.4097ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿPseudoÿR2ÿÿÿÿÿÿÿÿÿ=ÿÿÿÿÿ0.0000

    ------------------------------------------------------------------------------
    ÿtest_resultÿ|ÿÿÿÿÿÿÿÿRRRÿÿÿStd.ÿErr.ÿÿÿÿÿÿzÿÿÿÿP>|z|ÿÿÿÿÿ[95%ÿConf.ÿInterval]
    -------------+----------------------------------------------------------------
    0ÿÿÿÿÿÿÿÿÿÿÿÿ|ÿÿ(baseÿoutcome)
    -------------+----------------------------------------------------------------
    1ÿÿÿÿÿÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ1.121951ÿÿÿ.1703914ÿÿÿÿÿ0.76ÿÿÿ0.449ÿÿÿÿÿ.8331083ÿÿÿÿ1.510937
    -------------+----------------------------------------------------------------
    2ÿÿÿÿÿÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ.9634146ÿÿÿ.1518818ÿÿÿÿ-0.24ÿÿÿ0.813ÿÿÿÿÿ.7073295ÿÿÿÿ1.312214
    -------------+----------------------------------------------------------------
    3ÿÿÿÿÿÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ1.195122ÿÿÿ.1788663ÿÿÿÿÿ1.19ÿÿÿ0.234ÿÿÿÿÿ.8912885ÿÿÿÿÿ1.60253
    -------------+----------------------------------------------------------------
    4ÿÿÿÿÿÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ1.109756ÿÿÿ.1689752ÿÿÿÿÿ0.68ÿÿÿ0.494ÿÿÿÿÿ.8234187ÿÿÿÿ1.495665
    -------------+----------------------------------------------------------------
    5ÿÿÿÿÿÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ1.134146ÿÿÿ.1718066ÿÿÿÿÿ0.83ÿÿÿ0.406ÿÿÿÿÿÿÿÿ.8428ÿÿÿÿ1.526208
    -------------+----------------------------------------------------------------
    6ÿÿÿÿÿÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ1.243902ÿÿÿ.1844965ÿÿÿÿÿ1.47ÿÿÿ0.141ÿÿÿÿÿ.9301125ÿÿÿÿ1.663555
    -------------+----------------------------------------------------------------
    7ÿÿÿÿÿÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ1.146341ÿÿÿ.1732206ÿÿÿÿÿ0.90ÿÿÿ0.366ÿÿÿÿÿ.8524938ÿÿÿÿ1.541476
    -------------+----------------------------------------------------------------
    8ÿÿÿÿÿÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ1.146341ÿÿÿ.1732206ÿÿÿÿÿ0.90ÿÿÿ0.366ÿÿÿÿÿ.8524938ÿÿÿÿ1.541476
    -------------+----------------------------------------------------------------
    9ÿÿÿÿÿÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ.9512195ÿÿÿÿ.150448ÿÿÿÿ-0.32ÿÿÿ0.752ÿÿÿÿÿÿ.697673ÿÿÿÿ1.296909
    -------------+----------------------------------------------------------------
    10ÿÿÿÿÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ1.182927ÿÿÿ.1774564ÿÿÿÿÿ1.12ÿÿÿ0.263ÿÿÿÿÿÿ.881587ÿÿÿÿ1.587269
    ------------------------------------------------------------------------------

    .ÿ
    .ÿexit

    endÿofÿdo-file


    .
    Attached Files

    Comment


    • #3
      include a weight?
      Doug Hemken
      SSCC, Univ. of Wisc.-Madison

      Comment


      • #4
        Sorry, for my late response, that worked perfectly, thanks Joseph Coveney

        Comment

        Working...
        X