connectivity网易轻博客

1ElectronicSupplementaryMaterialGeneticAssessmentofEnvironmentalFeaturesthatInfluenceDeerDispersal:ImplicationsforPrion-InfectedPopulationsAmyC.
Kelly,NohraE.
Mateus-Pinilla,WilliamBrown,MarilynO.
Ruiz,MarlisR.
Douglas,MichaelE.
Douglas,PaulShelton,TomBeissel,JanNovakofskiMicrosatelliteMarkersThefollowingmicrosatelliteswereemployedinthisstudy:BM1225,BM4107,CSN3,(Bishopetal.
1994),IGF-1(Kirkpatrick1992),OBCAM(Friesetal.
1993),OarFcb304(Buchananetal.
1993),RT20,RT23,RT27(Wilsonetal.
1997)andSrcrsp-10(Bhebheetal.
1994).
Welabeledforwardprimerswithfluorescentdyes(NED,HEX,FAM)andseparatedmicrosatellitefragmentsonanABI3730XLcapillarysequencer(AppliedBiosystems,FosterCity,CA).
WevisualizedmicrosatellitegenotypeswithGeneMapper(v.
4.
0;AppliedBiosystems,FosterCity,CA).
WeusedMicro-checker(v.
2.
2.
3;VanOosterhoutetal.
2004)toevaluategenotypingerrorsusingexpectedallelefrequenciesderivedunderHardy-Weinbergequilibrium(HWE).
FSTSurfaceProjectionWeusedtheSingleSpeciesGeneticDivergenceoptionwithintheGeneticLandscapesGIS(GeographicInformationSystem)ToolboxtoprojectasurfacefrompairwiseFSTvaluescalculatedbetweenall31studysites.
TheprogramfirstassociatedpairwiseFST2valueswithmidpointsbetweenallstudysitesandanetworkofnearestneighbors.
Spatialinterpolationwasthenperformedusinganinversedistanceweightedinterpolationalgorithmtoestimategeneticdistancesalongagridoverlaidonthestudyarea.
GeneticdistancesforallpointsacrossthegridwereinterpolatedsuchthatmidpointFSTvaluesthatwerespatiallycloserinfluencedtheestimatemoresothanthosethatweredistant.
Moredetailsontheinterpolationprocedurearedescribedinhttp://www.
werc.
usgs.
gov/productdetails.
aspxid=4017.
FRAGSTATSmetricsTheConnectanceIndex(CONNECT)measuresfunctionalconnectivity,meaningthatgridcellsinthedatathatdepictthetargetvariablearenotliterallyadjacent,buttheyareconsideredadjacent(orconnected)withinagiventhresholddistance.
Inthiscase,adjacencywasdefinedascellswithin100mofeachother.
Theuser-defined100mthresholdwasusedtoaccountforpotentialimprecisionofdataclassificationsatfinespatialresolutionsandtoprovideamorerealistic(i.
e.
,functional)depictionofhowdeermightinteractwiththelandscape.
Themetricitselfisapercentage,witharangeof0to100.
Morespecifically,itmeasuresthepercentageoftargetvariableadjacencies(connectionsorjoins)relativetoallpossibleadjacencies.
FormoreinformationontheConnectanceIndexsee:http://www.
umass.
edu/landeco/research/fragstats/documents/Metrics/Connectivity%20Metrics/Metrics/C122%20-%20CONNECT.
htm3ThePatchCohesionIndex(COHESION)isasecondmeasureofconnectivityofalandscapevariable.
Thismetrictakesintoaccountphysicaladjacency(withoutathreshold)incombinationwiththesizeandshapeofthepatches.
Takingforestasanexample,ahigherCOHESIONvaluewouldoccurinalandscapewithlargerandcompactpatchescomparedtoonewithsmallorconvolutedpatches.
FormoreinformationonthePatchCohesionIndexsee:http://www.
umass.
edu/landeco/research/fragstats/documents/Metrics/Connectivity%20Metrics/Metrics/C121%20-%20COHESION.
htmTheClumpinessIndex(CLUMPY)isametricindicatinghowcontiguousordispersedaretheadjacentpatchesofalandscapevariable.
AhighervalueofCLUMPYwouldoccurifseveralpatcheswerelocatedclosetogetherratherthanbeingmoreuniformlydistributed.
FormoreinformationontheClumpinessIndex(CLUMPY)seehttp://www.
umass.
edu/landeco/research/fragstats/documents/Metrics/Contagion%20-%20Interspersion%20Metrics/Metrics/C115%20-%20CLUMPY.
htmThePerimeter-AreaFractalDimension(PAFRAC)isashapemetricdeterminedacrossarangeofspatialscales.
PARFRACislowforpatcheswithsimpleperimetersandincreasesforpatchshapeswithhighlyconvolutedperimeters.
FormoreinformationonthePerimeter-AreaFractalDimensionIndex(PAFRAC),seehttp://www.
umass.
edu/landeco/research/fragstats/documents/Metrics/Shape%20Metrics/Metrics/C23%20-%20PAFRAC.
htm.
Multivariatelinearregressionanalysis4DescriptionandsourceoflandscapevariablesincludedinmultivariateregressionanalysisarelistedinTableS1.
Topreventoverlyinfluentialobservationsfrombiasingourmodels,weusedleveragescores,Cook'sDvalues,andstandardizedinfluencevaluestoidentifyoutliers(Kieetal.
2002;ChatterjeeandHadi2009;Anlaufetal.
2011).
Leveragescoresidentifyobservationsthatresultinlargechangesinregressionlinefitupontheirdeletion.
Wecalculatedleverage(pi)accordingtoChatterjeeandHadi(1986)andconsideredobservationsoverlyinfluentialwhenpi>2p/N(p=numberofindependentvariablesinthemodel;N=numberofobservations).
Cook'sDvalueswerecalculatedaccordingtoCook(1977)andcomparedtoanFdistributionwithα=0.
05and(N-p)degreesoffreedom.
AllCook'sDvalues>thecriticalFvaluewereconsideredoverlyinfluentialandremovedfromthemodel(Cook1977).
LeveragescoresandCook'sDallowedustodeterminetheeffectsofoutliersontheoverallmodel,butstandardizedinfluencevalues(DFFITS)allowedustoexaminetheinfluenceofeachobservationonitspredictedvalue.
WecalculatedDFFITSaccordingtoChatterjeeandHadi(1986)andeliminatedobservationsyieldingvalues>2)/(Np(ChatterjeeandHadi1986).
Usingthesethreecriteria,weidentifiedthirteenobservationsoutof465(2.
8%)thatwereoutliersandafterstringentlyevaluatingtheirbasis(Motulsky2010),weomittedthemduringfurtheranalyses.
Themajorityoftheoutliersremoved(7/13)involvedstudysitesthathadrelativelylowsamplesizes.
Threeofthirteenoutliersinvolvedpairwisecomparisonswithstudysite27,thoughtheremainingtenoutliersappearedtoinvolvestudysitesthatwererandomlydistributedgeographically.
AsingleoutlierhadthehighestFSTvalueobserved,thoughtheremainingoutliersdidnotexhibitunusuallyhighorlowFSTvaluesascomparedtotherestofthe5dataset.
WecomparedvaluesofdependentvariablesofoutlierstovaluesfortherestofthedatabyexaminingboxplotsandplottingdependentvariablesagainstFSTvalues(datanotshown).
Trendsinthedistributionofvaluesfordependentvariablewerenotapparentinoutliersascomparedtotherestofthedata.
Whentwoormorelandscapevariableswerehighlycorrelated(Pearson'srP>0.
7),thepredictorwiththelowestpartialcorrelationinthefullmodelwasremoved.
RemovinglandscapevariableswithrP>0.
7(n=7)resultedinagenerallackofcollinearityamongpredictorsasdeterminedbyvarianceinflationfactors.
CorrelatedpredictorsthatwereremovedfromthemodelarelistedinTableS2.
Weusedvarianceinflationfactors(VIF)toevaluatetheincreaseinvarianceforestimatedregressioncoefficientsresultingfromcollinearpredictors,withVIF>10indicativeofhighmulticollinearity(Kutneretal.
2004).
Afterremovinghighlycorrelatedvariables,wecalculatedvarianceinflationfactorsforindependentvariablesandfoundthatthevarianceofestimatedregressioncoefficientswasnotsubstantiallyincreasedbycollinearpredictorsasVIFvaluesforallpredictorswere0.
7thatweresubsequentlyremovedfromthemodel.
VariableCorrelateDirectionofCorrelationVariableRemoved*%GrasslandSlope+%GrasslandForestCONNECTDevelopedCONNECT+ForestCONNECT%GrasslandGrasslandCONNECT-%GrasslandForestCONNECTGrasslandCONNECT+ForestCONNECTForestCONNECTWaterCONNECT+ForestCONNECTAgricultureCLUMPY%Agriculture-AgricultureCLUMPY%RiparianSlope+Slope%GrasslandForestCLUMPY-%GrasslandForestCONNECTDistance-ForestCONNECTSlopeGrasslandCOHESION+SlopeDevelopedCONNECTGrasslandCONNECT+DevelopedCONNECTGrasslandPAFRACSlope+SlopeDevelopedCONNECTWaterCONNECT+DevelopedCONNECT%GrasslandAgricultureCLUMPY-%Grassland%GrasslandAgriculturePAFRAC+%GrasslandDistanceDevelopedCONNECT-DevelopedCONNECTForestCONNECTWaterCONNECT+ForestCONNECT%AgricultureAgricultureCOHESION+AgricultureCOHESIONWaterCOHESIONWaterCLUMPY+WaterCLUMPY*thepredictorwiththelowestpartialcorrelationinthefullmodelwasremoved.
10TableS3.
Percentsignificant(P<0.
05)localr,rangeoflocalr,andmeanlocalrforfive,15and25nearestneighborsingroupsofwhite-taileddeerinnorthernIllinois(NIL),DuPageCounty(DuP),andWisconsin(WI).
GroupNumberofNearestNeighbors51525%P<0.
051MaxrMeanr%P<0.
051MaxrMeanr%P<0.
051MaxrMeanrAdultMales5.
70.
160.
134.
40.
110.
087.
90.
080.
06MaleYearlings7.
00.
280.
1711.
60.
180.
0914.
10.
120.
06MaleFawns9.
30.
190.
1511.
30.
090.
078.
20.
060.
05AdultMalesandFemaleYearlings6.
40.
270.
147.
60.
120.
088.
10.
090.
06AdultFemales14.
70.
320.
1618.
80.
240.
0920.
50.
150.
07FemaleYearlings5.
70.
160.
124.
80.
110.
074.
80.
070.
05FemaleFawns17.
10.
240.
1415.
20.
130.
0919.
50.
090.
06AdultFemalesandFawns16.
00.
310.
1622.
80.
230.
1024.
50.
190.
081NumberofautocorrelationcoefficientsthatweresignificantatP<0.
05dividedbythetotalnumberautocorrelationcoefficientscalculatedforeachgroup*100.
Includingonlysignificantlocalrvalues.