法国大学讲义：从spss到stata操作研究

Ateliersméthodologiques

DépartementdeBiostatistiqueetunitéEpitraumac

SPSS-STATAQuelquesélémentsdecomparaisonM.Dramaix13novembre2007Lesfichiers–lesfenêtresSPSS:approchesimilairedecelledesautreslogicielsWINDOWS(EXCEL)Données:éditablesdirectementdanslafenêtreExportationetimportationdenombreuxtypesdefichiers:EXCEl,DBF,STATA,SAS…STATA:plusproched’EPIINFOPasdesauvegardepossibledel’outputàposterioriAl’écran,longueurdel’outputlisiblelimitéeDonnées:visiblesavecDATABROWSERDonnées:éditablesavecDATAEDITORExportation–importationpluslimité:fichiersASCII,SASXPORT13novembre2007LesfenêtresdeSPSS13novembre2007LesfenêtresdeSTATA13novembre2007SPSSDATAVIEW13novembre2007STATADATABROWSER13novembre2007Lescommandes(I)DansSPSSPlussimpleLaplupartdirectementaccessiblesparlesmenusCommandesfacilementidentifiablesCommandesnonaccessibles:utilisationdanslaprogrammation,parex:DOREPEAT13novembre2007Lescommandes(II)DansSTATA

PluscomplexeMaislogiquedansl’écriturepouroptions,sélection…Passimplederetrouvertouteslescommandesparlesmenus13novembre2007Lescommandes(III)DansSTATADanscertainscas,plusieurscommandesnécessairespourunecommandeSPSSPlusgrandenécessitédeconnaîtrelathéoriepourtrouverlescommandesadéquates

moinsderisqued’erreurdanslechoixdes

méthodesPasd’abréviations,

fautes«orthographe»,neconfondpasmajusculesetminuscules13novembre2007Exemple1:testt-SPSST-TESTGROUPS=PPN(01)/MISSING=ANALYSIS/VARIABLES=PoidM/CRITERIA=CI(.95).13novembre2007Exemple1:testt-STATArobvarPoidM,by(PPN)ttestPoidM,by(PPN)ttestPoidM,by(PPN)unequal13novembre2007Exemple2–Régression(I)13novembre2007Exemple2–Régression

STATA(suite)13novembre2007Exemple2–Régression(2)REGRESSION/DESCRIPTIVESMEANSTDDEVCORRSIGN/MISSINGLISTWISE/STATISTICSCOEFFOUTSRANOVACOLLINTOL/CRITERIA=PIN(.05)POUT(.10)/NOORIGIN/DEPENDENTPnaiss/METHOD=ENTERPOIDMKAGEMTabacProbur/SCATTERPLOT=(*ZRESID,*ZPRED)/RESIDUALSHIST(ZRESID)NORM(ZRESID).regressPnaissPOIDMKAGEMTabacProbur,betapredictfitpredictsdres,rstandardqnormsdrestwoway(scattersdresfit)estatvifCommandeSPSS(paste)CommandesSTATA13novembre2007Exemple3–procédurepasàpasstepwise,pe(0.05):regressPOIENFAGEMETCIVPERBGBAUDEBMIxi:stepwise,pe(0.05):regressPOIENFAGEMi.ETUMETCIVPERBGBAUDEBMIi.grpariteREGRESSION/MISSINGLISTWISE/STATISTICSCOEFFOUTSRANOVA/CRITERIA=PIN(.05)POUT(.10)/NOORIGIN/DEPENDENTPOIENF/METHOD=STEPWISEAGEMETCIVPERBGBAUDEBMI.CommandeSPSSCommandeSTATACommandeSTATAavecdécompositionenindicateursstepwise,pe(0.05):regressPOIENFAGEMETCIVPERBGBAUDEBMICommandeSTATAstepwise,pe(0.05):regressPOIENFAGEMETCIVPERBGBAUDEBMIxi:stepwise,pe(0.05):regressPOIENFAGEMi.ETUMETCIVPERBGBAUDEBMIi.grpariteCommandeSTATAavecdécompositionenindicateursxi:stepwise,pe(0.05):regressPOIENFAGEMi.ETUMETCIVPERBGBAUDEBMIi.grpariteCommandeSTATAavecdécompositionenindicateursxi:stepwise,pe(0.05):regressPOIENFAGEMi.ETUMETCIVPERBGBAUDEBMIi.grparite13novembre2007LesméthodesdisponiblesPlusétenduesdansSTATAPlusdeméthodesappliquéesenépidémiologieEx:analysedestables,tailleéchantillon…Sortiesdifférentesaveclesdeuxlogiciels:plusoumoinsdestatistiquesavecl’unoul’autreEx:mesuresassociationdanslestables,ANOVA,procédurespasàpas…plusdétaillédansSPSSEx:testd’égalitédesvariances,mesuresépidémiologiquesdanslestables,courbesROC…plusdétaillédansSTATA13novembre2007Méthodes13novembre2007Sortietestt-SPSS13novembre2007Sortietestt–STATA(I).robvarPoidM,by(PPN)|SummaryofPoidsmèreavant|(pounds)PPN|MeanStd.Dev.Freq.------------+------------------------------------>=2500g|133.331.724016130<2500g|122.1355926.55927559------------+------------------------------------Total|129.8148130.57938189W0=1.12854604df(1,187)Pr>F=0.28945567W50=0.50995689df(1,187)Pr>F=0.47604792W10=0.74660652df(1,187)Pr>F=0.3886597213novembre2007Sortietestt–STATA(II).ttestPoidM,by(PPN)Two-samplettestwithequalvariances------------------------------------------------------------------------------Group|ObsMeanStd.Err.Std.Dev.[95%Conf.Interval]---------+-------------------------------------------------------------------->=2500g|130133.32.7823831.72402127.795138.805<2500g|59122.13563.45772326.55928115.2142129.057---------+--------------------------------------------------------------------combined|189129.81482.22432330.57938125.427134.2027---------+--------------------------------------------------------------------diff|11.164414.7432971.80715720.52166------------------------------------------------------------------------------diff=mean(>=2500g)-mean(<2500g)t=2.3537Ho:diff=0degreesoffreedom=187

