Figure and table goals! - lizzieinvancouver/heattolerance GitHub Wiki

Table of varieties we worked on

Give: Variety name // number of individuals (n) // how many of those flowered? // mean BB date (EL 7) // mean flowering date (skip, unless we have a lot of plants that did not go into experiment that flowered in greenhouse) // in experiment? (yes/no column)

... as a start we could make this in a csv file or Excel file. Focus only on the ones we would have included in experiment (see list in heattol_vars.xls).

We want a figure showing flowering by treatment ....

But first we need to figure out which data are best to show this. So:

  1. days to 25% and 50% flowering

  2. percflow (try mean across two stems, with na.rm=TRUE)

  3. bagbuds (try mean across two stems, with na.rm=TRUE)

  4. Soil moisture (last or first figure?)

SKIPPING: perflowering as measured by vFcount and vFexpec (try mean across two stems, with na.rm=TRUE) ... need to complete code in Flowgraphs that Lizzie started to get these estimates into percents....not the best measures so do not present. Also skipping caps. Flower buds were different sizes, so flowering data better seen in percent flowering and bagged buds. Could consider adding back in if needed for comparison to bagged buds and/or combine the two into total successful buds.

For each of these make a simple dotplot or boxplot (see ggplot) of each by treatment.

As a bonus! Try to make a plot of each data point versus day where day zero is the date the plant entered the chamber. Color code by variety. And facet_wrap by Treatment!

N leaves and stem length

For only the chamber data show:

  1. Calculate average leaves gained (or lost) from number when individual entered chamber; do the same for stem length. Then plot the results as dotplot or boxplot by treatment.

  2. n leaves and stem length by time (again: Try to make a plot of each data point versus day where day zero is the date the plant entered the chamber. Color code by variety and facet_wrap by Treatment).

Soil moisture

Dotplot and boxplots by treatment.

Update as of 20 April 2017

See chambgraphs.R notes ...

(1) Update figure noted as (a) with:

Extend y axis to -1.

Calculate an Anova on data by treat ... (install the car package!)

myobject <- lm(capfall_mean~Treat, data=XXX)

Anova(myobject)

Try to extract means from model and plot on ggplot .... OR figure out how to plot means on geom_point:

(2) Repeat as above for buds fallen into bag.

(3) Repeat as above for mean(percent flowering) and max(percent flowering) and min(percent flowering) ... so you'll need to calculate each of these at the individual level before plotting or doing any stats.

(4) Repeat as above for stem length change.

(5) Repeat as above for leaf number change.

(6) And! Repeat as above for 50% flowering.

Figure out how to visualize differences between stems in flowering, caps, stem, leaves etc.