document: the kepler mission @xcite , with the discovery of over 4100 planetary candidates in 3200 systems , has spawned a revolution in our understanding of planet occurrence rates around stars of all types . 
   one of kepler s profound discoveries is that small planets ( @xmath8 ) are nearly ubiquitous ( e.g. , * ? ? ? 
   * ; * ? ? ? 
   * ; * ? ? ? 
   * ; * ? ? ? 
   * ; * ? ? ? 
   * ) and , in particular , some of the most common planets have sizes between earth - sized and neptune - sized  a planet type not found in our own solar system . 
   indeed , it is within this group of super - earths to mini - neptunes that there is a transition from `` rocky '' planets to `` non - rocky planets '' ; the transition is near a planet radius of @xmath9 and is very sharp  occurring within @xmath10 of this transition radius @xcite .    unless an intra - system comparison of planetary radii is performed where only the relative planetary sizes are important @xcite , having accurate ( as well as precise ) planetary radii is crucial to our comprehension of the distribution of planetary structures . 
   in particular , understanding the radii of the planets to within @xmath11 is necessary if we are to understand the relative occurrence rates of `` rocky '' to `` non - rocky '' planets , and the relationship between radius , mass , and bulk density .. while there has been a systematic follow - up observation program to obtain spectroscopy and high resolution imaging , only approximately half of the kepler candidate stars have been observed ( mostly as a result of the brightness distribution of the candidate stars ) . 
   those stars that have been observed have been done mostly to eliminate false positives , to determine the stellar parameters of host stars , and to search for nearby stars that may be blended in the kepler photometric apertures . 
   stars that are identified as possible binary or triple stars are noted on the kepler community follow - up observation program website , and are often handled in individual papers ( e.g. , * ? ? ? 
   * ; * ? ? ? 
   * ) . the false positive assessment of an koi ( or all of the kois ) can take into account the likelihood of stellar companions ( e.g. , * ? ? ? 
   * ; * ? ? ? 
   * ) , and a false positive probability will likely be included in future koi lists . but presently , the current production of the planetary candidate koi list and the associated parameters are derived assuming that _ all _ of the koi host stars are single . 
   that is , the kepler pipeline treats each kepler candidate host star as a single star ( e.g. , * ? ? ? 
   * ; * ? ? ? 
   * ; * ? ? ? 
   thus , statistical studies based upon the kepler candidate lists are also assuming that all the stars in the sample set are single stars . 
   the exact fraction of multiple stars in the kepler candidate list is not yet determined , but it is certainly not zero . 
   recent work suggests that a non - negligible fraction ( @xmath12 ) of the kepler host stars may be multiple stars @xcite , although other work may indicate that ( giant ) planet formation may be suppressed in multiple star systems @xcite . 
   the presence of a stellar companion does not necessarily invalidate a planetary candidate , but it does change the observed transit depths and , as a result , the planetary radii . thus , assuming all of the stars in the kepler candidate list are single can introduce a systematic uncertainty into the planetary radii and occurrence rate distributions . 
   this has already been discussed for the occurrence rate of hot jupiters in the kepler sample where it was found that @xmath13 of hot jupiters were classified as smaller planets because of the unaccounted effects of transit dilution from stellar companions @xcite .    in this paper , we explore the effects of undetected stellar gravitationally bound companions on the observed transit depths and the resulting derived planetary radii for the entire kepler candidate sample . 
   we do not consider the dilution effects of line - of - sight background stars , rather only potential bound companions , as companions within @xmath14 are most likely bound companions ( e.g. , * ? ? ? 
   * ; * ? ? ? 
   * ) , and most stars beyond @xmath14 are either in the kepler input catalog @xcite or in the ukirt survey of the kepler field and , thus , are already accounted for with regards to flux dilution in the kepler project transit fitting pipeline . 
   within 1 , the density of blended background stars is fairly low , ranging between @xmath15 stars/@xmath16 @xcite . 
   thus , within a radius of 1 , we expect to find a blended background ( line - of - sight ) star only @xmath17% of the time . therefore , the primary contaminant within 1  of the host stars are bound companions . 
   we present here probabilistic uncertainties of the planetary radii based upon expected stellar multiplicity rates and stellar companion sizes . 
   we show that , in the absence of any spectroscopic or high resolution imaging observations to vet companions , the observed planetary radii will be systematically too small . 
   however , if a candidate host star is observed with high resolution imaging ( hri ) or with radial velocity ( rv ) spectroscopy to screen the star for companions , the underestimate of the true planet radius is significantly reduced . 
   while imaging and radial velocity vetting is effective for the kepler candidate host stars , it will be even more effective for the k2 and tess candidates which will be , on average , 10 times closer than the kepler candidate host stars . 
   the planetary radii are not directly observed ; rather , the transit depth is the observable which is then related to the planet size . 
   the observed depth ( @xmath18 ) of a planetary transit is defined as the fractional difference in the measured out - of - transit flux ( @xmath19 ) and the measured in - transit flux @xmath20 : @xmath21 if there are @xmath22 stars within a system , then the total out - of - transit flux in the system is given by @xmath23 and if the planet transits the @xmath24 star in the system , then the in  transit flux can be defined as @xmath25 where @xmath26 is the flux of the star with the transiting planet , @xmath27 is the radius of the planet , and @xmath28 is the radius of the star being transited . substituting into equation  ( [ eq - single - flux ] ) , the generalized transit depth equation ( in the absence of limb darkening or star spots ) becomes @xmath29    for a single star , @xmath30 and the transit depth expression simplifies to just the square of the size ratio between the planet and the star . 
   however , for a multiple star system , the relationship between the observed transit depth and the true planetary radius depends upon the brightness ratio of the transited star to the total brightness of the system _ and _ on the stellar radius which changes depending on which star the planet is transiting : @xmath31    the kepler planetary candidates parameters are estimated assuming the star is a single star @xcite , and , therefore , may incorrectly report the planet radius if the stellar host is really a multiple star system . the extra flux contributed by the companion stars 
   will dilute the observed transit depth , and the derived planet radius depends on the size of star presumed to be transited . 
   the ratio of the true planet radius , @xmath32 , to the observed planet radius assuming a single star with no companions , @xmath33 , can be described as : @xmath34 where @xmath35 is the radius of the ( assumed single ) primary star , and @xmath26 and @xmath28 are the brightness and the radius , respectively , of the star being transited by the planet . 
   this ratio reduces to unity in the case of a single star ( @xmath36 and @xmath37 ) . for a multiple star system where the planet orbits the primary star ( @xmath38 ) , 
   the planet size is underestimated only by the flux dilution factor : @xmath39 however , if the planet orbits one of the companion stars and not the primary star , then the ratio of the primary star radius ( @xmath35 ) to the radius of the companion star being transited ( @xmath28 ) affects the observed planetary radius , in addition to the flux dilution factor .    ) 
   are plotted as a function of companion - to - primary brightness ratios ( _ bottom axis _ ) and mass ratios ( _ top axis _ ) for possible binary systems ( _ top plot _ ) or triple ( _ bottom plot _ ) systems . 
   this figure is an example for the g - dwarf koi-299 ; similar calculations have been made for every koi . in each plot , the dark blue stars represent the correction factors if the planet orbits the primary star ( equation [ eq - primary - ratio ] ) ; the red circles represent the correction factors if the planet orbits the secondary star , and the light blue triangles represent the correction factors if the planet orbits the tertiary star ( equation [ eq - full - ratio ] ) . 
   the lines are third order polynomials fit to the distributions . 
   unity is marked with a horizontal dashed line . ]    , but for the m - dwarf koi-1085 , demonstrating that the details of the derived correction factors are dependent upon the koi properties . ] 
   to explore the possible effects of the undetected stellar companions on the derived planetary parameters , we first assess what companions are possible for each koi . 
   for this work , we have downloaded the cumulative kepler candidate list and stellar parameters table from the nasa exoplanet archive . 
   the cumulative list is updated with each new release of the koi lists ; as a result , the details of any one star and planet may have changed since the analysis for this paper was done . 
   however , the overall results of the paper presented here should remain largely unchanged .    for the koi lists , 
   the stellar parameters for each koi were determined by fitting photometric colors and spectroscopically derived parameters ( where available ) to the dartmouth stellar evolution database @xcite . 
   the planet parameters were then derived based upon the transit curve fitting and the associated stellar parameters . 
   other stars listed in the kepler input catalog or ukirt imaging that may be blended with the koi host stars were accounted for in the transit fitting , but , in general , as mentioned above , each planetary host star was assumed to be a single star . 
   we have restricted the range of possible bound stellar companions to each koi host star by utilizing the same dartmouth isochrones used to determine the stellar parameters . 
   possible gravitationally bound companions are assumed to lie along on the same isochrone as the primary star . 
   for each koi host star , we found the single best fit isochrone ( characterized by mass , metallicity , and age ) by minimizing the chi - square fit to the stellar parameters ( effective temperature , surface gravity , radius , and metallicity ) listed in the koi table . 
   we did not try to re - derive stellar parameters or independently find the best isochrone fit for the star ; we simply identified the appropriate dartmouth isochrone as used in the determination of the stellar parameters @xcite . 
   we note that there exists an additional uncertainty based upon the isochrone finding . in this work 
   , we did not try to re - derive the stellar parameters of the host stars , but rather , we simply find the appropriate isochrone that matches the koi stellar parameters . 
   thus , any errors in the stellar parameters derivations in the koi list are propagated here . 
   this is likely only a significant source of uncertainty for nearly equal brightness companions . 
   once an isochrone was identified for a given star , all stars along an isochrone with ( absolute ) kepler magnitudes fainter than the ( absolute ) kepler magnitude of the host star were considered to be viable companions ; i.e. , the primary host star was assumed to be the brightest star in the system . 
   the fainter companions listed within that particular isochrone were then used to establish the range of possible planetary radii corrections ( equation [ eq - full - ratio ] ) assuming the host star is actually a binary or triple star . 
   higher order ( e.g. , quadruple ) stellar multiples are not considered here as they represent only @xmath40% of the stellar population @xcite . 
   we have considered six specific multiplicity scenarios :    1 . 
   single star ( @xmath41 ) 2 . 
   binary star 
   planet orbits primary star 3 .   binary star 
   planet orbits secondary star 4 . 
   triple star 
   planet orbits primary star 5 . 
   triple star  planet orbits secondary star 6 . 
   triple star 
   planet orbits tertiary star .    based upon the brightness and size differences between the primary star and the putative secondary or tertiary companions , we have calculated for each koi the possible factor by which the planetary radii are underestimated ( @xmath0 ) . if the star is single , the correction factor is unity , and if , in a multiple star system , the planet orbits the primary star , only flux dilution affects the observed transit depth and the derived planetary radius ( eq .  [ eq - primary - ratio ] ) .    for the scenarios where the planets orbits the secondary or tertiary star , the planet size correction factors ( eq . 
   [ eq - full - ratio ] ) were determined only for stellar companions where the stellar companion could physically account for the observed transit depth . 
   if more than 100% of the stellar companion light had to be eclipsed in order to produce the observed transit in the presence of the flux dilution , then that star ( and all subsequent stars on the isochrone with lower mass ) was not considered viable as a potential source of the transit . 
   for example , for an observed 1% transit , no binary companions can be fainter than the primary star by 5 magnitudes or more ; an eclipse of such a secondary star would need to be more than 100% deep . 
   the stellar brightness limits were calculated independently for each planet within a koi system so as to not assume that all planets within a system necessarily orbited the same star . 
   figures [ fig - koi299 ] and [ fig - koi1085 ] show representative correction factors ( @xmath0 ) for koi-299 ( a g - dwarf with a super - earth sized @xmath42 planet ) and for koi-1085 ( an m - dwarf with an earth - sized @xmath43 planet ) . 
   the planet radius correction factors ( @xmath0 ) are shown as a function of the companion  to  primary brightness ratio ( bottom x - axis of plots ) and the companion  to  primary mass ratio ( top x - axis of plots ) and are determined for the koi assuming it is a binary - star system ( top plot ) or a triple - star system ( bottom plot ) .    the amplitude of the correction factor ( @xmath0 ) varies strongly depending on the particular system and which star the planet may orbit . if the planet orbits the primary star , then the largest the correction factors are for equal brightness companions ( @xmath44 for a binary system and @xmath45 for a triple system ) with an asymptotic approach to unity as the companion stars become fainter and fainter 
   . if the planet orbits the secondary or tertiary star , the planet radius correction factor can be significantly larger  ranging from @xmath46 for binary systems and @xmath47 or more for triple systems  depending on the size and brightness of the secondary or tertiary star . 
   it is important to recognize the full range of the possible correction factors , but in order to have a better understanding of the statistical correction any given koi ( or the koi list as a whole ) may need , we must understand the mean correction for any one multiplicity scenario and convert these into a single mean correction factor for each star . to do this 
   , we must take into account the probability the star may be a multiple star , the distribution of mass ratios if the star is a multiple , the probability that the planet orbits any one star if the stellar system has multiple stars , and whether or not the star has been vetted ( and how well it is has been vetted ) for stellar companions .    in order to calculate an average correction factor for each multiplicity scenario , 
   we have fitted the individual scenario correction factors as a function of mass ratio with a 3@xmath48-order polynomial ( see fig . 
   [ fig - koi299 ] and [ fig - koi1085 ] ) . because the isochrones are not evenly sampled in mass , taking a mean straight from the isochrone points would skew the results ; the polynomial parameterization of the correction factor as a function of the mass ratio enables a more robust determination of the mean correction factor for each multiplicity scenario . 
   if the companion  to  primary mass ratio distribution was uniform across all mass ratios , then a straight mean of the correction values determined from each polynomial curve would yield the average correction factor for each multiplicity scenario . 
   however , the mass ratio distribution is likely not uniform , and we have adopted the form displayed in figure 16 of @xcite . that distribution is a nearly - flat frequency distribution across all mass ratios with a @xmath49 enhancement for nearly equal mass companion stars ( @xmath50 ) . 
   this distribution is in contrast to the gaussian distribution shown in @xcite ; however , the more recent results of @xcite incorporate more stars , a broader breadth of stellar properties , and multiple companion detection techniques . 
   the mass ratio distribution is convolved with the polynomial curves fitted for each multiplicity scenario , and a weighted mean for each multiplicity scenario was calculated for every koi . 
   for example , in the case of koi-299 ( fig . 
   [ fig - koi299 ] ) , the single star mean correction factor is 1.0 ( by definition ) . for the binary star cases , the average scenario correction factors are 1.14 
   ( planet orbits primary ) and 2.28 ( planet orbits secondary ) ; for the triple stars cases , the correction factors are 1.16 ( planet orbits primary ) , 2.75 ( planet orbits secondary ) , and 4.61 ( planet orbits tertiary ) . for koi-1085 ( fig . [ fig - koi1085 ] ) , the weighted mean correction factors are 1.18 , 1.56 , 1.24 , 1.61 , and 2.29 , respectively . 
   to turn these individual scenario correction factors into an overall single mean correction factor @xmath3 per koi , the six scenario corrections are convolved with the probability that a koi will be a single star , a binary star , or a triple star . 
   the multiplicity rate of the kepler stars is still unclear @xcite , and , indeed , there may be some contradictory evidence for the the exact value for the multiplicity rates of the koi host stars ( e.g. , * ? ? ? 
   * ; * ? ? ? 
   * ) , but the multiplicity rates appear to be near @xmath51 , similar to the general field population . in the absence of a more definitive estimate , we have chosen to utilize the multiplicity fractions from @xcite 
   : a 54% single star fraction , a 34% binary star fraction , and a 12% triple star fraction @xcite . 
   we have grouped all higher order multiples ( @xmath52 ) into the single category of `` triples '' , given the relatively rarity of the quadruple and higher order stellar systems . for the scenarios 
   where there are multiple stars in a system , we have assumed that the planets are equally likely to orbit any one of the stars ( 50% for binaries , 33.3% for triples ) . 
   the final mean correction factors @xmath3 per koi are displayed in figure  [ fig - mean - factor ] ; the median value of the correction factor and the dispersion around that median is @xmath53 . 
   this median correction factor implies that assuming a star in the koi list is single , in the absence of any ( observational ) companion vetting , yields a statistical bias on the derived planetary radii where the radii are underestimated , _ on average _ , by a factor @xmath54 , and the mass density of the planets are overestimated by a factor of @xmath55 .    from figure [ fig - mean - factor ] , it is clear that the mean correction factor @xmath3 depends upon the stellar temperature of the host star . as most of the stars in the koi list are dwarfs , the lower temperature stars are typically lower mass stars and , thus , have a smaller range of possible stellar companions . 
   thus , an average value for the correction factor 1.5 represents the sample as a whole , but a more accurate value for the correction factor can be derived for a given star , with a temperature between @xmath56k , using the fitted @xmath57-order polynomial : @xmath58 where @xmath59 .    in the absence of any specific knowledge of the stellar properties ( other than the effective temperature ) and in the absence of any radial velocity or high resolution imaging to assess the specific companion properties of a given koi , ( see section [ sub - vetting ] ) , the above parameterization ( equation [ eq - factor - unvetted ] ) can be used to derive a mean radii correction factor @xmath3 for a given star . for g - dwarfs and hotter stars , 
   the correction factor is near @xmath60 . as the stellar temperature ( mass ) of the primary decreases to the range of m - dwarfs 
   , the correction factor can be as low as @xmath61 .    ) to the quoted radii uncertainties ( @xmath62)from the cumulative koi list ( see equation [ eq - total - unc ] ) . 
   for the red histogram , it is assumed that the kois are single as is the case in the published koi list ; for the blue histogram , it is assumed that each koi has been vetted with radial velocity ( rv ) and high resolution imaging ( see section [ sub - vetting ] ) . 
   the vertical dashed lines represent the median values of the distributions : @xmath63 for the unvetted kois and @xmath64 for the vetted kois ( see section [ sub - vetting ] ) . ] 
   the mean correction factor is useful for understanding how strongly the planetary radii may be underestimated , but an additional uncertainty term derived from the mean radius correction factor is potentially more useful as it can be added in quadrature to the formal planetary radii uncertainties . 
   the formal uncertainties , presented in the koi list , are derived from the uncertainties in the transit fitting and the uncertainty in the knowledge of the stellar radius , and they are calculated assuming the kois are single stars . 
   we can estimate an additional planet radius uncertainty term based upon the mean radii correction factor as @xmath66 where @xmath27 is the observed radius of the planet . 
   adding in quadrature to the reported uncertainty , a more complete uncertainty on the planetary radius can be reported as @xmath67 where @xmath68 is the uncertainty of the planetary radius as presented in the koi list .    the distribution of the ratio of the more complete koi radius uncertainties ( @xmath69 ) to the reported koi radius uncertainties ( @xmath68 ) is shown in figure  [ fig - ratio - unc ] . including the possibility that a koi may be a multiple star increases the planetary radii uncertainties . while the distribution has a long tail dependent upon the specific system , the planetary radii uncertainties are underestimated as reported in the koi list , _ on average _ , by a factor of 1.7 . 
   the above analysis has assumed that the kois have undergone no companion vetting , as is the assumption in the current koi list . in reality , the kepler project has funded a substantial ground - based follow - up observation program which includes radial velocity vetting and high resolution imaging . in this section , 
   we explore the effectiveness of the observational vetting . 
   the observational vetting reduces the fraction of undetected companions . 
   if there is no vetting or all stars are assumed to be single , as is the case for the published koi list , then the fraction of undetected companions is 100% and the mean correction factors @xmath3 are as presented above . 
   if every stellar companion is detected and accounted for in the planetary parameter derivations , then the fraction of undetected companions is 0% , and the mean correction factors are unity . 
   reality is somewhere in between these two extremes .    to explore the effectiveness of the observational vetting on reducing the radii corrections factors ( and the associated radii uncertainties ) , we have assumed that every koi has been vetted equally , and all companions within the reach of the observations have been detected and accounted . 
   thus , the corrections factors depend only on the fraction of companions stars that remain out of the reach of vetting and undetected .    in this simulation 
   , we have assumed that all companions with orbital periods of 2 years or less and all companions with angular separations of @xmath70 or greater have been detected . 
   this , of course , will not quite be true as random orbital phase effects , inclination effects , companion mass distribution , stellar rotation effects , etc . 
   will diminish the efficiency of the observations to detect companions . 
   we recognize the simplicity of these assumptions ; however , the purpose of this section is to assess the usefulness of observational vetting on reducing the uncertainties of the planetary radii estimates , not to explore fully the sensitivities and completeness of the vetting . 
   typical follow - up observations include stellar spectroscopy , a few radial velocity measurements , and high resolution imaging . 
   the radial velocity observations usually include @xmath71 measurements over the span of @xmath72 months and are typically sufficient to identify potential stellar companions with orbital periods of @xmath73 years or less . 
   while determining full orbits and stellar masses for any stellar companions detected typically requires more intensive observing , we have estimated that 3 measurements spanning @xmath74 months is sufficient to enable the detection of an rv trend for orbital periods of @xmath75 years or less and mark the star as needing more detailed observations . 
   the amplitude of the rv signature , and hence the ability to detect companions , does depend upon the masses of the primary and companion stars ; massive stars with low mass companions will display relatively low rv signatures . 
   however , rv vetting for the kepler program has been done at a level of @xmath76 m / s , which is sufficient to detect ( at @xmath77 ) a late - type m - dwarf companion in a two - year orbit around a mid b - dwarf primary . 
   indeed , the rv vetting is made even more effective by searching for companions via spectral signatures @xcite . 
   the high resolution imaging via adaptive optics , `` lucky imaging '' , and/or speckle observations typically has resolutions of @xmath78 ( e.g. , * ? ? ? 
   * ; * ? ? ? 
   * ; * ? ? ? 
   * ; * ? ? ? 
   * ; * ? ? ? 
   * ; * ? ? ? 
   * ; * ? ? ? 
   * ; * ? ? ? 
   * ; * ? ? ? 
   * ; * ? ? ? 
   * ; * ? ? ? 
   * ; * ? ? ? 
   * ; * ? ? ? 
   based upon monte carlo simulations in which we have averaged over random orbital inclinations and eccentricities , we have calculated the fraction of time within its orbit a companion will be detectable via high resolution imaging . with typical high resolution imaging of 0.05 , we have estimated that @xmath79 of the stellar companions will be detected at one full - width half - maximum ( fwhm=0.05 ) of the image resolution and beyond and @xmath80% at @xmath81 fwhm ( 0.1 ) of the image resolution and beyond .    to determine what fraction of possible stellar companions would be detected in such a scenario , we have used the nearly log - normal orbital period distribution from @xcite . to convert the high resolution imaging limits into period - limits , we have estimated the distance to each koi by determining a distance modulus from the observed kepler magnitude and the absolute kepler magnitude associated with the fitted isochrone . 
   the median distance to the kois was found to be @xmath82 pc , corresponding to @xmath83 au for 0.1   imaging . using the isochrone stellar mass , the semi - major axis detection limits were converted to orbital period limits ( assuming circular orbits ) .    combining the 2-year radial velocity limit and the @xmath70 imaging limit 
   , we were able to estimate the fraction of undetected companions for each individual koi ( see figure  [ fig - logp ] ) . 
   the distribution of the fraction of undetected companions ranges from @xmath84% and , on average , the ground - based observations leave @xmath85 of the possible companions undetected for the kois ( see figure  [ fig - logp ] ) . 
   the mean correction factors @xmath3 are only applicable to the undetected companions . for the stars that are vetted with radial velocity and/or high resolution imaging , the intrinsic stellar companion rate for the kois of 46% @xcite 
   is reduced by the unvetted companion fraction for each koi . 
   that is , we assume that companion stars detected in the vetting have been accounted for in the planetary radii determinations , and the unvetted companion fraction is the relevant companion rate for determining the correction factors . in the koi-299 example ( fig . 
   [ fig - logp ] ) , the undetected companion rate used to calculate the mean radii correction factor is @xmath86 . this lower fraction of undetected companions in turn reduces the mean correction factors for the vetted stars which are displayed in figure  [ fig - mean - factor ] ( blue points ) .    instead of a mean correction factor of @xmath87 , 
   the average correction factor is @xmath88 if the stars are vetted with radial velocity and high resolution imaging . 
   the mean correction factor still changes as a function of the primary star effective temperature but the dependence is much more shallow with coefficients for equation  [ eq - factor - unvetted ] of @xmath89 ( see figure  [ fig - mean - factor ] ) . 
   the above analysis has concentrated on the kepler mission and the associated koi list , but the same effects will apply to all transit surveys including k2 @xcite and tess @xcite . if the planetary host stars from k2 and tess are also assumed to be single with no observational vetting , the planetary radii will be underestimated by the same amount as the kepler kois ( fig .  [ fig - mean - factor ] and eq . 
   [ eq - factor - unvetted ] ) . 
   many k2 targets and nearly all of tess targets will be stars that are typically @xmath90 magnitudes brighter than the stars observed by kepler , and therefore , k2 and tess targets will be @xmath91 times closer than the kepler targets . 
   the effectiveness of the radial velocity vetting will remain mostly unaffected by the brighter and closer stars , but the effectiveness of the high resolution imaging will be significantly enhanced . instead of probing the stars to within @xmath92au , the imaging will be able to detect companion stars within @xmath93 au of the stars .    as a result , 
   the fraction of undetected companions will decrease significantly . even for the kepler stars that undergo vetting via radial velocity and high resolution 
   imaging , @xmath94% of the companions remain undetected . 
   but for the stars that are 10 times closer that fraction decreases to @xmath95% ( see figure  [ fig - logp ] ) . 
   this has the strong benefit of greatly reducing the mean correction factors for the stars that are observed by k2 and tess and are vetted for companions with radial velocity and high resolution imaging . 
   the mean correction factor for vetted k2/tess - like stars is only @xmath96 . 
   the correction factor has a much flatter dependence on the primary star effective temperature , because the majority of the possible stellar companions are detected by the vetting . 
   the coefficients for equation [ eq - factor - unvetted ] become @xmath97 . 
   the mean radii correction factors for vetted k2/tess planetary host stars correspond to a correction to the planetary radii uncertainties of only @xmath75% , in comparison to a correction of @xmath98% if the k2/tess stars remain unvetted .    for k2 and tess , where the number of candidate planetary systems may outnumber the kois by an order of magnitude ( or more ) , single epoch high resolution imaging may prove to be the most important observational vetting performed . while the imaging will not reach the innermost stellar companions , radial velocity observations require multiple visits over a baseline comparable to the orbital periods an observer is trying to sample . 
   in contrast , the high resolution imaging requires a single visit ( or perhaps one per filter on a single night ) and will sample the majority of the expected stellar companion period distribution . 
   understanding the occurrence rates of the earth - sized planets is one of the primary goals of the kepler mission and one of the uses of the koi list @xcite . 
   it has been shown that the transition from  rocky  to  non - rocky  planets occurs near a radius of @xmath99 and the transition is very sharp @xcite . 
   however , the amplitude of the uncertainties resulting from undetected companions may be large enough to push planets across this boundary and affect our knowledge of the fraction of earth - sized planets . 
   we have explored the possible effects of undetected companions on the derived occurrence rates . 
   the planetary radii can not simply be multiplied by a mean correction factor @xmath0 , as that factor is only a measure of the statistical uncertainty of the planetary radius resulting from assuming the stars are single and only a fraction of the stars are truly multiples . 
   instead a monte carlo simulation has been performed to assign randomly the effect of unseen companions on the kois . 
   the simulation was performed 10,000 times for each koi . for each realization of the simulation 
   , we have randomly assigned the star to be single , binary , or triple star via the 54% , the 34% and the 12% fractions @xcite . 
   if the koi is assigned to be a single star , the mean correction factors for the planets in that system are unity : @xmath100 . 
   if the koi star is a multiple star system , we have randomly assigned the stellar companion masses according to the masses available from the fitted isochrones and using the mass ratio distribution of @xcite . 
   finally , the planets are randomly assigned to the primary or to the companion stars ( i.e. , 50% fractions for binary stars and 33.3% fractions for triple stars ) . 
   once the details for the system are set for a particular realization , the final correction factor for the planets are determined from the polynomial fits for the individual multiplicity scenarios ( e.g. , fig . 
   [ fig - koi299 ] and [ fig - koi1085 ] ) .    for each set of the simulations , 
   we compiled the fraction of planets within the following planet - radii bins : @xmath101 ; @xmath102 ; @xmath103 corresponding to earth - sized , super - earth / mini - neptune - sized , and neptune - to - jupiter - sized planets . 
   the raw fractions directly from the koi - list , for these three categories of planets , are 33.3% , 46.0% , and 20.7% . 
   note that these are the raw fractions and are not corrected for completeness or detectability as must be done for a true occurrence rate calculation ; these fractions are necessary for comparing how unseen companions affect the determination of fractions . 
   finally , we repeated the simulations , but using the undetected multiple star fractions after vetting with radial velocity and high resolution imaging had been performed , thus , effectively increasing the fraction of stars with correction factors of unity . 
   the distributions of the change in the fractions of planets in each planet category , compared to the raw koi fractions , are shown in figure  [ fig - occurrence - rates ] . 
   if the occurrence rates utilize the assumed - single koi list ( i.e. , unvetted ) , then the earth - sized planet fraction may be overestimated by as much as @xmath7% and the giant - planet fraction may be underestimated by as much as 30% . 
   interestingly , the fraction of super - earth / mini - neptune planets does not change substantially ; this is a result of smaller planets moving into this bin , and larger planets moving out of the bin . 
   in contrast , if all of the kois undergo vetting via radial velocity and high resolution imaging , the fractional changes to these bin fractions are much smaller : @xmath104% for the earth - sized planets and @xmath105% for the neptune / jupiter - sizes planets . 
   we present an exploration of the effect of undetected companions on the measured radii of planets in the kepler sample . 
   we find that if stars are assumed to be single ( as they are in the current kepler objects of interest list ) and no companion vetting with radial velocity and/or high resolution imaging is performed , the planetary radii are underestimated , _ on average _ , by a factor of @xmath106 , corresponding to an overestimation of the planet bulk density by a factor of @xmath107 . 
   because lower mass stars will have a smaller range of stellar companion masses than higher mass stars , the planet radius mean correction factor has been quantified as a function of stellar effective temperature .    if the kois are vetted with radial velocity observations and high resolution imaging , the planetary radius mean correction necessary to account for undetected companions is reduced significantly to a factor of @xmath108 . 
   the benefit of radial velocity and imaging vetting is even more powerful for missions like k2 and tess , where the targets are , on average , ten times closer than the kepler objects of interest . 
   with vetting , the planetary radii for k2 and tess targets will only be underestimated , on average , by 10% . 
   given the large number of candidates expected to be produced by k2 and tess , single epoch high resolution imaging may be the most effective and efficient means of reducing the mean planetary radius correction factor . 
   finally , we explored the effects of undetected companions on the occurrence rate calculations for earth - sized , super - earth / mini - neptune - sized , and neptune - sized and larger planets . 
   we find that if the kepler objects of interest are all assumed to be single ( as they currently are in the koi list ) , then the fraction of earth - sized planets may be overestimated by as much as 15 - 20% and the fraction of large planets may be underestimated by as much as 30%    the particular radial velocity observations or high resolution imaging vetting that any one koi may ( or may not ) have undergone differs from star to star . 
   companion vetting simulations presented here show that a full understanding and characterization of the planetary companions is dependent upon also understanding the presence of stellar companions , but is also dependent upon understanding the limits of those observations . for a final occurrence rate determination of earth - sized planets and , more importantly , an uncertainty on that occurrence rate , 
   the stellar companion detections ( or lack thereof ) must be taken into account . 
   the authors would like to thank ji wang , tim morton , and gerard van belle for useful discussions during the writing of this paper . 
   this research has made use of the nasa exoplanet archive , which is operated by the california institute of technology , under contract with the national aeronautics and space administration under the exoplanet exploration program . 
   portions of this work were performed at the california institute of technology under contract with the national aeronautics and space administration .
summary: we present a study on the effect of undetected stellar companions on the derived planetary radii for the kepler objects of interest ( kois ) . 
   the current production of the koi list assumes that the each koi is a single star . 
   not accounting for stellar multiplicity statistically biases the planets towards smaller radii . the bias towards smaller radii depends on the properties of the companion stars and whether the planets orbit the primary or the companion stars . 
   defining a planetary radius correction factor @xmath0 , we find that if the kois are assumed to be single , then , _ on average _ , the planetary radii may be underestimated by a factor of @xmath1 . 
   if typical radial velocity and high resolution imaging observations are performed and no companions are detected , this factor reduces to @xmath2 . 
   the correction factor @xmath3 is dependent upon the primary star properties and ranges from @xmath4 for a and f stars to @xmath5 for k and m stars . for missions like k2 and tess where the stars may be closer than the stars in the kepler target sample , observational vetting ( primary imaging ) reduces the radius correction factor to @xmath6 . 
   finally , we show that if the stellar multiplicity rates are not accounted for correctly , occurrence rate calculations for earth - sized planets may overestimate the frequency of small planets by as much as @xmath7% .