Title stata.com probit — Probit regression

1 Probit regressionSyntaxMenuDescriptionOptionsRe marks and examplesStored resultsMethods and formulasReferencesAlso seeSyntaxprobitdepvar[indepvars] [if] [in] [weight] [,options]optionsDescriptionModelnoconst antsuppress constant termoffset(varname)includevarnamein model with coefficient constrained to 1asisretain perfect predictor variablesconstraints(constraints)apply specified linear constraintscollinearkeep collinear variablesSE/Robustvce(vcetype)vcetypemay beoim,robust,clusterclustvar,bootstrap, orjackknifeReportinglevel(#)set confidence level; default islevel(95)nocnsreportdo not display constraintsdisplayoptionscontrol column formats, row spacing, line width, display of omittedvariables and base and empty cells, and factor-variable labelingMaximizationmaximizeoptionscontr ol the maximization process; seldom usednocoefdo not display the coefficient table; seldom usedcoeflegenddisplay legend instead of statisticsindepvarsmay contain factor variables; see[U] Factor contain time-series operators.

2 See[U] Time-series ,by,fp,jackknife,mfp,mi estimate,nestreg,rolling,statsby,stepwis e, andsvyare allowed;see[U] Prefix (bootstrap)andvce(jackknife)are not allowed with themi estimateprefix; see [MI]mi are not allowed with thebootstrapprefix; see [R] (),nocoef, and weights are not allowed with thesvyprefix; see [SVY] ,iweights, andpweights are allowed; see[U] not appear in the dialog [U] 20 Estimation and postestimation commandsfor more capabilities of estimation Probit Probit regressionMenuStatistics>Binary outcomes> Probit regressionDescriptionprobitfits a maximum-likelihood Probit estimating on grouped data, see thebprobitcommand described in [R] auxiliary commands may be run afterprobit,logit, orlogistic.

3 See [R]logisticpostestimationfor a description of these [R]logisticfor a list of related estimation Model noconstant,offset(varname),constraints(c onstraints),collinear; see [R]estimation that all specified variables and observations be retained in the maximization option is typically not specified and may introduce numerical instability. Normallyprobitdrops variables that perfectly predict success or failure in the dependent variable along with theirassociated observations. In those cases, the effective coefficient on the dropped variables is infinity(negative infinity) for variables that completely determine a success (failure).

4 Dropping the variableand perfectly predicted observations has no effect on the likelihood or estimates of the remainingcoefficients and increases the numerical stability of the optimization process. Specifying this optionforces retention of perfect predictor variables and their associated observations. SE/Robust vce(vcetype)specifies the type of standard error reported, which includes types that are derivedfrom asymptotic theory (oim), that are robust to some kinds of misspecification (robust), thatallow for intragroup correlation (clusterclustvar), and that use bootstrap or jackknife methods(bootstrap,jackknife); see [R]vceoption.

5 Reporting level(#); see [R]estimation ; see [R]estimation :noomitted,vsquish,noemptycells,baseleve ls,allbaselevels,nofvla-bel,fvwrap(#),fv wrapon(style),cformat(%fmt),pformat(%fmt ),sformat(%fmt), andnolstretch; see [R]estimation options. Maximization maximizeoptions:difficult,technique(algo rithmspec),iterate(#),[no]log,trace,grad ient,showstep,hessian,showtolerance,tole rance(#),ltolerance(#),nrtolerance(#),no nrtolerance, andfrom(initspecs); see [R]maximize. These options areseldom Probit regression 3 The following options are available withprobitbut are not shown in the dialog box:nocoefspecifies that the coefficient table not be displayed.

6 This option is sometimes used byprogrammers but is of no use ; see [R]estimation and are presented under the following headings:Robust standard errorsModel identificationprobitfits maximum likelihood models with dichotomous dependent (left-hand-side) variablescoded as 0/1 (more precisely, coded as 0 and not 0).Example 1We have data on the make, weight, and mileage rating of 22 foreign and 52 domestic wish to fit a Probit model explaining whether a car is foreign based on its weight and is an overview of our data.

7 Use (1978 Automobile Data). keep make mpg weight foreign. describeContains data from : 74 1978 Automobile Datavars: 4 13 Apr 2013 17:45size: 1,702 (_dta has notes)storage display valuevariable name type format label variable labelmake str18 %-18s Make and Modelmpg int % Mileage (mpg)weight int % Weight (lbs.)foreign byte % origin Car typeSorted by: foreignNote: dataset has changed since last saved.

8 Inspect foreignforeign: Car type Number of ObservationsTotal Integers Nonintegers# Negative - - -# Zero 52 52 -# Positive 22 22 -## # Total 74 74 -# # Missing -0 1 74(2 unique values)foreign is labeled and all values are documented in the Probit Probit regressionTheforeignvariable takes on two unique values, 0 and 1.

9 The value 0 denotes a domestic car,and 1 denotes a foreign model that we wish to fit isPr(foreign= 1) = ( 0+ 1weight+ 2mpg)where is the cumulative normal fit this model, we type. Probit foreign weight mpgIteration 0: log likelihood = 1: log likelihood = (output omitted)Iteration 5: log likelihood = regression Number of obs = 74LR chi2(2) = > chi2 = likelihood = Pseudo R2 = Std. Err. z P>|z| [95% Conf.]

10 Interval] .0005661 .0515689 find that heavier cars are less likely to be foreign and that cars yielding better gas mileage arealso less likely to be foreign, at least holding the weight of the car [R]maximizefor an explanation of the noteStata interprets a value of 0 as a negative outcome (failure) and treats all other values (exceptmissing) as positive outcomes (successes). Thus if your dependent variable takes on the values 0 and1, then 0 is interpreted as failure and 1 as success.