Assignments

Tuesday, 7 May

1. Write a script to manipulate studyOneCleanData
   • Start with the ascii table of data
   • Import into a dataframe
   • Mean and SD of columns
   • Add a new column for average of high and low temperature
   • Graph Tempo vs. Mid-temperature
   • Save graph as pdf
   • Save the new data.frame as R object
   • Export the new table as tab-delimited ascii
2. Write a script to manipulate study3Opera
   • Start with the ascii table of data
   • Import into a dataframe
   • Mean and SD of columns
   • Mean authority for males vs. females
   • Mean tessitura for males vs. females
   • Graph tessitura vs. age, different colored points for males and females
   • Save graph as pdf

Wednesday, 8 May

1. Regression: —– data
   • Possible terms
     1. simple model: log(AGB) vs. log(dbh)
     2. second order term: log(dbh) squared
     3. rainfall and elevation
   • Graph
     1. log(AGB) vs. log(dbh)
     2. add curve of best fit
     3. overlay curves for high rainfall and low rainfall
   • Character variable (factors)
     1. use ForestType in the model
        • log(volume) ~ log(dbh) + ForestType
     2. compare to 3 independent models (3 forest types)
        • log(volume[dry]) ~ log(dbh[dry])
        • etc.
     2. Save a graph

Thursday, 9 May

1. Write a function to do regression, draw graph, add line
   1. Arguments x, y
   2. Save graph
   3. Loop through several x variables
2. Write a function with a loop

Friday, 10 May

1. Use lmer for regression of logagb on logdbh with species and forest type as factors
   1. Include squared term for logdbh
   2. Add locality as a group effect (does it change the fixed effect)
   3. Test forest type as fixed effect and as group effect
   4. Graph points and lines
   5. Compare alternative models
2. Use lmer for regression of pup Wt on momcat with year as a mixed effect
   1. Variable intercept, slope, or both
   2. Graph all points
   3. Use xyplot for groups (Lattice)
   4. Overlay lines of all random effects
   5. Compare alternative models
3. Use lmer for regression of log(ht) on log(dbh) with species as a mixed effect