#ReadMe Dryad File ##Data from: Intraspecific variation of a dominant grass and local adaptation in reciprocal garden communities along a US Great PlainsŐ precipitation gradient: Implications for grassland restoration with climate change 1. Contains 4 data files. Cover data 032312.csv date=date of measurement GDD =growing degree days as calcuated in the paper %GDD=percentage of GDD elapsed Site=4 reciprocal garden sites, Colby KS , Hays KS , Manhattan KS, Carbondale IL Block within site, 4 blocks per ecotype at each site Ecotype= Central KS (CKS), Eastern KS (EKS), and Southern Illinois (SIL) Time=2 time periods per season bbs, grass, forb, bare=5 of cover occupied by big bluestem (bbs), grass (other), forb, and bare ground, n=4 1m2 quadrats per plot, based on each quadrat having 81 wire intersections at which was determined cover MAP =Mean annual precipitation (followed by the year) MSP=Mean seasonal precipitation (followed by the year) GDD=cumulative growing degree days (followed by the year) PET=potential evaporation Aridity= MAP/MAT Gas Exchange data 040412.csv Broken up into two files-Time 1 and Time 2 date=date of measurement GDD =growing degree days as calcuated in the paper %GDD=percentage of GDD elapsed Site=4 reciprocal garden sites, Colby KS , Hays KS , Manhattan KS, Carbondale IL Block within site, 4 blocks per ecotype at each site Ecotype= Central KS (CKS), Eastern KS (EKS), and Southern Illinois (SIL) Photo=instantantous photosynthetic rate as described in paper Cond= Conductance as described in paper Ci=internal CO2 as described in paper Trmmol=transpiration rate as described in paper intr. WU=intrinsic water use efficiency as described in paper SPAD Seeded plots Data 031013.csv Site=4 reciprocal garden sites, Colby KS , Hays KS , Manhattan KS, Carbondale IL %GDD=percentage of GDD elapsed Ecotype= Central KS (CKS), Eastern KS (EKS), and Southern Illinois (SIL) Block within site, 4 blocks per ecotype at each site SPAD=spad reading based on the average of two measurements per leaf Tiller count data 2009 050411 revised.csv Site=4 reciprocal garden sites, Colby KS , Hays KS , Manhattan KS, Carbondale IL Block within site, 4 blocks per ecotype at each site Ecotype= Central KS (CKS), Eastern KS (EKS), and Southern Illinois (SIL) tiller counts=number of tlllers in 625cm2, 4 quadrats per plot, 4 replicate plots per ecotype 2. Below contains the SAS code used for the statistical analyses of each data file. Establishment: Tiller counts: FILE NAME Tiller count 2009 050411 revised.csv Proc IMPORT Datafile="C:\...\Tiller count data 2009 050411 revised.xls" out=Tiller2009 DBMS=Excel Replace; Sheet = 'DataSAS'; Getnames = Yes; Run; proc glimmix data=Tiller2009 order=data;*method=laplace nobound ; class Ecotype Planting_Site Block; model Tiller_counts = Planting_Site|Ecotype/ dist=negbin ddfm=kr; random intercept Ecotype / subject=Block(Planting_Site); nloptions tech=nrridg; lsmeans Planting_Site Ecotype/ ilink pdiff adjust=tukey plots=meanplot(ilink join); lsmeans Planting_Site*Ecotype/ ilink slicediff=Planting_Site adjust=bon plots=meanplot(ilink sliceby=Ecotype join); lsmeans Planting_Site*Ecotype/ ilink slicediff=Ecotype adjust=bon plots=meanplot(ilink sliceby=Ecotype join); run; Cover: FILE NAME: Cover data 032312.csv Proc IMPORT Datafile="C:\...\Cover data 032312.xls" out=CoverT1 DBMS=Excel Replace; Sheet = 'T1-GDD SAS'; Getnames = Yes; Run; Proc IMPORT Datafile="C:\Users\Nora\Documents\Consulting\Johnson Tetreault 2011\Cover data 032312.xls" out=CoverT2 DBMS=Excel Replace; Sheet = 'T2-GDD SAS'; Getnames = Yes; Run; Data Cover; set CoverT1 CoverT2; output; run; proc glimmix data=Cover order=data plots=studentpanel ;*nobound ; class Site Block Ecotype Time; model bbs = Site|Ecotype|_GDD/dist=Normal ddfm=kr; random intercept Ecotype/ subject=Block(Site) ; lsmeans Site*Ecotype /cl at _GDD=0.20 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ) ; lsmeans Site*Ecotype /cl at _GDD=0.35 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ) ; lsmeans Site*Ecotype /cl at _GDD=0.50 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ) ; lsmeans Site*Ecotype /cl at _GDD=0.65 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ); lsmeans Site*Ecotype /cl at _GDD=0.75 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ); lsmeans Ecotype/ cl at _GDD=0.35 adjust=bon; lsmeans Ecotype/ cl at _GDD=0.50 adjust=bon; lsmeans Ecotype/ cl at _GDD=0.65 adjust=bon; ods output lsmeans=lsmeans Slicediffs=Slicediffs; *output out=diagnostics student=student predicted=predicted; nloptions tech=nrridg; run; SPAD: FILE NAME SPAD seeded plots Data 031013.csv Proc IMPORT Datafile="C:\...\SPAD Seeded plots Data 031013.xls" out=SPAD DBMS=Excel Replace; Sheet = 'Sheet1'; Getnames = Yes; Run; proc glimmix data=SPAD order=data plots=studentpanel nobound ; GDD2 = Percent_GDD * Percent_GDD; class Site Block Ecotype; model SPAD = Site|Ecotype|Percent_GDD GDD2 GDD2*Site GDD2*Ecotype GDD2*Site*Ecotype/dist=Normal solution;*ddfm=kr; random intercept Ecotype / subject=block(Site) ; *lsmeans Site*Ecotype /at Percent_GDD=0.30 at GDD2=0.09 slicediff=Site adjust=bon ;*plot=meanplot(join sliceby=Ecotype ) ; *lsmeans Site*Ecotype /at Percent_GDD=0.40 at GDD2=0.16 slicediff=Site adjust=bon ;*plot=meanplot(join sliceby=Ecotype ) ; *lsmeans Site*Ecotype /at Percent_GDD=0.50 at GDD2=0.25 slicediff=Site adjust=bon ;*plot=meanplot(join sliceby=Ecotype ) ; *lsmeans Site*Ecotype /at Percent_GDD=0.60 at GDD2=0.36 slicediff=Site adjust=bon ; *plot=meanplot(join sliceby=Ecotype ); *lsmeans Site*Ecotype /at Percent_GDD=0.70 at GDD2=0.49 slicediff=Site adjust=bon ; * plot=meanplot(join sliceby=Ecotype ); *lsmeans Site*Ecotype /at Percent_GDD=0.80 at GDD2=0.64 slicediff=Site adjust=bon ;*plot=meanplot(join sliceby=Ecotype ) ; *lsmeans Site*Ecotype /at Percent_GDD=0.90 at GDD2=0.81 slicediff=Site adjust=bon ;*plot=meanplot(join sliceby=Ecotype ) ; lsmeans Site/ at Percent_GDD=0.50 at GDD2=0.25 adjust=bon; lsmeans Site/ at Percent_GDD=0.60 at GDD2=0.36 adjust=bon; lsmeans Site/ at Percent_GDD=0.70 at GDD2=0.49 adjust=bon; ods output lsmeans=lsmeans diffs=diffs; nloptions tech=nrridg; run; ALL OF THE FOLLOWING ARE IN THE FILE NAME: Gas Exchange data 040412.csv Photosynthesis: Proc IMPORT Datafile="C:\...\Gas Exchange data 040412.xls" out=Gas DBMS=Excel Replace; Sheet = 'Sheet1'; Getnames = Yes; Run; proc glimmix data=Gas order=data plots=studentpanel;*method=mspl;*nobound ; GDD2 = _GDD *_GDD; class Site Block Ecotype ; model Photo = Site|Ecotype|_GDD GDD2 GDD2*Site GDD2*Ecotype GDD2*Site*Ecotype/dist=Normal ;*ddfm=kr; *; random intercept Ecotype/ subject=Block(Site) ; *lsmeans Site*Ecotype /cl at _GDD=0.30 at GDD2=0.09 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ) ; *lsmeans Site*Ecotype /cl at _GDD=0.40 at GDD2=0.16 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ) ; *lsmeans Site*Ecotype /cl at _GDD=0.50 at GDD2= 0.25 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ); *lsmeans Site*Ecotype /cl at _GDD=0.60 at GDD2=0.36 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ); *lsmeans Site*Ecotype /cl at _GDD=0.70 at GDD2=0.49 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ); *lsmeans Site*Ecotype /cl at _GDD=0.80 at GDD2=0.64 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ); lsmeans Site/cl at _GDD=0.40 at GDD2=0.16 pdiff adjust=bon ; lsmeans Site/cl at _GDD=0.50 at GDD2=0.25 pdiff adjust=bon ; lsmeans Site/cl at _GDD=0.60 at GDD2=0.36 pdiff adjust=bon ; lsmeans Site/cl at _GDD=0.70 at GDD2=0.49 pdiff adjust=bon ; ods output lsmeans=lsmeans Slicediffs=Slicediffs; *output out=diagnostics student=student predicted=predicted; nloptions tech=nrridg; run; Water Use: Proc IMPORT Datafile="C:\...\Gas Exchange data 040412.xls" out=Gas DBMS=Excel Replace; Sheet = 'Sheet1'; Getnames = Yes; Run; proc glimmix data=Gas order=data plots=studentpanel method=mspl;*nobound ; GDD2 = _GDD *_GDD; class Site Block Ecotype ; model intr__WUE = Site|Ecotype|_GDD GDD2 GDD2*Site GDD2*Ecotype GDD2*Site*Ecotype/dist=Normal ddfm=kr; *GDD2 GDD2*Site GDD2*Ecotype GDD2*Site*Ecotype; random intercept Ecotype/ subject=Block(Site) ; random _residual_ / group=site; /*lsmeans Site*Ecotype / at _GDD=0.30 at GDD2=0.09 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ) ; lsmeans Site*Ecotype / at _GDD=0.40 at GDD2=0.16 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ) ; lsmeans Site*Ecotype / at _GDD=0.55 at GDD2=0.3025 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ); lsmeans Site*Ecotype / at _GDD=0.70 at GDD2=0.49 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ); lsmeans Site*Ecotype / at _GDD=0.80 at GDD2=0.64 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ); */ lsmeans Site/ at _GDD=0.40 at GDD2=0.16 pdiff adjust=bon plot=meanplot(join ); lsmeans Site/ at _GDD=0.50 at GDD2=0.25 pdiff adjust=bon plot=meanplot(join ); lsmeans Site/ at _GDD=0.60 at GDD2=0.36 pdiff adjust=bon plot=meanplot(join ); lsmeans Site/ at _GDD=0.70 at GDD2=0.49 pdiff adjust=bon plot=meanplot(join ); ods output lsmeans=lsmeans Slicediffs=Slicediffs; *output out=diagnostics student=student predicted=predicted; nloptions tech=nrridg; run; Transpiration: Proc IMPORT Datafile="C:\...\Gas Exchange data 040412.xls" out=Gas DBMS=Excel Replace; Sheet = 'Sheet1'; Getnames = Yes; Run; proc glimmix data=Gas order=data plots=studentpanel method=mspl;*nobound ; GDD2 = _GDD *_GDD; class Site Block Ecotype ; model Trmmol = Site|Ecotype|_GDD GDD2 GDD2*Site GDD2*Ecotype GDD2*Site*Ecotype/dist=Normal ddfm=kr; *GDD2 GDD2*Site GDD2*Ecotype GDD2*Site*Ecotype; random intercept Ecotype/ subject=Block(Site) ; /*lsmeans Site*Ecotype / at _GDD=0.30 at GDD2=0.09 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ) ; lsmeans Site*Ecotype / at _GDD=0.40 at GDD2=0.16 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ) ; lsmeans Site*Ecotype / at _GDD=0.55 at GDD2=0.3025 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ); lsmeans Site*Ecotype / at _GDD=0.70 at GDD2=0.49 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ); lsmeans Site*Ecotype / at _GDD=0.80 at GDD2=0.64 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ); */ lsmeans Site/ at _GDD=0.40 at GDD2=0.16 pdiff adjust=bon plot=meanplot(join ); lsmeans Site/ at _GDD=0.50 at GDD2=0.25 pdiff adjust=bon plot=meanplot(join ); lsmeans Site/ at _GDD=0.60 at GDD2=0.36 pdiff adjust=bon plot=meanplot(join ); lsmeans Site/ at _GDD=0.70 at GDD2=0.49 pdiff adjust=bon plot=meanplot(join ); ods output lsmeans=lsmeans Slicediffs=Slicediffs; *output out=diagnostics student=student predicted=predicted; nloptions tech=nrridg; run; Internal CO2: Proc IMPORT Datafile="C:\...\Gas Exchange data 040412.xls" out=Gas DBMS=Excel Replace; Sheet = 'Sheet1'; Getnames = Yes; Run; proc glimmix data=Gas order=data plots=studentpanel ;*method=mspl;*nobound ; GDD2 = _GDD *_GDD; class Site Block Ecotype ; model Ci = Site|Ecotype|_GDD GDD2 GDD2*Site GDD2*Ecotype GDD2*Site*Ecotype/dist=Normal ddfm=kr; *GDD2 GDD2*Site GDD2*Ecotype GDD2*Site*Ecotype; random /*intercept*/ Ecotype/ subject=Block(Site) ; random _residual_/group=Site; /*lsmeans Site*Ecotype / cl at _GDD=0.30 at GDD2=0.09 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ) ; lsmeans Site*Ecotype / cl at _GDD=0.40 at GDD2=0.16 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ) ; lsmeans Site*Ecotype / cl at _GDD=0.55 at GDD2=0.3025 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ); lsmeans Site*Ecotype / cl at _GDD=0.70 at GDD2=0.49 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ); lsmeans Site*Ecotype / cl at _GDD=0.80 at GDD2=0.64 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ); */ lsmeans Site/ cl at _GDD=0.40 at GDD2=0.16 pdiff adjust=bon plot=meanplot(join ); lsmeans Site/ cl at _GDD=0.50 at GDD2=0.25 pdiff adjust=bon plot=meanplot(join ); lsmeans Site/cl at _GDD=0.60 at GDD2=0.36 pdiff adjust=bon plot=meanplot(join ); lsmeans Site/ cl at _GDD=0.70 at GDD2=0.49 pdiff adjust=bon plot=meanplot(join ); ods output lsmeans=lsmeans Slicediffs=Slicediffs; *output out=diagnostics student=student predicted=predicted; nloptions tech=nrridg; run; Conductance: Proc IMPORT Datafile="C:\...\Gas Exchange data 040412.xls" out=Gas DBMS=Excel Replace; Sheet = 'Sheet1'; Getnames = Yes; Run; proc glimmix data=Gas order=data plots=studentpanel ;*method=mspl;*nobound ; GDD2 = _GDD *_GDD; class Site Block Ecotype ; model Cond = Site|Ecotype|_GDD GDD2 GDD2*Site GDD2*Ecotype GDD2*Site*Ecotype/dist=Normal ddfm=kr; random intercept Ecotype/ subject=Block(Site) ; *lsmeans Site*Ecotype / cl at _GDD=0.30 at GDD2=0.09 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ) ; *lsmeans Site*Ecotype / cl at _GDD=0.40 at GDD2=0.16 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ) ; *lsmeans Site*Ecotype / cl at _GDD=0.55 at GDD2=0.3025 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ); *lsmeans Site*Ecotype / cl at _GDD=0.70 at GDD2=0.49 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ); *lsmeans Site*Ecotype / cl at _GDD=0.80 at GDD2=0.64 slicediff=Site adjust=bon plot=meanplot(join sliceby=Ecotype ); lsmeans Site/cl at _GDD=0.40 at GDD2=0.16 pdiff adjust=bon ; lsmeans Site/cl at _GDD=0.50 at GDD2=0.25 pdiff adjust=bon ; lsmeans Site/cl at _GDD=0.60 at GDD2=0.36 pdiff adjust=bon ; lsmeans Site/cl at _GDD=0.70 at GDD2=0.49 pdiff adjust=bon ; ods output lsmeans=lsmeans Slicediffs=Slicediffs; *output out=diagnostics student=student predicted=predicted; nloptions tech=nrridg; run;