Appendix 1 from PWR (Davies et al. 'Phylogenetically weighted regression: a method for modeling non-stationarity on evolutionary trees')

R code by Master Regetz
Questions? Email j.davies@ubc.ca

Files:

arnell.R - code to run  Arnell PWR analysis

data/arnell.phy - phylogeny for Arnell data

data/arnell.txt - Arnell flowering and seed size data

pwr-crossvalidation.R - PWR cross-validation (CV), includes f(x)s and CV of Arnell data 

pwr-functions.R - defines core PWR functions

pwr-plots.R - defines PWR plot helper functions

pwr-treesims.R - generates simple tree/data to evaluate the performance of PWR across various tree topologies (stemminess and imbalance)

simple.R - generates simple tree/data to illustrate the value of PWR over conventional global models like PGLS