Treatment and intention-to-treat propensity score analysis to evaluate the impact of video-assisted thoracic surgery on 90-day mortality after anatomical resection for lung cancer
Abstract
GEVATS is supported by Ethicon and Spanish Society of Thoracic Surgery.
Full text
Treatment and intenon‐to‐treat propensity score analysis to evaluate the impact of video‐ assisted thoracic surgery on 90‐day mortality aer anatomical resecon for lung cancer Authors JoseLuisRecuero‐Díaz¹†,IñigoRoyo‐Crespo¹,DavidGómezde‐Antonio²,SergiCall³,Borja Aguinagalde⁴,MaríaTeresaGómez‐Hernández⁵,JorgeHernández‐Ferrández⁶,DavidSánchez‐ Lorente⁷,JulioSesma‐Romero⁸,EduardoRivo⁹,NicolásMoreno‐Mata¹⁰,andRaulEmbun¹,¹¹*† onbehalfoftheSpanishGroupofVideo‐assistedThoracicSurgery(GEVATS) †Theseauthorscontributedequallytothiswork. Affiliaons 1. DepartmentofThoracicSurgery,HospitalUniversitarioMiguelServet,HospitalClínico UniversitarioLozanoBlesa,IISAragón,Zaragoza,Spain 2. DepartmentofThoracicSurgery,HospitalUniversitarioPuertadeHierroMajadahonda, Madrid,Spain 3. DepartmentofThoracicSurgery,HospitalUniversitariMútuaTerrasa,Universitatde Barcelona,Terrasa,Barcelona,Spain 4. DepartmentofThoracicSurgery,HospitalUniversitariodeDonosa,SanSebasán‐ Donosa,Spain 5. DepartmentofThoracicSurgery,HospitalUniversitariodeSalamanca,IBSAL, UniversidaddeSalamanca,Salamanca,Spain 6. DepartmentofThoracicSurgery,HospitalUniversitarioSagratCor,Barcelona,Spain 7. DepartmentofThoracicSurgery,HospitalClínicdeBarcelona,InstutoRespiratorio, UniversityofBarcelona,Barcelona,Spain 8. DepartmentofThoracicSurgery,HospitalGeneralUniversitarioAlicante,Alicante, Spain 9. DepartmentofThoracicSurgery,HospitalUniversitarioSanagodeCompostela, SanagodeCompostela,Spain 10. DepartmentofThoracicSurgery,HospitalRamónyCajal,Madrid,Spain 11. DepartmentofSurgery,FacultyofMedicine,UniversityofZaragoza,Zaragoza,Spain *Corresponding author This is a pre-copyedited, author-produced version of an article accepted for publication in European Journal of Cardiothoracic Surgery following peer review. The version of recordJose Luis Recuero-Díaz, Iñigo Royo-Crespo, David Gómez de-Antonio, Sergi Call, Borja Aguinagalde, María Teresa Gómez-Hernández, Jorge Hernández-Ferrández, David Sánchez-Lorente, Julio Sesma-Romero, Eduardo Rivo, Nicolás Moreno-Mata, Raul Embun, on behalf of the Spanish Group of Video-assisted Thoracic Surgery (GEVATS), Treatment and intention-to-treat propensity score analysis to evaluate the impact of video-assisted thoracic surgery on 90-day mortality after anatomical resection for lung cancer, European Journal of Cardio-Thoracic Surgery, Volume 62, Issue 3, September 2022, ezac122 is available online at: https://doi.org/10.1093/ejcts/ezac122.
RaulEmbun DepartmentofThoracicSurgery HospitalUniversitarioMiguelServet PaseoIsabellaCatólica1,Zaragoza50009,Spain Tel:+34‐635492179 Email:raulem[email protected] Abstract OBJECTIVES:Theaimofthisstudywastoknowthetreatmenteffectofvideo‐assistedthoracic surgery(VATS)on90‐daymortalityaeranatomicallungreseconbasedonanaonwide cohort. METHODS:Thisisamulcentreprospecvecohortof2721anatomicalreseconsforlung cancerfromDecember2016toMarch2018.Treatmentandintenon‐to‐treat(ITT)analyses wereperformedaerinverseprobabilityscoreweighnganddifferentpropensityscore matchingalgorithms.Covariatebalancewasassessedbystandardizedmeandifferences.The esmatorsreportedweretheaveragetreatmenteffect,theaveragetreatmenteffectonthe treatedandoddsraosaercondionallogiscmodelswith95%confidenceintervals.The unconfoundednessassumponwasevaluatedbysensivityanalysisforaveragetreatment effect(c‐dependence)andaveragetreatmenteffectonthetreated(Γ). RESULTS:VATSwastheinialapproachin1911paents(70.2%),though273cases(14.3%)had tobeconvertedtothoracotomy.Ninety‐daymortalityrateswere:treatmentanalysis(VATS 1.16%vsopen3.9%,P<0.001),ITTanalysis(VATS1.78%vsopen3.36%,P=0.012).Aer inverseprobabilityscoreweighngandpropensityscorematching,inthetreatmentanalysis, VATSmeantabsoluteriskreduconsbetween2.25%and2.96%andrelaveriskreducons between65%and70%(OR=0.34,95%confidenceinterval0.15–0.79,allP‐values<0.004). However,alltheesmatorsturnedouttobenon‐significantintheITTanalyses.Ahigh sensivitytounobservableconfounderswasproved(c‐dependence0.135,Γ=1.5). CONCLUSIONS:VATScanreducetheriskof90‐daymortalityaeranatomicallungresecon. However,theimplicaonsofconversiontothoracotomy,comparingITTversustreatment analysis,andthepotenalimpactofhiddenbiasshoulddeservefurtheraenoninthefuture. Keywords:Video‐assistedthoracicsurgery,90‐Daymortality,Intenon‐to‐treatanalysis, Thoracicsurgery,Anatomicallungresecon,Lungcancer Abbreviaons ATE‐Averagetreatmenteffect ATT‐Averagetreatmenteffectonthetreated CIs‐Confidenceintervals GEVATS‐SpanishGroupofVideo‐assistedThoracicSurgery ITT‐Intenon‐to‐treat IPSW‐Inverseprobabilityscoreweighng PSM‐Propensityscorematching
SMD‐Standardizedmeandifferences VATS‐Video‐assistedthoracicsurgery Introducon Althoughtherearemulpleretrospecveseriesthathaveshownalowerrateofcomplicaons andpostoperavestayinpaentsoperatedonbyvideo‐assistedthoracicsurgery(VATS),there aremorediscrepanciesontheimpactthatthesurgicalapproachcouldhaveonpostoperave mortality.Inthisregard,arecentrandomizedclinicaltrial,designedtoevaluateshort‐termand oncologicefficacyofVATS,failedtodemonstratedifferencesinpostoperavemortality betweenVATSandthoracotomy[1]. Regardingpostoperavemortality,recentserieshaveshownthat90‐daymortalitycould double30‐dayorin‐hospitalmortalityaerlungresecon[2–4].Althoughmostseries comparingVATSwiththoracotomydonotmenon90‐daymortality,someimportantstudies havenotshownsignificantdifferencesineitherin‐hospital,30‐dayor90‐daymortality[5–7].A recentpublicaonfromtheESTSdatabasedidshowsignificantdifferencesinmortalityat discharge(VATS1%vsthoracotomy1.9%,P=0.020),notmenoningmortalityat90days[8]. Fewstudieshaveperformedanintenon‐to‐treat(ITT)analysiscomparingVATSand thoracotomyforlungresecon,andevenlesshavecomparedtheresultsobtaineddepending onthestrategyofanalysisperformed[9,10].However,sinceanon‐negligibleproporonof paentsundergoingVATSmustbefinallyconvertedtothoracotomy,ITTanalysisseemstobe themostappropriatestrategytoevaluatetreatmenteffectsrelatedtoVATSinarealscenario [11]. InviewoftheconflicngresultsandtheunderreportedITTanalysisintheliterature,the objecveofourstudywastodeterminetheimpactofthesurgicalapproachon90‐day mortality,comparingtreatmentandITTanalysisinpaentswhounderwentandanatomical lungreseconforlungcancerinthenaonalcohortbytheSpanishGroupofVideo‐assisted ThoracicSurgery(GEVATS)[12]. Materials and Methods Ethical statement Thisprojectwasapprovedbyallthelocalethicscommieesandinformedconsentwas obtainedfromtherecruitedpaentstousetheirclinicaldataforscienficpurposes(Approval byEthicsCommieeofAragonHealthResearchInstuteon20May2015PI15/0072). Data source TheGEVATSoftheSpanishSocietyofThoracicSurgeryisaprospecvevoluntarymulcentre observaonalstudywithatotalof33ThoracicSurgeryDepartmentparcipang.TheCentres werenotselectedbasedontheirexperienceinVATSoranyothercriterion.Thecohort includedpaentswhounderwentananatomicallungreseconfrom20December2016to20 March2018.TheGEVATSobjecvesweretoknowtheimpactofsurgicalapproachonshort‐ andlong‐termoutcomes.ThemethodoftheGEVATS,includingsamplesizejusficaonand theauditprocessperformed,wasrecentlypublished[12].
PaentallocaonintoVATSorthoracotomydependedonclinicalpracce,experienceand beliefsofeachoftheparcipangsurgeonsanddepartments. Inthisprospecveobservaonalcohortstudy,wetriedtospecificallyelucidatetheimpactof thesurgicalapproachon90‐daymortalityaeranatomicallungreseconforlungcancer.For thispurpose,thosepaentswithadiagnosisotherthanlungcarcinomaandthosewho underwentpneumonectomyorextendedlungreseconwereexcluded.Extendedresecon wasconsideredincaseofchestwall,diaphragmorsleeveresecon.Ourmanuscriptis reportedaccordingtotheSTROBErecommendaonsandESTSStascalPrimerforpropensity scoreanalysis. Stascal analysis Adoubletypeofanalysiswascarriedout,treatmentandITT,dependingontheneedof conversionfromVATStoopenthroughouttheprocedure.Theassociaonofthesurgical approach(openversusVATS)asthetreatmentvariable,withbaselineoncologicandsurgical variablesthatcouldinfluenceontheoutcomevariable(90‐daymortality)andtheapproachto bechosen,wasanalysedbytwo‐tailedstascalhypothesistesng,usingMann–Whitneyand Chi‐squaretestsandstandardizedmeandifferences(SMD).ThosevariableswithaP‐valueof lessthan0.2and/orstandardizeddifferencesgreaterthan0.1werethecovariatesusedto buildthepropensityscoretocorrectforseleconbias. Missingdataweredealtbycasewisedeleonanalysiswhenlessthan5%ofpaentshad incompleteregistries. Thepropensityscorewasesmatedbyalogitmodelandtheoverlapassumponwasassessed ondensityplotsfortreatmentandITTanalysis.Thetreatmenteffectswereevaluatedbasedon theinverseprobabilityscoreweighng(IPSW)andthepropensityscorematching(PSM) throughthenearest‐neighbormethodwithandwithoutreplacement,usingdifferentcalliper widths(0.035,0.05and0.1)andmatchingraos(1:1,1:2and1:3). CovariatebalancedwasassessedbySMD,beforeandaerweighngormatching.SMDless than0.1or0.05wereconsideredgoodorexcellent,respecvely,toexcluderesidualimbalance [14].BalanceofcovariateswasdisplayedondotplotsforIPSWandPSM,separately. Thetreatmenteffectswereesmatedbyweightedmeanandmatchingoutcomemodels reporngtheaveragetreatmenteffect(ATE),basedonthedifferenceinpotenaloutcome means,andtheaveragetreatmenteffectonthetreated(ATT).ATEandATTwerereportedas absoluteriskreducons(IPSWandPSM)andrelaveriskreducons(IPSW).Inaddion,in caseofPSM1:1withoutreplacementandcalliper0.035,acondionallogiscfixed‐effects regressionmodelwasconducted.Inthiscase,treatmenteffectswerereportedasoddsraos. TheimpactofsurgeonexperienceinVATSprocedures(≤50versus>50cases),surgeonseniority (residentversusfaculty<10yearsversusfaculty10–20yearsversusfaculty>20years),surgical volumeandVATSratebydepartment(discretevariables)wereusedtoadjusttheoddsraos reportedbythecondionallogiscmodelspreviouslydescribed.Tocomputesurgicalvolume andVATSratebyinstuonthroughoutthe15‐monthrecruitmentperiod,weconsideredthe reportssubmiedbytheHeadsoftheAdministraveDepartmentsfromeachInstuon. Thesereportswereusedintheauditprocesswepreviouslypublished[12].95%confidence intervals(CIs)werecalculatedfromrobuststandarderrorsandP‐valueslessthan0.05were consideredstascallysignificant.
ThecondionalparaldependencemethodproposedbyMastenetal.wasusedtoevaluate thesensivityofconclusionsabouttheATE.BoundsontheATEgivenasetofc‐dependence values(between0and1)andthebreakdownpoint(maximumvalueofthec‐dependence parameterunderwhichtheconclusionsllholds)werereported.Inaddion,theimpactof hiddenbiasontheATTesmatoraerPSM(1:1,calliper0.035,noreplacement)wasassessed withtheboundingapproachproposedbyRosenbaum.TheΓparameterandcorrespondingP‐ valueswereusedtomeasurethesensivityofATTtounobservableconfounders. TheTreatmentEffectsSuiteinStata/MP16.0andtheStatapackagesStddiff,Psmatch2, Calipmatch,TesensivityandMhboundswereusedforthestascalanalysis.TablauDesktop 2020.3.1wasusedforplotrepresentaonofthecovariatebalance. Results Atotalof3533paentswererecruited,including1917VATScases(54.3%).Aerexclusionof paentswithadiagnosisdifferenttolungcancer(448paents,12.7%),pneumonectomy(236 paents,6.7%)andextendedresecon(165,4.7%),2721paents(77%oftheenreGEVATS cohort)mettheinclusioncriteria.Thetypesofreseconsincludedwere:2444lobectomies (90%),111bilobectomies(4%)and166anatomicalsegmentectomies(6%). VATSwastheinialapproachin1911paents(70.2%),represenngthetreatmentgroupin theITTanalysis.However,273cases(14.3%)hadtobeconvertedtothoracotomyand, therefore,theVATSarminthetreatmentanalysisconsistedof1638paents(60.2%).The unadjustedanalysisshowedanimportantassociaonbetweenin‐hospitalmortalityand90‐day mortalitywiththesurgicalapproachperformed(Table1).Theoutcomevariable(90‐day mortality)wasmissingin11cases(0.4%)thatwerenotconsideredintheanalysisoftreatment effects.Thepercentagesofmissingvalueswerenegligibleforalltheconfoundersexceptfor DLCO.ThemainanalysisexcludedDLCOasconfounder,so14variableswereusedtobuildthe propensityscoreinacomplete‐caseanalysissinceonly1.8%ofthepaentshadmissingvalues insomeofthecovariatesincluded(Table2).Thepropensityscoredistribuontoprovethe overlapassumponisshowninadensityplot(Figure1).Only9paentsintheVATSgroup (0.6%)hadapropensityscorehigherthanthemaximumpropensityscorevalueintheopen group,whileonly2paentsintheopengroup(0.2%)hadapropensityscorelowerthanthe minimumvalueintheVATSgroup. AerIPSW,theraosofVATStothoracotomywere1339:1324paentsinthetreatment analysisand1347:1317paentsintheITTanalysis.Covariatebalancewasproved(Figure2). PSM1:1withoutreplacementandacalliperwidthof0.035yielded872matchedpairsinthe treatmentanalysisand705matchedpairsintheITTanalysis.Mostofthecovariatesreached anexcellentbalanceaermatching(Table3).Theproporonsofcasesmatchedto1control were54.2%(treatmentanalysis)and37.6%(ITTanalysis). AerPSMwithreplacement,2664paents(1608VATSand1056open)inthetreatment analysisand2663paents(1877VATSand786open)intheITTanalysiswerematchedto1,2 or3counterparts.Independentlyofthematchingraosandcalliperwidths,mostofthe covariatesreachedanexcellentbalance(Figure3). Basedonthetreatmentanalysis,VATSwasconsistentlyassociatedtoalower90‐daymortality rateaccordingtotheATEandATTesmators,whateverthepropensityscoretechniqueor
algorithmused.However,inthecaseofITTanalysis,theslllowermortalityaerVATSwas associatedtonon‐significantATEandATT,intermsofabsoluteandrelaveriskreducons,in allthecases(Table4). InthecondionallogiscmodelaerPSM1:1(noreplacementandcalliperwidth0.035),the resultsobtainedaertreatmentanalysis(OR=0.38;95%CI0.20–0.73;P=0.004)andITT analysis(OR=0.71;95%CI0.38–1.33;P=0.283)wereconsistentwiththeATEandATT esmators.Aeradjusngthefixed‐effectmodelsbysurgeonexperience(numberofVATS proceduresandseniority)anddepartmentexperience(surgicalvolumeandVATSrate)theodds raosweresimilartothoseinthenon‐adjustedmatchedtreatmentanalysis(OR=0.34;95%CI 0.15–0.79;P=0.012)andlower,thoughsllnon‐significant,intheITTanalysis(OR=0.41;95% CI0.13–1.25;P=0.117). Inthesensivityanalysis,theclosetozeroc‐dependencevaluesunderwhichtheconclusions sllheld(treatmentanalysis0.135andITTanalysis0.09)meantahighsensivityoftheATE esmatortohiddenbias(Figure4).Inaddion,themaximumΓvalueatwhichsignificantP‐ valueswereobtained(Γ=1.5)showedthattheATTesmaonofthisstudyissensivetobias (unobservedvariables)abletoincreasetheoddsofreceivingVATSinmorethan50%. Inthemainanalysis,wedidnotconsiderDLCObecauseofthehighrateofmissingvalues (15.2%).However,inasecondaryanalysisincludingDLCOasacovariatetobuildthepropensity score,theATEandATTesmatorsaerIPSWwereequivalenttothoseobtainedinthemain analysis.Moreover,thesensivityanalysisshowedanevenhigherpotenalimpactof unmeasuredorunobservableconfounders(Γ=1.15). Discussion OurmainfindingwasthatVATScanreduce90‐daymortalityaerananatomicallungresecon forlungcancer.However,thisbeneficialeffectdecreasedandturnedouttobenon‐significant whenanITTanalysiswasperformed.Noteworthy,thesefindingswereconsistentalongallthe propensityscorealgorithmscarriedout.Finally,thehighsensivityofouranalysistohidden biashighlightsthatthepotenalimpactofunmeasuredandunobservableconfoundersin observaonalstudiescomparingVATSandopensurgerycouldbedeterminant. Theresultsofthefewstudiescomparingtheimpactofthesurgicalapproachon90‐day mortalityleadtocontradictoryconclusions.Limitedstascalpowersecondarytosmall samplesizesandsurprisinglylowmortalityratescouldpartlyexplainnon‐conclusivefindings [1,6,7].However,whenanalysingtheresultsobtainedbylargeinstuonalcohorts,significant stascaldifferencescouldonlycorrespondtoclinicallytrivialtreatmenteffect[15]. Inthethoracicsurgeryliterature,usefulrecommendaonshavebeenmadetoimprove reporngbasedonthepropensityscoreanalysis,sincethisanalycapproachisnotas straighorwardasregressionanalysis[16].ConsideringtherecommendaonsbytheEuropean SocietyofThoracicSurgeonsinitsstascalprimerreport,wecarriedoutthe2advised methods,IPSWandPSMwithcalliper[13].Inaddion,toassesstheconsistencyofourresults, weputintopraccedifferentmatchingalgorithms(treatedtocontrolraos,calliperwidths andreplacementversusnoreplacementmatching). Theuseof2propensityscorestechniquescomparingoutcomesaerVATSoropenlung reseconisanexcepon[10].Pagèsetal.publishedtheresultsbyaFrenchnaonwidestudy
evaluangtheimpactofsurgicalapproachonshort‐andlong‐termoutcomesaerlobectomy forlungcancer.Inthislargecohort(n=24,811)recruitedthroughoutan8‐yearperiod,only 4.9%ofthecaseswereoperatedonbyVATS.Itisnoteworthythatthesignificancelevel obtainedinafewoutcomesdifferedbetweenthepropensityscorestrategiesused(IPSW versusPSM).However,regarding30‐daymortality,bothIPSWandPSMfailedtodemonstratea significantreduconaerVATS(PSMOR=0.89,95%CI0.45–1.81;IPSWOR=0.74,95%CI 0.37–1.45). TheimportanceofanITTstrategywhenevaluangtheimpactofthesurgicalapproachwas highlightedinoneofthefirstmeta‐analysespublishedcomparingVATSandthoracotomy[17]. However,theinformaonaboutconversionfromVATStoopensurgeryisnotalwaysregistered, evenincaseofnaonalandinternaonalregistries[6,8,10,18,19].Consequently,ITTanalysis isnotawidepracceinthethoracicsurgeryliteratureand,therefore,thebenefitconferredto VATSwhenevaluangtreatmenteffectscouldbemisleading.Infact,thehighertherateof conversionfromVATStoopenthehigherthepotenaldiscrepancybetweentreatmentandITT analyses.Inthisregard,despiteahigherrateofVATSinourcohort(60.2%)thantheone reportedbymanynaonwideandmulcentreregistries,theproporonofconversionweare reporng(14.3%)isnotablysuperiorcomparedtosomerecognizedseriesrangingfrom2%to 9%[11,20,21].Riskofconversioninarecentmetanalysiswas9.6%(95%CI:6.6–13.9%)[22]. Wecannotconcludeaboutthereasonsforourhigherconversionrate.However,naonwide representaveness,dataauditthatinvolvedalltheparcipangcentres,orjustamaerof disposionofsurgeons,morepronetostarngtheprocedurebyVATSinourcountry,couldbe someofthehypotheses.Inaddion,althoughanyonecouldexpectthathigherVATSratesare paralleltolowerconversionratesbecauseofexperienceacquision,thisreasoningcouldbe conflicngdependingonthescenarioandthetypeofcohortinqueson,forexamplesingle centreversusmulcentre.Inthisregard,ourrelavelyhighconversionratecouldsimply representaconsequenceofourdirectlyproporonalhighVATSrate. AlthoughsomestudiesarguethatconversionfromVATStoopendoesnotentailasurgical failure,otherserieshaveassociatedconversiontoworseoutcomesintermsofpostoperave morbidity[20,23].Ourstudywouldsupportadetrimentaleffectassociatedtoconversion. However,moredetailedandspecificanalysescomparingconvertedVATStostraight thoracotomycouldgenerateinsighulknowledgeabouttheimplicaonsofconversionitself andthemostappropriatedisposiontowardsstarngtheprocedurebyVATS. ITTanalysisisclaimedasthegoldstandardinrandomizedcontroltrialstoreflectapragmac clinicalscenario,tomaintainprognoscbalancegeneratedfromtheoriginalrandomizaonand togiveanunbiasedesmateoftreatmenteffect.However,weconsideredimportantatthis stagewhenmostofthestudiesdealingwithVATSaresllreporngtheirresultsbasedona treatmentanalysis,toincludethislessconservavestrategyofanalysissothatwecouldsll compareourfindingswiththecoetaneousliterature.Toourknowledge,thisisthefirst manuscriptthatevaluatetheimpactofthesurgicalapproachonshort‐termmortalityaer lungresecon,comparingatreatmentandITTanalysis. Limitaons Themaindrawbackinmulcentrevoluntaryregistriesisrelatedtoseleconbiasanddata quality.Althoughthedetailsofourauditwerepreviouslyreported,wecannotrejectresidual biasatthisstage[12].
Aspreviouslypublished,tocalculateoursamplesize,weconsidered2%absoluterisk differencein90‐daymortalityasrelevantinpracce(4%vs2%),andanexpected25%ofcases byVATS[12].Becauseoftheunderesmated90‐daymortalityratedifferencebetweenVATS andOpen(1.16%vs3.9%)andthehigher‐than‐expectedVATSrate(60%),thestascalpower ofour'posive'conclusioninthetreatmentanalysis(99%)ishigherthantheconvenonal 80%.However,intheITTanalysis,thelowerdifferenceinmortalitybetweenbothapproaches (1.78%vs3.36%)andthelowerproporonofcontrolpaents(29%)couldbethereasonsfora limitedpower(69%)ofour'negave'results. Missingvaluesareanothercornerstoneinobservaonalstudies.AlthoughDLCOwastheonly covariatewithamissingrateofhigherthan2%(415/2721paents15%),thisrespiratory parameterhasbeenprovedasoneofthemostdeterminantriskfactorsaerlungresecon, eveninourowncohort[24–26].Despiteequivalentresultsinaposthocanalysisbasedonthe subgroupofpaentswiththisvaluepresent,confoundingbiascannotbeexcluded. DespitethenaonalrepresentavenessoftheGEVATScohort,sinceitwouldhaveincluded 50%oftheanatomicallungreseconsperformedinSpainoverthe15‐monthrecruitment period,thehighauditedrateofVATSregistered(60%)couldcompromisetheexportaonof ourresultstoothernaonalorlocalscenarioswithanimportantdiscrepancyinVATS implementaon[12].However,theequipoiseinourstudyintermsoftheproporonof paentsineachgroupcomparedtootherseriesmakesitlesslikelythateitherouropengroup representsastronglyselectedcolleconofchallengingcasesorourVATSgroupastrongly selectedcolleconoffavourablecases[6,10,27]. PSMisconsideredtojeopardizethegeneralizaonofitsconclusionswhenmatchedsample sizeislowerthan50%oftheoriginalsample.InourstudybycarryingoutPSM1:1without replacementandcalliperwidth0.035,theproporonofmatchedsamplewas64%and52% aertreatmentandITTanalysis,respecvely.Nevertheless,theresultsobtainedaersuch PSMalgorithmswereequivalenttothosereachedaermakinguseofallthesamplewith otherpropensityscoretechniques,whichseemstopreserveourconclusionatamoregeneral seng. Finally,hiddenbiasisaubiquitousproblemineveryobservaonalstudy,evenincaseofa rigorouspropensityscoreanalysis,andsensivitytesttriestogetabeerunderstandingofthis obstacletoobtainreliableesmators.However,thistypeofanalysisisexceponalinthoracic surgery[10].ThestudypublishedbyPagèsetal.fromtheEpithordatabasereportedahigh sensivitytounobservableconfounderswhenesmangpostoperavedeath(Γ=1.6).Our resultswereverysimilar(Γ=1.5),evenaerincludingDLCOasacovariateinthepropensity scoreanalysis(Γ=1.2,thecloserto1themoresensiveouresmaonstohiddenbias). Consequently,thereliabilityofouresmaonsshouldbecauouslyadmiedunlmore robustevidence. Conclusion Havinginmindtheexclusioncriteriapreviouslyreferred(pneumonectomyandextended resecon),VATSwasconsistentlyshowntoreduce90‐daymortalityby65–70%aeran anatomicallungreseconforlungcancerintheGEVATScohort.However,drawnontheITT analysisincludingpaentsconvertedtoopensurgeryintotheVATSgroup,thisbenefit decreasedtotheextentofobtainingstascallynon‐significant,thoughpotenally
underpowered,differences.UnlupcomingevidencecomparingconvertedVATStostraight thoracotomy,aneffortshouldbemadetoopmizethebeneficialeffectofVATSonshort‐term mortalityaerananatomicalreseconforlungcancer.Intheanalysisofhiddenbias,we demonstratedthatthetreatmenteffectesmatorswerehighlysensivetounobservable confounders,whichcouldonlybeovercomeinasufficientlylargemulcentrerandomizedtrial. Acknowledgements Wethankallthelocalresearcherswhocontributedwithpaentrecruitmentanddataentry. Coordinator:RaulEmbun(UniversityHospitalMiguelServet). Scienfic commiee:DavidGómezdeAntonio(UniversityHospitalPuertadeHierro Majadahonda,Madrid);SergiCall(MútuaTerrasaUniversityHospital,UniversityofBarcelona, Terrasa,Barcelona);NicolásMoreno‐Mata(RamónyCajalUniversityHospital,Madrid); MarceloF.Jiménez(SalamancaUniversityHospital,UniversityofSalamanca,IBSAL, Salamanca);MiguelCongregado(VirgenMacarenaUniversityHospital,Seville);andSergio Bolufer‐Nadal(GeneralUniversityHospitalofAlicante). Local researchers:JoseLuisRecueroandIñigoRoyoCrespo(UniversityHospitalMiguelServet andLozanoBlesa);BorjaAguinagaldeandIkerLópezSanz(DonosaUniversityHospital,San Sebasán‐Donosa);SergioAmor‐AlonsoandFranciscoJavierMoradiellos‐Díez(Quironsalud MadridUniversityHospital,Madrid);MiguelJesúsArrarás(ValencianOncologyInstute Foundaon,Valencia);AnaIsabelBlancoOrozco(UniversityHospitalVirgendelRocío,Sevilla); MarcBoadaandDavidSánchez(HospitalClinicofBarcelona,RespiratoryInstute,Universityof Barcelona,Barcelona);AlbertoCabañeroSánchez(RamónyCajalUniversityHospital,Madrid); IsabelCalVázquezandRamónMoreno‐Balsalobre(UniversityHospitalLaPrincesa,Madrid); ÁngelCillerueloRamos(UniversityClinicalHospital,Valladolid);SilvanaCrowleyCarrasco (UniversityHospitalPuertadeHierroMajadahonda,Madrid);ElenaFernández‐Marnand FlorennoHernandoTrancho(HospitalClínicoSanCarlos,Madrid);SanagoGarcía‐Barajasand CiprianoLópezGarcía(UniversityHospitalofBadajoz,Badajoz);MariaDoloresGarcía‐Jiménez (UniversityHospitalofAlbacete,Albacete);JoseMaríaGarcía‐PrimandEduardoRivo(Sanago deCompostelaUniversityHospital,SanagodeCompostela);JoseAlbertoGarcia‐Salcedo(12 deOctubreUniversityHospital,Madrid);JuanJoséGelbenzu‐ZazpeandMaríaElenaRamírez‐ Gil(ComplejoHospitalarioofNavarre,Pamplona);CarlosFernandoGiraldo‐OspinaandRoberto MongilPoce(HospitalRegionalUniversitario,Málaga);MaríaTeresaGómezHernández (SalamancaUniversityHospital,UniversityofSalamanca,IBSAL,Salamanca);JorgeHenández andJuanJoséFiblaAlfara(SagratCorUniversityHospital,Barcelona);JenniferD.IllanaWolf (PuertadelMarHospital,Cádiz);AlbertoJaureguiAbularach(Valld'HebronUniversityHospital, Barcelona);UnaiJiménezandRafaelRojo‐Marcos(UniversityHospitalCruces,Bilbao);NéstorJ. Marnez‐Hernández(LaRiberaUniversityHospital,Alcira,Valencia);ElisabethMarnez‐Téllez andJuanCarlosTrujilloReyes(SantaCreuandSantPauHospital,AutonomousUniversityof Barcelona,Barcelona);LucíaMillaCollado(HospitalArnaudeVilanova,Lleida);SergioB. MorenoMerino(VirgenMacarenaUniversityHospital,Seville);CarmeObiols(MútuaTerrasa UniversityHospital,UniversityofBarcelona,Terrasa,Barcelona);FlorencioQuero‐Valenzuela (VirgendelasNievesHospital,Granada);RicardRamos‐Izquierdo(UniversityHospitalof Bellvitge,HospitaletdeLlobregat,Barcelona);AlbertoRodríguez‐Fuster(HospitaldelMar, HospitaldelMarMedicalResearchInstute,Barcelona);LauraSanchezMoreno(University
11. DecaluweH,PetersenRH,HansenH,PiwkowskiC,AugusnF,BrunelliAetal.;ESTS MinimallyInvasiveThoracicSurgeryInterestGroup(MITIG).Majorintraoperave complicaonsduringvideo‐assistedthoracoscopicanatomicallungresecons:an intenon‐to‐treatanalysis.EurJCardiothoracSurg2015;48:588–99. 12. EmbunR,Royo‐CrespoI,RecueroDíazJL,BoluferS,CallS,CongregadoMetal.Spanish Video‐AssistedThoracicSurgeryGroup:method,auding,andinialresultsfroma naonalprospecvecohortofpaentsreceivinganatomicallungresecons.Arch Bronconeumol2020;56:718–24. 13. BenedeoU,HeadSJ,AngeliniGD,BlackstoneEH.Stascalprimer:propensityscore matchinganditsalternaves.EurJCardiothoracSurg2018;53:1112–7. 14. AusnPC.Balancediagnoscsforcomparingthedistribuonofbaselinecovariates betweentreatmentgroupsinpropensity‐scorematchedsamples.StatMed2009 10;28:3083–107. 15. AbdelsaarZM,AllenMS,ShenKR,CassiviSD,NicholsFC,WigleDAetal.Variaonin hospitaladoponratesofvideo‐assistedthoracoscopiclobectomyforlungcancerand theeffectonoutcomes.AnnThoracSurg2017;103:454–60. 16. AusnPC.Propensity‐scorematchinginthecardiovascularsurgeryliteraturefrom 2004to2006:asystemacreviewandsuggesonsforimprovement.JThorac CardiovascSurg2007;134:1128–35. 17. YanTD,BlackD,BannonPG,McCaughanBC.Systemacreviewandmeta‐analysisof randomizedandnonrandomizedtrialsonsafetyandefficacyofvideo‐assistedthoracic surgerylobectomyforearly‐stagenon‐small‐celllungcancer.JClinOncol 2009;27:2553–62. 18. PuriV,GaissertHA,WormuthDW,GroganEL,BurfeindWR,ChangACetal.Thoracic surgeonsparcipantsperformingthoracoscopiclobectomy.AnnThoracSurg 2019;107:202–8. 19. PaulS,AltorkiNK,ShengS,LeePC,HarpoleDH,OnaisMWetal.Thoracoscopic lobectomyisassociatedwithlowermorbiditythanopenlobectomy:apropensity‐ matchedanalysisfromtheSTSdatabase.JThoracCardiovascSurg2010;139:366–78. 20. BongiolaS,GonfioA,ViggianoD,BorgianniS,PoliL,CrisciRetal.;ItalianVATS Group.RiskfactorsandimpactofconversionfromVATStoopenlobectomy:analysis fromanaonaldatabase.SurgEndosc20191;33:3953–62. 21. LaursenLØ,PetersenRH,HansenHJ,JensenTK,RavnJ,KongeL.Video‐assisted thoracoscopicsurgerylobectomyforlungcancerisassociatedwithalower30‐day morbiditycomparedwithlobectomybythoracotomy.EurJCardiothoracSurg 2016;49:870–5. 22. PowerAD,MerriRE,Abdel‐RasoulM,Moffa‐BruceSD,D'SouzaDM,KneuertzPJ. Esmangtheriskofconversionfromvideo‐assistedthoracoscopiclungsurgeryto thoracotomy—asystemacreviewandmeta‐analysis.JThoracDis2021;13:812–23.
23. FourdrainA,DeDominicisF,IquilleJ,LafieS,MerluscaG,Wie‐PfisterAetal. Intraoperaveconversionduringvideo‐assistedthoracoscopydoesnotconstutea treatmentfailure.EurJCardiothoracSurg2019;55:660–5. 24. BrunelliA,DineshP,Woodcock‐ShawJ,LilechildD,PompiliC.Ninety‐daymortality aervideo‐assistedthoracoscopiclobectomy:incidenceandriskfactors.AnnThorac Surg2017;104:1020–6. 25. GómezdeAntonioD,CrowleyCarrascoS,RomeroRománA,RoyuelaA,SánchezCalle Á,ObiolsFornellCetal.Surgicalriskfollowinganatomiclungreseconinthoracic surgery:aprediconmodelderivedfromaSpanishMulcenterDatabase.Arch Bronconeumol2021Feb24:S0300‐2896(21)00070‐3. 26. AguinagaldeB,InsausA,LopezI,SanchezL,BoluferS,EmbunR.VATSlobectomy morbidityandmortalityislowerinpaentswiththesameppoDLCO:analysisofthe databaseoftheSpanishVideo‐AssistedThoracicSurgeryGroup.ArchBronconeumol 2021;57:750–6. 27. XuJ,NiH,WuY,CaoJ,HanX,LiuLetal.Perioperavecomparisonofvideo‐assisted thoracicsurgeryandopenlobectomyforpT1‐stagenon‐smallcelllungcancerpaents inChina:amul‐centerpropensityscore‐matchedanalysis.TranslLungCancerRes 2021;10:402–14. Note:Thisistheauthor'sacceptedmanuscriptversion.Thefinalpublishedarclemaycontain addionaleditsandformangchangesmadeduringtheproduconprocess.Pleasecitethe publishedversionwhenreferencingthiswork.
Figure_1
Figure_2
Figure_3
Figure_4