Statalist


[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: st: RE: large numbers in comb(n,k) function: no success


From   Rafal Raciborski <rraciborski@gmail.com>
To   statalist@hsphsun2.harvard.edu
Subject   Re: st: RE: large numbers in comb(n,k) function: no success
Date   Thu, 5 Feb 2009 08:35:58 -0500

if you have access to Mathematica, you can calculate it but i am not
sure whether it will be helpful, for example, here is the result of
Binomial[180000,2000]:

152517439466441991164450127611194712847529937178533379587654534576248355254531\
769463839975669536718026413455843661312478666027392290336784287127984185915677\
238417670099802220636148426827215980792048571828138587881611667807724251676306\
654768363896274293226399517765318241706996915730654898637853073245523683217375\
950822175172528516274560751954053705336587900275599167367308285904209194721513\
247599977155722439110621279586153323123082788657252422899363927426477266208262\
676657081727196330248289305705422508108488376257444519013568170479068807446854\
256042505091737800373037224230302972622534710029802423146794325436715084107554\
907669851532044765094178406173295161695402752290403306211576760486114495312573\
678706358328233678575601429860877597992474257572256588548203975776900481716015\
201395354794241422669083128479686848264285509119702994065807226662357733305323\
497019115334535067753531538834549941130380917843586839879431655611622687835901\
932962676219190543314861254357447527596180338193319460354827342679652323452864\
780869203251910906085617835068411610781193069151533335506137452259790161680556\
709001307889715765390821277247198952613051892777420993900174151647923955954877\
473821389380644572565985226209107513991808109507382999675333424049988222238061\
922084130658255046913412563524495791299065879762834881501097758743456772414528\
281620389469674744158598189543324307398208175823088821684280703640009112503578\
833718508569441473685559980015042894027519606687120226875687221666472383131416\
840206460392016154757500906478213204581220571520470533488611555555723865094114\
775500863861968818872527785408296672674648690178007773974784748687870380547005\
885169993748428660636279603581878181405887320974912477346188950309047531301336\
130808970258008914380114075940500941216744784661927629280657806660366715344085\
462721782722841194348336152260765636876719104960889677325080549140717356587125\
819615647407352948197882835083950719919385215903516065649976234797486196121534\
790133974077235261128004556133123982761898621019008510473531918506333849668198\
277848278193416474649094360151613920852261398511921063502162677437631079046262\
071343600972522365624764307241076106019736548181353404876489609171862585124340\
593840544095534179039319482987808400094819872307397538440873075419965995484273\
895032084258234451748864893602126612916598790215427568127292534466182119015403\
220497055132110395022420519243843452125793846659879961641473014972827803346239\
269524946359890255996467873189635205294579129983920925514618978015052702303557\
835953984459577577259720432899556138826115541147780099933630107426448995358028\
094466550200179216708570733289129770992844411631239564263412812563361061079504\
374937635943519812966135378299316698303483445149353560332643585399901440154127\
391434721230174059922765497121814061567005211484863001004753256945975002183068\
865898862634623602592188120585623652626664975582932381384270298757444769901234\
170705860342967242610418204201422896698983688946372740625489936003287546759570\
059563330031662274710757905680361783140568219413403110337651089787762956845609\
777577466524911502616860033139748971763652792094702593430891493977709122071287\
365235908574608074978689267077437023307965270233390163831271512737473602797842\
984014827947665665619034749476593816992735620900345687912472167979768740590015\
882648863012701770475077145673127100729584929121042178008507473898790966560437\
303174294104691462837341155804450068703353340870831180088160990885536821542281\
735838189426301662188931583955889628508360023793853230415344103073313781799584\
338027228620847402754500105983204104632999325341809211303368052056326898714510\
859135654350695878104833779790479780294450497500860936367588322263541138247098\
257071598308592664030789058334025704430419211203756309309087045328159389208796\
310810006015649579924183208080571661820121684951978351278876758914869677673749\
805979296476865157046744168326883083631798447427958510273314684110890599658784\
177014163928161063356454929484784585760330945364281475395083706045433798562140\
877159777757701195963753853603960459538680876482640305145270088713199046533845\
151410694974436811863488488027592143724435503664143466015626902203058178608783\
577316360626991560895377980946207909133329414162474835164977602635652627248791\
674911683250749576948526920021718711661263237709968613409497281608110999731268\
731954986629182650410149374568504322453534096171877773720383107656185355263585\
712779042569007515896753016716138889265857078489903422593565102186841458945038\
208397178140611738611019670264616858203234665622143738687430051733046683335240\
636043341724656510885455380682397990066030346894991918010378827217501574301556\
194985444382326457199476414805537635058836145017677973323359841057486688872547\
671748165984793007424427545106052945017031356635413767098609836861449607728127\
0258773556800




On Thu, Feb 5, 2009 at 8:09 AM, Maarten buis <maartenbuis@yahoo.co.uk> wrote:
>
> --- On Thu, 5/2/09, Inna Becher wrote:
>> I have to calculate  probabilities for observations to be
>> sampled. That's why I use exp with lnfactorial...
>
> Looks like you want to calculate a likelihood function. Is that right?
> If so, then you can use the fact that you almost always maximize the log likelihood, which will alow you to avoid taking the exp of lnfactorial.
>
> -- Maarten
>
> -----------------------------------------
> Maarten L. Buis
> Department of Social Research Methodology
> Vrije Universiteit Amsterdam
> Boelelaan 1081
> 1081 HV Amsterdam
> The Netherlands
>
> visiting address:
> Buitenveldertselaan 3 (Metropolitan), room N515
>
> +31 20 5986715
>
> http://home.fsw.vu.nl/m.buis/
> -----------------------------------------
>
>
>
>
>
>
>
> *
> *   For searches and help try:
> *   http://www.stata.com/help.cgi?search
> *   http://www.stata.com/support/statalist/faq
> *   http://www.ats.ucla.edu/stat/stata/
>
*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/



© Copyright 1996–2014 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index