Categories

file for Ian’s phenology analysis

Ian wants to quantify the overlap of flowering time between all pairs of plants in his experiment. This following R script reads his file with the flowering schedule of all 31 plants in his experiment and writes a file called ianPhenPairs.csv that has the flowering schedule of every possible pair combination (one per line). Note that there is a separate record for each plant as a sire and dam.

# script ian.phenology.r

pp <- read.csv("https://echinaceaproject.org/wp-content/uploads/ian.phenology.final.csv")

str(pp)

p <- merge(pp,pp, by = NULL)

str(p)

31^2 == dim(p)[1]

write.csv(p, “ianPhenPairs.csv”, row.names = FALSE)

twitterpinterestmail

Leave a Reply

You can use these HTML tags

<a href="" title=""> <abbr title=""> <acronym title=""> <b> <blockquote cite=""> <cite> <code> <del datetime=""> <em> <i> <q cite=""> <s> <strike> <strong>