Ha:diff<0Ha:diff!=0Ha:diff>0

Pr(T<t)=0.9902Pr(|T|>|t|)=0.0196Pr(T>t)=0.009813novembre2007SortiesTables-SPSS13novembre2007SortiesTables–STATA(I).tabulateTabacPPN,chi2columnexactexpectedlrchi2V+--------------------+|Key||--------------------||frequency||expectedfrequency||columnpercentage|+--------------------+Tabagisme|durant|PPNgrossesse|>=2500g<2500g|Total-----------+----------------------+----------Non|8629|115|79.135.9|115.0|66.1549.15|60.85-----------+----------------------+----------Oui|4430|74|50.923.1|74.0|33.8550.85|39.15-----------+----------------------+----------Total|13059|189|130.059.0|189.0|100.00100.00|100.00Pearsonchi2(1)=4.9237Pr=0.026

likelihood-ratiochi2(1)=4.8674Pr=0.027Cramér'sV=0.1614Fisher'sexact=0.0361-sidedFisher'sexact=0.02013novembre2007.ccPPNTabacProportion|ExposedUnexposed|TotalExposed-----------------+------------------------+------------------------Cases|3029|590.5085Controls|4486|1300.3385-----------------+------------------------+------------------------Total|74115|1890.3915|||Pointestimate|[95%Conf.Interval]|------------------------+------------------------Oddsratio|2.021944|1.0290923.965864(exact)

Attr.frac.ex.|.5054264|.0282695.7478481(exact)Attr.frac.pop|.2569965|

+-------------------------------------------------chi2(1)=4.92Pr>chi2=0.0265SortiesTables–STATA(II)

Tablespour“épidémiologistes”13novembre2007SortiesAnalyseStratifiée–SPSS(I)13novembre2007SortiesAnalyseStratifiée–SPSS(II)13novembre2007SortiesAnalysesStratifiée-STATA.ccPPNTabac,by(Ethnie)pool

Ethnie|OR[95%Conf.Interval]M-HWeight-----------------+-------------------------------------------------Blanc|5.7575761.65757425.13881.375(exact)Noir|3.3.486538523.45437.7692308(exact)Autres|1.25.2734955.2782292.089552(exact)-----------------+-------------------------------------------------Crude|2.0219441.0290923.965